Jpp 21.0.0-rc.1
the software that should make you happy
Loading...
Searching...
No Matches
JCALIBRATE::JK40Parameters Struct Reference

Fit parameters for two-fold coincidence rate due to K40. More...

#include <JFitK40.hh>

Inheritance diagram for JCALIBRATE::JK40Parameters:
JCALIBRATE::JK40Parameters_t JCALIBRATE::JModel_t JCALIBRATE::JModel

Public Member Functions

 JK40Parameters ()
 Default constructor.
 
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_tgetGradient (const double ct) const
 Get gradient.
 
const JK40Parameters_tgetK40Parameters () 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 JK40ParametersgetInstance ()
 Get default values.
 

Public Attributes

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
 

Detailed Description

Fit parameters for two-fold coincidence rate due to K40.

Definition at line 721 of file JFitK40.hh.

Constructor & Destructor Documentation

◆ JK40Parameters()

JCALIBRATE::JK40Parameters::JK40Parameters ( )
inline

Default constructor.

Definition at line 727 of file JFitK40.hh.

728 {}

Member Function Documentation

◆ getInstance()

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

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.

740 {
741 static JK40Parameters parameters;
742
743 parameters.R .set(18.473257);
744 parameters.p1.set( 3.030307);
745 parameters.p2.set(-0.965429);
746 parameters.p3.set( 1.290367);
747 parameters.p4.set( 0.405618);
748 parameters.cc.set( 0.0);
749 parameters.bc.set( 0.0);
750
751 return parameters;
752 }
JK40Parameters()
Default constructor.
Definition JFitK40.hh:727

◆ getN()

size_t JCALIBRATE::JK40Parameters::getN ( ) const
inline

Get number of fit parameters.

Returns
number of parameters

Definition at line 760 of file JFitK40.hh.

761 {
762 return ((R .isFree() ? 1 : 0) +
763 (p1.isFree() ? 1 : 0) +
764 (p2.isFree() ? 1 : 0) +
765 (p3.isFree() ? 1 : 0) +
766 (p4.isFree() ? 1 : 0) +
767 (cc.isFree() ? 1 : 0) +
768 (bc.isFree() ? 1 : 0));
769 }
bool isFree() const
Check if parameter is free.
Definition JFitK40.hh:244
JParameter_t bc
constant background
Definition JFitK40.hh:714
JParameter_t R
maximal coincidence rate [Hz]
Definition JFitK40.hh:708
JParameter_t p1
1st order angle dependence coincidence rate
Definition JFitK40.hh:709
JParameter_t p2
2nd order angle dependence coincidence rate
Definition JFitK40.hh:710
JParameter_t p3
3rd order angle dependence coincidence rate
Definition JFitK40.hh:711
JParameter_t p4
4th order angle dependence coincidence rate
Definition JFitK40.hh:712
JParameter_t cc
fraction of signal correlated background
Definition JFitK40.hh:713

◆ getIndex()

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

Get index of parameter.

Parameters
ppointer to data member
Returns
index

Definition at line 778 of file JFitK40.hh.

779 {
780 if (!(this->*p).isFree()) {
781 return INVALID_INDEX;
782 }
783
784 int N = 0;
785
786 if (p == &JK40Parameters::R) { return N; } if (R .isFree()) { ++N; }
787 if (p == &JK40Parameters::p1) { return N; } if (p1.isFree()) { ++N; }
788 if (p == &JK40Parameters::p2) { return N; } if (p2.isFree()) { ++N; }
789 if (p == &JK40Parameters::p3) { return N; } if (p3.isFree()) { ++N; }
790 if (p == &JK40Parameters::p4) { return N; } if (p4.isFree()) { ++N; }
791 if (p == &JK40Parameters::cc) { return N; } if (cc.isFree()) { ++N; }
792 if (p == &JK40Parameters::bc) { return N; } if (bc.isFree()) { ++N; }
793
794 return INVALID_INDEX;
795 }
static const int INVALID_INDEX
invalid index
Definition JFitK40.hh:60

