Auxiliary methods for light properties of deep-sea water. More...

Classes
class	JAbsorptionLengthOfPureWater
	Absorption length of pure water. More...

class	JAbstractLED
	Interface for emission profile from LED. More...

class	JAbstractMedium
	Medium interface. More...

class	JAbstractPDF
	Probability Density Functions of the time response of a PMT with an implementation for the JDispersionInterface interface. More...

class	JAbstractPMT
	PMT interface. More...

class	JACoeffInterface
	Interface for calculation of ionization constant. More...

class	JACoeffSource
	Implementation for calculation of ionization constant. More...

class	JCDFTable
	Multi-dimensional CDF table for arrival time of Cherenkov light. More...

class	JCDFTable1D
	Custom class for CDF table in 1 dimension. More...

class	JCDFTable2D
	Custom class for CDF table in 2 dimensions. More...

class	JDIS
	Deep-inelastic muon-nucleon scattering. More...

class	JDispersion
	Implementation of dispersion for water in deep sea. More...

class	JDispersionInterface
	Light dispersion inteface. More...

class	JDISSource
	Implementation for calculation of inverse interaction length and shower energy due to deep-inelastic muon-nucleon scattering. 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...

struct	JK40Rates
	Auxiliary class for K40 rates. 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...

struct	JMobley
	Auxiliary data structure for scattering lengths of deep-sea water. More...

class	JNPETable
	Custom class for integrated values of the PDF of the arrival time of Cherenkov light. 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	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	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	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	JPDFTransformer_t
	Transformer for the 1D probability density function (PDF) of the time response of a PMT to a muon. 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	JPetzhold
	Measurement of light scattering in water. 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
	Implementation for calculation of inverse interaction length and shower energy. More...

struct	JRadiationSource_t
	Auxiliary data structure for handling member methods of class JRadiation. More...

struct	JSter
	Struct for the Sternheimer coefficients. More...

struct	JSterCoefficient
	Auxiliary data structure Ster coefficients as a function of proton number. More...

Typedefs
typedef JTOOLS::JGridPolint1Function1D_t	JGeantFunction1D_t

typedef double	JRateL0_t
	Type definition of singles rate [Hz].

typedef std::vector< double >	JRateL1_t
	Type definition of count rate as a function of multiplicty [Hz] The multiples rate start counting at two-fold coincidences.

typedef size_t	multiplicity_type
	Type definition of multiplicity.

typedef JTOOLS::JElement2D< double, double >	JElement2D_t

typedef JTOOLS::JElement3D< double, double >	JElement3D_t

typedef JTOOLS::JRange< double >	JTimeRange
	Type definition for time range (unit [ns]).

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_ELONGATED_EMSHOWER = 113 , 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
const double	getSpeedOfLight ()
	Get speed of light.

const double	getInverseSpeedOfLight ()
	Get inverse speed of light.

double	getIndexOfRefraction ()
	Get average index of refraction of water corresponding to group velocity.

double	getIndexOfRefractionPhase ()
	Get average index of refraction of water corresponding to phase velocity.

double	getTanThetaC ()
	Get average tangent of Cherenkov angle of water corresponding to group velocity.

double	getCosThetaC ()
	Get average cosine of Cherenkov angle of water corresponding to group velocity.

double	getSinThetaC ()
	Get average sine of Cherenkov angle of water corresponding to group velocity.

