VGaisserHillasParameter.h

Go to the documentation of this file.

#ifndef _evt_VGaisserHillasParameter_h_
#define _evt_VGaisserHillasParameter_h_


namespace evt {

  class VGaisserHillasParameter {

  public:
    typedef int IsClonableTag;

    VGaisserHillasParameter() { }
    VGaisserHillasParameter(const double nMax, const double xMax) : fXMax(xMax), fNMax(nMax) { }
    virtual ~VGaisserHillasParameter() { }

    virtual VGaisserHillasParameter* Create() const = 0;
    virtual VGaisserHillasParameter* Clone() const = 0;

    double GetXMax() const { return fXMax; }
    double GetXMaxError() const { return fXMaxError; }
    double GetNMax() const { return fNMax; }
    double GetNMaxError() const { return fNMaxError; }
    double GetNMaxXMaxCorrelation() const { return fRhoNMaxXMax; }

    double GetChiSquare() const { return fChiSqr; }
    unsigned int GetNdof() const { return fNdof; }

    bool IsdEdXProfile() const { return fIsEnergyDeposit; }

    void SetXMax(const double xMax, const double error)
    { fXMax = xMax; fXMaxError = error; }

    void SetNMax(const double nMax, const double error, const bool isEnergyDeposit = false)
    { fNMax = nMax; fNMaxError = error; fIsEnergyDeposit = isEnergyDeposit; }

    void SetNMaxXMaxCorrelation(const double rho)
    { fRhoNMaxXMax = rho; }

    void SetChiSquare(const double chi, const unsigned int ndof)
    { fChiSqr = chi; fNdof = ndof; }

    virtual double operator()(const double depth) const
    { return this->Eval(depth); }

    virtual double Eval(const double depth) const = 0;
    virtual double GetIntegral() const = 0;
    virtual double GetIntegralError() const = 0;

  protected:
    double fXMax = 0;
    double fXMaxError = 0;
    double fNMax = 0;
    double fNMaxError = 0;
    double fRhoNMaxXMax = 0;
    double fChiSqr = 0;
    unsigned int fNdof = 0;

    double fGammaIntegral = 0;
    double fGammaError = 0;

    bool fIsEnergyDeposit = false;

  };

}


#endif