1 #ifndef __JPHYSICS__JGEANZ__
2 #define __JPHYSICS__JGEANZ__
15 namespace JPP {
using namespace JPHYSICS; }
61 const double a =
a0 +
a1 * log(E);
62 const double y =
pow(z,a-1.0) *
exp(-z/
b) / (
pow(
b,a) * std::tgamma(a));
100 const double a =
a0 +
a1 * log(E);
101 const double x = z /
b;
124 const double eps = 1.0e-3)
const
131 const double Q = P * (1.0 - eps);
133 for (
int i = 100; i != 0; --i) {
135 const double z = 0.5 * (zmin + zmax);
138 if (fabs(p-Q) < p*eps) {
148 return 0.5 * (zmin + zmax);
164 const double a =
a0 +
a1 * log(E);
Q(UTCMax_s-UTCMin_s)-livetime_s
Auxiliary methods for mathematics.
Function object for longitudinal profile of EM-shower.
double Gamma(const double a, const double x)
Incomplete gamma function.
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable STRING $argv[2] set_array QUANTILES set_variable FORMULA *[0] exp(-0.5 *(x-[1])*(x-[1])/([2]*[2]))" set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"` source JAcoustics.sh -- typeset -A TRIPODS get_tripods $WORKDIR/tripod.txt TRIPODS XMEAN
JGeanz(const double __a0, const double __a1, const double __b)
constructor
static const JGeanz geanz(1.85, 0.62, 0.54)
Function object for longitudinal EM-shower profile.
T pow(const T &x, const double y)
Power .
static double getMinimalShowerSize()
Get minimal shower size.
double getIntegral(const double E, const double z) const
Integral of PDF (starting from 0).
double getMaximum(const double E) const
Get depth of shower maximum.
double getLength(const double E, const double P, const double eps=1.0e-3) const
Get shower length for a given integrated probability.
double getProbability(const double E, const double z) const
Probability Density Function.
double operator()(const double E, const double z) const
Probability Density Function.