Parametrisation of time-over-threshold distribution. More...

#include <JFitToT.hh>

Inheritance diagram for JCALIBRATE::JFitToT:

Public Member Functions
	JFitToT (const JPMTParameters &parameters, const JRange_t &range)
	Constructor. More...

void	fitInit (TH1 &h1, const double gradientThreshold=-0.005)
	Fit initializer. More...

TFitResultPtr	operator() (TH1 &h1, const double minWeight, const std::string &option)
	Fit histogram. More...

Double_t	getValue (Double_t x, Double_t par)
	Get rate as a function of the fit parameters. More...

const double	getInitGain () const
	Retrieve initial gain. More...

const double	getInitGainSpread () const
	Retrieve initial gainspread. More...

const Double_t *	getModelParameters () const
	Get model parameters. More...

Double_t *	getModelParameters ()
	Get model parameters. More...

void	setModelParameter (const int i, const Double_t value)
	Set model parameter. More...

void	setModelParameters (const Double_t *data)
	Set model parameters. More...

const Double_t	getModelParameter (const int i) const
	Get model parameter. More...

JFitParameter_t	getModelParameter (Double_t JFitToTParameters::*p) const
	Get model parameter. More...

Static Public Member Functions
static Int_t	getNumberOfModelParameters ()
	Get number of model parameters. More...

Public Attributes
Double_t	gain
	PMT gain. More...

Double_t	gainSpread
	PMT gain spread. More...

Static Private Member Functions
static JPMTAnalogueSignalProcessor &	getInstance ()
	Get unique instance of PMT analogue signal processor. More...

Private Attributes
JRange_t	fitRange
	fit-range [ns] More...

double	WfitRange
	Cumulative weight within fit-range. More...

double	init_gain
	Initial input value for gain. More...

double	init_gainSpread
	Initial input value for gainSpread. More...

const int	NPE
	True number of photo-electrons. More...

Detailed Description

Parametrisation of time-over-threshold distribution.

Note that for use in ROOT fit operations, the member method JFitToT::getValue is static.

Definition at line 173 of file JFitToT.hh.

Constructor & Destructor Documentation

JCALIBRATE::JFitToT::JFitToT	(	const JPMTParameters &	parameters,
		const JRange_t &	range
	)

inline

Constructor.

Parameters

parameters	parameters
range	abscissa range

Definition at line 184 of file JFitToT.hh.

                                                                      :
 
       JFitToTParameters(parameters),
 
       TF1              (FITTOT_FNAME.c_str(),
                         this,
                         &JFitToT::getValue,
                         range.getLowerLimit(),
                         range.getUpperLimit(),
                         JFitToT::getNumberOfModelParameters()),
 
       fitRange         (range),
       WfitRange        (1.0),
       
       NPE              (1)
     {
       getInstance().setPMTParameters(parameters);
 
       init_gain       = getInstance().gain;
       init_gainSpread = getInstance().gainSpread; 
     }

Member Function Documentation

void JCALIBRATE::JFitToT::fitInit	(	TH1 &	h1,
		const double	gradientThreshold = `-0.005`
	)

inline

Fit initializer.

Parameters

h1	ROOT 1D-histogram
gradientThreshold	treshold gradient for search of maximum

