39#ifndef HEP_PHYSICAL_CONSTANTS_H
40#define HEP_PHYSICAL_CONSTANTS_H
49static const double pi = 3.14159265358979323846;
50static const double twopi = 2*pi;
51static const double halfpi = pi/2;
52static const double pi2 = pi*pi;
57static const double Avogadro = 6.02214179e+23/mole;
63static const double c_light = 2.99792458e+8 * m/s;
64static const double c_squared = c_light * c_light;
71static const double h_Planck = 6.62606896e-34 * joule*s;
72static const double hbar_Planck = h_Planck/twopi;
73static const double hbarc = hbar_Planck * c_light;
74static const double hbarc_squared = hbarc * hbarc;
79static const double electron_charge = - eplus;
80static const double e_squared = eplus * eplus;
87static const double electron_mass_c2 = 0.510998910 * MeV;
88static const double proton_mass_c2 = 938.272013 * MeV;
89static const double neutron_mass_c2 = 939.56536 * MeV;
90static const double amu_c2 = 931.494028 * MeV;
91static const double amu = amu_c2/c_squared;
97static const double mu0 = 4*pi*1.e-7 * henry/m;
98static const double epsilon0 = 1./(c_squared*mu0);
103static const double elm_coupling = e_squared/(4*pi*epsilon0);
104static const double fine_structure_const = elm_coupling/hbarc;
105static const double classic_electr_radius = elm_coupling/electron_mass_c2;
106static const double electron_Compton_length = hbarc/electron_mass_c2;
107static const double Bohr_radius = electron_Compton_length/fine_structure_const;
109static const double alpha_rcl2 = fine_structure_const
110 *classic_electr_radius
111 *classic_electr_radius;
113static const double twopi_mc2_rcl2 = twopi*electron_mass_c2
114 *classic_electr_radius
115 *classic_electr_radius;
119static const double k_Boltzmann = 8.617343e-11 * MeV/kelvin;
124static const double STP_Temperature = 273.15*kelvin;
125static const double STP_Pressure = 1.*atmosphere;
126static const double kGasThreshold = 10.*mg/cm3;
131static const double universe_mean_density = 1.e-25*g/cm3;