NeutralParticleTrack.h

Go to the documentation of this file.

#ifndef __NeutralParticleTrack_H
#define __NeutralParticleTrack_H

#include "jbltools/kinfit/TrackFitObject.h"

#include <cassert>

// Class NeutralParticleTrack

class NeutralParticleTrack: public TrackFitObject {
  public:
    NeutralParticleTrack (const char *name_,
                          double pt,
                          double phi0,
                          double theta,
                          double dca,
                          double z0,
                          double mass,
                          double sstart,
                          double sstop
                          );
    NeutralParticleTrack (const char *name_,
                          double pt,
                          double phi0,
                          double theta,
                          double dca,
                          double z0,
                          double mass,
                          double sstart
                          );
    NeutralParticleTrack (const char *name_,
                          const ThreeVector& vertex,
                          const ThreeVector& momentum,
                          double mass
                         );
    virtual ~NeutralParticleTrack();

    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 double getMass () 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 void addToGlobCov(double *glcov, int idim) 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 bool fixVertexParam (int ivertex,    
                                 bool fix=true  
                                 ) { return fixParam (6+ivertex, fix); }
    
    static const double parfact[NPARMAX];
  
  protected:
    
    virtual void initCov(const float cov_[15]);
    
    void updateCache() const;
  
    enum {NPAR = 8};
    
    // Cache variables
    mutable double pt;            
    mutable double phi0;          
    mutable double theta;         
    mutable double dca;           
    mutable double z0;            
    mutable double mass;          
    
    mutable double sphi0;              
    mutable double cphi0;              
    mutable double px;                 
    mutable double py;                 
    mutable double x0;                 
    mutable double y0;                 
    
    mutable double cottheta;           
    mutable double sintheta;           
    mutable double sin2theta;          
    mutable double momentum;           
    mutable double momderfact;         
    mutable double energy;             
    mutable double beta;               
  
};

#endif // __NeutralParticleTrack_H

---