1 #ifndef __JPHYSICS__JK40RATES___
2 #define __JPHYSICS__JK40RATES___
14 namespace JPP {
using namespace JPHYSICS; }
79 if (rates_Hz.
rateL1.empty()) {
83 rates_Hz.
rateL1.push_back(500.0);
84 rates_Hz.
rateL1.push_back( 50.0);
85 rates_Hz.
rateL1.push_back( 5.0);
175 for (JRateL1_t::iterator i =
rateL1.begin(); i !=
rateL1.end(); ++i) {
191 const double rateL0 =
object.rateL0;
193 if (in >>
object.rateL0) {
195 object.rateL1.clear();
197 for (
double x; in >> x; ) {
198 object.rateL1.push_back(x);
219 out <<
object.rateL0;
221 for (JRateL1_t::const_iterator i =
object.
rateL1.begin(); i !=
object.rateL1.end(); ++i) {
JK40Rates()
Default constructor.
Auxiliary methods for mathematics.
JK40Rates(const JRateL0_t rateL0_Hz, const JRateL1_t &rateL1_Hz=JRateL1_t())
Constructor.
double binomial(const int n, const int k)
Binomial function.
static const multiplicity_type LOWER_L1_MULTIPLICITY
Lower L1 multiplicity.
multiplicity_type getLowerL1Multiplicity() const
Get lower multiplicty.
friend std::istream & operator>>(std::istream &in, JK40Rates &object)
Read K40 rates from input.
JRateL1_t rateL1
multiples rates [Hz]
friend std::ostream & operator<<(std::ostream &out, const JK40Rates &object)
Write K40 rates to output.
T pow(const T &x, const double y)
Power .
double getMultiplesRate(const multiplicity_type M) const
Get multiples rate.
then usage $script[distance] fi case set_variable R
double getSinglesRate() const
Get singles rate.
size_t multiplicity_type
Type definition of multiplicity.
multiplicity_type getUpperL1Multiplicity() const
Get upper multiplicty.
void correct(const double QE)
Correct rates for global efficiency,.
std::vector< double > JRateL1_t
Type definition of count rate as a function of multiplicty [Hz] The multiples rate start counting at ...
double JRateL0_t
Type definition of singles rate [Hz].
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID CHECK_EXIT_CODE typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
static JK40Rates & getInstance()
Get reference to unique instance of this class object.
JRateL0_t rateL0
singles rate [Hz]
Auxiliary class for K40 rates.