Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Classes | Typedefs | Enumerations | Functions | Variables
JPHYSICS Namespace Reference

Auxiliary classes and methods for calculation of PDF and muon energy loss. More...

Classes

class  JAbstractMedium
 Medium interface. More...
 
class  JAbstractPMT
 PMT interface. More...
 
class  JCDFTable
 Multi-dimensional CDF table for arrival time of Cherenkov light. More...
 
class  JDispersion
 Implementation of dispersion for water in deep sea. More...
 
class  JDispersionInterface
 Light dispersion inteface. More...
 
class  JGeane
 Interface for muon energy loss. More...
 
class  JGeane_t
 Function object for the energy loss of the muon. More...
 
class  JGeaneWater
 Function object for energy dependent energy loss of the muon. More...
 
class  JGeant
 Function object for the probability density function of photon emission from EM-shower as a function of the index of refraction and the cosine of the emission angle. More...
 
class  JGeant_t
 Base class for the probability density function of photon emission from EM-shower as a function of the index of refraction and the cosine of the emission angle. More...
 
class  JGeanx
 Probability density function of photon emission from EM-shower as a function of cosine of the emission angle. More...
 
class  JGeanz
 Function object for longitudinal profile of EM-shower. More...
 
class  JAbstractLED
 Interface for emission profile from LED. More...
 
class  JLED
 Probability Density Functions of the time response of a PMT. More...
 
class  JLED_C
 Probability Density Functions of the time response of a PMT (C-like interface) More...
 
class  JNPETable
 Custom class for integrated values of the PDF of the arrival time of Cherenkov light. More...
 
class  JNPETable1D
 Custom class for the integrated values of the PDF of the arrival time of Cherenkov light. More...
 
class  JPDF
 Low level interface for the calculation of the Probability Density Functions (PDFs) of the arrival time of Cherenkov light from a muon or an EM-shower on a photo-multiplier tube (PMT). More...
 
class  JAbstractPDF
 Probability Density Functions of the time response of a PMT with an implementation for the JDispersionInterface interface. More...
 
class  JPDF_C
 Probability Density Functions of the time response of a PMT with an implementation of the JAbstractPMT and JAbstractMedium interfaces via C-like methods. More...
 
class  JPDFTable
 Multi-dimensional PDF table for arrival time of Cherenkov light. More...
 
class  JAbsorptionLengthOfPureWater
 Absorption length of pure water. More...
 
class  JPDFTransformer_t
 Transformer for the 1D probability density function (PDF) of the time response of a PMT to a muon. More...
 
class  JPD0Transformer_t
 Transformer for the 1D probability density function (PDF) of the time response of a PMT due to a point source. More...
 
class  JPDGTransformer_t
 Transformer for the 2D probability density function (PDF) of the time response of a PMT due to an EM shower. More...
 
class  JPDFTransformer
 Template definition of transformer of the probability density function (PDF) of the time response of a PMT. More...
 
class  JPDFTransformer< 2, JArgument_t >
 Template specialisation of transformer of the 2D probability density function (PDF) of the time response of a PMT due to a bright point. More...
 
class  JPDFTransformer< 3, JArgument_t >
 Template specialisation of transformer of the 3D probability density function (PDF) of the time response of a PMT due to a muon. More...
 
class  JPDFTransformer< 4, JArgument_t >
 Template specialisation of transformer of the 4D probability density function (PDF) of the time response of a PMT due to an EM shower. More...
 
class  JPDFTransformer< 5, JArgument_t >
 Template specialisation of transformer of the 5D probability density function (PDF) of the time response of a PMT due to an EM shower. More...
 
class  JRadiation
 Auxiliary class for the calculation of the muon radiative cross sections. More...
 
class  JRadiationFunction
 Fast implementation of class JRadiation. More...
 
class  JRadiationInterface
 Interface for calculation of inverse interaction length and shower energy. More...
 
class  JRadiationSource
 Source of radiation. More...
 
class  JCDFTable1D
 Custom class for CDF table in 1 dimension. More...
 
class  JCDFTable2D
 Custom class for CDF table in 2 dimensions. More...
 

Typedefs

typedef
JTOOLS::JGridPolint1Function1D_t 
JGeantFunction1D_t
 
