JPHYSICS::JK40Rates Struct Reference

Auxiliary class for K40 rates. More...

#include <JK40Rates.hh>

Inheritance diagram for JPHYSICS::JK40Rates:

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

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

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

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

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

JPHYSICS::JK40Rates::JK40Rates ( )

inline

Default constructor.

Definition at line 45 of file JK40Rates.hh.

                :
      rateL0(0.0),
      rateL1()
    {}

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

double JPHYSICS::JK40Rates::getSinglesRate ( ) const

inline

Get singles rate.

Returns

rate [Hz]

Definition at line 71 of file JK40Rates.hh.

    {
      return rateL0;
    }

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

inline

Get multiples rate.

Returns

rate [Hz]

Definition at line 82 of file JK40Rates.hh.

    {
      return rateL1;
    }

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;
    }

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;
    }

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;
    }

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 Function Documentation

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;
    }

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

const multiplicity_type JPHYSICS::JK40Rates::LOWER_L1_MULTIPLICITY = 2

static

Lower L1 multiplicity.

Definition at line 217 of file JK40Rates.hh.

JRateL0_t JPHYSICS::JK40Rates::rateL0

protected

singles rate [Hz]

Definition at line 221 of file JK40Rates.hh.

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