double	getKappaC ()
	Get average R-dependence of arrival time of Cherenkov light (a.k.a.

double	geanc ()
	Equivalent muon track length per unit shower energy.

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.

double	getMinimalWavelength ()
	Get minimal wavelength for PDF evaluations.

double	getMaximalWavelength ()
	Get maximal wavelength for PDF evaluations.

double	cherenkov (const double lambda, const double n)
	Number of Cherenkov photons per unit track length and per unit wavelength.

double	getDeltaRayTmin ()
	Get minimum delta-ray kinetic energy.

double	getDeltaRayTmax (const double E, const double M)
	Get maximum delta-ray kinetic energy for given lepton energy and mass.

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 given energy and given mass.

template<class JFormFactor_t >
double	getDeltaRays (const double E, const double M, const double Tmin, const double Tmax, const double Z, const double A, const JFormFactor_t &F, const int N=1000000)
	Get equivalent EM-shower energy due to delta-rays per unit track length for an ionising particle with given energy and given mass and for a given form factor.

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.

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	getDeltaRaysFromMuonFit (const double E)
	Equivalent EM-shower energy due to delta-rays per unit muon track length.

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.

double	getDeltaRaysFromTauFit (const double E)
	Equivalent EM-shower energy due to delta-rays per unit tau track length.

double	getDeltaRayProbability (const double x)
	Emission profile of photons from delta-rays.

const double	getRayleighCrossSection (const double n, const double lambda)
	Rayleigh cross section.

const double	getRayleighScatteringLength (const double n, const double lambda)
	Rayleigh scattering length.

std::string	getLabel (const JPDFType_t pdf)
	Get PDF label.

int	getPDFType (const std::string &file_name)
	Get PDF type.

std::string	getFilename (const std::string &file_name, const JPDFType_t pdf)
	Get PDF file name.

bool	is_muon (const int pdf)
	Test if given PDF type corresponds to Cherenkov light from muon.

bool	is_bremsstrahlung (const int pdf)
	Test if given PDF type corresponds to Cherenkov light from Bremsstrahlung.

bool	is_deltarays (const int pdf)
	Test if given PDF type corresponds to Cherenkov light from delta-rays.

bool	is_scattered (const int pdf)
	Test if given PDF type corresponds to scattered light.

double	getThetaMCS (const double E, const double x, const double X0, const double M, const double Q)
	Get multiple Coulomb scattering angle.

double	getThetaMCS (const double E, const double x)
	Get multiple Coulomb scattering angle for muon.

double	henyey_greenstein (const double g, const double x)
	Auxiliary method to describe light scattering in water (Henyey-Greenstein).

double	henyey_greenstein (const double x)
	Auxiliary method to describe light scattering in water (Heneyey-Greenstein).

double	rayleigh (const double a, const double x)
	Auxiliary method to describe light scattering in water (Rayleigh).

double	rayleigh (const double x)
	Auxiliary method to describe light scattering in water (Rayleigh).

double	f4 (const double x)
	Model specific function to describe light scattering in water (f4).

double	p00075 (const double x)
	Model specific function to describe light scattering in water (p00075).

double	getKineticEnergy (const double E, const double m)
	Get kinetic energy of particle with given energy and mass.

Variables
static const double	C = 0.299792458
	Physics constants.

static const double	C_INVERSE = 1.0/C
	Inverse speed of light in vacuum [ns/m].

static const double	AVOGADRO = 6.0221415e23
	Avogadro's number.

static const double	NUCLEON_MOLAR_MASS = 1.0
	nucleon molar mass [g/mol]

static const double	H = 4.13566733e-15
	Planck constant [eV s].

static const double	HBAR = H/(2*PI)
	Planck constant [eV s].

static const double	HBARC = HBARC1.0e9
	Planck constant [eV m].

static const double	ALPHA_ELECTRO_MAGNETIC = 1.0/137.036
	Electro-Magnetic coupling constant.

static const double	THETA_MCS = 13.6e-3
	Multiple Coulomb scattering constant [GeV].

static const double	R_EARTH_KM = 6371
	Geophysics constants.

static const double	DENSITY_EARTH = 5.51
	Average density of the Earth [gr/cm³].

static const double	DENSITY_SEA_WATER = 1.038
	Fixed environment values.

static const double	DENSITY_ROCK = 2.65
	Density of rock [g/cm^3].

static const double	SALINITY_SEA_WATER = 0.035
	Salinity of sea water.

static const double	INDEX_OF_REFRACTION_WATER = 1.3800851282
	Average index of refraction of water corresponding to the group velocity.

static const double	X0_WATER_M = 0.36
	Radiation length pure water [m].

static const double	TAN_THETA_C_WATER = sqrt((INDEX_OF_REFRACTION_WATER - 1.0) * (INDEX_OF_REFRACTION_WATER + 1.0))
	Derived quantities of optical medium.

static const double	COS_THETA_C_WATER = 1.0 / INDEX_OF_REFRACTION_WATER
	Average cosine corresponding to the group velocity.

static const double	SIN_THETA_C_WATER = TAN_THETA_C_WATER * COS_THETA_C_WATER
	Average sine corresponding to the group velocity.

static const double	KAPPA_WATER = 0.96
	Average R-dependence of arrival time of Cherenkov light.

static const double	MASS_PHOTON = 0.0
	Particle masses.

static const double	MASS_ELECTRON_NEUTRINO = 0.0
	electron neutrino mass [GeV]

static const double	MASS_MUON_NEUTRINO = 0.0
	muon neutrino mass [GeV]

static const double	MASS_TAU_NEUTRINO = 0.0
	tau neutrino mass [GeV]

static const double	MASS_ELECTRON = 0.510998946e-3
	electron mass [GeV]

static const double	MASS_MUON = 0.1056583745
	muon mass [GeV]

static const double	MASS_TAU = 1.77682
	tau mass [GeV]

static const double	MASS_NEUTRAL_PION = 0.1349766
	pi_0 mass [GeV]

static const double	MASS_CHARGED_PION = 0.13957018
	pi^+/- mass [GeV]

static const double	MASS_NEUTRAL_KAON = 0.497614
	K_0 mass [GeV].

static const double	MASS_CHARGED_KAON = 0.493677
	K^+/- mass [GeV].

static const double	MASS_NEUTRAL_RHO = 0.77526
	rho_0 mass [GeV]

static const double	MASS_CHARGED_RHO = 0.77511
	rho^+/- mass [GeV]

static const double	MASS_NEUTRAL_D = 1.86483
	D_0 mass [GeV].

static const double	MASS_CHARGED_D = 1.86965
	D^+/- mass [GeV].

static const double	MASS_CHARGED_D_S = 1.96834
	D_s^+/- mass [GeV].

static const double	MASS_PROTON = 0.9382720813
	proton mass [GeV]

static const double	MASS_NEUTRON = 0.9395654133
	neutron mass [GeV]

static const double	MASS_DELTA_1232 = 1.232
	Delta (1232) mass [GeV].

static const double	MASS_LAMBDA = 1.115683
	Lambda mass [GeV].

static const double	MASS_NEUTRAL_SIGMA = 1.192642
	Sigma_0 mass [GeV].

static const double	MASS_CHARGED_SIGMA = 1.18937
	Sigma^+/- mass [GeV].

static const double	MASS_NEUTRAL_XI = 1.31486
	Xi_0 mass [GeV].

static const double	MASS_CHARGED_XI = 1.32171
	Xi^+/- mass [GeV].

static const double	MASS_CHARGED_OMEGA = 1.67245
	Omega^+/- mass [GeV].

static const double	MASS_CHARGED_LAMBDA_C = 2.28646
	Lambda_c^+/- mass [GeV].

static const double	MASS_DOUBLYCHARGED_SIGMA_C = 2.45397
	Sigma_c^++/– mass [GeV]

static const double	MASS_CHARGED_SIGMA_C = 2.4529
	Sigma_c^+/- mass [GeV].

static const double	MASS_NEUTRAL_SIGMA_C = 2.45375
	Sigma_c_0 mass [GeV].

static const double	MASS_CHARGED_XI_C = 2.46793
	Xi_c^+/- mass [GeV].

static const double	MASS_NEUTRAL_XI_C = 2.47091
	Xi_c_0 mass [GeV].

static const double	MASS_NEUTRAL_OMEGA_C = 2.6952
	Omega_c_0 mass [GeV].

static const double	MASS_NEUTRAL_B = 5.27958
	B_0 mass [GeV].

static const double	MASS_CHARGED_B = 5.27926
	B^+/- mass [GeV].

static const double	MASS_NEUTRAL_B_S = 5.36677
	B_s^0 mass [GeV].

static const double	MASS_NEUTRAL_LAMBDA_B = 5.6194
	Lambda_b^0 mass [GeV].

static const double	MASS_NEUTRAL_XI_B = 5.7878
	Xi_b^0 mass [GeV].

static const double	MASS_CHARGED_XI_B = 5.7911
	Xi_b^+/- mass [GeV].

static const double	MASS_CHARGED_OMEGA_B = 6.071
	Omega_b^+/- mass [GeV].

static const double	MASS_CHARGED_B_C = 6.2756
	B_c^+/- mass [GeV].

static const JGeaneWater	gWater
	Function object for energy loss of muon in sea water.

static const JGeant	geant (geanx, 0.0001)
	Function object for the number of photons from EM-shower as a function of emission angle.

static const JGeanx	geanx (0.35, -5.40)
	Function object for the number of photons from EM-shower as a function of emission angle.

static const JGeanz	geanz (1.85, 0.62, 0.54)
	Function object for longitudinal EM-shower profile.

static JSterCoefficient	getSterCoefficient
	Function object for Ster coefficients.

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

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

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

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

static double	MODULE_RADIUS_M = 0.216
	Radius of optical module [m].

static const double	DELTARAY_TMIN = 0.000915499
	Minimum allowed delta-ray kinetic energy [GeV].

static const double	DELTARAY_TMAX = 1.0e10
	Maximum allowed delta-ray kinetic energy [GeV].

static const JPetzhold	petzhold
	Function object for measurement of light scattering in water.

static const JAbsorptionLengthOfPureWater	getAbsorptionLengthOfPureWater
	Function object for absorption length of pure water.

static const JRadiationSource_t	EErad_t = { &JRadiation::TotalCrossSectionEErad, &JRadiation::EfromEErad, &JRadiation::ThetaRMSfromEErad }

static const JRadiationSource_t	Brems_t = { &JRadiation::TotalCrossSectionBrems, &JRadiation::EfromBrems, &JRadiation::ThetaRMSfromBrems }

static const JRadiationSource_t	GNrad_t = { &JRadiation::TotalCrossSectionGNrad, &JRadiation::EfromGNrad, &JRadiation::ThetaRMSfromGNrad }

Detailed Description

Auxiliary methods for light properties of deep-sea water.

Author: mdejong; vcarretero

Typedef Documentation

◆ JGeantFunction1D_t

typedef JTOOLS::JGridPolint1Function1D_t JPHYSICS::JGeantFunction1D_t

Definition at line 21 of file JGeant_t.hh.

◆ JRateL0_t

typedef double JPHYSICS::JRateL0_t

Type definition of singles rate [Hz].

Definition at line 21 of file JK40Rates.hh.

◆ JRateL1_t

typedef std::vector<double> JPHYSICS::JRateL1_t

Type definition of count rate as a function of multiplicty [Hz] The multiples rate start counting at two-fold coincidences.

Definition at line 27 of file JK40Rates.hh.

◆ multiplicity_type

typedef size_t JPHYSICS::multiplicity_type

Type definition of multiplicity.

Definition at line 33 of file JK40Rates.hh.

◆ JElement2D_t

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

Definition at line 27 of file JLED.hh.

◆ JElement3D_t

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

Definition at line 28 of file JLED.hh.

◆ JTimeRange

typedef JTOOLS::JRange<double> JPHYSICS::JTimeRange

Type definition for time range (unit [ns]).

Definition at line 19 of file JPhysics/JTimeRange.hh.

Enumeration Type Documentation

◆ JPDFType_t

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
DIRECT_LIGHT_FROM_BRIGHT_POINT	direct light from bright point
SCATTERED_LIGHT_FROM_BRIGHT_POINT	scattered light from bright point
LIGHT_FROM_ELONGATED_EMSHOWER	light from elongated 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
LIGHT_FROM_BRIGTH_POINT	direct and scattered light from brigth point

Definition at line 24 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
 
    //DIRECT_LIGHT_FROM_SHOWER           =   15,      //!< direct + scattered light from HADRONIC shower
 
    DIRECT_LIGHT_FROM_BRIGHT_POINT     =   23,      //!< direct    light from bright point
    SCATTERED_LIGHT_FROM_BRIGHT_POINT  =   24,      //!< scattered light from bright point
 
    LIGHT_FROM_ELONGATED_EMSHOWER      =  113,      //!< light from elongated 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
    LIGHT_FROM_BRIGTH_POINT            = 1023       //!< direct and scattered light from brigth point
  };

