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
	Function object for the 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 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
	Template definition of transformer of the Probability Density Functions of the time response of a PMT. More...

class	JPDFTransformer< 1, JArgument_t >
	Template specialisation of transformer of the 1D Probability Density Functions of the time response of a PMT. More...

class	JPDFTransformer< 2, JArgument_t >
	Template specialisation of transformer of the 2D Probability Density Functions of the time response of a PMT. More...

class	JPDFTransformer< 3, JArgument_t >
	Template specialisation of transformer of the 3D Probability Density Functions of the time response of a PMT. More...

class	JPDFTransformer< 4, JArgument_t >
	Template specialisation of transformer of the 4D Probability Density Functions of the time response of a PMT. More...

class	JPDFTransformer< 5, JArgument_t >
	Template specialisation of transformer of the 5D Probability Density Functions of the time response of a PMT. 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::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, LIGHT_FROM_MUON = 1001, LIGHT_FROM_EMSHOWERS = 1003, LIGHT_FROM_DELTARAYS = 1005, LIGHT_FROM_EMSHOWER = 1013 }
	PDF types. More...

Functions
double	geanc ()
	Equivalent muon track length per unit shower energy. 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 JGeane	gWater (2.67e-1 JTOOLS::DENSITY_SEA_WATER, 3.4e-4 JTOOLS::DENSITY_SEA_WATER)
	Function object for Energy loss of muon in sea water. More...

static const JGeane	gRock (2.67e-1 0.9 JTOOLS::DENSITY_ROCK, 3.4e-4 1.2 JTOOLS::DENSITY_ROCK)
	Function object for Energy loss of muon in rock. 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

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

enum JPHYSICS::JPDFType_t

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
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

Definition at line 27 of file JPDFTypes.hh.

                   {
 
     DIRECT_LIGHT_FROM_MUON         =    1,      //!< direct    light from muon
     SCATTERED_LIGHT_FROM_MUON      =    2,      //!< scattered light from muon
 
     DIRECT_LIGHT_FROM_EMSHOWERS    =    3,      //!< direct    light from EM showers
     SCATTERED_LIGHT_FROM_EMSHOWERS =    4,      //!< scattered light from EM showers
 
     DIRECT_LIGHT_FROM_DELTARAYS    =    5,      //!< direct    light from delta-rays
     SCATTERED_LIGHT_FROM_DELTARAYS =    6,      //!< scattered light from delta-rays
 
     SCATTERED_LIGHT_FROM_MUON_5D   =   12,      //!< scattered light from muon
 
     DIRECT_LIGHT_FROM_EMSHOWER     =   13,      //!< direct    light from EM shower
     SCATTERED_LIGHT_FROM_EMSHOWER  =   14,      //!< scattered light from EM shower
 
     LIGHT_FROM_MUON                = 1001,      //!< direct and scattered light from muon
     LIGHT_FROM_EMSHOWERS           = 1003,      //!< direct and scattered light from EM showers
     LIGHT_FROM_DELTARAYS           = 1005,      //!< direct and scattered light from delta-rays
     LIGHT_FROM_EMSHOWER            = 1013       //!< direct and scattered light from EM shower
   };

Function Documentation

double JPHYSICS::geanc ( )

inline

Equivalent muon track length per unit shower energy.

Returns: equivalent muon track length [m/Gev]

Definition at line 23 of file JGeane.hh.

   {
     return 4.7;                   // dx/dE [m/GeV]
   }

double JPHYSICS::getMinimalWavelength ( )

inline

Get minimal wavelength for PDF evaluations.

Returns: wavelength of light [nm]

Definition at line 38 of file JPDFToolkit.hh.

   {
     return 300.0;
   }

double JPHYSICS::getMaximalWavelength ( )

inline

Get maximal wavelength for PDF evaluations.

Returns: wavelength of light [nm]

Definition at line 49 of file JPDFToolkit.hh.

   {
     return 700.0;
   }

double JPHYSICS::cherenkov	(	const double	lambda,
		const double	n
	)

inline

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

Parameters

lambda	wavelength of light [nm]
n	index of refraction

Returns: number of photons per unit track length and per unit wavelength [m^-1 nm^-1]

Definition at line 62 of file JPDFToolkit.hh.

   {
     using namespace JTOOLS;
 
     const double x = n*lambda;
     
     return 1.0e9 * 2 * PI * ALPHA_ELECTRO_MAGNETIC * (n*n - 1.0) / (x*x);
   }

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