typedef JTOOLS::JElement2D
< double, double > 
JElement2D_t
 
typedef JTOOLS::JElement3D
< double, double > 
JElement3D_t
 

Enumerations

enum  JPDFType_t {
  DIRECT_LIGHT_FROM_MUON = 1, SCATTERED_LIGHT_FROM_MUON = 2, DIRECT_LIGHT_FROM_EMSHOWERS = 3, SCATTERED_LIGHT_FROM_EMSHOWERS = 4,
  DIRECT_LIGHT_FROM_DELTARAYS = 5, SCATTERED_LIGHT_FROM_DELTARAYS = 6, SCATTERED_LIGHT_FROM_MUON_5D = 12, DIRECT_LIGHT_FROM_EMSHOWER = 13,
  SCATTERED_LIGHT_FROM_EMSHOWER = 14, DIRECT_LIGHT_FROM_BRIGHT_POINT = 23, SCATTERED_LIGHT_FROM_BRIGHT_POINT = 24, LIGHT_FROM_MUON = 1001,
  LIGHT_FROM_EMSHOWERS = 1003, LIGHT_FROM_DELTARAYS = 1005, LIGHT_FROM_EMSHOWER = 1013, LIGHT_FROM_BRIGTH_POINT = 1023
}
 PDF types. More...
 

Functions

double geanc ()
 Equivalent muon track length per unit shower energy. More...
 
static const JGeane_t gRock (2.67e-1 *0.9 *DENSITY_ROCK, 3.40e-4 *1.2 *DENSITY_ROCK)
 Function object for energy loss of muon in rock. More...
 
double getMinimalWavelength ()
 Get minimal wavelength for PDF evaluations. More...
 
double getMaximalWavelength ()
 Get maximal wavelength for PDF evaluations. More...
 
double cherenkov (const double lambda, const double n)
 Number of Cherenkov photons per unit track length and per unit wavelength. More...
 
double getDeltaRaysFromMuon (const double E)
 Equivalent EM-shower energy due to delta-rays per unit muon track length. More...
 
double getDeltaRaysFromTau (const double E)
 Equivalent EM-shower energy due to delta-rays per unit tau track length. More...
 
const double getRayleighCrossSection (const double n, const double lambda)
 Rayleigh cross section. More...
 
const double getRayleighScatteringLength (const double n, const double lambda)
 Rayleigh scattering length. More...
 
std::string getLabel (const JPDFType_t pdf)
 Get PDF label. More...
 
int getPDFType (const std::string &file_name)
 Get PDF type. More...
 
bool hasWildCard (const std::string &file_name)
 Check wild card. More...
 
std::string getFilename (const std::string &file_name, const JPDFType_t pdf)
 Get PDF file name. More...
 
bool is_muon (const int pdf)
 Test if given PDF type corresponds to Cherenkov light from muon. More...
 
bool is_bremsstrahlung (const int pdf)
 Test if given PDF type corresponds to Cherenkov light from Bremsstrahlung. More...
 
bool is_deltarays (const int pdf)
 Test if given PDF type corresponds to Cherenkov light from delta-rays. More...
 
bool is_scattered (const int pdf)
 Test if given PDF type corresponds to scattered light. More...
 
double getThetaMCS (const double E, const double x, const double X0, const double M, const double Q)
 Get multiple Coulomb scattering angle. More...
 
double getThetaMCS (const double E, const double x)
 Get multiple Coulomb scattering angle. More...
 

Variables

static const JGeaneWater gWater
 Function object for energy loss of muon in sea water. More...
 
static const JGeant geant (geanx, 0.0001)
 Function object for the number of photons from EM-shower as a function of emission angle. More...
 
static const JGeanx geanx (0.35,-5.40)
 Function object for the number of photons from EM-shower as a function of emission angle. More...
 
static const JGeanz geanz (1.85, 0.62, 0.54)
 Function object for longitudinal EM-shower profile. More...
 
static const JCCnu cc_nu
 Function object for charged current neutrino cross section [cm^2] as a function of neutrino energy [GeV]. More...
 
static const JNCnu nc_nu
 Function object for neutral current neutrino cross section [cm^2] as a function of neutrino energy [GeV]. More...
 