Function Documentation

◆ getSpeedOfLight()

const double JPHYSICS::getSpeedOfLight ( )

inline

Get speed of light.

Returns: speed of light [m/ns]

Definition at line 112 of file JPhysics/JConstants.hh.

  {
    return C;
  }

◆ getInverseSpeedOfLight()

const double JPHYSICS::getInverseSpeedOfLight ( )

inline

Get inverse speed of light.

Returns: inverse speed of light [ns/m]

Definition at line 123 of file JPhysics/JConstants.hh.

  { 
    return C_INVERSE; 
  }

◆ getIndexOfRefraction()

double JPHYSICS::getIndexOfRefraction ( )

inline

Get average index of refraction of water corresponding to group velocity.

Returns: index of refraction

Definition at line 134 of file JPhysics/JConstants.hh.

  {
    return INDEX_OF_REFRACTION_WATER;
  }

◆ getIndexOfRefractionPhase()

double JPHYSICS::getIndexOfRefractionPhase ( )

inline

Get average index of refraction of water corresponding to phase velocity.

Returns: index of refraction

Definition at line 145 of file JPhysics/JConstants.hh.

  {
    return 1.35;
  }

◆ getTanThetaC()

double JPHYSICS::getTanThetaC ( )

inline

Get average tangent of Cherenkov angle of water corresponding to group velocity.

Returns: tan(theta_C)

Definition at line 156 of file JPhysics/JConstants.hh.

  {
    return TAN_THETA_C_WATER;
  }

◆ getCosThetaC()

double JPHYSICS::getCosThetaC ( )

inline

Get average cosine of Cherenkov angle of water corresponding to group velocity.

Returns: cos(theta_C)

Definition at line 167 of file JPhysics/JConstants.hh.

  {
    return COS_THETA_C_WATER;
  }

◆ getSinThetaC()

double JPHYSICS::getSinThetaC ( )

inline

Get average sine of Cherenkov angle of water corresponding to group velocity.

Returns: sin(theta_C)

Definition at line 178 of file JPhysics/JConstants.hh.

  {
    return SIN_THETA_C_WATER;
  }

◆ getKappaC()

double JPHYSICS::getKappaC ( )

inline

Get average R-dependence of arrival time of Cherenkov light (a.k.a.

kappa).

Returns: kappa

Definition at line 189 of file JPhysics/JConstants.hh.

  {
    return KAPPA_WATER;
  }

◆ geanc()

double JPHYSICS::geanc ( )

inline

Equivalent muon track length per unit shower energy.

See ANTARES internal note ANTARES-SOFT-2002-015, J. Brunner.

Returns: equivalent muon track length [m/Gev]

Definition at line 28 of file JGeane.hh.

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

◆ gRock()

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.

◆ getMinimalWavelength()

double JPHYSICS::getMinimalWavelength ( )

inline

Get minimal wavelength for PDF evaluations.

Returns: wavelength of light [nm]

Definition at line 35 of file JPDFToolkit.hh.

  {
    return 300.0;
  }

◆ getMaximalWavelength()

double JPHYSICS::getMaximalWavelength ( )

inline

Get maximal wavelength for PDF evaluations.

Returns: wavelength of light [nm]

Definition at line 46 of file JPDFToolkit.hh.

  {
    return 700.0;
  }

◆ cherenkov()

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 59 of file JPDFToolkit.hh.

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

◆ getDeltaRayTmin()

double JPHYSICS::getDeltaRayTmin ( )

inline

Get minimum delta-ray kinetic energy.

Returns: minimum delta-ray kinetic energy [GeV]

Definition at line 73 of file JPDFToolkit.hh.

  {
    const double Emin = MASS_ELECTRON / getSinThetaC(); // delta-ray Cherenkov threshold [GeV]
 
    return getKineticEnergy(Emin, MASS_ELECTRON);
  }