E	muon energy [GeV]

Returns: equivalent energy loss [GeV/m]

Definition at line 81 of file JPDFToolkit.hh.

   {
     static const double a =  3.186e-01;
     static const double b =  3.384e-01;
     static const double c = -2.759e-02;
     static const double d =  1.630e-03;
     static const double Emin = 0.13078; // [GeV]
 
     if (E > Emin) {
 
       const double x = log10(E);                      //
       const double y = a + x*(b + x*(c + x*(d)));     // [MeV g^-1 cm^2]
 
       return y * JTOOLS::DENSITY_SEA_WATER * 1.0e-1;  // [GeV/m]
     }
 
     return 0.0;
   }

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

E	tau energy [GeV]

Returns: equivalent energy loss [GeV/m]

Definition at line 109 of file JPDFToolkit.hh.

   {
     static const double a = -2.374e-01;
     static const double b =  5.143e-01;
     static const double c = -4.213e-02;
     static const double d =  1.804e-03;
     static const double Emin = 2.19500; // [GeV]
     
     if (E > Emin) {
 
       const double x = log10(E);                      //
       const double y = a + x*(b + x*(c + x*(d)));     // [MeV g^-1 cm^2]
 
       return y * JTOOLS::DENSITY_SEA_WATER * 1.0e-1;  // [GeV/m]
     }
 
     return 0.0;
   }

const double JPHYSICS::getRayleighCrossSection	(	const double	n,
		const double	lambda
	)

inline

Rayleigh cross section.

Parameters

n	index of refraction
lambda	wavelength of light [nm]

Returns: cross section [cm^2]

Definition at line 136 of file JPDFToolkit.hh.

   {
     using JTOOLS::PI;
 
     static const double d = 0.36;                    // size of H2O molecule [nm]
     static const double U = PI*PI*PI*PI*PI*2.0/3.0;
     static const double V = d*d*d*d*d*d;
 
     const double W     = (n*n - 1.0) / (n*n + 2.0);
     const double sigma = 1.0e-14 * U*V*W*W / (lambda*lambda*lambda*lambda);   // [cm^2]
 
     return sigma;
   }

const double JPHYSICS::getRayleighScatteringLength	(	const double	n,
		const double	lambda
	)

inline

Rayleigh scattering length.

Parameters

n	index of refraction
lambda	wavelength of light [nm]

Returns: scattering length [m]

Definition at line 159 of file JPDFToolkit.hh.

   {
     using namespace JTOOLS;
 
     static const double amu = 18.01528; // H20 mass in Atomic mass units
 
     const double sigma = getRayleighCrossSection(n, lambda);
     const double ls    = 1.0e-2 / (DENSITY_SEA_WATER * AVOGADRO * sigma / amu);     // [m]
 
     return ls;
   }

std::string JPHYSICS::getLabel ( const JPDFType_t pdf )

inline

Get PDF label.

Parameters

pdf PDF type

Returns: PDF label

Definition at line 59 of file JPDFTypes.hh.

   {
     std::ostringstream os;
 
     os << pdf;
 
     return os.str();
   }

int JPHYSICS::getPDFType ( const std::string & file_name )

inline

Get PDF type.

Parameters

file_name file name

Returns: PDF type (-1 in case of error)

Definition at line 75 of file JPDFTypes.hh.

   {
     using namespace std;
 
     static const char* digits = "0123456789";
     
     int type = -1;
     
     string buffer = JEEP::getFilename(file_name);
     
     string::size_type pos = buffer.find_first_of(digits);
 
     if (pos != string::npos) {
 
       string::size_type len = buffer.substr(pos).find_first_not_of(digits);
 
       istringstream(buffer.substr(pos, len)) >> type;
     }
 
     return type;
   }

bool JPHYSICS::hasWildCard ( const std::string & file_name )

inline

Check wild card.

Parameters

file_name file name

Returns: true if wild card present; else false

Definition at line 104 of file JPDFTypes.hh.

   {
     return (file_name.find(WILD_CARD) != std::string::npos);
   }

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_name	input file name
pdf	PDF type

Returns: output file name

Definition at line 120 of file JPDFTypes.hh.

   {
     using namespace std;
 
     string buffer = file_name;
 
     string::size_type pos = buffer.find(WILD_CARD);
     
     if (pos == string::npos) {    
       throw JException(string("Method getFilename(): Missing wild card character \'") + WILD_CARD + "\'.");
     }
 
     return buffer.replace(pos, 1, getLabel(pdf));
   }

