JCALIBRATE::JModel_t Struct Reference

Fit model. More...

#include <JFitK40.hh>

Inheritance diagram for JCALIBRATE::JModel_t:

Public Member Functions
size_t	getN () const
	Get number of fit parameters.
int	getIndex (JParameter_t JK40Parameters::*p) const
	Get index of parameter.
double	getValue (const double ct) const
	Get K40 coincidence rate as a function of cosine angle between PMT axes.
const JK40Parameters_t &	getGradient (const double ct) const
	Get gradient.
const JK40Parameters_t &	getK40Parameters () const
	Get K40 parameters.
void	setK40Parameters (const JK40Parameters_t &parameters)
	Set K40 parameters.
void	reset ()
	Reset.
void	print (std::ostream &out) const
	Print model parameters to output stream conform include files.

Static Public Member Functions
static const JK40Parameters &	getInstance ()
	Get default values.

Public Attributes
JPMTParameters_t	parameters [NUMBER_OF_PMTS]
JParameter_t	R
	maximal coincidence rate [Hz]
JParameter_t	p1
	1st order angle dependence coincidence rate
JParameter_t	p2
	2nd order angle dependence coincidence rate
JParameter_t	p3
	3rd order angle dependence coincidence rate
JParameter_t	p4
	4th order angle dependence coincidence rate
JParameter_t	cc
	fraction of signal correlated background
JParameter_t	bc
	constant background

Private Attributes
JK40Parameters_t	gradient

Friends
std::ostream &	operator<< (std::ostream &out, const JModel_t &object)
	Write model parameters to output stream.

Detailed Description

Fit model.

Definition at line 842 of file JFitK40.hh.

Member Function Documentation

getInstance()

static const JK40Parameters & JCALIBRATE::JK40Parameters::getInstance ( )

inlinestaticinherited

Get default values.

Values obtained with $JPP_DIR/examples/JCalibrate/JOMGsim.sh type B (see $JPP_DIR/examples/JCalibrate/README.md). If you change these values, you may also want to change the corresponding values in JK40DefaultSimulator.hh.

Returns

parameters

Definition at line 739 of file JFitK40.hh.

    {
      static JK40Parameters parameters;
 
      parameters.R .set(18.473257);
      parameters.p1.set( 3.030307);
      parameters.p2.set(-0.965429);
      parameters.p3.set( 1.290367);
      parameters.p4.set( 0.405618);
      parameters.cc.set( 0.0);
      parameters.bc.set( 0.0);
 
      return parameters;
    }

getN()

size_t JCALIBRATE::JK40Parameters::getN ( ) const

inlineinherited

Get number of fit parameters.

Returns

number of parameters

Definition at line 760 of file JFitK40.hh.

    {
      return ((R .isFree() ? 1 : 0) +
              (p1.isFree() ? 1 : 0) +
              (p2.isFree() ? 1 : 0) +
              (p3.isFree() ? 1 : 0) +
              (p4.isFree() ? 1 : 0) +
              (cc.isFree() ? 1 : 0) +
              (bc.isFree() ? 1 : 0));
    }

getIndex()

int JCALIBRATE::JK40Parameters::getIndex ( JParameter_t JK40Parameters::* p ) const

inlineinherited

Get index of parameter.

Parameters

p	pointer to data member

Returns

index

Definition at line 778 of file JFitK40.hh.

    {
      if (!(this->*p).isFree()) {
        return INVALID_INDEX;
      }
 
      int N = 0;
 
      if (p == &JK40Parameters::R)  { return N; } if (R .isFree()) { ++N; }
      if (p == &JK40Parameters::p1) { return N; } if (p1.isFree()) { ++N; }
      if (p == &JK40Parameters::p2) { return N; } if (p2.isFree()) { ++N; }
      if (p == &JK40Parameters::p3) { return N; } if (p3.isFree()) { ++N; }
      if (p == &JK40Parameters::p4) { return N; } if (p4.isFree()) { ++N; }
      if (p == &JK40Parameters::cc) { return N; } if (cc.isFree()) { ++N; }
      if (p == &JK40Parameters::bc) { return N; } if (bc.isFree()) { ++N; }
 
      return INVALID_INDEX;
    }

getValue()

double JCALIBRATE::JK40Parameters::getValue ( const double ct ) const

inlineinherited

Get K40 coincidence rate as a function of cosine angle between PMT axes.

Parameters

ct	cosine angle between PMT axes

Returns

rate [Hz]

Definition at line 804 of file JFitK40.hh.

    {
      return R * exp(ct*(p1+ct*(p2+ct*(p3+ct*p4))) - (p1+p2+p3+p4));
    }

getGradient()

const JK40Parameters_t & JCALIBRATE::JK40Parameters::getGradient ( const double ct ) const

inlineinherited

Get gradient.

Parameters

ct	cosine angle between PMT axes

Returns

gradient