static const JCCnubar cc_nubar
 Function object for charged current anti-neutrino cross section [cm^2] as a function of neutrino energy [GeV]. More...
 
static const JNCnubar nc_nubar
 Function object for neutral current anti-neutrino cross section [cm^2] as a function of neutrino energy [GeV]. More...
 
static double MODULE_RADIUS_M = 0.25
 Radius of optical module [m]. More...
 
static const
JAbsorptionLengthOfPureWater 
getAbsorptionLengthOfPureWater
 Function object for absorption length of pure water. More...
 
static const char WILD_CARD = '%'
 wild card character for file name substitution More...
 

Detailed Description

Auxiliary classes and methods for calculation of PDF and muon energy loss.

Author
mdejong

Typedef Documentation

Definition at line 21 of file JGeant_t.hh.

typedef JTOOLS::JElement2D<double, double> JPHYSICS::JElement2D_t

Definition at line 26 of file JLED.hh.

typedef JTOOLS::JElement3D<double, double> JPHYSICS::JElement3D_t

Definition at line 27 of file JLED.hh.

Enumeration Type Documentation

PDF types.

Enumerator
DIRECT_LIGHT_FROM_MUON 

direct light from muon

SCATTERED_LIGHT_FROM_MUON 

scattered light from muon

DIRECT_LIGHT_FROM_EMSHOWERS 

direct light from EM showers

SCATTERED_LIGHT_FROM_EMSHOWERS 

scattered light from EM showers

DIRECT_LIGHT_FROM_DELTARAYS 

direct light from delta-rays

SCATTERED_LIGHT_FROM_DELTARAYS 

scattered light from delta-rays

SCATTERED_LIGHT_FROM_MUON_5D 

scattered light from muon

DIRECT_LIGHT_FROM_EMSHOWER 

direct light from EM shower

SCATTERED_LIGHT_FROM_EMSHOWER 

scattered light from EM shower

DIRECT_LIGHT_FROM_BRIGHT_POINT 

direct light from bright point

SCATTERED_LIGHT_FROM_BRIGHT_POINT 

scattered light from bright point

LIGHT_FROM_MUON 

direct and scattered light from muon

LIGHT_FROM_EMSHOWERS 

direct and scattered light from EM showers

LIGHT_FROM_DELTARAYS 

direct and scattered light from delta-rays

LIGHT_FROM_EMSHOWER 

direct and scattered light from EM shower

LIGHT_FROM_BRIGTH_POINT 

direct and scattered light from brigth point

Definition at line 27 of file JPDFTypes.hh.

27  {
28 
29  DIRECT_LIGHT_FROM_MUON = 1, //!< direct light from muon
30  SCATTERED_LIGHT_FROM_MUON = 2, //!< scattered light from muon
31 
32  DIRECT_LIGHT_FROM_EMSHOWERS = 3, //!< direct light from EM showers
33  SCATTERED_LIGHT_FROM_EMSHOWERS = 4, //!< scattered light from EM showers
34 
35  DIRECT_LIGHT_FROM_DELTARAYS = 5, //!< direct light from delta-rays
36  SCATTERED_LIGHT_FROM_DELTARAYS = 6, //!< scattered light from delta-rays
37 
38  SCATTERED_LIGHT_FROM_MUON_5D = 12, //!< scattered light from muon
39 
40  DIRECT_LIGHT_FROM_EMSHOWER = 13, //!< direct light from EM shower
41  SCATTERED_LIGHT_FROM_EMSHOWER = 14, //!< scattered light from EM shower
42 
43  //DIRECT_LIGHT_FROM_SHOWER = 15, //!< direct + scattered light from HADRONIC shower
44 
45  DIRECT_LIGHT_FROM_BRIGHT_POINT = 23, //!< direct light from bright point
46  SCATTERED_LIGHT_FROM_BRIGHT_POINT = 24, //!< scattered light from bright point
47 
48  LIGHT_FROM_MUON = 1001, //!< direct and scattered light from muon
49  LIGHT_FROM_EMSHOWERS = 1003, //!< direct and scattered light from EM showers
50  LIGHT_FROM_DELTARAYS = 1005, //!< direct and scattered light from delta-rays
51  LIGHT_FROM_EMSHOWER = 1013, //!< direct and scattered light from EM shower
52  LIGHT_FROM_BRIGTH_POINT = 1023 //!< direct and scattered light from brigth point
53  };
direct and scattered light from brigth point
Definition: JPDFTypes.hh:52
direct and scattered light from EM shower
Definition: JPDFTypes.hh:51
direct and scattered light from muon
Definition: JPDFTypes.hh:48
scattered light from muon
Definition: JPDFTypes.hh:38
scattered light from EM shower
Definition: JPDFTypes.hh:41
direct light from EM showers
Definition: JPDFTypes.hh:32
direct light from bright point
Definition: JPDFTypes.hh:45
direct light from muon
Definition: JPDFTypes.hh:29
scattered light from muon
Definition: JPDFTypes.hh:30
scattered light from delta-rays
Definition: JPDFTypes.hh:36
direct light from EM shower
Definition: JPDFTypes.hh:40
scattered light from EM showers
Definition: JPDFTypes.hh:33
scattered light from bright point
Definition: JPDFTypes.hh:46
direct light from delta-rays
Definition: JPDFTypes.hh:35
direct and scattered light from delta-rays
Definition: JPDFTypes.hh:50
direct and scattered light from EM showers
Definition: JPDFTypes.hh:49

