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.

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

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 794 of file JFitK40.hh.

Member Function Documentation

getInstance()

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

inlinestaticinherited

Get default values.

The default parameter values are set to those obtained from a designated simulation of K40 decays (see http://wiki.km3net.de/index.php/OMGsim_simulations_for_K40_fit).
If you change these values, you may also want to change the corresponding values in JK40DefaultSimulator.hh.

Returns

parameters

Definition at line 695 of file JFitK40.hh.

    {
      static JK40Parameters parameters;
 
      parameters.R .set(18.460546);
      parameters.p1.set( 3.0767);
      parameters.p2.set(-1.2078);
      parameters.p3.set( 0.9905);
      parameters.p4.set( 0.9379);
      parameters.cc.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 715 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));
    }

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 732 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; }
 
      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 757 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 769 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; }
 
      return gradient;
    }

getK40Parameters()

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

inlineinherited

Get K40 parameters.

Returns

K40 parameters

Definition at line 633 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 644 of file JFitK40.hh.

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

reset()

void JCALIBRATE::JK40Parameters_t::reset ( )

inlineinherited

Reset.

Definition at line 653 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);
    }

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 807 of file JFitK40.hh.

    {
      using namespace std;
 
      out << "Rate [Hz] " << FIXED(12,6) << object.R  << endl;
      out << "p1        " << FIXED(12,6) << object.p1 << endl;
      out << "p2        " << FIXED(12,6) << object.p2 << endl;
      out << "p3        " << FIXED(12,6) << object.p3 << endl;
      out << "p4        " << FIXED(12,6) << object.p4 << endl;
      out << "cc        " << FIXED(12,6) << object.cc << endl;
 
      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 797 of file JFitK40.hh.

gradient

JK40Parameters_t JCALIBRATE::JK40Parameters::gradient

mutableprivateinherited

Definition at line 787 of file JFitK40.hh.

R

JParameter_t JCALIBRATE::JK40Parameters_t::R

inherited

maximal coincidence rate [Hz]

Definition at line 664 of file JFitK40.hh.

p1

JParameter_t JCALIBRATE::JK40Parameters_t::p1

inherited

1st order angle dependence coincidence rate

Definition at line 665 of file JFitK40.hh.

p2

JParameter_t JCALIBRATE::JK40Parameters_t::p2

inherited

2nd order angle dependence coincidence rate

Definition at line 666 of file JFitK40.hh.

p3

JParameter_t JCALIBRATE::JK40Parameters_t::p3

inherited

3rd order angle dependence coincidence rate

Definition at line 667 of file JFitK40.hh.

p4

JParameter_t JCALIBRATE::JK40Parameters_t::p4

inherited

4th order angle dependence coincidence rate

Definition at line 668 of file JFitK40.hh.

cc

JParameter_t JCALIBRATE::JK40Parameters_t::cc

inherited

fraction of signal correlated background

Definition at line 669 of file JFitK40.hh.

---

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

• JFitK40.hh