1 #ifndef __JPHYSICS__JPDFTOOLKIT__
2 #define __JPHYSICS__JPDFTOOLKIT__
16 namespace JPP {
using namespace JPHYSICS; }
53 const double x = n*lambda;
69 static const double a = 3.186e-01;
70 static const double b = 3.384e-01;
71 static const double c = -2.759e-02;
72 static const double d = 1.630e-03;
73 static const double Emin = 0.13078;
77 const double x = log10(E);
78 const double y = a + x*(b + x*(c + x*(
d)));
97 static const double a = -2.374e-01;
98 static const double b = 5.143e-01;
99 static const double c = -4.213e-02;
100 static const double d = 1.804e-03;
101 static const double Emin = 2.19500;
105 const double x = log10(E);
106 const double y = a + x*(b + x*(c + x*(
d)));
125 static const double d = 0.36;
126 static const double U =
PI*
PI*
PI*
PI*
PI*2.0/3.0;
127 static const double V = d*d*d*d*d*
d;
129 const double W = (n*n - 1.0) / (n*n + 2.0);
130 const double sigma = 1.0e-14 * U*V*W*W / (lambda*lambda*lambda*lambda);
146 static const double amu = 18.01528;
double getDeltaRaysFromTau(const double E)
Equivalent EM-shower energy due to delta-rays per unit tau track length.
static const double AVOGADRO
Avogadro's number [gr^-1].
static const double DENSITY_SEA_WATER
Fixed environment values.
double getMinimalWavelength()
Get minimal wavelength for PDF evaluations.
const double getRayleighScatteringLength(const double n, const double lambda)
Rayleigh scattering length.
const double getRayleighCrossSection(const double n, const double lambda)
Rayleigh cross section.
double getMaximalWavelength()
Get maximal wavelength for PDF evaluations.
static const double PI
Mathematical constants.
double cherenkov(const double lambda, const double n)
Number of Cherenkov photons per unit track length and per unit wavelength.
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
static const double ALPHA_ELECTRO_MAGNETIC
Electro-Magnetic coupling constant.
double getDeltaRaysFromMuon(const double E)
Equivalent EM-shower energy due to delta-rays per unit muon track length.