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 double xmin, const double xmax)
	Constructor. More...
TFitResultPtr	operator() (TH1 &h1, const std::string &option)
	Fit histogram. More...
const Double_t *	getModelParameters () const
	Get model parameters. More...
Double_t *	getModelParameters ()
	Get model parameters. 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 Double_t	getValue (const Double_t *x, const Double_t *par)
	Get rate as a function of the fit parameters. More...
static Int_t	getNumberOfModelParameters ()
	Get number of model parameters. More...

Public Attributes
Double_t	W
	normalisation More...
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
double	Xmin
double	Xmax

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

Constructor & Destructor Documentation

JCALIBRATE::JFitToT::JFitToT ( const JPMTParameters &  parameters, const double  xmin, const double  xmax  )

inline

Constructor.

Parameters

parameters	parameters
xmin	minimal x
xmax	maximal x

Definition at line 164 of file JFitToT.hh.

                                        :

      JFitToTParameters(parameters),

      TF1("fittot",
          JFitToT::getValue,
          xmin, xmax,
          JFitToT::getNumberOfModelParameters()),

      Xmin(xmin),
      Xmax(xmax)
    {
      getInstance().setPMTParameters(parameters);
    }

Member Function Documentation

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

Definition at line 191 of file JFitToT.hh.

    {
      using namespace std;
      using namespace JPP;

      Double_t x0 = numeric_limits<Double_t>::max();
      Double_t y0 = 0.0;
      this->W     = 0.0;
      
      for (Int_t i = 1; i <= h1.GetXaxis()->GetNbins(); ++i) {

        const Double_t x = h1.GetBinCenter (i);
        const Double_t y = h1.GetBinContent(i);

        if (x >= this->Xmin && 
            x <= this->Xmax) {
          this->W += y;
        }

        if (x >= 0.5*TIME_OVER_THRESHOLD_NS && 
            x <= 2.0*TIME_OVER_THRESHOLD_NS) {

          if (y >= y0) {
            x0 = x;
            y0 = y;
          }
        }
      }

      double xmin = Xmin + x0 - TIME_OVER_THRESHOLD_NS;
      double xmax = Xmax + x0 - TIME_OVER_THRESHOLD_NS;
      
      if (xmin < FITTOT_TOT_MIN_NS) {
        xmin = FITTOT_TOT_MIN_NS;
      }

      this->SetRange(xmin, xmax);

      this->gain = getInstance().getNPE(x0);

      setParameter(*this, this->getModelParameter(&JFitToT::gain));
      setParameter(*this, this->getModelParameter(&JFitToT::gainSpread));

      setParLimits(*this, this->getModelParameter(&JFitToT::gain),       FITTOT_GAIN_MIN,       FITTOT_GAIN_MAX);
      setParLimits(*this, this->getModelParameter(&JFitToT::gainSpread), FITTOT_GAINSPREAD_MIN, FITTOT_GAINSPREAD_MAX);

      this->SetParameters(this->getModelParameters());

      getInstance().gain       = this->gain;
      getInstance().gainSpread = this->gainSpread;

      const TFitResultPtr result = h1.Fit(this, option.c_str());

      this->setModelParameters(this->GetParameters());

      return result;
    }

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

inlinestatic

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

    {
      static const int NPE = 1;
      
      const Double_t W         = par[0];
      getInstance().gain       = par[1];
      getInstance().gainSpread = par[2];
      
      const Double_t tot_ns = x[0];
      const Double_t npe    = getInstance().getNPE(tot_ns);
      
      return W * getInstance().getProbability(npe, NPE) * getInstance().getDerivative(npe);
    }

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

inlinestaticprivate

Get unique instance of PMT analogue signal processor.

Returns

reference to PMT analogue signal processor.

Definition at line 279 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 73 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 84 of file JFitToT.hh.

    {
      return &this->W;
    }

Double_t* JCALIBRATE::JFitToTParameters::getModelParameters ( )

inlineinherited

Get model parameters.

Returns

pointer to parameters

Definition at line 95 of file JFitToT.hh.

    {
      return &this->W;
    }

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

inlineinherited

Set model parameters.

Parameters

data	pointer to parameters

Definition at line 106 of file JFitToT.hh.

    {
      for (Int_t i = 0; i != getNumberOfModelParameters(); ++i) {
        getModelParameters()[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 120 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 132 of file JFitToT.hh.

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

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

Member Data Documentation

double JCALIBRATE::JFitToT::Xmin

private

Definition at line 286 of file JFitToT.hh.

double JCALIBRATE::JFitToT::Xmax

private

Definition at line 287 of file JFitToT.hh.

Double_t JCALIBRATE::JFitToTParameters::W

inherited

normalisation

Definition at line 142 of file JFitToT.hh.

Double_t JCALIBRATE::JFitToTParameters::gain

inherited

PMT gain.

Definition at line 143 of file JFitToT.hh.

Double_t JCALIBRATE::JFitToTParameters::gainSpread

inherited

PMT gain spread.

Definition at line 144 of file JFitToT.hh.

---

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

• JFitToT.hh