Definition at line 213 of file JFitToT.hh.

     {
       using namespace std;
       using namespace JPP;
       
       // Find gain and gainspread initial values
       const Int_t N = h1.GetXaxis()->GetNbins();
 
       double max    = h1.GetBinContent(N);
       double mode   = h1.GetBinCenter(N);
 
       for (int i = h1.GetBinCenter(N-1);
                i > h1.GetBinCenter(h1.FindBin(getInstance().mean_ns));
              --i) {
 
         const double x = h1.GetBinCenter (i);
         const double y = h1.GetBinContent(i);
 
         const double gradient = ( (h1.GetBinContent(i-1) - h1.GetBinContent(i+1)) /
                                   (h1.GetBinCenter (i+1) - h1.GetBinCenter (i-1)) );
 
         if (y > max) {
 
           mode = x;
           max  = y;
           
         } else if (gradient < gradientThreshold) {
 
           break;
         }
       }
 
       init_gain         = getInstance().getNPE(mode);
       init_gainSpread   = ( init_gain / 3.0 > FITTOT_GAINSPREAD_MIN ?
                             init_gain / 3.0 : FITTOT_GAINSPREAD_MIN ); // Range rule of thumb
 
       
       // Set parameter names and initial values
       const Int_t Ngain        = this->getModelParameter(&JFitToT::gain);
       const Int_t NgainSpread  = this->getModelParameter(&JFitToT::gainSpread);
       
       this->SetParName(Ngain,       MAKE_CSTRING("gain"));
       this->SetParName(NgainSpread, MAKE_CSTRING("gainSpread"));
       
       getInstance().gain       = init_gain;
       getInstance().gainSpread = init_gainSpread;
       
       setModelParameter(Ngain,       init_gain);
       setModelParameter(NgainSpread, init_gainSpread);
       
       this->SetParameters(this->getModelParameters());
 
 
       // Set parameter ranges
       setParLimits(*this, Ngain,       FITTOT_GAIN_MIN,       FITTOT_GAIN_MAX);
       setParLimits(*this, NgainSpread, FITTOT_GAINSPREAD_MIN, FITTOT_GAINSPREAD_MAX);
 
 
       // Set fit range and compute corresponding histogram weight
       if (!fitRange.is_valid()) {
         fitRange.setRange(mode - FITTOT_TOT_LMARGIN,
                           mode + FITTOT_TOT_RMARGIN);
       }
 
       fitRange.setRange(std::max(h1.GetXaxis()->GetXmin(), fitRange.getLowerLimit()),
                         std::min(h1.GetXaxis()->GetXmax(), fitRange.getUpperLimit()));
 
       WfitRange = h1.Integral(h1.FindBin(fitRange.getLowerLimit()),
                               h1.FindBin(fitRange.getUpperLimit()));
       
       this->SetRange(fitRange.getLowerLimit(), fitRange.getUpperLimit());
     }

TFitResultPtr JCALIBRATE::JFitToT::operator()	(	TH1 &	h1,
		const double	minWeight,
		const std::string &	option
	)

inline

Fit histogram.

Note that the PMT parameters which are part of the model are reset before the fit according the status of each PMT and the obtained fit parameters are copied back to the model parameters after the fit.

Parameters

h1	ROOT 1D-histogram
minWeight	minimum histogram weight in fit-range
option	fit option

Definition at line 297 of file JFitToT.hh.

     {
       using namespace std;
       using namespace JPP;
 
       fitInit(h1);
 
       if (WfitRange >  minWeight &&
           init_gain >= FITTOT_GAIN_MIN && init_gain <= FITTOT_GAIN_MAX) {
 
         const TFitResultPtr result = h1.Fit(this, option.c_str());
         
         this->setModelParameters(this->GetParameters());
 
         const double ToTcheck = getInstance().getTimeOverThreshold(gain);
 
         return (fitRange(ToTcheck) ? result : TFitResultPtr());
         
       } else {
 
         return TFitResultPtr();
       }
     }

Double_t JCALIBRATE::JFitToT::getValue	(	Double_t *	x,
		Double_t *	par
	)

inline

Get rate as a function of the fit parameters.

Parameters

x	pointer to abscissa values
par	pointer to parameter values

Returns: rate distribution at specified time-over-threshold [Hz/ns]

Definition at line 329 of file JFitToT.hh.

     {
       const Double_t tot_ns = x[0];
 
       getInstance().gain        = par[0];
       getInstance().gainSpread  = par[1];
 
       // Compute total weight (= histogram weight within fit range
       //                         plus pdf weight outside fit range)
       const double Wtot = (WfitRange +
                            1.0 - getInstance().getIntegralOfTimeOverThresholdProbability(fitRange.getLowerLimit(),
                                                                                          fitRange.getUpperLimit(),
                                                                                          NPE));
       
       return Wtot * getInstance().getTimeOverThresholdProbability(tot_ns, NPE);
     }

const double JCALIBRATE::JFitToT::getInitGain ( ) const

