fit::HistoPoissonLikelihoodRatio< T > Class Template Reference

#include <HistoPoissonLikelihoodRatio.h>

Public Member Functions
T &	function ()
const T &	function () const
	HistoPoissonLikelihoodRatio ()
	HistoPoissonLikelihoodRatio (T &t, TH1 *histo, double rangeMin, double rangeMax)
size_t	numberOfBins () const
double	operator() () const
void	setHistos (TH1 *histo)

Private Attributes
std::vector< double >	cont_
double	deltaX_
size_t	nBins_
double	rangeMax_
double	rangeMin_
T *	t_
double	xMax_
double	xMin_

Detailed Description

template<typename T>
class fit::HistoPoissonLikelihoodRatio< T >

Definition at line 11 of file HistoPoissonLikelihoodRatio.h.

Constructor & Destructor Documentation

template<typename T>

fit::HistoPoissonLikelihoodRatio< T >::HistoPoissonLikelihoodRatio ( )

inline

Definition at line 13 of file HistoPoissonLikelihoodRatio.h.

{}

template<typename T>

fit::HistoPoissonLikelihoodRatio< T >::HistoPoissonLikelihoodRatio ( T &  t, TH1 *  histo, double  rangeMin, double  rangeMax  )

inline

Definition at line 14 of file HistoPoissonLikelihoodRatio.h.

References fit::HistoPoissonLikelihoodRatio< T >::cont_, fit::HistoPoissonLikelihoodRatio< T >::deltaX_, mps_fire::i, fit::HistoPoissonLikelihoodRatio< T >::nBins_, fit::HistoPoissonLikelihoodRatio< T >::xMax_, and fit::HistoPoissonLikelihoodRatio< T >::xMin_.

        : t_(&t), rangeMin_(rangeMin), rangeMax_(rangeMax) {
      nBins_ = histo->GetNbinsX();
      xMin_ = histo->GetXaxis()->GetXmin();
      xMax_ = histo->GetXaxis()->GetXmax();
      deltaX_ = (xMax_ - xMin_) / nBins_;
      for (size_t i = 0; i < nBins_; ++i) {
        cont_.push_back(histo->GetBinContent(i + 1));
      }
    }

Member Function Documentation

template<typename T>

T& fit::HistoPoissonLikelihoodRatio< T >::function ( )

inline

Definition at line 52 of file HistoPoissonLikelihoodRatio.h.

References fit::HistoPoissonLikelihoodRatio< T >::t_.

{ return *t_; }

template<typename T>

const T& fit::HistoPoissonLikelihoodRatio< T >::function ( ) const

inline

Definition at line 53 of file HistoPoissonLikelihoodRatio.h.

References fit::HistoPoissonLikelihoodRatio< T >::t_.

{ return *t_; }

template<typename T>

size_t fit::HistoPoissonLikelihoodRatio< T >::numberOfBins ( ) const

inline

Definition at line 43 of file HistoPoissonLikelihoodRatio.h.

References fit::HistoPoissonLikelihoodRatio< T >::deltaX_, mps_fire::i, fit::HistoPoissonLikelihoodRatio< T >::nBins_, fit::HistoPoissonLikelihoodRatio< T >::rangeMax_, fit::HistoPoissonLikelihoodRatio< T >::rangeMin_, and fit::HistoPoissonLikelihoodRatio< T >::xMin_.

Referenced by fit::RootMinuitResultPrinter< HistoPoissonLikelihoodRatio< T > >::print().

                                {
      size_t fullBins = 0;
      for (size_t i = 0; i < nBins_; ++i) {
        double x = xMin_ + (i + .5) * deltaX_;
        if ((x > rangeMin_) && (x < rangeMax_))
          fullBins++;
      }
      return fullBins;
    }

template<typename T>

double fit::HistoPoissonLikelihoodRatio< T >::operator() ( ) const

inline

Definition at line 24 of file HistoPoissonLikelihoodRatio.h.

References fit::HistoPoissonLikelihoodRatio< T >::cont_, fit::HistoPoissonLikelihoodRatio< T >::deltaX_, mps_fire::i, log, fit::HistoPoissonLikelihoodRatio< T >::nBins_, fit::HistoPoissonLikelihoodRatio< T >::rangeMax_, fit::HistoPoissonLikelihoodRatio< T >::rangeMin_, and fit::HistoPoissonLikelihoodRatio< T >::xMin_.

                              {
      double chi2lambda = 0;
      for (size_t i = 0; i < nBins_; ++i) {
        double x = xMin_ + (i + .5) * deltaX_;
        if ((x > rangeMin_) && (x < rangeMax_)) {
          double nu = (*t_)(x);
          if (nu > 0 && cont_[i] > 0)
            chi2lambda += nu - cont_[i] + cont_[i] * log(cont_[i] / nu);
        }
      }
      chi2lambda *= 2;
      return chi2lambda;
    }

