BCHistogramFitter.h

Go to the documentation of this file.

#ifndef __BCMODELHISTOGRAMFITTER__H
#define __BCMODELHISTOGRAMFITTER__H

/*
 * Copyright (C) 2008-2011, Daniel Kollar and Kevin Kroeninger.
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */

// ---------------------------------------------------------

#include <vector>

#include "BCModel.h"

class BCHistogramFitterToyDataInterface;

// ROOT classes
class TH1D;
class TF1;

// ---------------------------------------------------------

class BCHistogramFitter : public BCModel
{
   public:

      /* @{ */

      BCHistogramFitter();

      BCHistogramFitter(TH1D * hist, TF1 * func);

      ~BCHistogramFitter();

      /* @} */

      /* @{ */

      TH1D * GetHistogram()
         { return fHistogram; };

      TH1D * GetHistogramExpected()
         { return fHistogramExpected; };


      TF1 * GetFitFunction()
         { return fFitFunction; };

      TGraph * GetErrorBand()
         { return fErrorBand; };

      TGraph * GetGraphFitFunction()
         { return fGraphFitFunction; };

      /* @} */

      /* @{ */

      int SetHistogram(TH1D * hist);

      int SetHistogramExpected(const std::vector <double>& parameters);

      int SetFitFunction(TF1 * func);

      void SetFlagIntegration(bool flag)
         { fFlagIntegration = flag; };

      /* @} */
      /* @{ */

//    virtual double LogAPrioriProbability(std::vector <double> parameters);

      virtual double LogLikelihood(std::vector <double> parameters);

//    void PrintHistogram(const char * options = "", const char * filename = "");

      double FitFunction(std::vector <double> x, std::vector <double> parameters);

      int Fit();

      int Fit(TH1D * hist, TF1 * func);

      void DrawFit(const char * options = "HIST", bool flaglegend = false);

      int CalculatePValueFast(std::vector<double> par, double &pvalue, int nIterations = 100000);

      int CalculatePValueFast(std::vector<double> par, BCHistogramFitterToyDataInterface* callback,
            double &pvalue,  int nIterations = 100000);

      int CalculatePValueLikelihood(std::vector<double> par, double &pvalue);

      int CalculatePValueLeastSquares(std::vector<double> par, double &pvalue, bool weightExpect=true);

      int CalculatePValueKolmogorov(std::vector<double> par, double &pvalue);


      double CDF(const std::vector<double>& parameters, int index, bool lower=false);

      /* @} */

   private:

      TH1D * fHistogram;

      TF1 * fFitFunction;

      bool fFlagIntegration;

      TGraph * fErrorBand;

      TGraph * fGraphFitFunction;

      TH1D * fHistogramExpected;

};

// ---------------------------------------------------------

class BCHistogramFitterToyDataInterface
{
public:
   virtual void operator()(const std::vector<double>& expectation, const std::vector<int>& toyData)=0;

   virtual ~BCHistogramFitterToyDataInterface(){}
};

#endif