HelixClass_double.h

#ifndef HelixClass_double_H
#define HelixClass_double_H 1
#include <vector>
#include "LineClass.h"
class HelixClass_double;

class HelixClass_double {
 public:

    HelixClass_double();
    ~HelixClass_double();
    void Initialize_VP(double * pos, double * mom, double q, double B);

    void Initialize_BZ(double xCentre, double yCentre, double radius, 
                                   double bZ, double phi0, double B, double signPz,
                                   double zBegin);
    void Initialize_Canonical(double phi0, double d0, double z0, double omega, 
                              double tanlambda, double B);
    const double * getMomentum();

    const double * getReferencePoint();

    double getPhi0();

    double getD0();

    double getZ0();

    double getOmega();

    double getTanLambda();

    double getPXY();


    double getXC();

    double getYC();


    double getRadius();   


    double getPointInXY(double x0, double y0, double ax, double ay, 
                              double * ref , double * point);

    double getPointInZ(double zLine, double * ref, double * point);

    double getDistanceToPoint(double * xPoint, double * Distance);

    double getDistanceToPoint(const double* xPoint, double distCut);
    double getDistanceToPoint(const std::vector<double>& xPoint, double distCut);

    double getPointOnCircle(double Radius, double * ref, double * point);

    double getDistanceToHelix(HelixClass_double * helix, double * pos, double * mom);

    void setHelixEdges(double * xStart, double * xEnd);

    double * getStartingPoint() {return _xStart;}

    double * getEndPoint() {return _xEnd;}
    
    double getBz();

    double getPhiZ();

    void getExtrapolatedMomentum(double * pos, double * momentum);

    double getCharge();

 private:    
    double _momentum[3]; // momentum @ ref point 
    double _referencePoint[3]; // coordinates @ ref point
    double _phi0; // phi0 in canonical parameterization 
    double _d0;   // d0 in canonical parameterisation
    double _z0;   // z0 in canonical parameterisation
    double _omega; // signed curvuture in canonical parameterisation
    double _tanLambda; // TanLambda 
    double _pxy; // Transverse momentum
    double _charge; // Particle Charge
    double _bField; // Magnetic field (assumed to point to Z>0)
    double _radius; // radius of circle in XY plane
    double _xCentre; // X of circle centre
    double _yCentre; // Y of circle centre
    double _phiRefPoint; // Phi w.r.t. (X0,Y0) of circle @ ref point
    double _phiAtPCA; // Phi w.r.t. (X0,Y0) of circle @ PCA 
    double _xAtPCA; // X @ PCA
    double _yAtPCA; // Y @ PCA
    double _pxAtPCA; // PX @ PCA
    double _pyAtPCA; // PY @ PCA
    double _phiMomRefPoint; // Phi of Momentum vector @ ref point
    double _const_pi; // PI
    double _const_2pi; // 2*PI
    double _const_pi2; // PI/2    
    double _FCT; // 2.99792458E-4
    double _xStart[3]; // Starting point of track segment
    double _xEnd[3]; // Ending point of track segment

    double _bZ;
    double _phiZ;

};


#endif