1 #ifndef __JPHYSICS__JPDFTOOLKIT__
2 #define __JPHYSICS__JPDFTOOLKIT__
19 namespace JPP {
using namespace JPHYSICS; }
62 const double x = n*lambda;
98 const double gamma = E / M;
99 const double beta = sqrt((1.0 + 1.0/gamma) * (1.0 - 1.0/gamma));
102 (1.0 + 2.0*gamma*ratio + ratio*ratio) );
133 const double K = 0.307075;
135 const double gamma = E / M;
136 const double beta = sqrt((1.0 + 1.0/gamma) * (1 - 1.0/gamma));
138 const double a = 0.25/(E*
E);
139 const double b = beta*beta/Tmax;
140 const double c = 1.0;
142 const double W = 0.5 * K * (Z/
A) * (1.0/(beta*beta));
144 const double sT = Tmax + Tmin;
145 const double dT = Tmax - Tmin;
146 const double rT = Tmax / Tmin;
148 const double weight = W * (a*sT*dT - b*dT + c*log(rT));
175 template<
class JFormFactor_t>
182 const JFormFactor_t&
F,
183 const int N = 1000000)
189 const double K = 0.307075;
191 const double gamma = E / M;
192 const double beta = sqrt((1.0 + 1.0/gamma) * (1 - 1.0/gamma));
194 const double W = 0.5 * K * (Z/
A) * (1.0/(beta*beta));
196 const double xmin = log(Tmin);
197 const double xmax = log(Tmax);
198 const double dx = (xmax -
xmin) / ((
double)
N);
202 for (
double x = xmin;
x <=
xmax;
x += dx) {
204 const double T =
exp(
x);
205 const double y = W *
F(T) * dx;
236 static const double Z = 10.0;
237 static const double A = 18.0;
263 static const double Z = 10.0;
264 static const double A = 18.0;
288 static const double a = 3.195e-01;
289 static const double b = 3.373e-01;
290 static const double c = -2.731e-02;
291 static const double d = 1.610e-03;
292 static const double Emin = 0.13078;
296 const double x =
log10(E);
297 const double y = a + x*(b + x*(c + x*(
d)));
325 static const double Z = 10.0;
326 static const double A = 18.0;
350 static const double a = -2.178e-01;
351 static const double b = 4.995e-01;
352 static const double c = -3.906e-02;
353 static const double d = 1.615e-03;
354 static const double Emin = 2.19500;
358 const double x =
log10(E);
359 const double y = a + x*(b + x*(c + x*(
d)));
381 return 0.188 *
exp(-1.25 *
pow(fabs(x - 0.90), 1.30));
395 static const double d = 0.36;
396 static const double U =
PI*
PI*
PI*
PI*
PI*2.0/3.0;
397 static const double V = d*d*d*d*d*
d;
399 const double W = (n*n - 1.0) / (n*n + 2.0);
400 const double sigma = 1.0e-14 * U*V*W*W / (lambda*lambda*lambda*lambda);
416 static const double amu = 18.01528;
static const uint32_t K[64]
then usage $script< input file >[option[primary[working directory]]] nWhere option can be E
double getDeltaRayTmax(const double E, const double M)
Get maximum delta-ray kinetic energy for given lepton energy and mass.
double getDeltaRayTmin()
Get minimum delta-ray kinetic energy.
static const double MASS_MUON
muon mass [GeV]
static const double DELTARAY_TMAX
Maximum allowed delta-ray kinetic energy [GeV].
static const double AVOGADRO
Avogadro's number.
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)
double getDeltaRaysFromTauFit(const double E)
Equivalent EM-shower energy due to delta-rays per unit tau track length.
const double getRayleighCrossSection(const double n, const double lambda)
Rayleigh cross section.
double getDeltaRaysFromMuon(const double E, const JRange< double > T_GeV=JRange< double >(DELTARAY_TMIN, DELTARAY_TMAX))
Equivalent EM-shower energy due to delta-rays per unit muon track length.
double getMaximalWavelength()
Get maximal wavelength for PDF evaluations.
set_variable E_E log10(E_{fit}/E_{#mu})"
static const double MASS_TAU
tau mass [GeV]
do set_variable OUTPUT_DIRECTORY $WORKDIR T
T pow(const T &x, const double y)
Power .
static const double PI
Mathematical constants.
double getDeltaRays(const double E, const double M, const double Tmin, const double Tmax, const double Z, const double A)
Get equivalent EM-shower energy due to delta-rays per unit track length for an ionising particle with...
static const double MASS_ELECTRON
electron mass [GeV]
double getKineticEnergy(const Trk &trk)
Get track kinetic energy.
double getDeltaRaysFromMuonFit(const double E)
Equivalent EM-shower energy due to delta-rays per unit muon track length.
$WORKDIR ev_configure_dqsimulator txt echo process $DQ_SIMULATOR $i $SOURCE_HOST[$index] csh c(setenv ROOTSYS $ROOTSYS &&source $JPP_DIR/setenv.csh $JPP_DIR &&($DQ_SIMULATOR\-u\$NAME\$\-H\$SERVER\$\-M\$LOGGER\$\-d $DEBUG</dev/null > &/dev/null &))'
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
then usage $script< input file >[option[primary[working directory]]] nWhere option can be N
Auxiliary class to define a range between two values.
static const double ALPHA_ELECTRO_MAGNETIC
Electro-Magnetic coupling constant.
double getDeltaRayProbability(const double x)
Emission profile of photons from delta-rays.
double getDeltaRaysFromTau(const double E, const JRange< double > T_GeV=JRange< double >(DELTARAY_TMIN, DELTARAY_TMAX))
Equivalent EM-shower energy due to delta-rays per unit tau track length.
do set_variable MODULE getModule a $WORKDIR detector_a datx L $STRING JEditDetector a $WORKDIR detector_a datx M $MODULE setz o $WORKDIR detector_a datx JEditDetector a $WORKDIR detector_b datx M $MODULE setz o $WORKDIR detector_b datx done echo Output stored at $WORKDIR detector_a datx and $WORKDIR tripod_a txt JDrawDetector2D a $WORKDIR detector_a datx a $WORKDIR detector_b datx L BL o detector $FORMAT $BATCH JDrawDetector2D T $WORKDIR tripod_a txt T $WORKDIR tripod_b txt L BL o tripod $FORMAT $BATCH JCompareDetector a $WORKDIR detector_a datx b $WORKDIR detector_b datx o $WORKDIR abc root &dev null for KEY in X Y Z
double getDeltaRaysFromElectron(const double E, const JRange< double > T_GeV=JRange< double >(DELTARAY_TMIN, DELTARAY_TMAX))
Equivalent EM-shower energy due to delta-rays per unit electron track length.
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
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
double getSinThetaC()
Get average sine of Cherenkov angle of water corresponding to group velocity.
static const double DELTARAY_TMIN
Minimum allowed delta-ray kinetic energy [GeV].