◆ getValue()

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

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

Parameters
ctcosine angle between PMT axes
Returns
rate [Hz]

Definition at line 804 of file JFitK40.hh.

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

◆ getGradient()

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

Get gradient.

Parameters
ctcosine angle between PMT axes
Returns
gradient

Definition at line 816 of file JFitK40.hh.

817 {
818 gradient.reset();
819
820 const double rate = getValue(ct);
821 const double ct2 = ct * ct;
822
823 if (R .isFree()) { gradient.R = rate / R; }
824 if (p1.isFree()) { gradient.p1 = rate * ct - rate; }
825 if (p2.isFree()) { gradient.p2 = rate * ct2 - rate; }
826 if (p3.isFree()) { gradient.p3 = rate * ct2 * ct - rate; }
827 if (p4.isFree()) { gradient.p4 = rate * ct2 * ct2 - rate; }
828 if (cc.isFree()) { gradient.cc = rate; }
829 if (bc.isFree()) { gradient.bc = 1.0; }
830
831 return gradient;
832 }
JK40Parameters_t gradient
Definition JFitK40.hh:835
double getValue(const double ct) const
Get K40 coincidence rate as a function of cosine angle between PMT axes.
Definition JFitK40.hh:804

◆ getK40Parameters()

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

Get K40 parameters.

Returns
K40 parameters

Definition at line 625 of file JFitK40.hh.

626 {
627 return static_cast<const JK40Parameters_t&>(*this);
628 }
JK40Parameters_t()
Default constructor.
Definition JFitK40.hh:614

◆ setK40Parameters()

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

Set K40 parameters.

Parameters
parametersK40 parameters

Definition at line 636 of file JFitK40.hh.

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

◆ reset()

void JCALIBRATE::JK40Parameters_t::reset ( )
inlineinherited

Reset.

Definition at line 645 of file JFitK40.hh.

646 {
647 R .set(0.0);
648 p1.set(0.0);
649 p2.set(0.0);
650 p3.set(0.0);
651 p4.set(0.0);
652 cc.set(0.0);
653 bc.set(0.0);
654 }
void set()
Set current value.
Definition JFitK40.hh:275

◆ print()

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

Print model parameters to output stream conform include files.

Parameters
outoutput stream

Definition at line 662 of file JFitK40.hh.

663 {
664 using namespace std;
665
666 out << "JFitK40.hh" << endl;
667 out << "parameters.R .set(" << FIXED(9,6) << this->R () << ");" << endl;
668 out << "parameters.p1.set(" << FIXED(9,6) << this->p1() << ");" << endl;
669 out << "parameters.p2.set(" << FIXED(9,6) << this->p2() << ");" << endl;
670 out << "parameters.p3.set(" << FIXED(9,6) << this->p3() << ");" << endl;
671 out << "parameters.p4.set(" << FIXED(9,6) << this->p4() << ");" << endl;
672 out << "cc " << FIXED(9,6) << this->cc() << endl;
673 out << "bc " << FIXED(9,6) << this->bc() << endl;
674 out << endl;
675
676 out << "JK40DefaultSimulator.hh" << endl;
677 out << "static constexpr double p1 = " << FIXED(9,6) << this->p1() << ";" << endl;
678 out << "static constexpr double p2 = " << FIXED(9,6) << this->p2() << ";" << endl;
679 out << "static constexpr double p3 = " << FIXED(9,6) << this->p3() << ";" << endl;
680 out << "static constexpr double p4 = " << FIXED(9,6) << this->p4() << ";" << endl;
681 out << endl;
682 }
Auxiliary data structure for floating point format specification.
Definition JManip.hh:448

Member Data Documentation

◆ gradient

JK40Parameters_t JCALIBRATE::JK40Parameters::gradient
mutableprivate

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: