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;
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)));
126 return 0.188 *
exp(-1.25 *
pow(fabs(x - 0.90), 1.30));
140 static const double d = 0.36;
141 static const double U =
PI*
PI*
PI*
PI*
PI*2.0/3.0;
142 static const double V = d*d*d*d*d*
d;
144 const double W = (n*n - 1.0) / (n*n + 2.0);
145 const double sigma = 1.0e-14 * U*V*W*W / (lambda*lambda*lambda*lambda);
161 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.
V(JDAQEvent-JTriggerReprocessor)*1.0/(JDAQEvent+1.0e-10)
const double getRayleighCrossSection(const double n, const double lambda)
Rayleigh cross section.
double getMaximalWavelength()
Get maximal wavelength for PDF evaluations.
set_variable E_E log10(E_{fit}/E_{#mu})"
T pow(const T &x, const double y)
Power .
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.
double getDeltaRayProbability(const double x)
Emission profile of photons from delta-rays.
$WORKDIR ev_configure_domsimulator txt echo process $DOM_SIMULATOR $i $SOURCE_HOST[$index] csh c(setenv ROOTSYS $ROOTSYS &&source $JPP_DIR/setenv.csh $JPP_DIR &&($DOM_SIMULATOR\-u\$NAME\$\-H\$SERVER\$\-M\$LOGGER\$\-d $DEBUG</dev/null > &/dev/null &))'
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 JAcousticsToolkit.sh typeset -A TRIPODS get_tripods $WORKDIR/tripod.txt TRIPODS XMEAN