JPHYSICS::JK40Rates Struct Reference

Auxiliary class for K40 rates. More...

#include <JK40Rates.hh>

Inheritance diagram for JPHYSICS::JK40Rates:

Public Member Functions
	JK40Rates ()
	Default constructor.
	JK40Rates (const JRateL0_t rateL0_Hz, const JRateL1_t &rateL1_Hz=JRateL1_t())
	Constructor.
double	getSinglesRate () const
	Get singles rate.
const JRateL1_t &	getMultiplesRates () const
	Get multiples rate.
double	getMultiplesRate (const multiplicity_type M) const
	Get multiples rate at given multiplicity.
multiplicity_type	getLowerL1Multiplicity () const
	Get lower multiplicty.
multiplicity_type	getUpperL1Multiplicity () const
	Get upper multiplicty.
void	correct (const double QE)
	Correct rates for global efficiency,.

Static Public Attributes
static const multiplicity_type	LOWER_L1_MULTIPLICITY = 2
	Lower L1 multiplicity.

Protected Attributes
JRateL0_t	rateL0
	singles rate [Hz]
JRateL1_t	rateL1
	multiples rates [Hz]

Friends
std::istream &	operator>> (std::istream &in, JK40Rates &object)
	Read K40 rates from input.
std::ostream &	operator<< (std::ostream &out, const JK40Rates &object)
	Write K40 rates to output.

Detailed Description

Auxiliary class for K40 rates.

The singles rate refers to the counting rate of a PMT and the multiples rate to the rate of genuine coincidences due to K40 decays.

Definition at line 41 of file JK40Rates.hh.

Constructor & Destructor Documentation

JK40Rates() [1/2]

JPHYSICS::JK40Rates::JK40Rates ( )

inline

Default constructor.

Definition at line 45 of file JK40Rates.hh.

                :
      rateL0(0.0),
      rateL1()
    {}

JK40Rates() [2/2]

JPHYSICS::JK40Rates::JK40Rates ( const JRateL0_t rateL0_Hz, const JRateL1_t & rateL1_Hz = JRateL1_t() )

inline

Constructor.

The multiples rates start counting at two-fold coincidences.

Parameters

rateL0_Hz	singles rate [Hz]
rateL1_Hz	multiples rates [Hz]

Definition at line 59 of file JK40Rates.hh.

                                                        :
      rateL0(rateL0_Hz),
      rateL1(rateL1_Hz)
    {}

Member Function Documentation

getSinglesRate()

double JPHYSICS::JK40Rates::getSinglesRate ( ) const

inline

Get singles rate.

Returns

rate [Hz]

Definition at line 71 of file JK40Rates.hh.

    {
      return rateL0;
    }

getMultiplesRates()

const JRateL1_t & JPHYSICS::JK40Rates::getMultiplesRates ( ) const

inline

Get multiples rate.

Returns

rate [Hz]

Definition at line 82 of file JK40Rates.hh.

    {
      return rateL1;
    }

getMultiplesRate()

double JPHYSICS::JK40Rates::getMultiplesRate ( const multiplicity_type M ) const

inline

Get multiples rate at given multiplicity.

Parameters

M	multiplicity (M >= JK40Rates::LOWER_L1_MULTIPLICITY)

Returns

rate [Hz]

Definition at line 94 of file JK40Rates.hh.

    {
      if (M >= LOWER_L1_MULTIPLICITY && M - LOWER_L1_MULTIPLICITY < (multiplicity_type) rateL1.size())
        return rateL1[M - LOWER_L1_MULTIPLICITY];
      else
        return 0.0;
    }

getLowerL1Multiplicity()

multiplicity_type JPHYSICS::JK40Rates::getLowerL1Multiplicity ( ) const

inline

Get lower multiplicty.

Returns

lower multiplicity

Definition at line 108 of file JK40Rates.hh.

    {
      return LOWER_L1_MULTIPLICITY;
    }

getUpperL1Multiplicity()

multiplicity_type JPHYSICS::JK40Rates::getUpperL1Multiplicity ( ) const

inline

Get upper multiplicty.

Returns

upper multiplicity

Definition at line 119 of file JK40Rates.hh.

    {
      return rateL1.size() + 1;
    }

correct()

void JPHYSICS::JK40Rates::correct ( const double QE )

inline

Correct rates for global efficiency,.

Parameters

QE	global efficiency

Definition at line 130 of file JK40Rates.hh.

    {
      if (QE > 0.0) {
      
        rateL0 /= QE;
 
        JRateL1_t buffer = rateL1;
        
        for (multiplicity_type M = getUpperL1Multiplicity(); M >= getLowerL1Multiplicity(); --M) {
 
          // determine contribution from higher multiplicities
          
          double R = 0.0;
          
          for (multiplicity_type i = M + 1; i <= getUpperL1Multiplicity(); ++i) {
            R += buffer[i - LOWER_L1_MULTIPLICITY] * JMATH::binomial(i, M) * pow(QE, M) * pow(1.0 - QE, i - M);
          }
 
          if (getMultiplesRate(M) > R)
            buffer[M - LOWER_L1_MULTIPLICITY] = (getMultiplesRate(M) - R) / pow(QE, M);
          else
            buffer[M - LOWER_L1_MULTIPLICITY] =  0.0;
        }
 
        rateL1 = buffer;
        
      } else {
 
        rateL0 = 0.0;
 
        for (JRateL1_t::iterator i = rateL1.begin(); i != rateL1.end(); ++i) {
          *i = 0.0;
        }
      }
    }

Friends And Related Symbol Documentation

operator>>

std::istream & operator>> ( std::istream & in, JK40Rates & object )

friend

Read K40 rates from input.

Parameters

in	input stream
object	K40 rates

Returns

input stream

Definition at line 174 of file JK40Rates.hh.

    {
      const double rateL0 = object.rateL0;
 
      if (in >> object.rateL0) {
 
        object.rateL1.clear();
 
        for (double x; in >> x; ) {
          object.rateL1.push_back(x);
        }
 
      } else {
 
        object.rateL0 = rateL0;
      }
 
      return in;
    }

operator<<

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

friend

Write K40 rates to output.

Parameters

out	output stream
object	K40 rates

Returns

output stream

Definition at line 202 of file JK40Rates.hh.

    {
      out << object.rateL0;
 
      for (JRateL1_t::const_iterator i = object.rateL1.begin(); i != object.rateL1.end(); ++i) {
        out << ' ' << *i;
      }
 
      return out;
    }

Member Data Documentation

LOWER_L1_MULTIPLICITY

const multiplicity_type JPHYSICS::JK40Rates::LOWER_L1_MULTIPLICITY = 2

static

Lower L1 multiplicity.

Definition at line 217 of file JK40Rates.hh.

rateL0

JRateL0_t JPHYSICS::JK40Rates::rateL0

protected

singles rate [Hz]

Definition at line 221 of file JK40Rates.hh.

rateL1

JRateL1_t JPHYSICS::JK40Rates::rateL1

protected

multiples rates [Hz]

Definition at line 222 of file JK40Rates.hh.

---

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

• JK40Rates.hh