◆ getDeltaRayTmax()

double JPHYSICS::getDeltaRayTmax	(	const double	E,
		const double	M )

inline

Get maximum delta-ray kinetic energy for given lepton energy and mass.

This formula is taken from reference https://pdg.lbl.gov/2020/reviews/rpp2020-rev-passage-particles-matter.pdf

N.B.: This function should not be used for electrons.
For electrons use 0.5 * getKineticEnergy(E, MASS_ELECTRON) instead.

Parameters

E	particle energy [GeV]
M	particle mass [GeV]

Returns: maximum delta-ray kinetic energy [GeV]

Definition at line 93 of file JPDFToolkit.hh.

  {
    const double ratio = MASS_ELECTRON / M;
 
    const double gamma = E / M;
    const double beta  = sqrt((1.0 + 1.0/gamma) * (1.0 - 1.0/gamma));    
 
    return ( (2.0*MASS_ELECTRON*beta*beta*gamma*gamma) /
             (1.0 + 2.0*gamma*ratio + ratio*ratio) );
  }                          

◆ getDeltaRays() [1/2]

double JPHYSICS::getDeltaRays	(	const double	E,
		const double	M,
		const double	Tmin,
		const double	Tmax,
		const double	Z,
		const double	A )

inline

Get equivalent EM-shower energy due to delta-rays per unit track length
for an ionising particle with given energy and given mass.

N.B: For this function a form factor $F(T) = \frac{1}{2E^2}T^2 - \frac{\beta^2}{T_{\rm max}} + 1$ is assumed, where
:

corresponds to the delta-ray kinetic energy,
$T_{\rm max}$ to the maximum delta-ray kinetic energy,
to the ionising particle's energy and
$\beta$ to the corresponding particle velocity, relative to the speed of light.

Parameters

E	particle energy [GeV]
M	particle mass [GeV]
Tmin	minimum delta-ray kinetic energy [GeV]
Tmax	maximum delta-ray kinetic energy [GeV]
Z	atomic number [unit]
A	atomic mass [g/mol]

Returns: equivalent energy loss [GeV/m]

Definition at line 124 of file JPDFToolkit.hh.

  {      
    if (Tmax > Tmin) {
 
      const double K = 0.307075;                            // [MeV mol^-1 cm^2]
 
      const double gamma = E / M;
      const double beta  = sqrt((1.0 + 1.0/gamma) * (1 - 1.0/gamma));
      
      const double a = 0.25/(E*E);
      const double b = beta*beta/Tmax;
      const double c = 1.0;
 
      const double W = 0.5 * K * (Z/A) * (1.0/(beta*beta)); // [MeV g^-1 cm^2]
 
      const double sT = Tmax + Tmin;
      const double dT = Tmax - Tmin;
      const double rT = Tmax / Tmin;
        
      const double weight = W * (a*sT*dT - b*dT + c*log(rT));
 
      return weight * DENSITY_SEA_WATER * 1.0e-1;           // [GeV m^-1]
        
    } else {
 
      return 0.0;
    }
  }

◆ getDeltaRays() [2/2]

template<class JFormFactor_t >

double JPHYSICS::getDeltaRays	(	const double	E,
		const double	M,
		const double	Tmin,
		const double	Tmax,
		const double	Z,
		const double	A,
		const JFormFactor_t &	F,
		const int	N = 1000000 )

inline

Get equivalent EM-shower energy due to delta-rays per unit track length
for an ionising particle with given energy and given mass and for a given form factor.

The template parameter corresponds to a class which contains an operator()(const double)
to compute the form factor corresponding to a given delta-ray kinetic energy.

Parameters

E	particle energy [GeV]
M	particle mass [GeV]
Tmin	minimum delta-ray kinetic energy [GeV]
Tmax	maximum delta-ray kinetic energy [GeV]
Z	atomic number [unit]
A	atomic mass [g/mol]
F	form factor functor
N	number of points for numeric integration

Returns: equivalent energy loss [GeV/m]

Definition at line 176 of file JPDFToolkit.hh.

  {
    using namespace JPP;
 
    if (Tmax > Tmin) {
 
      const double K = 0.307075;                            // [MeV mol^-1 cm^2]
 
      const double gamma = E / M;
      const double beta  = sqrt((1.0 + 1.0/gamma) * (1 - 1.0/gamma));
 
      const double W = 0.5 * K * (Z/A) * (1.0/(beta*beta)); // [MeV g^-1 cm^2]
 
      const double xmin = log(Tmin);
      const double xmax = log(Tmax);
      const double dx   = (xmax - xmin) / ((double) N);
        
      double weight = 0.0;
 
      for (double x = xmin; x <= xmax; x += dx) {
 
        const double T = exp(x);
        const double y = W * F(T) * dx;                     // [MeV g^-1 cm^-2]
          
        weight += y;
      }
 
      return weight * DENSITY_SEA_WATER * 1.0e-1;           // [GeV m^-1]
        
    } else {
 
      return 0.0;
    }
  }

◆ getDeltaRaysFromElectron()

double JPHYSICS::getDeltaRaysFromElectron	(	const double	E,
		const JRange< double >	T_GeV = JRange<double>(DELTARAY_TMIN, DELTARAY_TMAX) )

inline

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

N.B.: This function assumes a medium with atomic number Z=10 and atomic mass A=18
and a form factor $F(T) = \frac{1}{2E^2}T^2 - \frac{\beta^2}{T_{\rm max}} + 1$ , where
:

corresponds to the delta-ray kinetic energy,
$T_{\rm max}$ to the maximum delta-ray kinetic energy,
to the ionising particle's energy and
$\beta$ to the corresponding particle velocity, relative to the speed of light.

Parameters

E	electron energy [GeV]
T_GeV	allowed kinetic energy range of delta-rays [GeV]

Returns: equivalent energy loss [GeV/m]

Definition at line 233 of file JPDFToolkit.hh.

  {
    static const double Z = 10.0; // atomic number [unit]
    static const double A = 18.0; // atomic mass   [g/mol]
    
    const double Tmin = T_GeV.constrain(getDeltaRayTmin());
    const double Tmax = T_GeV.constrain(0.5 * getKineticEnergy(E, MASS_ELECTRON));
    
    return getDeltaRays(E, MASS_ELECTRON, Tmin, Tmax, Z, A);
  }

◆ getDeltaRaysFromMuon()

double JPHYSICS::getDeltaRaysFromMuon	(	const double	E,
		const JRange< double >	T_GeV = JRange<double>(DELTARAY_TMIN, DELTARAY_TMAX) )

inline

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

