00001 00002 #ifndef _CONSTANTS_ 00003 #define _CONSTANTS_ 00004 #include <math.h> 00005 namespace MuELoss { 00006 00007 static const double a = 1/137.036; // fine structure constant 00008 static const double a_2 = a*a; // fine structure constant^2 00009 static const double a_3 = a_2*a; // fine structure constant^3 00010 static const double a_4 = a_3*a; // fine structure constant^4 00011 static const double Me = 0.51099906; // electron mass (MeV) 00012 static const double Me_2 = Me*Me; // electron mass^2 (MeV^2) 00013 static const double Mm = 105.658389; // muon mass (MeV) 00014 static const double Mm_2 = Mm*Mm; // muon mass^2 (MeV^2) 00015 static const double Mm_Me = Mm/Me; // muon mass / electron mass 00016 static const double Mm_Me_2 = Mm_Me*Mm_Me; // (muon mass / electron mass)^2 00017 static const double Me_Mm = Me/Mm; // electron mass / muon mass 00018 static const double Me_Mm_2 = Me_Mm*Me_Mm; // (electron mass / muon mass)^2 00019 static const double Mp = 938.27231; // proton mass (MeV) 00020 static const double Le = 3.8616; // Compton wavelength of the electron * 10^11 cm 00021 static const double Le_2 = Le*Le; // Compton wavelength of the electron^2 * 10^22 cm^2 00022 static const double e = 2.718; 00023 static const double sqrt_e = sqrt(e); // sqrt(2.718...) 00024 static const double pi = 3.141; 00025 static const double Na = 6.023; // Avogadro's number*10^(-23) 00026 static const double R = 189; // value of radiation logarithm ? 00027 static const double m1_2 = 0.54e6; // m1^2 from photonuclear diff. cross section (in MeV^2) 00028 static const double m2_2 = 1.80e6; // m2^2 from photonuclear diff. cross section (in MeV^2) 00029 00030 00031 static const double megaelectronvolt = 1.0; 00032 00033 static const double electronvolt = 1.e-06 *megaelectronvolt; 00034 static const double kiloelectronvolt = 1.e+03 *electronvolt; 00035 static const double gigaelectronvolt = 1.e+09 *electronvolt; 00036 static const double teraelectronvolt = 1.e+12 *electronvolt; 00037 static const double petaelectronvolt = 1.e+15 *electronvolt; 00038 00039 static const double eV = electronvolt; 00040 static const double keV = kiloelectronvolt; 00041 static const double MeV = megaelectronvolt; 00042 static const double GeV = gigaelectronvolt; 00043 static const double TeV = teraelectronvolt; 00044 static const double PeV = petaelectronvolt; 00045 00046 } // namespace MuELoss 00047 00048 #endif // _CONSTANTS_ 00049
1.3.9.1