1 #ifndef __JSIRENE__JPYTHIA__
2 #define __JSIRENE__JPYTHIA__
14 namespace JPP {
using namespace JSIRENE; }
38 using namespace JAANET;
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);
static const JPythia pythia
Function object for relative light yield as a function of GEANT particle code.
fi JEventTimesliceWriter a
static double getEnergy(const double Es)
Get equivalent EM-energy for given pion energy.
then print_variable DETECTOR INPUT_FILE INTERMEDIATE_FILE check_input_file $DETECTOR $INPUT_FILE check_output_file $INTERMEDIATE_FILE $OUTPUT_FILE JMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Auxiliary class to determine EM-equivalent energy as a function of PDG particle code and energy...
JPythia()
Default constructor.
double operator()(const int type, const double Es) const
Get EM-equivalent energy.