N.B.: This function assumes a medium with atomic number Z=10 and atomic mass A=18
and a form factor $F(T) = \frac{1}{2E^2}T^2 - \frac{\beta^2}{T_{\rm max}} + 1$ , where
:

corresponds to the delta-ray kinetic energy,
$T_{\rm max}$ to the maximum delta-ray kinetic energy,
to the ionising particle's energy and
$\beta$ to the corresponding particle velocity, relative to the speed of light.

Parameters

E	muon energy [GeV]
T_GeV	allowed kinetic energy range of delta-rays [GeV]

Returns: equivalent energy loss [GeV/m]

Definition at line 260 of file JPDFToolkit.hh.

  {
    static const double Z = 10.0; // atomic number [unit]
    static const double A = 18.0; // atomic mass   [g/mol]
    
    const double Tmin = T_GeV.constrain(getDeltaRayTmin());
    const double Tmax = T_GeV.constrain(getDeltaRayTmax(E, MASS_MUON));
    
    return getDeltaRays(E, MASS_MUON, Tmin, Tmax, Z, A);
  }

◆ getDeltaRaysFromMuonFit()

double JPHYSICS::getDeltaRaysFromMuonFit ( const double E )

inline

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

N.B.: This function assumes a medium with atomic number Z=10 and atomic mass A=18
and a form factor $F(T) = \frac{1}{2E^2}T^2 - \frac{\beta^2}{T_{\rm max}} + 1$ , where
:

corresponds to the delta-ray kinetic energy,
$T_{\rm max}$ to the maximum delta-ray kinetic energy,
to the ionising particle's energy and
$\beta$ to the corresponding particle velocity, relative to the speed of light.

Parameters

E	muon energy [GeV]

Returns: equivalent energy loss [GeV/m]

Definition at line 286 of file JPDFToolkit.hh.

  {
    static const double a =  3.195e-01;
    static const double b =  3.373e-01;
    static const double c = -2.731e-02;
    static const double d =  1.610e-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]
      
      if (y > 0.0) {
        return y * DENSITY_SEA_WATER * 1.0e-1;        // [GeV/m]
      }
    }
    
    return 0.0;
  }

◆ getDeltaRaysFromTau()

double JPHYSICS::getDeltaRaysFromTau	(	const double	E,
		const JRange< double >	T_GeV = JRange<double>(DELTARAY_TMIN, DELTARAY_TMAX) )

inline

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

N.B.: This function assumes a medium with atomic number Z=10 and atomic mass A=18
and a form factor $F(T) = \frac{1}{2E^2}T^2 - \frac{\beta^2}{T_{\rm max}} + 1$ , where
:

corresponds to the delta-ray kinetic energy,
$T_{\rm max}$ to the maximum delta-ray kinetic energy,
to the ionising particle's energy and
$\beta$ to the corresponding particle velocity, relative to the speed of light.

Parameters

E	tau energy [GeV]
T_GeV	allowed kinetic energy range of delta-rays [GeV]

Returns: equivalent energy loss [GeV/m]

Definition at line 322 of file JPDFToolkit.hh.

  {
    static const double Z = 10.0; // atomic number [unit]
    static const double A = 18.0; // atomic mass   [g/mol]
    
    const double Tmin = T_GeV.constrain(getDeltaRayTmin());
    const double Tmax = T_GeV.constrain(getDeltaRayTmax(E, MASS_TAU));
    
    return getDeltaRays(E, MASS_TAU, Tmin, Tmax, Z, A);
  }

◆ getDeltaRaysFromTauFit()

double JPHYSICS::getDeltaRaysFromTauFit ( const double E )

inline

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

N.B.: This function assumes a medium with atomic number Z=10 and atomic mass A=18
and a form factor $F(T) = \frac{1}{2E^2}T^2 - \frac{\beta^2}{T_{\rm max}} + 1$ , where
:

corresponds to the delta-ray kinetic energy,
$T_{\rm max}$ to the maximum delta-ray kinetic energy,
to the ionising particle's energy and
$\beta$ to the corresponding particle velocity, relative to the speed of light.

Parameters

E	tau energy [GeV]

Returns: equivalent energy loss [GeV/m]

Definition at line 348 of file JPDFToolkit.hh.

  {
    static const double a = -2.178e-01;
    static const double b =  4.995e-01;
    static const double c = -3.906e-02;
    static const double d =  1.615e-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]
   
      if (y > 0) {
        return y * DENSITY_SEA_WATER * 1.0e-1;        // [GeV/m]
      }
    }
   
    return 0.0;
  }

◆ getDeltaRayProbability()

double JPHYSICS::getDeltaRayProbability ( const double x )

inline

Emission profile of photons from delta-rays.

Profile is taken from reference ANTARES-SOFT-2002-015, J. Brunner (fig. 3).

Parameters

x	cosine emission angle

Returns: probability

Definition at line 378 of file JPDFToolkit.hh.

  {
    //return 1.0 / (4.0 * PI);
    return 0.188 * exp(-1.25 * pow(fabs(x - 0.90), 1.30));
  }

◆ getRayleighCrossSection()

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 392 of file JPDFToolkit.hh.

  {
    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;
  }

◆ getRayleighScatteringLength()

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 413 of file JPDFToolkit.hh.

  {
    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;
  }

◆ getLabel()

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

inline

Get PDF label.

Parameters

pdf PDF type

Returns: PDF label

Definition at line 61 of file JPDFTypes.hh.

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

◆ getPDFType()

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 77 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;
  }

◆ getFilename()

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 JEEP::FILENAME_WILDCARD 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 110 of file JPDFTypes.hh.

  {
    return JEEP::setWildCard(file_name, getLabel(pdf));
  }

◆ is_muon()

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 123 of file JPDFTypes.hh.

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

◆ is_bremsstrahlung()

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 137 of file JPDFTypes.hh.

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

◆ is_deltarays()

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 151 of file JPDFTypes.hh.

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

◆ is_scattered()

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 165 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);
  }

◆ getThetaMCS() [1/2]

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 40 of file JPhysicsSupportkit.hh.

  {
    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;
  }

◆ getThetaMCS() [2/2]

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

inline

Get multiple Coulomb scattering angle for muon.

Parameters

E	Energy [GeV]
x	distance [m]

Returns: angle [rad]

Definition at line 65 of file JPhysicsSupportkit.hh.

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

◆ henyey_greenstein() [1/2]

double JPHYSICS::henyey_greenstein	(	const double	g,
		const double	x )

inline

Auxiliary method to describe light scattering in water (Henyey-Greenstein).

Parameters

g	angular dependence parameter
x	cosine scattering angle

Returns: probability

Definition at line 78 of file JPhysicsSupportkit.hh.

  {
    const double a0 = (1.0 - g*g) / (4*PI);
    const double y  =  1.0 + g*g - 2.0*g*x;
    
    return a0 / (y*sqrt(y));
  }