bool JPHYSICS::is_muon ( const int pdf )

inline

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

Parameters

pdf PDF type

Returns: true if PDF corresponds to muon; else false

Definition at line 143 of file JPDFTypes.hh.

   {
     return (pdf == DIRECT_LIGHT_FROM_MUON           ||
             pdf == SCATTERED_LIGHT_FROM_MUON        ||
             pdf == LIGHT_FROM_MUON);
   }

bool JPHYSICS::is_bremsstrahlung ( const int pdf )

inline

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

Parameters

pdf PDF type

Returns: true if PDF corresponds to Bremsstrahlung; else false

Definition at line 157 of file JPDFTypes.hh.

   {
     return (pdf == DIRECT_LIGHT_FROM_EMSHOWERS      ||
             pdf == SCATTERED_LIGHT_FROM_EMSHOWERS   ||
             pdf == LIGHT_FROM_EMSHOWERS);
   }

bool JPHYSICS::is_deltarays ( const int pdf )

inline

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

Parameters

pdf PDF type

Returns: true if PDF corresponds to delta-rays; else false

Definition at line 171 of file JPDFTypes.hh.

   {
     return (pdf == DIRECT_LIGHT_FROM_DELTARAYS      ||
             pdf == SCATTERED_LIGHT_FROM_DELTARAYS   ||
             pdf == LIGHT_FROM_DELTARAYS);
   }

bool JPHYSICS::is_scattered ( const int pdf )

inline

Test if given PDF type corresponds to scattered light.

Parameters

pdf PDF type

Returns: true if PDF corresponds to scattered light; else false

Definition at line 185 of file JPDFTypes.hh.

   {
     return (pdf == SCATTERED_LIGHT_FROM_MUON        ||
             pdf == SCATTERED_LIGHT_FROM_EMSHOWERS   ||
             pdf == SCATTERED_LIGHT_FROM_DELTARAYS   ||
             pdf == SCATTERED_LIGHT_FROM_MUON_5D     ||
             pdf == SCATTERED_LIGHT_FROM_EMSHOWER);
   }

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

E	Energy [GeV]
x	distance [m]
X0	radiation length [m]
M	mass [GeV]
Q	charge [unit]

Returns: angle [rad]

Definition at line 32 of file JPhysicsToolkit.hh.

   {
     using JTOOLS::THETA_MCS;
     
     if (E > M) {
 
       const double p    = sqrt((E + M) * (E - M));
       const double beta = p / E;
       
       return THETA_MCS * Q * sqrt(x/X0) * (1.0 + 0.038*log(x/X0)) / (beta*p);
     }
 
     return 0.0;
   }

double JPHYSICS::getThetaMCS	(	const double	E,
		const double	x
	)

inline

Get multiple Coulomb scattering angle.

Parameters

E	Energy [GeV]
x	distance [m]

Returns: angle [rad]

Definition at line 59 of file JPhysicsToolkit.hh.

   {
     using namespace JTOOLS;
 
     return getThetaMCS(E, x, X0_WATER_M, MASS_MUON, 1.0);
   }

Variable Documentation

const JGeane JPHYSICS::gWater(2.67e-1 *JTOOLS::DENSITY_SEA_WATER, 3.4e-4 *JTOOLS::DENSITY_SEA_WATER)

static

Function object for Energy loss of muon in sea water.

Approximate values of energy loss parameters taken from reference: Proceedings of ICRC 2001, "Precise parametrizations of muon energy losses in water", S. Klimushin, E. Bugaev and I. Sokalski.

const JGeane JPHYSICS::gRock(2.67e-1 *0.9 *JTOOLS::DENSITY_ROCK, 3.4e-4 *1.2 *JTOOLS::DENSITY_ROCK)

static

Function object for Energy loss of muon in rock.

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 39 of file JPDF.hh.

const JAbsorptionLengthOfPureWater JPHYSICS::getAbsorptionLengthOfPureWater

static

Function object for absorption length of pure water.

Definition at line 329 of file JPDFToolkit.hh.

const char JPHYSICS::WILD_CARD = '%'

static

wild card character for file name substitution

Definition at line 50 of file JPDFTypes.hh.

Classes

Typedefs

Enumerations

Functions

Variables

Detailed Description

Typedef Documentation

Enumeration Type Documentation

Function Documentation

Variable Documentation