LeptonFitObject.h

Go to the documentation of this file.

#ifndef __LEPTONFITOBJECT_H
#define __LEPTONFITOBJECT_H

#include "ParticleFitObject.h"
#include "EVENT/Track.h"
#include "lcio.h"

using namespace lcio ;

// Class LeptonFitObject
// Class for leptons with (q/pt, theta, phi) in kinematic fits

class LeptonFitObject : public ParticleFitObject {
  public:

    LeptonFitObject(double ptinv, double theta, double phi,
                  double Dptinv, double Dtheta, double Dphi,
                  double m = 0);

    LeptonFitObject(double ptinv, double theta, double phi,
                  double Dptinv, double Dtheta, double Dphi,
                  double Rhoptinvtheta, double Rhoptinvphi, double Rhothetaphi,
                  double m = 0);

    LeptonFitObject(Track* track, double Bfield, double m = 0);

    LeptonFitObject(const TrackState* trackstate, double Bfield, double m = 0);
                 
    LeptonFitObject (const LeptonFitObject& rhs              
                   );
    LeptonFitObject& operator= (const LeptonFitObject& rhs   
                             );

    virtual ~LeptonFitObject();

    virtual LeptonFitObject *copy() const;

    virtual LeptonFitObject& assign (const BaseFitObject& source   
                                   );

    virtual const char *getParamName (int ilocal     
                                      ) const;

    virtual bool   updateParams (double p[],   
                                 int idim      
                                 );

    // these depend on actual parametrisation!
    virtual double getDPx(int ilocal) const;
    virtual double getDPy(int ilocal) const;
    virtual double getDPz(int ilocal) const;
    virtual double getDE(int ilocal) const;

    //  virtual double getChi2() const;

    double getFirstDerivative_Meta_Local( int iMeta, int ilocal , int metaSet ) const; // first derivative of intermediate variable iMeta wrt local parameter ilocal
    double getSecondDerivative_Meta_Local( int iMeta, int ilocal , int jlocal, int metaSet ) const; // second derivative of intermediate variable iMeta wrt local parameters ilocal and jlocal
    virtual int getNPar() const {return NPAR;}

  protected:

    //  virtual void initCov();

    void updateCache() const;

    //    mutable bool cachevalid;

    mutable double ctheta, stheta, stheta2, cphi, sphi, cottheta,
      p2, p, e, e2, pt, pt2, pt3, px, py, pz, dpdptinv, dpdtheta, dptdptinv,
      dpxdptinv, dpydptinv, dpzdptinv, dpxdtheta, dpydtheta, dpzdtheta, dpxdphi, dpydphi, dpzdphi,
      chi2, dEdptinv, dEdtheta, dEdp, qsign, ptinv2;

    static bool adjustPtinvThetaPhi (double& m, double &ptinv, double& theta, double& phi);

    enum {NPAR=3};

};



#endif // __LEPTONFITOBJECT_H