◆ henyey_greenstein() [2/2]

double JPHYSICS::henyey_greenstein ( const double x )

inline

Auxiliary method to describe light scattering in water (Heneyey-Greenstein).

Parameters

x	cosine scattering angle

Returns: probability

Definition at line 95 of file JPhysicsSupportkit.hh.

  {
    static const double g = 0.924;
    
    return henyey_greenstein(g, x);
  }

◆ rayleigh() [1/2]

double JPHYSICS::rayleigh	(	const double	a,
		const double	x )

inline

Auxiliary method to describe light scattering in water (Rayleigh).

Parameters

a	angular dependence parameter
x	cosine scattering angle

Returns: probability

Definition at line 110 of file JPhysicsSupportkit.hh.

  {
    const double a0 = 1.0 / (1.0 + a/3.0) / (4*PI);
 
    return a0 * (1.0 + a*x*x); 
  }

◆ rayleigh() [2/2]

double JPHYSICS::rayleigh ( const double x )

inline

Auxiliary method to describe light scattering in water (Rayleigh).

Parameters

x	cosine scattering angle

Returns: probability

Definition at line 126 of file JPhysicsSupportkit.hh.

  {
    return rayleigh(0.835, x);
  }

◆ f4()

double JPHYSICS::f4 ( const double x )

inline

Model specific function to describe light scattering in water (f4).

Parameters

x	cosine scattering angle

Returns: probability

Definition at line 138 of file JPhysicsSupportkit.hh.

  {
    static const double g1 = 0.77;
    static const double g2 = 0.00;
    static const double f  = 1.00;
 
    return f * henyey_greenstein(g1,x)  +  (1.0 - f) * henyey_greenstein(g2,x);
  }

◆ p00075()

double JPHYSICS::p00075 ( const double x )

inline

Model specific function to describe light scattering in water (p00075).

Parameters

x	cosine scattering angle

Returns: probability

Definition at line 154 of file JPhysicsSupportkit.hh.

  {
    static const double g = 0.924;
    static const double f = 0.17;
 
    return f * rayleigh(x)  +  (1.0 - f) * henyey_greenstein(g,x);
  }

◆ getKineticEnergy()

double JPHYSICS::getKineticEnergy	(	const double	E,
		const double	m )

inline

Get kinetic energy of particle with given energy and mass.

Parameters

E	energy [GeV]
m	mass [GeV]

Returns: kinetic energy [GeV]

Definition at line 24 of file JPhysicsToolkit.hh.

  {
    if (E > m) {
      return sqrt((E - m) * (E + m));
    } else {
      return 0.0;
    }
  }

Variable Documentation

◆ C

const double JPHYSICS::C = 0.299792458

static

Physics constants.

Speed of light in vacuum [m/ns]

Definition at line 25 of file JPhysics/JConstants.hh.

◆ C_INVERSE

const double JPHYSICS::C_INVERSE = 1.0/C

static

Inverse speed of light in vacuum [ns/m].

Definition at line 26 of file JPhysics/JConstants.hh.

◆ AVOGADRO

const double JPHYSICS::AVOGADRO = 6.0221415e23

static

Avogadro's number.

Definition at line 27 of file JPhysics/JConstants.hh.

◆ NUCLEON_MOLAR_MASS

const double JPHYSICS::NUCLEON_MOLAR_MASS = 1.0

static

nucleon molar mass [g/mol]

Definition at line 28 of file JPhysics/JConstants.hh.

◆ H

const double JPHYSICS::H = 4.13566733e-15

static

Planck constant [eV s].

Definition at line 29 of file JPhysics/JConstants.hh.

◆ HBAR

const double JPHYSICS::HBAR = H/(2*PI)

static

Planck constant [eV s].

Definition at line 30 of file JPhysics/JConstants.hh.

◆ HBARC

const double JPHYSICS::HBARC = HBAR*C*1.0e9

static

Planck constant [eV m].

Definition at line 31 of file JPhysics/JConstants.hh.

◆ ALPHA_ELECTRO_MAGNETIC

const double JPHYSICS::ALPHA_ELECTRO_MAGNETIC = 1.0/137.036

static

Electro-Magnetic coupling constant.

Definition at line 32 of file JPhysics/JConstants.hh.

◆ THETA_MCS

const double JPHYSICS::THETA_MCS = 13.6e-3

static

Multiple Coulomb scattering constant [GeV].

Definition at line 33 of file JPhysics/JConstants.hh.

◆ R_EARTH_KM

const double JPHYSICS::R_EARTH_KM = 6371

static

Geophysics constants.

Radius of the Earth [m]

Definition at line 38 of file JPhysics/JConstants.hh.

◆ DENSITY_EARTH

const double JPHYSICS::DENSITY_EARTH = 5.51

static

Average density of the Earth [gr/cm³].

Definition at line 39 of file JPhysics/JConstants.hh.

◆ DENSITY_SEA_WATER

const double JPHYSICS::DENSITY_SEA_WATER = 1.038

static

Fixed environment values.

Density of sea water [g/cm^3]

Definition at line 44 of file JPhysics/JConstants.hh.

◆ DENSITY_ROCK

const double JPHYSICS::DENSITY_ROCK = 2.65

static

Density of rock [g/cm^3].

Definition at line 45 of file JPhysics/JConstants.hh.

◆ SALINITY_SEA_WATER

const double JPHYSICS::SALINITY_SEA_WATER = 0.035

static

Salinity of sea water.

Definition at line 46 of file JPhysics/JConstants.hh.

◆ INDEX_OF_REFRACTION_WATER

const double JPHYSICS::INDEX_OF_REFRACTION_WATER = 1.3800851282

static

Average index of refraction of water corresponding to the group velocity.

Definition at line 47 of file JPhysics/JConstants.hh.

◆ X0_WATER_M

const double JPHYSICS::X0_WATER_M = 0.36

static

Radiation length pure water [m].

Definition at line 48 of file JPhysics/JConstants.hh.

◆ TAN_THETA_C_WATER

const double JPHYSICS::TAN_THETA_C_WATER = sqrt((INDEX_OF_REFRACTION_WATER - 1.0) * (INDEX_OF_REFRACTION_WATER + 1.0))

static

Derived quantities of optical medium.

Average tangent corresponding to the group velocity

Definition at line 53 of file JPhysics/JConstants.hh.

◆ COS_THETA_C_WATER

const double JPHYSICS::COS_THETA_C_WATER = 1.0 / INDEX_OF_REFRACTION_WATER

static

Average cosine corresponding to the group velocity.

Definition at line 54 of file JPhysics/JConstants.hh.