template<typename T>

void fit::HistoPoissonLikelihoodRatio< T >::setHistos ( TH1 *  histo )

inline

Definition at line 37 of file HistoPoissonLikelihoodRatio.h.

References fit::HistoPoissonLikelihoodRatio< T >::deltaX_, fit::HistoPoissonLikelihoodRatio< T >::nBins_, fit::HistoPoissonLikelihoodRatio< T >::xMax_, and fit::HistoPoissonLikelihoodRatio< T >::xMin_.

                               {
      nBins_ = histo->GetNbinsX();
      xMin_ = histo->GetXaxis()->GetXmin();
      xMax_ = histo->GetXaxis()->GetXmax();
      deltaX_ = (xMax_ - xMin_) / nBins_;
    }

Member Data Documentation

template<typename T>

std::vector<double> fit::HistoPoissonLikelihoodRatio< T >::cont_

private

Definition at line 60 of file HistoPoissonLikelihoodRatio.h.

Referenced by fit::HistoPoissonLikelihoodRatio< T >::HistoPoissonLikelihoodRatio(), and fit::HistoPoissonLikelihoodRatio< T >::operator()().

template<typename T>

double fit::HistoPoissonLikelihoodRatio< T >::deltaX_

private

Definition at line 59 of file HistoPoissonLikelihoodRatio.h.

Referenced by fit::HistoPoissonLikelihoodRatio< T >::HistoPoissonLikelihoodRatio(), fit::HistoPoissonLikelihoodRatio< T >::numberOfBins(), fit::HistoPoissonLikelihoodRatio< T >::operator()(), and fit::HistoPoissonLikelihoodRatio< T >::setHistos().

template<typename T>

size_t fit::HistoPoissonLikelihoodRatio< T >::nBins_

private

Definition at line 58 of file HistoPoissonLikelihoodRatio.h.

Referenced by fit::HistoPoissonLikelihoodRatio< T >::HistoPoissonLikelihoodRatio(), fit::HistoPoissonLikelihoodRatio< T >::numberOfBins(), fit::HistoPoissonLikelihoodRatio< T >::operator()(), and fit::HistoPoissonLikelihoodRatio< T >::setHistos().

template<typename T>

double fit::HistoPoissonLikelihoodRatio< T >::rangeMax_

private

Definition at line 57 of file HistoPoissonLikelihoodRatio.h.

Referenced by fit::HistoPoissonLikelihoodRatio< T >::numberOfBins(), and fit::HistoPoissonLikelihoodRatio< T >::operator()().

template<typename T>

double fit::HistoPoissonLikelihoodRatio< T >::rangeMin_

private

Definition at line 57 of file HistoPoissonLikelihoodRatio.h.

Referenced by fit::HistoPoissonLikelihoodRatio< T >::numberOfBins(), and fit::HistoPoissonLikelihoodRatio< T >::operator()().

template<typename T>

T* fit::HistoPoissonLikelihoodRatio< T >::t_

private

Definition at line 56 of file HistoPoissonLikelihoodRatio.h.

Referenced by fit::HistoPoissonLikelihoodRatio< T >::function().

template<typename T>

double fit::HistoPoissonLikelihoodRatio< T >::xMax_

private

Definition at line 59 of file HistoPoissonLikelihoodRatio.h.

Referenced by fit::HistoPoissonLikelihoodRatio< T >::HistoPoissonLikelihoodRatio(), and fit::HistoPoissonLikelihoodRatio< T >::setHistos().

template<typename T>

double fit::HistoPoissonLikelihoodRatio< T >::xMin_

private

Definition at line 59 of file HistoPoissonLikelihoodRatio.h.

Referenced by fit::HistoPoissonLikelihoodRatio< T >::HistoPoissonLikelihoodRatio(), fit::HistoPoissonLikelihoodRatio< T >::numberOfBins(), fit::HistoPoissonLikelihoodRatio< T >::operator()(), and fit::HistoPoissonLikelihoodRatio< T >::setHistos().