1 #ifndef __JFIT__JSHOWERNPE__
2 #define __JFIT__JSHOWERNPE__
17 namespace JPP {
using namespace JFIT; }
81 inline double getH0(
const size_t M)
const
83 return this->
getY1(M);
93 inline double getH1(
const double E_GeV)
const
95 using namespace JTOOLS;
96 using namespace JPHYSICS;
98 return E_GeV * this->
getY();
109 double getP(
const double E_GeV,
const bool hit)
const
122 double getChi2(
const double E_GeV,
const bool hit)
const
double getChi2(const double E_GeV, const bool hit) const
Get chi2 for observing a hit or not as a function of shower energy.
JShowerNPE(const JK40 &k40, const double y)
Constructor.
double getY0() const
Get expectation value for number of single hits.
Auxiliary class for converting various rates to expectation values of the number of hits within a giv...
double getP(const double E_GeV, const bool hit) const
Get probability for observing a hit or not as a function of shower energy.
double getH0() const
Expected number of photo-electrons for random background hypothesis.
double getY1(const size_t M) const
Get expectation value for number of multiple hits.
double getY() const
Get light yield due to bremsstrahlung.
double getP(size_t M) const
Get probability to observe a hit with given multiplicity.
double __y
light yield [npe/GeV]
double getP(const double expval, bool hit)
Get Poisson probability to observe a hit or not for given expectation value for the number of hits...
JShowerNPE()
Default constructor.
double getH1(const double E_GeV) const
Expected number of photo-electrons for EM shower hypothesis as a function of shower energy...
Auxiliary class for handling EM shower light yield.
double getChi2(const double P)
Get chi2 corresponding to given probability.
double getH0(const size_t M) const
Expected number of photo-electrons for random background hypothesis.