ChargedParticleTrack.h

Go to the documentation of this file.

#ifndef __CHARGEDPARTICLETRACK_H
#define __CHARGEDPARTICLETRACK_H

#include "jbltools/kinfit/TrackFitObject.h"

#include <cassert>

// Class ChargedParticleTrack

class ChargedParticleTrack: public TrackFitObject {
  public:
    ChargedParticleTrack (const char *name_,
                          double kappa,
                          double phi0,
                          double theta,
                          double dca,
                          double z0,
                          double mass,
                          double charge_,
                          double sstart,
                          double sstop
                          );
    ChargedParticleTrack (const char *name_,
                          double kappa,
                          double phi0,
                          double theta,
                          double dca,
                          double z0,
                          double mass,
                          double charge_,
                          double sstart
                          );
    ChargedParticleTrack (const char *name_,
                          const ThreeVector& vertex,
                          const ThreeVector& momentum,
                          double mass_,
                          double charge_ = 1
                         );
    ChargedParticleTrack (const char *name_,
                          const float par_[5],
                          const float cov_[15],
                          double mass_,
                          double charge_ = 1,
                          double sstart = 0,
                          double sstop = 0
                         );
    virtual ~ChargedParticleTrack();

    virtual int getNPar() const;
    
    virtual bool   setParam (int ilocal,         
                             double par_,        
                             bool measured_,     
                             bool fixed_ = false 
                            );  
    virtual bool   setParam (int ilocal,         
                             double par_        
                            );  

    virtual bool setParameters (int ivertex,                  
                                const ThreeVector& vertex,    
                                const FourVector& momentum,   
                                double charge_                
                               );
                                
    virtual const char *getParamName (int ilocal     
                                     ) const;
 
    virtual void getTrajectoryPointEx (double s,  
                                       ThreeVector& p 
                                      ) const; 
                                    
    virtual void getVertexEx (int ivertex, 
                              ThreeVector& p
                             ) const; 

    virtual void setVertex (int ivertex,         
                             const TwoVector& p  
                           ); 
    
    virtual void getTrajectoryDerivativeEx (double s, 
                                            int ilocal, 
                                            ThreeVector& p
                                           ) const;
    virtual void getVertexDerivativeEx (int ivertex, 
                                        int ilocal, 
                                        ThreeVector& p
                                       ) const;

    virtual void getMomentumAtTrajectoryEx (double s, 
                                            FourVector& p
                                           ) const; 
    virtual void getMomentumEx (int ivertex, 
                                FourVector& p 
                               ) const; 

    virtual double getCharge () const; 

    virtual void getMomentumDerivativeAtTrajectoryEx (double s, 
                                                      int ilocal, 
                                                      FourVector& p
                                                     ) const; 
    virtual void getMomentumDerivativeEx (int ivertex, 
                                          int ilocal, 
                                          FourVector& p
                                         ) const;
    virtual double getArcLength (int i
                                ) const;
       
    virtual double getDChi2DParam(int ilocal) const {assert (0);};
    virtual double getD2Chi2DParam2(int ilocal1, int ilocal2) const {assert (0);};
    
    virtual void addToGlobalChi2DerMatrix (int idim, double *M) const {assert (0);};
    virtual void addToGlobalDerMatrix (int idim, double c, double *M) const {assert (0);};

    virtual JBLHelix getTangentialHelix (double s     
                                        );
    virtual double getMass() const;
    
    static const double parfact[NPARMAX];
    
    virtual bool fixVertexParam (int ivertex,    
                                 bool fix=true  
                                 ) { return fixParam (6+ivertex, fix); }
  
  protected:
    void updateCache() const;
    
    virtual void initCov(const float cov_[15]);
    
    double getNormalS (double s) const;
  
    enum {NPAR = 8};
    
    // Cache variables
    mutable double kappa;         
    mutable double phi0;          
    mutable double theta;         
    mutable double dca;           
    mutable double z0;            
    mutable double mass;          
    mutable double r;             
    mutable double sphi0;         
    mutable double cphi0;         
    mutable double dcamir;        
    mutable double cottheta;      
    mutable double sintheta;      
    mutable double sin2theta;     
    mutable double cBq;           
    mutable double pt;            
    mutable double momentum;      
    mutable double momderfact;    
    mutable double energy;        
    mutable double beta;          
  
    double charge;                


};

#endif // __CHARGEDPARTICLETRACK_H

---