NeutrinoFitObject.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 __NEUTRINOFITOBJECT_H
#define __NEUTRINOFITOBJECT_H
 
#include "analysis/OrcaKinFit/ParticleFitObject.h"
 
#include <cmath>
 
namespace Belle2 {
  namespace OrcaKinFit {
 
    class NeutrinoFitObject : public ParticleFitObject {
    public:
      NeutrinoFitObject(double E, double theta, double phi,
                        double DE = 1, double Dtheta = 0.1, double Dphi = 0.1);
 
      NeutrinoFitObject(const NeutrinoFitObject& rhs               
                       );
      NeutrinoFitObject& operator= (const NeutrinoFitObject& rhs   
                                   );
 
      virtual ~NeutrinoFitObject();
 
      virtual NeutrinoFitObject* copy() const override;
 
      virtual NeutrinoFitObject& assign(const BaseFitObject& source    
                                       ) override;
 
      virtual const char* getParamName(int ilocal      
                                      ) const override;
 
      virtual bool   updateParams(double p[],    
                                  int idim      
                                 ) override;
 
      // 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 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 int getNPar() const override {return NPAR;}
 
    protected:
      void updateCache() const override;
 
      enum {NPAR = 3};
 
      mutable double ctheta, stheta, cphi, sphi,
              pt, px, py, pz, dptdE,
              dpxdE, dpydE, dpxdtheta, dpydtheta,
              chi2;
 
    };
 
  }// end OrcaKinFit namespace
} // end Belle2 namespace
 
 
#endif // __NEUTRINOFITOBJECT_H