JetFitObject.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * See https://github.com/tferber/OrcaKinfit, forked from                 *
 * https://github.com/iLCSoft/MarlinKinfit                                *
 *                                                                        *
 * Further information about the fit engine and the user interface        *
 * provided in MarlinKinfit can be found at                               *
 * https://www.desy.de/~blist/kinfit/doc/html/                            *
 * and in the LCNotes LC-TOOL-2009-001 and LC-TOOL-2009-004 available     *
 * from http://www-flc.desy.de/lcnotes/                                   *
 *                                                                        *
 * Adopted by: Torben Ferber (torben.ferber@desy.de) (TF)                 *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#ifndef __JETFITOBJECT_H
#define __JETFITOBJECT_H
 
#include "analysis/OrcaKinFit/ParticleFitObject.h"
 
namespace Belle2 {
  namespace OrcaKinFit {
 
// 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 override;
 
      virtual JetFitObject& assign(const BaseFitObject& source    
                                  ) override;
 
      virtual const char* getParamName(int ilocal      
                                      ) const override;
 
      virtual bool   updateParams(double p[],    
                                  int idim      
                                 ) override;
 
      virtual int getNPar() const override {return NPAR;}
 
      // these depend on actual parametrisation!
 
      virtual double getDPx(int ilocal) const override;
      virtual double getDPy(int ilocal) const override;
      virtual double getDPz(int ilocal) const override;
      virtual double getDE(int ilocal) const override;
 
      virtual double getCov(int
                            ilocal,     
                            int jlocal     
                           ) const override;
 
      virtual double getError(int ilocal      
                             ) const override;
 
 
//    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 override;
      virtual double getSecondDerivative_Meta_Local(int iMeta, int ilocal , int jlocal , int metaSet) const override;
 
      //    virtual double getChi2() const;
 
    protected:
 
      enum {NPAR = 3};
 
      void updateCache() const override;
 
      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(const double& m, double& E, double& theta, double& phi);
 
      //    double calcChi2 () const;
 
    };
 
 
  }// end OrcaKinFit namespace
} // end Belle2 namespace
 
 
#endif // __JETFITOBJECT_H