JCALIBRATE::JFitToT Struct Reference

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

#include <JFitToT.hh>

Inheritance diagram for JCALIBRATE::JFitToT:

Public Member Functions
	JFitToT (const JPMTParameters &parameters, const JRange< double > &range)
	Constructor. More...
const JPMTAnalogueSignalProcessor &	getCPU () const
	Access method for the analogue signal processor. More...
TFitResultPtr	operator() (TH1 &h1, const std::string &option)
	Fit histogram. More...
Double_t	getValue (const Double_t *x, const Double_t *par)
	Get rate as a function of the fit parameters. 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...
Double_t	normalization

Private Attributes
JPMTAnalogueSignalProcessor	cpu

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 179 of file JFitToT.hh.

Constructor & Destructor Documentation

JCALIBRATE::JFitToT::JFitToT ( const JPMTParameters &  parameters, const JRange< double > &  range  )

inline

Constructor.

Parameters

parameters	parameters
range	abscissa range

Definition at line 190 of file JFitToT.hh.

                                                                           :

      TF1(FITTOT_FNAME.c_str(),
          this,
          &JFitToT::getValue,
          range.getLowerLimit(),
          range.getUpperLimit(),
          JFitToT::getNumberOfModelParameters()),

      JFitToTParameters(parameters),
      
      cpu(parameters)
    {}

Member Function Documentation

const JPMTAnalogueSignalProcessor& JCALIBRATE::JFitToT::getCPU ( ) const

inline

Access method for the analogue signal processor.

Returns

reference to analogue signal processor

Definition at line 210 of file JFitToT.hh.

    {
      return cpu;
    }

TFitResultPtr JCALIBRATE::JFitToT::operator() ( TH1 &  h1, 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
option	fit option

Returns

fit result

Definition at line 226 of file JFitToT.hh.

    {
      using namespace std;
      using namespace JPP;

      // Set initial gain and gain-spread
      const Int_t Ngain          = this->getModelParameter(&JFitToT::JFitToTParameters::gain);
      const Int_t NgainSpread    = this->getModelParameter(&JFitToT::JFitToTParameters::gainSpread);
      const Int_t Nnormalization = this->getModelParameter(&JFitToT::JFitToTParameters::normalization);      
      
      SetParName(Ngain,          MAKE_CSTRING(FITTOT_GAIN_PARNAME));
      SetParName(NgainSpread,    MAKE_CSTRING(FITTOT_GAINSPREAD_PARNAME));
      SetParName(Nnormalization, MAKE_CSTRING(FITTOT_NORMALIZATION_PARNAME));

      normalization = h1.Integral(h1.FindBin(this->GetXmin()),
                                  h1.FindBin(this->GetXmax()));

      SetParameters(getModelParameters());

      // Set parameter and fit ranges
      setParLimits(*this, Ngain,       FITTOT_GAIN_MIN,       FITTOT_GAIN_MAX);
      setParLimits(*this, NgainSpread, FITTOT_GAINSPREAD_MIN, FITTOT_GAINSPREAD_MAX);
      
      FixParameter(Nnormalization, normalization);

      // Fit histogram
      const TFitResultPtr result = h1.Fit(this, option.c_str());
        
      this->setModelParameters(this->GetParameters());
      
      cpu.gain       = getModelParameter(Ngain);
      cpu.gainSpread = getModelParameter(NgainSpread);
        
      return result;
    }

Double_t JCALIBRATE::JFitToT::getValue ( const Double_t *  x, const 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 270 of file JFitToT.hh.

    {
      const double tot_ns = x[0];

      // Set new parameter values
      cpu.gain            = par[0]; 
      cpu.gainSpread      = par[1];

      // Determine normalization factor
      const int    NPE    = 1;
      const double Whist  = par[2];
      const double Wpdf   = cpu.getIntegralOfTimeOverThresholdProbability(GetXmin(), GetXmax(), NPE);
      
      return (Wpdf > 0.0 ? (Whist / Wpdf) : Whist) * cpu.getTimeOverThresholdProbability(tot_ns, NPE);
    }

static Int_t JCALIBRATE::JFitToTParameters::getNumberOfModelParameters ( )

inlinestaticinherited

Get number of model parameters.

Returns

number of parameters

Definition at line 87 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 98 of file JFitToT.hh.

    {
      return &this->gain; 
    }

Double_t* JCALIBRATE::JFitToTParameters::getModelParameters ( )

inlineinherited

Get model parameters.

Returns

pointer to parameters

Definition at line 109 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 121 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 132 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 146 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 158 of file JFitToT.hh.

    {
      const Int_t i = &(this->*p) - getModelParameters();

      return JFitParameter_t(i, getModelParameter(i));
    }

Member Data Documentation

JPMTAnalogueSignalProcessor JCALIBRATE::JFitToT::cpu

private

Definition at line 288 of file JFitToT.hh.

Double_t JCALIBRATE::JFitToTParameters::gain

inherited

PMT gain.

Definition at line 168 of file JFitToT.hh.

Double_t JCALIBRATE::JFitToTParameters::gainSpread

inherited

PMT gain spread.

Definition at line 169 of file JFitToT.hh.

Double_t JCALIBRATE::JFitToTParameters::normalization

inherited

Definition at line 170 of file JFitToT.hh.

---

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

• JFitToT.hh