Go to the documentation of this file. 1 #ifndef __JCALIBRATE_JCALIBRATEK40__
2 #define __JCALIBRATE_JCALIBRATEK40__
30 static const char*
const WS_t =
"WS";
31 static const char*
const WB_t =
"WB";
33 static const char*
const _2S =
".2S";
34 static const char*
const _1B =
".1B";
35 static const char*
const _1L =
".1L";
36 static const char*
const _2R =
".2R";
37 static const char*
const _2F =
".2F";
53 const double epsilon = 1.0e-5) :
72 else if (max(first.first, first.second) == max(second.first, second.second))
73 return min(first.first, first.second) < min(second.first, second.second);
75 return max(first.first, first.second) < max(second.first, second.second);
107 inline double coincidenceP(
double E1,
double E2,
double ED,
int M_min,
int M_max)
110 double P = exp(-E1) * E1 * exp(-E2) * E2 * exp(-ED);
113 for (
int m = 2; m < M_min; m++) {
114 P = P * ( ED/(m-1) + E1/m + E2/m );
120 for (
int m = M_min+1; m <= M_max; m++) {
121 P = P * ( ED/(m-1) + E1/m + E2/m );
static const char *const _2F
Name extension for 2F rate.
static const char *const _1L
Name extension for 1D live time.
static const char *const _2S
Name extension for 2D counts.
static const char *const weights_hist_t
Histogram naming.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
double coincidenceP(double E1, double E2, double ED, int M_min, int M_max)
Coincidence probability of two PMTs.
const JPMT & getPMT(const int index) const
Get PMT.
static const char *const WB_t
Named bin for value of TMax_ns in JCalibrateK40.
Data structure for a composite optical module.
static const char *const W1_overall_t
Named bin for duration of the run.
JPairwiseComparator(const JModule &module, const double epsilon=1.0e-5)
Constructor.
JCombinatorics::pair_type pair_type
double getDot(const pair_type &pair) const
Get cosine of space angle between PMT axes.
static const char *const _1B
Name extension for 1D background.
Auxiliary class to sort pairs of PMT addresses within optical module.
bool operator()(const pair_type &first, const pair_type &second) const
Comparison of two pairs of PMT addresses.
static const char *const WS_t
Named bin for time residual bin width.
static const char *const _2R
Name extension for 2D rate.
double getDot(const JFirst_t &first, const JSecond_t &second)
Get dot product of objects.