ISRPhotonFitObject.h

 
#ifndef __ISRPHOTONFITOBJECT_H
#define __ISRPHOTONFITOBJECT_H

#include "ParticleFitObject.h"



class ISRPhotonFitObject : public ParticleFitObject {
  public:
    ISRPhotonFitObject(double px, double py, double pz,                   
                        double b_, double PzMaxB_, double PzMinB_ = 0.);  
    
    ISRPhotonFitObject (const ISRPhotonFitObject& rhs              
                   );
    ISRPhotonFitObject& operator= (const ISRPhotonFitObject& rhs   
                             );
                             
    virtual ~ISRPhotonFitObject();
    
    virtual ISRPhotonFitObject *copy() const;
    
    virtual ISRPhotonFitObject& 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 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 int getNPar() const {return NPAR;}
  
  protected:
    
    enum {NPAR=3}; // well, it's actually 1...Daniel should update

    virtual double PgFromPz(double pz);
    
    void updateCache() const;
  
    mutable bool cachevalid;
    
    mutable double pt2, p2, p, pz,
                   dpx0, dpy0, dpz0, dE0, dpx1, dpy1, dpz1, dE1,
                   dpx2, dpy2, dpz2, dE2, d2pz22, d2E22,
                   chi2,                   
                   b, PzMinB, PzMaxB, dp2zFact;
};

#endif // __ISRPHOTONFITOBJECT_H