ChiZeroRegression.h

Go to the documentation of this file.

#ifndef FdProfileConstrainedGeometryFit_ChiZeroRegression_h
#define FdProfileConstrainedGeometryFit_ChiZeroRegression_h

#include <fevt/EyeRecData.h>

namespace FdProfileConstrainedGeometryFit {

class ChiZeroRegression {
 public:
  ChiZeroRegression(const fevt::EyeRecData& eyeRec)
    {PrepParams(eyeRec);}

  void
    operator()(const double chi0,
         double& chi2,
         double& rp,
         double& t0) const;

  void
   operator()(const double chi0,
              double& chi2,
              double& rp, double &rpErr,
              double& t0, double &t0Err) const;

  unsigned int GetNDof() const { return chii.size() - 2; }
  void SetUseLightFlux(bool use=true) {fUseLightFlux=use;}
  void SetRealAtm(bool realAtm, bool deex);

 protected:
  void
    PrepParams(const fevt::EyeRecData& eyeRec);

  void
    LinearFitErrors(const std::vector<double>& x,
                    const std::vector<double>& y,
                    const std::vector<double>& ey,
                    double& a0, double& a0err,
                    double& a1, double& a1err,
                    double& chi2) const;

  static void MinuitFitFunc(int& npar, double* gin, double& f, double* par, int iflag);

  void
    MinuitFitErrors(const double Chi0,
                    double& T0, double& T0err,
                    double& Rp, double& Rperr,
                    double& chi2) const;

  
  static std::vector<double> times;
  static std::vector<double> timeErrs;
  static std::vector<double> chii;
  static double fsdpTheta;
  static int fEyeId;
  static int fTelId;
  bool fUseLightFlux;
  bool fRealAtm;
  bool fDeex;
};
} //namespace

#endif