JetFitObject.h

Go to the documentation of this file.

#ifndef __JETFITOBJECT_H
#define __JETFITOBJECT_H

#include "ParticleFitObject.h"

// Class JetFitObject

class JetFitObject : public ParticleFitObject {
  public:
    JetFitObject(double E, double theta, double phi, 
                 double DE, double Dtheta, double Dphi, 
                 double m = 0);
                 
    JetFitObject (const JetFitObject& rhs              
                   );
    JetFitObject& operator= (const JetFitObject& rhs   
                             );

    virtual ~JetFitObject();
    
    virtual JetFitObject *copy() const;
    
    virtual JetFitObject& assign (const BaseFitObject& source   
                                 );
    
    virtual const char *getParamName (int ilocal     
                                     ) const;
 
    virtual bool   updateParams (double p[],   
                                 int idim      
                                );  

    virtual int getNPar() const {return NPAR;}
        
    // 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 getCov (int ilocal,    
                           int jlocal     
                           ) const ;

    virtual double getError (int ilocal     
                             ) const;

    
//    virtual void   addToDerivatives (double der[],      ///< Derivatives vector, length idim
//                                     int idim,          ///< Length of derivatives vector
//                                     double efact=0,    ///< Factor for dE/dx_i
//                                     double pxfact=0,   ///< Factor for dpx/dx_i
//                                     double pyfact=0,   ///< Factor for dpy/dx_i 
//                                     double pzfact=0    ///< Factor for dpz/dx_i
//                                     ) const;

    // daniel's new method
    // derivatives of intermediate variable wrt local variable
    virtual double getFirstDerivative_Meta_Local( int iMeta, int ilocal , int metaSet ) const;
    virtual double getSecondDerivative_Meta_Local( int iMeta, int ilocal , int jlocal , int metaSet ) const;

    //    virtual double getChi2() const;

  protected:
    
    enum {NPAR=3};
    
    void updateCache() const;

    mutable double ctheta, stheta, cphi, sphi,
                   p2, p, pt, px, py, pz, dpdE, dptdE, 
                   dpxdE, dpydE, dpzdE, dpxdtheta, dpydtheta,
                   chi2;
                   // d2pdE2, d2ptsE2;
                   
    static bool adjustEThetaPhi (double& m, double &E, double& theta, double& phi);
    
    //    double calcChi2 () const;
    
};



#endif // __JETFITOBJECT_H