QuadraticFitData.h

Go to the documentation of this file.

#ifndef _utl_QuadraticFitData_h_
#define _utl_QuadraticFitData_h_

#include <utl/QuadraticFitter-fwd.h>


namespace utl {

  class QuadraticFitData {
  public:
    double GetStart() const { return fStart; }
    double GetStop() const { return fStop; }
    double GetExtremePosition() const { return fExtremePosition; }
    double GetExtremePositionError() const { return fExtremePositionError; }
    double GetChi2() const { return fChi2; }
    int GetNdof() const { return fNdof; }
    const double* GetCoefficients() const { return fCoefficients; }
    double GetCoefficient(unsigned int i) const { return fCoefficients[i]; }

    double
    operator()(const double x)
      const
    {
      const double z = x - fExtremePosition;
      return fCoefficients[0] + z*(fCoefficients[1] + z*fCoefficients[2]);
    }

    void SetStart(const double start) { fStart = start; }
    void SetStop(const double stop) { fStop = stop; }
    void SetExtremePosition(const double pos, const double posErr = 0)
    { fExtremePosition = pos; fExtremePositionError = posErr; }
    void SetChi2Ndof(const double chi2, const int ndof)
    { fChi2 = chi2; fNdof = ndof; }
    void SetCoefficients(const double coeff[3])
    { for (int i = 0; i < 3; ++i) fCoefficients[i] = coeff[i]; }

  private:
    double fStart = 0;
    double fStop = 0;
    double fExtremePosition = 0;
    double fExtremePositionError = 0;
    double fChi2 = 0;
    int fNdof = 0;
    double fCoefficients[3] = { 0 };

    template<class Histogram, class ErrorPolicy>
    friend class utl::QuadraticFitter;
  };

}


#endif