◆ SIN_THETA_C_WATER

const double JPHYSICS::SIN_THETA_C_WATER = TAN_THETA_C_WATER * COS_THETA_C_WATER

static

Average sine corresponding to the group velocity.

Definition at line 55 of file JPhysics/JConstants.hh.

◆ KAPPA_WATER

const double JPHYSICS::KAPPA_WATER = 0.96

static

Average R-dependence of arrival time of Cherenkov light.

Definition at line 56 of file JPhysics/JConstants.hh.

◆ MASS_PHOTON

const double JPHYSICS::MASS_PHOTON = 0.0

static

Particle masses.

Note that the neutrino masses are set to zero. photon mass [GeV]

Definition at line 62 of file JPhysics/JConstants.hh.

◆ MASS_ELECTRON_NEUTRINO

const double JPHYSICS::MASS_ELECTRON_NEUTRINO = 0.0

static

electron neutrino mass [GeV]

Definition at line 64 of file JPhysics/JConstants.hh.

◆ MASS_MUON_NEUTRINO

const double JPHYSICS::MASS_MUON_NEUTRINO = 0.0

static

muon neutrino mass [GeV]

Definition at line 65 of file JPhysics/JConstants.hh.

◆ MASS_TAU_NEUTRINO

const double JPHYSICS::MASS_TAU_NEUTRINO = 0.0

static

tau neutrino mass [GeV]

Definition at line 66 of file JPhysics/JConstants.hh.

◆ MASS_ELECTRON

const double JPHYSICS::MASS_ELECTRON = 0.510998946e-3

static

electron mass [GeV]

Definition at line 68 of file JPhysics/JConstants.hh.

◆ MASS_MUON

const double JPHYSICS::MASS_MUON = 0.1056583745

static

muon mass [GeV]

Definition at line 69 of file JPhysics/JConstants.hh.

◆ MASS_TAU

const double JPHYSICS::MASS_TAU = 1.77682

static

tau mass [GeV]

Definition at line 70 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_PION

const double JPHYSICS::MASS_NEUTRAL_PION = 0.1349766

static

pi_0 mass [GeV]

Definition at line 72 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_PION

const double JPHYSICS::MASS_CHARGED_PION = 0.13957018

static

pi^+/- mass [GeV]

Definition at line 73 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_KAON

const double JPHYSICS::MASS_NEUTRAL_KAON = 0.497614

static

K_0 mass [GeV].

Definition at line 74 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_KAON

const double JPHYSICS::MASS_CHARGED_KAON = 0.493677

static

K^+/- mass [GeV].

Definition at line 75 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_RHO

const double JPHYSICS::MASS_NEUTRAL_RHO = 0.77526

static

rho_0 mass [GeV]

Definition at line 76 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_RHO

const double JPHYSICS::MASS_CHARGED_RHO = 0.77511

static

rho^+/- mass [GeV]

Definition at line 77 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_D

const double JPHYSICS::MASS_NEUTRAL_D = 1.86483

static

D_0 mass [GeV].

Definition at line 78 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_D

const double JPHYSICS::MASS_CHARGED_D = 1.86965

static

D^+/- mass [GeV].

Definition at line 79 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_D_S

const double JPHYSICS::MASS_CHARGED_D_S = 1.96834

static

D_s^+/- mass [GeV].

Definition at line 80 of file JPhysics/JConstants.hh.

◆ MASS_PROTON

const double JPHYSICS::MASS_PROTON = 0.9382720813

static

proton mass [GeV]

Definition at line 82 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRON

const double JPHYSICS::MASS_NEUTRON = 0.9395654133

static

neutron mass [GeV]

Definition at line 83 of file JPhysics/JConstants.hh.

◆ MASS_DELTA_1232

const double JPHYSICS::MASS_DELTA_1232 = 1.232

static

Delta (1232) mass [GeV].

Definition at line 84 of file JPhysics/JConstants.hh.

◆ MASS_LAMBDA

const double JPHYSICS::MASS_LAMBDA = 1.115683

static

Lambda mass [GeV].

Definition at line 85 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_SIGMA

const double JPHYSICS::MASS_NEUTRAL_SIGMA = 1.192642

static

Sigma_0 mass [GeV].

Definition at line 86 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_SIGMA

const double JPHYSICS::MASS_CHARGED_SIGMA = 1.18937

static

Sigma^+/- mass [GeV].

Definition at line 87 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_XI

const double JPHYSICS::MASS_NEUTRAL_XI = 1.31486

static

Xi_0 mass [GeV].

Definition at line 88 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_XI

const double JPHYSICS::MASS_CHARGED_XI = 1.32171

static

Xi^+/- mass [GeV].

Definition at line 89 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_OMEGA

const double JPHYSICS::MASS_CHARGED_OMEGA = 1.67245

static

Omega^+/- mass [GeV].

Definition at line 90 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_LAMBDA_C

const double JPHYSICS::MASS_CHARGED_LAMBDA_C = 2.28646

static

Lambda_c^+/- mass [GeV].

Definition at line 91 of file JPhysics/JConstants.hh.

◆ MASS_DOUBLYCHARGED_SIGMA_C

const double JPHYSICS::MASS_DOUBLYCHARGED_SIGMA_C = 2.45397

static

Sigma_c^++/– mass [GeV]

Definition at line 92 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_SIGMA_C

const double JPHYSICS::MASS_CHARGED_SIGMA_C = 2.4529

static

Sigma_c^+/- mass [GeV].

Definition at line 93 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_SIGMA_C

const double JPHYSICS::MASS_NEUTRAL_SIGMA_C = 2.45375

static

Sigma_c_0 mass [GeV].

Definition at line 94 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_XI_C

const double JPHYSICS::MASS_CHARGED_XI_C = 2.46793

static

Xi_c^+/- mass [GeV].

Definition at line 95 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_XI_C

const double JPHYSICS::MASS_NEUTRAL_XI_C = 2.47091

static

Xi_c_0 mass [GeV].

Definition at line 96 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_OMEGA_C

const double JPHYSICS::MASS_NEUTRAL_OMEGA_C = 2.6952

static

Omega_c_0 mass [GeV].

Definition at line 97 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_B

const double JPHYSICS::MASS_NEUTRAL_B = 5.27958

static

B_0 mass [GeV].

Definition at line 98 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_B

const double JPHYSICS::MASS_CHARGED_B = 5.27926

static

B^+/- mass [GeV].

Definition at line 99 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_B_S

const double JPHYSICS::MASS_NEUTRAL_B_S = 5.36677

static

B_s^0 mass [GeV].

Definition at line 100 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_LAMBDA_B

const double JPHYSICS::MASS_NEUTRAL_LAMBDA_B = 5.6194

static

Lambda_b^0 mass [GeV].

