19namespace JPP {
using namespace JFIT; }
117 inline double getH1(
const double E_GeV)
const
121 const double E = gWater.getE(E_GeV, this->
getZ());
123 if (E >= MASS_MUON * INDEX_OF_REFRACTION_WATER)
124 return (this->
getY1() +
125 this->
getYA() * getDeltaRaysFromMuon(E) +
139 double getP(
const double E_GeV,
const bool hit)
const
152 double getChi2(
const double E_GeV,
const bool hit)
const
Auxiliary classes and methods for linear and iterative data regression.
double getChi2(const double P)
Get chi2 corresponding to given probability.
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.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Auxiliary class for handling light yields due to K40 decays.
double getH0() const
Expected number of photo-electrons for random background hypothesis.
double y0
light yield due to random background [npe]
Auxiliary class for handling various light yields.
double getChi2(const double E_GeV, const bool hit) const
Get chi2 for observing a hit or not as a function of muon energy.
double yA
light yield due to delta-rays [npe*m/GeV]
double z
position along muon path [m]
JNPE()
Default constructor.
double getYB() const
Get light yield due to bremsstrahlung.
double getP(const double E_GeV, const bool hit) const
Get probability for observing a hit or not as a function of muon energy.
double y1
light yield due to minimum ionizing particle [npe]
JNPE(const double y0, const double y1, const double yA, const double yB, const double z)
Constructor.
double getY1() const
Get light yield due to minimum ionizing particle.
double yB
light yield due to bremsstrahlung [npe/GeV]
double getH1(const double E_GeV) const
Expected number of photo-electrons for muon hypothesis as a function of muon energy.
double getZ() const
Get position along muon path.
double getYA() const
Get light yield due to delta-rays.