SimpleHelix.h

#ifndef SimpleHelix_H
#define SimpleHelix_H 1

#include "Trajectory.h"

class SimpleHelix : public Trajectory {

public:

  virtual ~SimpleHelix() {} 
 
  SimpleHelix( double d0, double phi0, double omega,
               double z0, double tanLambda, 
               LCVector3D referencePoint, LCErrorMatrix* errors=0) ;
  
  virtual LCVector3D getPosition(double s, LCErrorMatrix* errors=0) const ;
  
  virtual LCVector3D getDirection(double s,  LCErrorMatrix* errors=0) const ;
  
  virtual LCErrorMatrix getCovarianceMatrix( double s) const ;

  virtual double getPathAt(const LCVector3D position ) const ;
  
  virtual double getRadius() const ;

  /*----------------------------------------------------------------------*/
  
  virtual double getIntersectionWithPlane( LCPlane3D p, bool& pointExists) const  ;
  
  virtual  double getIntersectionWithCylinder(const LCCylinder & cylinder,
                                              bool & pointExists) const ;

  virtual double getStart() const ;
  virtual double getEnd() const ;

  virtual bool setStart(double s);
  virtual bool setEnd(double s);
  virtual bool setStartEnd(double start, double end);

  virtual void printProperties();

protected:

  SimpleHelix() {} 

  virtual void init() ;

  virtual double getCentreX() const ;
  virtual double getCentreY() const ;
  virtual double getWindingLength() const ;
  virtual double getPitch();
  double _d0;
  double _phi0;
  double _omega;
  double _z0;
  double _tanLambda;

  //  double _Bz;

  double _helixStart;
  double _helixEnd;

  static const char* _names[];

  static const double _a; // = 2.99792458E-4;
  static const double _pi; // = 3.14159265358979323846;

  LCVector3D  _reference; 
  LCErrorMatrix* _errors;
}; // class 



// /** Physical trajectory describing a (charged) particle's  path in a B 
//  *  field and material. 
//  *  @author F.Gaede, DESY
//  *  @version $Id: SimpleHelix.h,v 1.7 2007-06-20 18:47:25 samson Exp $
//  */

// class PhysicalSimpleLine : public SimpleLine{
  
//   /** Particle's momentum at path length s. Implementations will have to  have knowledge
//    *  about the particle type, B-field and material.
//    */
//   virtual LCLorentzVector get4Momentum( double s ) const ;
// }


#endif /* ifndef SimpleLine_H */