JCALIBRATE::JFitK40_t< TF2 > Struct Template Reference

Template specialisation of two-fold coincidence rate due to K40 and other radioactive decays. More...

#include <JFitK40.hh>

Inheritance diagram for JCALIBRATE::JFitK40_t< TF2 >:

Public Member Functions
	JFitK40_t (const JModule &module, const Double_t xmin, const Double_t xmax, const Double_t ymin, const Double_t ymax, const bool option)
	Constructor. More...
TFitResultPtr	operator() (TH2 &h2, const std::string &option)
	Fit 2D-histogram. More...

Static Public Member Functions
static Double_t	getRate (const Double_t *x, const Double_t *data)
	Get K40 coincidence rate as a function of the fit parameters. More...

Detailed Description

template<>
struct JCALIBRATE::JFitK40_t< TF2 >

Template specialisation of two-fold coincidence rate due to K40 and other radioactive decays.

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

Definition at line 586 of file JFitK40.hh.

Constructor & Destructor Documentation

JCALIBRATE::JFitK40_t< TF2 >::JFitK40_t ( const JModule &  module, const Double_t  xmin, const Double_t  xmax, const Double_t  ymin, const Double_t  ymax, const bool  option  )

inline

Constructor.

Parameters

module	detector module
xmin	minimal x
xmax	maximal x
ymin	minimal y
ymax	maximal y
option	true = fix average t0; false = free average t0

Definition at line 600 of file JFitK40.hh.

                                     :
      JFitK40_t<>(module, option),
      TF2("f2",
          JFitK40_t<TF2>::getRate,
          xmin, xmax, ymin, ymax,
          getNumberOfModelParameters())
    {}

Member Function Documentation

static Double_t JCALIBRATE::JFitK40_t< TF2 >::getRate ( const Double_t *  x, const Double_t *  data  )

inlinestatic

Get K40 coincidence rate as a function of the fit parameters.

To speed up the calculation, it is assumed that the parameter values do not change unless also the x[0] value changes (i.e. the index of the PMT pair).

Parameters

x	pointer to data
data	pointer to parameter values

Returns

rate [Hz/ns]

Definition at line 624 of file JFitK40.hh.

    {
      static int       id    =  -1;
      static Double_t  t0;
      static Double_t  sigma;
      static Double_t  rate;
      
      const int      ix    = (int) x[0];
      const Double_t dt_ns = x[1];

      if (ix != id) {
        
        getInstance().setModelParameters(data);
        
        const pair_type& pair = getInstance().getPair(ix);

        t0     =  getInstance().getTimeOffset(pair);
        sigma  =  getInstance().getSigma     (pair);
        rate   =  getInstance().getValue     (pair);
        id     =  ix;
      }

      return getInstance().bg  +  rate * (getInstance().cc + TMath::Gaus(dt_ns, t0, sigma, kTRUE));
    }

TFitResultPtr JCALIBRATE::JFitK40_t< TF2 >::operator() ( TH2 &  h2, const std::string &  option  )

inline

Fit 2D-histogram.

Note that the PMT parameters are partially 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

h2	ROOT histogram
option	fit option

Definition at line 659 of file JFitK40.hh.

    {
      for (int pmt = 0; pmt != NUMBER_OF_PMTS; ++pmt) {

        parameters[pmt].QE  = JPMTParameters_t().QE;
        parameters[pmt].TTS = JPMTParameters_t().TTS;

        if (disable[pmt]) {

          parameters[pmt].QE = 0.0;

          fixParameter(*this, this->getModelParameter(pmt, &JPMTParameters_t::QE));
          fixParameter(*this, this->getModelParameter(pmt, &JPMTParameters_t::TTS));
          fixParameter(*this, this->getModelParameter(pmt, &JPMTParameters_t::t0));

        } else {

          // Note that setting limits to a parameter will also release it

          if (!isParameterFixed(*this, this->getModelParameter(pmt, &JPMTParameters_t::QE)))  {
            setParLimits(*this, this->getModelParameter(pmt, &JPMTParameters_t::QE),  FITK40_QE_MIN,            FITK40_QE_MAX);
          }
          if (!isParameterFixed(*this, this->getModelParameter(pmt, &JPMTParameters_t::TTS))) {
            setParLimits(*this, this->getModelParameter(pmt, &JPMTParameters_t::TTS), FITK40_TTS_MIN_NS,        FITK40_TTS_MAX_NS);
          }
          if (!isParameterFixed(*this, this->getModelParameter(pmt, &JPMTParameters_t::t0)))  {
            setParLimits(*this, this->getModelParameter(pmt, &JPMTParameters_t::t0),  h2.GetYaxis()->GetXmin(), h2.GetYaxis()->GetXmax());
          }
        }
      }

      if (index_of_average_t0 != -1) {

        this->setT0(index_of_average_t0, 0.0);

        fixParameter(*this, this->getModelParameter(index_of_average_t0, &JPMTParameters_t::t0));
      }

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

      getInstance() = *this;

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

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

      return result;
    }

---

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

• JFitK40.hh