Definition at line 101 of file JPhysics/JConstants.hh.

◆ MASS_NEUTRAL_XI_B

const double JPHYSICS::MASS_NEUTRAL_XI_B = 5.7878

static

Xi_b^0 mass [GeV].

Definition at line 102 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_XI_B

const double JPHYSICS::MASS_CHARGED_XI_B = 5.7911

static

Xi_b^+/- mass [GeV].

Definition at line 103 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_OMEGA_B

const double JPHYSICS::MASS_CHARGED_OMEGA_B = 6.071

static

Omega_b^+/- mass [GeV].

Definition at line 104 of file JPhysics/JConstants.hh.

◆ MASS_CHARGED_B_C

const double JPHYSICS::MASS_CHARGED_B_C = 6.2756

static

B_c^+/- mass [GeV].

Definition at line 105 of file JPhysics/JConstants.hh.

◆ gWater

const JGeaneWater JPHYSICS::gWater

static

Function object for energy loss of muon in sea water.

Definition at line 381 of file JGeane.hh.

◆ geant

const JGeant JPHYSICS::geant(geanx, 0.0001)	(	geanx	,
		0.	0001 )

static

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

◆ geanx

const JGeanx JPHYSICS::geanx(0.35, -5.40)	(	0.	35,
		-5.	40 )

static

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

◆ geanz

const JGeanz JPHYSICS::geanz(1.85, 0.62, 0.54)	(	1.	85,
		0.	62,
		0.	54 )

static

Function object for longitudinal EM-shower profile.

◆ getSterCoefficient

JSterCoefficient JPHYSICS::getSterCoefficient

static

Function object for Ster coefficients.

Definition at line 82 of file JIonization.hh.

◆ cc_nu

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.

◆ nc_nu

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.

◆ cc_nubar

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.

◆ nc_nubar

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.

◆ MODULE_RADIUS_M

double JPHYSICS::MODULE_RADIUS_M = 0.216

static

Radius of optical module [m].

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

Definition at line 40 of file JPDF.hh.

◆ DELTARAY_TMIN

const double JPHYSICS::DELTARAY_TMIN = 0.000915499

static

Minimum allowed delta-ray kinetic energy [GeV].

Definition at line 26 of file JPDFToolkit.hh.

◆ DELTARAY_TMAX

const double JPHYSICS::DELTARAY_TMAX = 1.0e10

static

Maximum allowed delta-ray kinetic energy [GeV].

Definition at line 27 of file JPDFToolkit.hh.

◆ petzhold

const JPetzhold JPHYSICS::petzhold

static

Function object for measurement of light scattering in water.

Definition at line 270 of file JPhysicsSupportkit.hh.

◆ getAbsorptionLengthOfPureWater

const JAbsorptionLengthOfPureWater JPHYSICS::getAbsorptionLengthOfPureWater

static

Function object for absorption length of pure water.

Definition at line 470 of file JPhysicsSupportkit.hh.

◆ EErad_t

const JRadiationSource_t JPHYSICS::EErad_t = { &JRadiation::TotalCrossSectionEErad, &JRadiation::EfromEErad, &JRadiation::ThetaRMSfromEErad }

static

Definition at line 512 of file JRadiation.hh.

512{ &JRadiation::TotalCrossSectionEErad, &JRadiation::EfromEErad, &JRadiation::ThetaRMSfromEErad };

◆ Brems_t

const JRadiationSource_t JPHYSICS::Brems_t = { &JRadiation::TotalCrossSectionBrems, &JRadiation::EfromBrems, &JRadiation::ThetaRMSfromBrems }

static

Definition at line 513 of file JRadiation.hh.

513{ &JRadiation::TotalCrossSectionBrems, &JRadiation::EfromBrems, &JRadiation::ThetaRMSfromBrems };

◆ GNrad_t

const JRadiationSource_t JPHYSICS::GNrad_t = { &JRadiation::TotalCrossSectionGNrad, &JRadiation::EfromGNrad, &JRadiation::ThetaRMSfromGNrad }

static

Definition at line 514 of file JRadiation.hh.

514{ &JRadiation::TotalCrossSectionGNrad, &JRadiation::EfromGNrad, &JRadiation::ThetaRMSfromGNrad };

Classes

Typedefs

Enumerations

Functions

Variables

Detailed Description

Typedef Documentation

◆ JGeantFunction1D_t

◆ JRateL0_t

◆ JRateL1_t

◆ multiplicity_type

◆ JElement2D_t

◆ JElement3D_t

◆ JTimeRange

Enumeration Type Documentation

◆ JPDFType_t

Function Documentation

◆ getSpeedOfLight()

◆ getInverseSpeedOfLight()

◆ getIndexOfRefraction()

◆ getIndexOfRefractionPhase()

◆ getTanThetaC()

◆ getCosThetaC()

◆ getSinThetaC()

◆ getKappaC()

◆ geanc()

◆ gRock()

◆ getMinimalWavelength()

◆ getMaximalWavelength()

◆ cherenkov()

◆ getDeltaRayTmin()

◆ getDeltaRayTmax()

◆ getDeltaRays() [1/2]

◆ getDeltaRays() [2/2]

◆ getDeltaRaysFromElectron()

◆ getDeltaRaysFromMuon()

◆ getDeltaRaysFromMuonFit()

◆ getDeltaRaysFromTau()

◆ getDeltaRaysFromTauFit()

◆ getDeltaRayProbability()

◆ getRayleighCrossSection()

◆ getRayleighScatteringLength()

◆ getLabel()

◆ getPDFType()

◆ getFilename()

◆ is_muon()

◆ is_bremsstrahlung()

◆ is_deltarays()

◆ is_scattered()

◆ getThetaMCS() [1/2]

◆ getThetaMCS() [2/2]

◆ henyey_greenstein() [1/2]

◆ henyey_greenstein() [2/2]

◆ rayleigh() [1/2]

◆ rayleigh() [2/2]

◆ f4()

◆ p00075()

◆ getKineticEnergy()

Variable Documentation

◆ C

◆ C_INVERSE

◆ AVOGADRO

◆ NUCLEON_MOLAR_MASS

◆ H

◆ HBAR

◆ HBARC

◆ ALPHA_ELECTRO_MAGNETIC

◆ THETA_MCS

◆ R_EARTH_KM

◆ DENSITY_EARTH

◆ DENSITY_SEA_WATER

◆ DENSITY_ROCK

◆ SALINITY_SEA_WATER

◆ INDEX_OF_REFRACTION_WATER

◆ X0_WATER_M

◆ TAN_THETA_C_WATER

◆ COS_THETA_C_WATER

◆ SIN_THETA_C_WATER

◆ KAPPA_WATER

◆ MASS_PHOTON