Jpp 21.0.0-rc.1-88-g0130508c4
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

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

Constructor & Destructor Documentation

◆ JK40Parameters()

JCALIBRATE::JK40Parameters::JK40Parameters ( )
inline

Default constructor.

Definition at line 730 of file JFitK40.hh.

731 {}

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

743 {
744 static JK40Parameters parameters;
745 parameters.R .set(18.430675);
746 parameters.p1.set( 2.919895);
747 parameters.p2.set(-0.831970);
748 parameters.p3.set( 1.407887);
749 parameters.p4.set( 0.170510);
750 parameters.cc.set( 0.0);
751 parameters.bc.set( 0.0);
752
753 return parameters;
754 }
JK40Parameters()
Default constructor.
Definition JFitK40.hh:730

◆ getN()

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

Get number of fit parameters.

Returns
number of parameters

Definition at line 762 of file JFitK40.hh.

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

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

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

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

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

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

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

◆ getK40Parameters()

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

Get K40 parameters.

Returns
K40 parameters

Definition at line 628 of file JFitK40.hh.

629 {
630 return static_cast<const JK40Parameters_t&>(*this);
631 }
JK40Parameters_t()
Default constructor.
Definition JFitK40.hh:617

◆ setK40Parameters()

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

Set K40 parameters.

Parameters
parametersK40 parameters

Definition at line 639 of file JFitK40.hh.

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

◆ reset()

void JCALIBRATE::JK40Parameters_t::reset ( )
inlineinherited

Reset.

Definition at line 648 of file JFitK40.hh.

649 {
650 R .set(0.0);
651 p1.set(0.0);
652 p2.set(0.0);
653 p3.set(0.0);
654 p4.set(0.0);
655 cc.set(0.0);
656 bc.set(0.0);
657 }
void set()
Set current value.
Definition JFitK40.hh:278

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

666 {
667 using namespace std;
668
669 out << "JFitK40.hh" << endl;
670 out << "parameters.R .set(" << FIXED(9,6) << this->R () << ");" << endl;
671 out << "parameters.p1.set(" << FIXED(9,6) << this->p1() << ");" << endl;
672 out << "parameters.p2.set(" << FIXED(9,6) << this->p2() << ");" << endl;
673 out << "parameters.p3.set(" << FIXED(9,6) << this->p3() << ");" << endl;
674 out << "parameters.p4.set(" << FIXED(9,6) << this->p4() << ");" << endl;
675 out << "cc " << FIXED(9,6) << this->cc() << endl;
676 out << "bc " << FIXED(9,6) << this->bc() << endl;
677 out << endl;
678
679 out << "JK40DefaultSimulator.hh" << endl;
680 out << "static constexpr double p1 = " << FIXED(9,6) << this->p1() << ";" << endl;
681 out << "static constexpr double p2 = " << FIXED(9,6) << this->p2() << ";" << endl;
682 out << "static constexpr double p3 = " << FIXED(9,6) << this->p3() << ";" << endl;
683 out << "static constexpr double p4 = " << FIXED(9,6) << this->p4() << ";" << endl;
684 out << endl;
685 }
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 837 of file JFitK40.hh.

◆ R

JParameter_t JCALIBRATE::JK40Parameters_t::R
inherited

maximal coincidence rate [Hz]

Definition at line 711 of file JFitK40.hh.

◆ p1

JParameter_t JCALIBRATE::JK40Parameters_t::p1
inherited

1st order angle dependence coincidence rate

Definition at line 712 of file JFitK40.hh.

◆ p2

JParameter_t JCALIBRATE::JK40Parameters_t::p2
inherited

2nd order angle dependence coincidence rate

Definition at line 713 of file JFitK40.hh.

◆ p3

JParameter_t JCALIBRATE::JK40Parameters_t::p3
inherited

3rd order angle dependence coincidence rate

Definition at line 714 of file JFitK40.hh.

◆ p4

JParameter_t JCALIBRATE::JK40Parameters_t::p4
inherited

4th order angle dependence coincidence rate

Definition at line 715 of file JFitK40.hh.

◆ cc

JParameter_t JCALIBRATE::JK40Parameters_t::cc
inherited

fraction of signal correlated background

Definition at line 716 of file JFitK40.hh.

◆ bc

JParameter_t JCALIBRATE::JK40Parameters_t::bc
inherited

constant background

Definition at line 717 of file JFitK40.hh.


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