Function Documentation

double JPHYSICS::geanc ( )
inline

Equivalent muon track length per unit shower energy.

Returns
equivalent muon track length [m/Gev]

Definition at line 31 of file JGeane.hh.

32  {
33  return 4.7; // dx/dE [m/GeV]
34  }
static const JGeane_t JPHYSICS::gRock ( 2.67e-1 *0.9 *  DENSITY_ROCK,
3.40e-4 *1.2 *  DENSITY_ROCK 
)
static

Function object for energy loss of muon in rock.

double JPHYSICS::getMinimalWavelength ( )
inline

Get minimal wavelength for PDF evaluations.

Returns
wavelength of light [nm]

Definition at line 34 of file JPDFToolkit.hh.

35  {
36  return 300.0;
37  }
double JPHYSICS::getMaximalWavelength ( )
inline

Get maximal wavelength for PDF evaluations.

Returns
wavelength of light [nm]

Definition at line 45 of file JPDFToolkit.hh.

46  {
47  return 700.0;
48  }
double JPHYSICS::cherenkov ( const double  lambda,
const double  n 
)
inline

Number of Cherenkov photons per unit track length and per unit wavelength.

Parameters
lambdawavelength of light [nm]
nindex of refraction
Returns
number of photons per unit track length and per unit wavelength [m^-1 nm^-1]

Definition at line 58 of file JPDFToolkit.hh.

60  {
61  using namespace JTOOLS;
62 
63  const double x = n*lambda;
64 
65  return 1.0e9 * 2 * PI * ALPHA_ELECTRO_MAGNETIC * (n*n - 1.0) / (x*x);
66  }
static const double PI
Constants.
Definition: JConstants.hh:20
static const double ALPHA_ELECTRO_MAGNETIC
Electro-Magnetic coupling constant.
Definition: JConstants.hh:28
const int n
Definition: JPolint.hh:628
double JPHYSICS::getDeltaRaysFromMuon ( const double  E)
inline

Equivalent EM-shower energy due to delta-rays per unit muon track length.

Internal parameters are obtained with application [script] JDeltaRays[.sh].

Parameters
Emuon energy [GeV]
Returns
equivalent energy loss [GeV/m]

Definition at line 77 of file JPDFToolkit.hh.

78  {
79  static const double a = 3.186e-01;
80  static const double b = 3.384e-01;
81  static const double c = -2.759e-02;
82  static const double d = 1.630e-03;
83  static const double Emin = 0.13078; // [GeV]
84 
85  if (E > Emin) {
86 
87  const double x = log10(E); //
88  const double y = a + x*(b + x*(c + x*(d))); // [MeV g^-1 cm^2]
89 
90  return y * JTOOLS::DENSITY_SEA_WATER * 1.0e-1; // [GeV/m]
91  }
92 
93  return 0.0;
94  }
static const double DENSITY_SEA_WATER
Fixed environment values.
Definition: JConstants.hh:34
fi JEventTimesliceWriter a
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
Definition: JPath.sh:52
then usage $script[input file[working directory[option]]] nWhere option can be E
Definition: JMuonPostfit.sh:37
double JPHYSICS::getDeltaRaysFromTau ( const double  E)
inline

