StableChargedParticle.h

Go to the documentation of this file.

#ifndef __StableChargedParticle_H
#define __StableChargedParticle_H

#include"FitParticle.h"


//  Class StableChargedParticle:

class StableChargedParticle : public FitParticle {
  public:
    StableChargedParticle();
    virtual ~StableChargedParticle();
    
    
    // these depend on actual parametrisation!
    
    virtual void   addToDerivatives (double der[],
                                     int idim,
                                     double pxfact=0, 
                                     double pyfact=0, 
                                     double pzfact=0, 
                                     double efact=0) const;

    // add to matrix der2 of size idim x idim
    // pxfact*d^2px/(dx_i dx_j) + pyfact...
    virtual void   addTo2ndDerivatives (double der2[],
                                        int idim,
                                        double pxfact, 
                                        double pyfact, 
                                        double pzfact, 
                                        double efact) const;
    
    virtual void setCovMatrix(double values[NMEAS*(NMEAS+1)/2]);
    
    virtual double getChi2() const;
    virtual double getDChi2DParam(int ilocal) const;
    virtual double getD2Chi2DParam2(int ilocal1, int ilocal2) const;
    
    virtual void addToGlobalChi2DerMatrix (int idim, double *M) const;
    virtual void addToGlobalDerMatrix (int idim, double c, double *M) const;
    
    virtual void addToGlobCov(double *cov, int idim) const; // Add covariance matrix elements to
                                         // global covariance matrix of size idim x idim

    virtual void invalidateCache() const;
  
  protected:
    void updateCache() const;
    
    
    void setVertex (FitThreeVector& v, FitThreeVector dv[NPAR], 
                    int n, double s) const;
    
    enum {NPAR = 6};
    enum {NMEAS = 5};
    
    double mass;
    double bfield;
  
    mutable bool cachevalid;
    
    mutable double chi2;
    mutable double r, sphi0, cphi0, dcamir, cottheta, cos2theta;
    mutable double resid[NMEAS];
    double cov[NMEAS][NMEAS];
    double covinv[NMEAS][NMEAS];
    mutable FitThreeVector v1;         
    mutable FitThreeVector dv1[NPAR];  
  
};



#endif // __StableChargedParticle_H

---