LCLine3D.h

#ifndef LCLine3D_H
#define LCLine3D_H 1

// #include "CLHEP/Vector/ThreeVector.h"
#include <LCGeometryTypes.h>
#include <LCPlane3D.h>

class LCLine3D {

public:

  LCLine3D();

  LCLine3D(const LCVector3D & point, const LCVector3D & direction) ;

  LCLine3D(const LCVector3D & point, 
           const LCVector3D & direction, 
           const LCVector3D & reference) ;

  LCLine3D(double d0, double phi0, double z0, double tanLambda) ;

  LCLine3D(double d0, double phi0, double z0, double tanLambda,
           const LCVector3D & reference) ;

  LCLine3D(const LCLine3D & line) ;

  ~LCLine3D() {}

  bool set(const LCVector3D & point, 
           const LCVector3D & direction, 
           const LCVector3D & reference) ;

  bool set(double d0, double phi0, double z0, double tanLambda,
           const LCVector3D & reference) ;

  LCLine3D & operator=(const LCLine3D & rhs) ;

  LCVector3D position(const double s = 0) const ;

  LCVector3D direction() const ;

  double distance(const LCVector3D & point) const ;  

  double projectPoint(const LCVector3D & point) const ;

  bool operator==(const LCLine3D & rhs) const ; 

  bool operator!=(const LCLine3D & rhs) const ; 

  double intersectionWithPlane(const LCPlane3D plane, bool& pointExists) const ;

protected:

  LCVector3D _point;
  LCVector3D _direction;
  LCVector3D _reference;
};

std::ostream & operator << (std::ostream &os, const LCLine3D &l) ;

#endif /* ifndef LCLine3D_H */