inline

Retrieve initial gain.

Returns: initial gain setting

Definition at line 352 of file JFitToT.hh.

                                      {
       return init_gain;
     }

const double JCALIBRATE::JFitToT::getInitGainSpread ( ) const

inline

Retrieve initial gainspread.

Returns: initial gainspread setting

Definition at line 362 of file JFitToT.hh.

                                            {
       return init_gainSpread;
     }

static JPMTAnalogueSignalProcessor& JCALIBRATE::JFitToT::getInstance ( )

inlinestaticprivate

Get unique instance of PMT analogue signal processor.

Returns: reference to PMT analogue signal processor.

Definition at line 374 of file JFitToT.hh.

     {
       static JPMTAnalogueSignalProcessor object;
 
       return object;
     }

static Int_t JCALIBRATE::JFitToTParameters::getNumberOfModelParameters ( )

inlinestaticinherited

Get number of model parameters.

Returns: number of parameters

Definition at line 82 of file JFitToT.hh.

     {
       return sizeof(JFitToTParameters) / sizeof(Double_t);
     }

const Double_t* JCALIBRATE::JFitToTParameters::getModelParameters ( ) const

inlineinherited

Get model parameters.

Returns: pointer to parameters

Definition at line 93 of file JFitToT.hh.

     {
       return &this->gain; 
     }

Double_t* JCALIBRATE::JFitToTParameters::getModelParameters ( )

inlineinherited

Get model parameters.

Returns: pointer to parameters

Definition at line 104 of file JFitToT.hh.

     {
       return &this->gain; 
     }

void JCALIBRATE::JFitToTParameters::setModelParameter	(	const int	i,
		const Double_t	value
	)

inlineinherited

Set model parameter.

Parameters

i	parameter index
value	parameter value

Definition at line 116 of file JFitToT.hh.

     {
       getModelParameters()[i] = value;
     }

void JCALIBRATE::JFitToTParameters::setModelParameters ( const Double_t * data )

inlineinherited

Set model parameters.

Parameters

data	pointer to parameters

Definition at line 127 of file JFitToT.hh.

     {
       for (Int_t i = 0; i != getNumberOfModelParameters(); ++i) {
         setModelParameter(i, data[i]);
       }
     }

const Double_t JCALIBRATE::JFitToTParameters::getModelParameter ( const int i ) const

inlineinherited

Get model parameter.

Parameters

i	parameter index

Returns: parameter value

Definition at line 141 of file JFitToT.hh.

     {
       return getModelParameters()[i];
     }

JFitParameter_t JCALIBRATE::JFitToTParameters::getModelParameter ( Double_t JFitToTParameters::* p ) const

inlineinherited

Get model parameter.

Parameters

p	pointer to data member

Returns: parameter index and value

Definition at line 153 of file JFitToT.hh.

     {
       const Int_t i = &(this->*p) - getModelParameters();
 
       return JFitParameter_t(i, getModelParameter(i));
     }

Member Data Documentation

JRange_t JCALIBRATE::JFitToT::fitRange

private

fit-range [ns]

Definition at line 381 of file JFitToT.hh.

double JCALIBRATE::JFitToT::WfitRange

private

Cumulative weight within fit-range.

Definition at line 382 of file JFitToT.hh.

double JCALIBRATE::JFitToT::init_gain

private

Initial input value for gain.

Definition at line 384 of file JFitToT.hh.

double JCALIBRATE::JFitToT::init_gainSpread

private

Initial input value for gainSpread.

Definition at line 385 of file JFitToT.hh.

const int JCALIBRATE::JFitToT::NPE

private

True number of photo-electrons.

Definition at line 387 of file JFitToT.hh.

Double_t JCALIBRATE::JFitToTParameters::gain

inherited

PMT gain.

Definition at line 163 of file JFitToT.hh.

Double_t JCALIBRATE::JFitToTParameters::gainSpread

inherited

PMT gain spread.

Definition at line 164 of file JFitToT.hh.

The documentation for this struct was generated from the following file:

JFitToT.hh

Public Member Functions

Static Public Member Functions

Public Attributes

Static Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation