1 #ifndef __JPHYSICS__JGEANX__
2 #define __JPHYSICS__JGEANX__
76 const double xmax)
const
91 const double x = fabs(ct - 1.0/
n);
106 const double xmax)
const
111 const double x = 1.0 /
n;
113 const double ai = 1.0 /
a;
114 const double bi = 1.0 /
b;
116 const double xl =
pow(fabs(
x -
xmin),
a) * -
b;
117 const double xr =
pow(fabs(
x -
xmax),
a) * -
b;
119 const double gp = tgamma(ai);
121 double yl =
Gamma(ai, xl) * gp;
122 double yr =
Gamma(ai, xr) * gp;
129 return ai *
pow(-bi,ai) * (yl + yr);
Auxiliary methods for mathematics.
Probability density function of photon emission from EM-shower as a function of cosine of the emissio...
double evaluate(const double ct) const
Functional dependence.
double evaluate(const double xmin, const double xmax) const
Integral.
double c
normalisation constant
JGeanx(const double __a, const double __b, const double __n=getIndexOfRefractionPhase())
Constructor.
double operator()(const double xmin, const double xmax) const
Integral number of photons from EM-shower between two emission angles.
double operator()(const double ct) const
Number of photons from EM-shower as a function of emission angle.
const double n
index of refraction
double Gamma(const double a, const double x)
Incomplete gamma function.
T pow(const T &x, const double y)
Power .
static const double PI
Mathematical constants.
Auxiliary methods for light properties of deep-sea water.
static const JGeanx geanx(0.35, -5.40)
Function object for the number of photons from EM-shower as a function of emission angle.
double getIndexOfRefractionPhase()
Get average index of refraction of water corresponding to phase velocity.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).