G/ChiZeroRegression.h

Go to the documentation of this file.

#ifndef FdProfileConstrainedGeometryFitPG_ChiZeroRegression_h
#define FdProfileConstrainedGeometryFitPG_ChiZeroRegression_h

#include <fevt/EyeRecData.h>
#include <utl/PhysicalConstants.h>

namespace FdProfileConstrainedGeometryFitPG {

class ChiZeroRegression {
 public:
  ChiZeroRegression(const fevt::Eye& eye, unsigned int& telId)
    {telId = PrepParams(eye);}

  void
    operator()(const double chi0,
         double& chi2,
         double& rp,
         double& t0,
         double Xmax=600*utl::g/utl::cm2) const;

  void
   operator()(const double chi0,
              double& chi2,
              double& rp, double &rpErr,
              double& t0, double &t0Err,
              double Xmax=600*utl::g/utl::cm2) const;

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

 protected:
  unsigned int
    PrepParams(const fevt::Eye& eye);

  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 std::vector<int> tels;
  static std::vector<int> pixids;
  static double fsdpTheta;
  static double fsdpPhi;
  static int fEyeId;
  static int fTelId;
  static double fXmax;
  bool fUseLightFlux;
  bool fRealAtm;
  bool fDeex;
  static bool fEmissionPointCorrection;
};
} //namespace

#endif