NeutrinoFitObject.h

#ifndef __NEUTRINOFITOBJECT_H
#define __NEUTRINOFITOBJECT_H

#include "ParticleFitObject.h"

#include <cmath>


// Class NeutrinoFitObject

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;
    
    virtual NeutrinoFitObject& 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:
    void updateCache() const;
    
    enum {NPAR=3};
  
    mutable double ctheta, stheta, cphi, sphi,
                   pt, px, py, pz, dptdE, 
                   dpxdE, dpydE, dpxdtheta, dpydtheta,
                   chi2;

};

#endif // __NEUTRINOFITOBJECT_H