Equivalent EM-shower energy due to delta-rays per unit tau track length.

Internal parameters are obtained with application [script] JDeltaRays[.sh].

Parameters
Etau energy [GeV]
Returns
equivalent energy loss [GeV/m]

Definition at line 105 of file JPDFToolkit.hh.

106  {
107  static const double a = -2.374e-01;
108  static const double b = 5.143e-01;
109  static const double c = -4.213e-02;
110  static const double d = 1.804e-03;
111  static const double Emin = 2.19500; // [GeV]
112 
113  if (E > Emin) {
114 
115  const double x = log10(E); //
116  const double y = a + x*(b + x*(c + x*(d))); // [MeV g^-1 cm^2]
117 
118  return y * JTOOLS::DENSITY_SEA_WATER * 1.0e-1; // [GeV/m]
119  }
120 
121  return 0.0;
122  }
static const double DENSITY_SEA_WATER
Fixed environment values.
Definition: JConstants.hh:34
fi JEventTimesliceWriter a
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
Definition: JPath.sh:52
then usage $script[input file[working directory[option]]] nWhere option can be E
Definition: JMuonPostfit.sh:37
const double JPHYSICS::getRayleighCrossSection ( const double  n,
const double  lambda 
)
inline

Rayleigh cross section.

Parameters
nindex of refraction
lambdawavelength of light [nm]
Returns
cross section [cm^2]

Definition at line 132 of file JPDFToolkit.hh.

134  {
135  using JTOOLS::PI;
136 
137  static const double d = 0.36; // size of H2O molecule [nm]
138  static const double U = PI*PI*PI*PI*PI*2.0/3.0;
139  static const double V = d*d*d*d*d*d;
140 
141  const double W = (n*n - 1.0) / (n*n + 2.0);
142  const double sigma = 1.0e-14 * U*V*W*W / (lambda*lambda*lambda*lambda); // [cm^2]
143 
144  return sigma;
145  }
static const double PI
Constants.
Definition: JConstants.hh:20
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
Definition: JPath.sh:52
alias put_queue eval echo n
Definition: qlib.csh:19
const double JPHYSICS::getRayleighScatteringLength ( const double  n,
const double  lambda 
)
inline

Rayleigh scattering length.

Parameters
nindex of refraction
lambdawavelength of light [nm]
Returns
scattering length [m]

Definition at line 155 of file JPDFToolkit.hh.

157  {
158  using namespace JTOOLS;
159 
160  static const double amu = 18.01528; // H20 mass in Atomic mass units
161 
162  const double sigma = getRayleighCrossSection(n, lambda);
163  const double ls = 1.0e-2 / (DENSITY_SEA_WATER * AVOGADRO * sigma / amu); // [m]
164 
165  return ls;
166  }
static const double DENSITY_SEA_WATER
Fixed environment values.
Definition: JConstants.hh:34
const double getRayleighCrossSection(const double n, const double lambda)
Rayleigh cross section.
Definition: JPDFToolkit.hh:132
const int n
Definition: JPolint.hh:628
static const double AVOGADRO
Avogadro&#39;s number [gr^-1].
Definition: JConstants.hh:24
std::string JPHYSICS::getLabel ( const JPDFType_t  pdf)
inline

Get PDF label.

Parameters
pdfPDF type
Returns
PDF label

Definition at line 65 of file JPDFTypes.hh.

66  {
67  std::ostringstream os;
68 
69  os << pdf;
70 
71  return os.str();
72  }
int JPHYSICS::getPDFType ( const std::string &  file_name)
inline

Get PDF type.

Parameters
file_namefile name
Returns
PDF type (-1 in case of error)

Definition at line 81 of file JPDFTypes.hh.

