1 #ifndef __JSIRENE__JPYTHIA__
2 #define __JSIRENE__JPYTHIA__
72 static const double a = 72.425;
73 static const double b = -49.417;
74 static const double c = 5.858;
75 static const double d = 207.252;
76 static const double e = 132.784;
77 static const double f = -10.277;
78 static const double g = -19.441;
79 static const double h = 58.598;
80 static const double i = 53.161;
81 static const double kref = 2.698;
83 static const double u = (
a - f) / kref;
85 const double x = log10(Es);
86 const double y = (e +
x*(d +
x*(c +
x*(b +
x*(
a +
x*
u))))) / (i +
x*(h +
x*(g +
x*(f +
x*
u))));
88 return pow(10.0,
y - kref);
Extensions to Evt data format.
@ TRACK_TYPE_CHARGED_PION_MINUS
@ TRACK_TYPE_CHARGED_PION_PLUS
@ TRACK_TYPE_NEUTRAL_ANTIPION
@ TRACK_TYPE_NEUTRAL_PION
@ TRACK_TYPE_ANTIELECTRON
T pow(const T &x, const double y)
Power .
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
static const JPythia pythia
Function object for relative light yield as a function of GEANT particle code.
Auxiliary class to determine EM-equivalent energy as a function of PDG particle code and energy.
double operator()(const int type, const double Es) const
Get EM-equivalent energy.
JPythia()
Default constructor.
static double getEnergy(const double Es)
Get equivalent EM-energy for given pion energy.