Definition at line 816 of file JFitK40.hh.

    {
      gradient.reset();
 
      const double rate  = getValue(ct);
      const double ct2   = ct * ct;
 
      if (R .isFree()) { gradient.R  = rate / R; }
      if (p1.isFree()) { gradient.p1 = rate * ct         -  rate; }
      if (p2.isFree()) { gradient.p2 = rate * ct2        -  rate; }
      if (p3.isFree()) { gradient.p3 = rate * ct2 * ct   -  rate; }
      if (p4.isFree()) { gradient.p4 = rate * ct2 * ct2  -  rate; }
      if (cc.isFree()) { gradient.cc = rate; }
      if (bc.isFree()) { gradient.bc = 1.0; }
 
      return gradient;
    }

getK40Parameters()

const JK40Parameters_t & JCALIBRATE::JK40Parameters_t::getK40Parameters ( ) const

inlineinherited

Get K40 parameters.

Returns

K40 parameters

Definition at line 625 of file JFitK40.hh.

    {
      return static_cast<const JK40Parameters_t&>(*this);
    }

setK40Parameters()

void JCALIBRATE::JK40Parameters_t::setK40Parameters ( const JK40Parameters_t & parameters )

inlineinherited

Set K40 parameters.

Parameters

parameters

K40 parameters

Definition at line 636 of file JFitK40.hh.

    {
      static_cast<JK40Parameters_t&>(*this) = parameters;
    }

reset()

void JCALIBRATE::JK40Parameters_t::reset ( )

inlineinherited

Reset.

Definition at line 645 of file JFitK40.hh.

    {
      R .set(0.0);
      p1.set(0.0);
      p2.set(0.0);
      p3.set(0.0);
      p4.set(0.0);
      cc.set(0.0);
      bc.set(0.0);
    }

print()

void JCALIBRATE::JK40Parameters_t::print ( std::ostream & out ) const

inlineinherited

Print model parameters to output stream conform include files.

Parameters

out	output stream

Definition at line 662 of file JFitK40.hh.

    {
      using namespace std;
 
      out << "JFitK40.hh" << endl;
      out << "parameters.R .set(" << FIXED(9,6) << this->R () << ");" << endl;
      out << "parameters.p1.set(" << FIXED(9,6) << this->p1() << ");" << endl;
      out << "parameters.p2.set(" << FIXED(9,6) << this->p2() << ");" << endl;
      out << "parameters.p3.set(" << FIXED(9,6) << this->p3() << ");" << endl;
      out << "parameters.p4.set(" << FIXED(9,6) << this->p4() << ");" << endl;
      out << "cc                " << FIXED(9,6) << this->cc() << endl;
      out << "bc                " << FIXED(9,6) << this->bc() << endl;
      out << endl;
 
      out << "JK40DefaultSimulator.hh" << endl;
      out << "static constexpr double p1  =  " << FIXED(9,6) << this->p1() << ";" << endl;
      out << "static constexpr double p2  =  " << FIXED(9,6) << this->p2() << ";" << endl;
      out << "static constexpr double p3  =  " << FIXED(9,6) << this->p3() << ";" << endl;
      out << "static constexpr double p4  =  " << FIXED(9,6) << this->p4() << ";" << endl;
      out << endl;
    }

Friends And Related Symbol Documentation

operator<<

std::ostream & operator<< ( std::ostream & out, const JModel_t & object )

friend

Write model parameters to output stream.

Parameters

out	output stream
object	model parameters

Returns

output stream

Definition at line 855 of file JFitK40.hh.

    {
      using namespace std;
 
      out << static_cast<const JK40Parameters&>(object);
 
      for (int i = 0; i != NUMBER_OF_PMTS; ++i) {
        out << "PMT[" << FILL(2,'0') << i << FILL() << "]." << object.parameters[i].status << endl << object.parameters[i];
      }
 
      return out;
    }

Member Data Documentation

parameters

JPMTParameters_t JCALIBRATE::JModel_t::parameters[NUMBER_OF_PMTS]

Definition at line 845 of file JFitK40.hh.

gradient

JK40Parameters_t JCALIBRATE::JK40Parameters::gradient

mutableprivateinherited

Definition at line 835 of file JFitK40.hh.

R

JParameter_t JCALIBRATE::JK40Parameters_t::R

inherited

maximal coincidence rate [Hz]

Definition at line 708 of file JFitK40.hh.

p1

JParameter_t JCALIBRATE::JK40Parameters_t::p1

inherited

1st order angle dependence coincidence rate

Definition at line 709 of file JFitK40.hh.

p2

JParameter_t JCALIBRATE::JK40Parameters_t::p2

inherited

2nd order angle dependence coincidence rate

Definition at line 710 of file JFitK40.hh.

p3

JParameter_t JCALIBRATE::JK40Parameters_t::p3

inherited

3rd order angle dependence coincidence rate

Definition at line 711 of file JFitK40.hh.

p4

JParameter_t JCALIBRATE::JK40Parameters_t::p4

inherited

4th order angle dependence coincidence rate

Definition at line 712 of file JFitK40.hh.

cc

JParameter_t JCALIBRATE::JK40Parameters_t::cc

inherited

fraction of signal correlated background

Definition at line 713 of file JFitK40.hh.

bc

JParameter_t JCALIBRATE::JK40Parameters_t::bc

inherited

constant background

Definition at line 714 of file JFitK40.hh.

---

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

• JFitK40.hh