82  {
83  using namespace std;
84 
85  static const char* digits = "0123456789";
86 
87  int type = -1;
88 
89  string buffer = JEEP::getFilename(file_name);
90 
91  string::size_type pos = buffer.find_first_of(digits);
92 
93  if (pos != string::npos) {
94 
95  string::size_type len = buffer.substr(pos).find_first_not_of(digits);
96 
97  istringstream(buffer.substr(pos, len)) >> type;
98  }
99 
100  return type;
101  }
std::string getFilename(const std::string &file_name)
Get file name part, i.e. part after last JEEP::PATHNAME_SEPARATOR if any.
Definition: JeepToolkit.hh:88
bool JPHYSICS::hasWildCard ( const std::string &  file_name)
inline

Check wild card.

Parameters
file_namefile name
Returns
true if wild card present; else false

Definition at line 110 of file JPDFTypes.hh.

111  {
112  return (file_name.find(WILD_CARD) != std::string::npos);
113  }
static const char WILD_CARD
Definition: JDAQTags.hh:34
std::string JPHYSICS::getFilename ( const std::string &  file_name,
const JPDFType_t  pdf 
)
inline

Get PDF file name.

The input file name should contain the wild card character WILD_CARD which will be replaced by the label corresponding to the given PDF type.

Parameters
file_nameinput file name
pdfPDF type
Returns
output file name

Definition at line 126 of file JPDFTypes.hh.

128  {
129  using namespace std;
130 
131  string buffer = file_name;
132 
133  string::size_type pos = buffer.find(WILD_CARD);
134 
135  if (pos == string::npos) {
136  throw JException(string("Method getFilename(): Missing wild card character \'") + WILD_CARD + "\'.");
137  }
138 
139  return buffer.replace(pos, 1, getLabel(pdf));
140  }
std::string getLabel(const JLocation &location)
Get module label for monitoring and other applications.
Definition: JLocation.hh:246
static const char WILD_CARD
Definition: JDAQTags.hh:34
bool JPHYSICS::is_muon ( const int  pdf)
inline

Test if given PDF type corresponds to Cherenkov light from muon.

Parameters
pdfPDF type
Returns
true if PDF corresponds to muon; else false

Definition at line 149 of file JPDFTypes.hh.

150  {
151  return (pdf == DIRECT_LIGHT_FROM_MUON ||
152  pdf == SCATTERED_LIGHT_FROM_MUON ||
153  pdf == LIGHT_FROM_MUON);
154  }
direct and scattered light from muon
Definition: JPDFTypes.hh:48
direct light from muon
Definition: JPDFTypes.hh:29
scattered light from muon
Definition: JPDFTypes.hh:30
bool JPHYSICS::is_bremsstrahlung ( const int  pdf)
inline

Test if given PDF type corresponds to Cherenkov light from Bremsstrahlung.

Parameters
pdfPDF type
Returns
true if PDF corresponds to Bremsstrahlung; else false

Definition at line 163 of file JPDFTypes.hh.

164  {
165  return (pdf == DIRECT_LIGHT_FROM_EMSHOWERS ||
167  pdf == LIGHT_FROM_EMSHOWERS);
168  }
direct light from EM showers
Definition: JPDFTypes.hh:32
scattered light from EM showers
Definition: JPDFTypes.hh:33
direct and scattered light from EM showers
Definition: JPDFTypes.hh:49
bool JPHYSICS::is_deltarays ( const int  pdf)
inline

Test if given PDF type corresponds to Cherenkov light from delta-rays.

Parameters
pdfPDF type
Returns
true if PDF corresponds to delta-rays; else false

Definition at line 177 of file JPDFTypes.hh.

178  {
179  return (pdf == DIRECT_LIGHT_FROM_DELTARAYS ||
181  pdf == LIGHT_FROM_DELTARAYS);
182  }
scattered light from delta-rays
Definition: JPDFTypes.hh:36
direct light from delta-rays
Definition: JPDFTypes.hh:35
direct and scattered light from delta-rays
Definition: JPDFTypes.hh:50
bool JPHYSICS::is_scattered ( const int  pdf)
inline

Test if given PDF type corresponds to scattered light.

Parameters
pdfPDF type
Returns
true if PDF corresponds to scattered light; else false

Definition at line 191 of file JPDFTypes.hh.

192  {
193  return (pdf == SCATTERED_LIGHT_FROM_MUON ||
198  }
scattered light from muon
Definition: JPDFTypes.hh:38
scattered light from EM shower
Definition: JPDFTypes.hh:41
scattered light from muon
Definition: JPDFTypes.hh:30
scattered light from delta-rays
Definition: JPDFTypes.hh:36
scattered light from EM showers
Definition: JPDFTypes.hh:33
double JPHYSICS::getThetaMCS ( const double  E,
const double  x,
const double  X0,
const double  M,
const double  Q 
)
inline

Get multiple Coulomb scattering angle.

The formula is taken from reference: Particle Data Book, formula 27.14.

Parameters
EEnergy [GeV]
xdistance [m]
X0radiation length [m]
Mmass [GeV]
Qcharge [unit]
Returns
angle [rad]

Definition at line 32 of file JPhysicsToolkit.hh.

37  {
38  using JTOOLS::THETA_MCS;
39 
40  if (E > M) {
41 
42  const double p = sqrt((E + M) * (E - M));
43  const double beta = p / E;
44 
45  return THETA_MCS * Q * sqrt(x/X0) * (1.0 + 0.038*log(x/X0)) / (beta*p);
46  }
47 
48  return 0.0;
49  }
static const double THETA_MCS
Multiple Coulomb scattering constant [GeV].
Definition: JConstants.hh:29
do $JPP JMEstimator M
Definition: JMEstimator.sh:37
then usage $script[input file[working directory[option]]] nWhere option can be E
Definition: JMuonPostfit.sh:37
double JPHYSICS::getThetaMCS ( const double  E,
const double  x 
)
inline

Get multiple Coulomb scattering angle.

Parameters
EEnergy [GeV]
xdistance [m]
Returns
angle [rad]

Definition at line 59 of file JPhysicsToolkit.hh.

60  {
61  using namespace JTOOLS;
62 
63  return getThetaMCS(E, x, X0_WATER_M, MASS_MUON, 1.0);
64  }
static const double MASS_MUON
muon mass [GeV]
Definition: JConstants.hh:59
static const double X0_WATER_M
Radiation length pure water [m].
Definition: JConstants.hh:38
double getThetaMCS(const double E, const double x, const double X0, const double M, const double Q)
Get multiple Coulomb scattering angle.
then usage $script[input file[working directory[option]]] nWhere option can be E
Definition: JMuonPostfit.sh:37

Variable Documentation

const JGeaneWater JPHYSICS::gWater
static

Function object for energy loss of muon in sea water.

Definition at line 328 of file JGeane.hh.

const JGeant JPHYSICS::geant(geanx, 0.0001)
static

Function object for the number of photons from EM-shower as a function of emission angle.

const JGeanx JPHYSICS::geanx(0.35,-5.40)
static

Function object for the number of photons from EM-shower as a function of emission angle.

const JGeanz JPHYSICS::geanz(1.85, 0.62, 0.54)
static

Function object for longitudinal EM-shower profile.

const JCCnu JPHYSICS::cc_nu
static

Function object for charged current neutrino cross section [cm^2] as a function of neutrino energy [GeV].

Definition at line 618 of file JNeutrino.hh.

const JNCnu JPHYSICS::nc_nu
static

Function object for neutral current neutrino cross section [cm^2] as a function of neutrino energy [GeV].

Definition at line 620 of file JNeutrino.hh.

const JCCnubar JPHYSICS::cc_nubar
static

Function object for charged current anti-neutrino cross section [cm^2] as a function of neutrino energy [GeV].

Definition at line 621 of file JNeutrino.hh.

const JNCnubar JPHYSICS::nc_nubar
static

Function object for neutral current anti-neutrino cross section [cm^2] as a function of neutrino energy [GeV].

Definition at line 622 of file JNeutrino.hh.

double JPHYSICS::MODULE_RADIUS_M = 0.25
static

Radius of optical module [m].

This parameter is used to implement shadowing of the PMT by the optical module.

Definition at line 42 of file JPDF.hh.

const JAbsorptionLengthOfPureWater JPHYSICS::getAbsorptionLengthOfPureWater
static

Function object for absorption length of pure water.

Definition at line 325 of file JPDFToolkit.hh.

const char JPHYSICS::WILD_CARD = '%'
static

wild card character for file name substitution

Definition at line 56 of file JPDFTypes.hh.