gambit is hosted by Hepforge, IPPP Durham
GAMBIT  v1.5.0-2191-ga4742ac
a Global And Modular Bsm Inference Tool
numerical_constants.hpp
Go to the documentation of this file.
1 // GAMBIT: Global and Modular BSM Inference Tool
2 // *********************************************
28 //
30 
31 #ifndef __numerical_constants_hpp__
32 #define __numerical_constants_hpp__
33 
34 #include <cmath>
35 
36 namespace Gambit
37 {
38 
39  const double pi = 3.141592653589793238462643383279502884197;
40  const double root2 = sqrt(2.0);
41  const double zeta3 = 1.2020569031595942855; // Riemann zeta function of 3
42  const double hbar = 6.582119514e-25; // GeV s (http://pdg.lbl.gov/2017/reviews/rpp2017-rev-phys-constants.pdf)
43  const double K2eV = 8.6173303e-5; // eV per K
44  const double eV2g = 1.782661907e-33; // eV per g
45  const double gev2cm = 197.327053e-16; // cm per GeV^-1
46  const double gev2cm2 = pow(197.327053e-16, 2.0); // cm^2 per GeV^-2
47  const double gev2pb = gev2cm2*1e36; // pb per GeV^-2
48  const double gev2tocm3s1 = 1.16733e-17; // cm^3 s^-1 per GeV^-2
49  const double s2cm = 2.99792458e10; // cm per s
50  const double m_planck = 1.220910e19; // Planck mass (GeV)
51  const double m_planck_red = m_planck/sqrt(8.0*pi); // reduced Planck mass (GeV)
52  const double atomic_mass_unit = 0.931494028; // atomic mass unit (GeV/c^2)
53  const double m_proton_amu = 1.00727646688; // proton mass (amu)
54  const double m_neutron_amu = 1.0086649156; // neutron mass (amu)
55  const double m_proton = m_proton_amu * atomic_mass_unit; // proton mass (GeV/c^2)
56  const double m_neutron = m_neutron_amu * atomic_mass_unit; // neutron mass (GeV/c^2)
57  const double m_electron = 0.5109989461e-3; // electron mass (GeV/c^2)
58  const double alpha_EM = 7.2973525664e-3; // fine structure constant
59 
61  const double c_SI = s2cm/100; // speed of light in m/s
62  const double eV_to_J = 1.6021766208e-19; // electron charge in C
63  const double Mpc_SI = 969394202136*pow(10,11)/pi; // Mpc in m
64 
65  const double GN_cgs = 6.67408e-8; // Newton's constant in cm^3.g^-1.s^-2
66  const double GN_SI = GN_cgs/1e3; // Newton's constant in m^3.kg^-1.s^-2
67 
68  const double kB_SI = 1.38064852e-23; // Boltzmann constant in Kg/K^4/s^3
69  const double kB_eV_over_K = kB_SI/eV_to_J; // Boltzmann constant in eV/K
70 
71  const double hP_SI = 6.626070040e-34; // Planck const. in Js
72  const double hP_eVs = hP_SI/eV_to_J; // Planck const. in eVs
73  const double hc_eVcm =hP_eVs*s2cm; // Planck const. x speed of light in eV cm
74 
75  const double sigmaB_SI = 2*pow(pi,5)*pow(kB_SI,4)/(15*pow(hP_SI,3)*c_SI*c_SI); // Stefan-Boltzman constant in W/m^2/K^4 = Kg/K^4/s^3
76 
77  const double Neff_SM = 3.045; // effective number of relativistic dof in the early Universe
78  // the value of 3.045 holds for 3 SM neutrinos in the absence of any non-standard particles
79  // or components. Value from de Salas, Pastor '16, arXiv:1606.06986
80 
85  constexpr double alpha_e_OS_thomson_limit = 0.0072973525664;
86 
91  constexpr double alpha_e_OS_MZ = 1. / 128.944;
92 
101  constexpr double delta_alpha_OS = 1. - alpha_e_OS_thomson_limit / alpha_e_OS_MZ;
102 
103  static const struct Mesons_masses
104  {
105  static constexpr double pi0 = 0.135; // neutral pion mass (GeV/c^2)
106  static constexpr double pi_plus = 0.1396; // charged pion mass (GeV/c^2)
107  static constexpr double pi_minus = 0.1396; // charged pion mass (GeV/c^2)
108  static constexpr double eta = 0.547; // eta mass (GeV/c^2)
109  static constexpr double rho0 = 0.775; // neutral rho meson mass (GeV/c^2)
110  static constexpr double rho_plus = 0.775; // charged rho meson mass (GeV/c^2)
111  static constexpr double rho_minus = 0.775; // charged rho meson mass (GeV/c^2)
112  static constexpr double omega = 0.7827; // omega meson mass (GeV/c^2)
113  static constexpr double kaon_plus = 0.4937; // charged kaon meson mass (GeV/c^2)
114  static constexpr double kaon_minus = 0.4937; // charged kaon meson mass (GeV/c^2)
115  static constexpr double kaon0 = 0.4976; // neutral kaon meson mass (GeV/c^2)
116  static constexpr double rho1450 = 1.465; // rho(1450) mass (GeV/c^2)
117  static constexpr double D_plus = 1.86962; // charged D meson mass (GeV/c^2)
118  static constexpr double D_s = 1.96847; // D_s meson mass (GeV/c^2)
119  static constexpr double B_plus = 5.27929; // charged B meson mass (GeV/c^2)
120  static constexpr double B_s = 5.36679; // B_s meson mass (GeV/c^2)
121  static constexpr double B_c = 6.2751; // B_c meson mass (GeV/c^2)
122  static constexpr double eta_prime = 0.95778; // eta prime meson mass (GeV/c^2)
123  static constexpr double eta_c = 2.9839; // eta_c meson mass (GeV/c^2)
124  static constexpr double Dstar_plus = 2.01026; // charged D* meson mass (GeV/c^2)
125  static constexpr double Dstar_s = 2.1122; // D*_s meson mass (GeV/c^2)
126  static constexpr double phi = 1.019461; // phi meson mass (GeV/c^2)
127  static constexpr double Jpsi = 3.0969; // Jpsi meson mass (GeV/c^2)
128  } meson_masses;
129 
130  static const struct Mesons_decay_constants
131  {
132  // PDG 2018
133  static constexpr double pi_plus = 0.1302; // (GeV)
134  static constexpr double pi0 = 0.1302; // (GeV)
135  static constexpr double K_plus = 0.1557; // (GeV)
136  static constexpr double D_plus = 0.2126; // (GeV)
137  static constexpr double D_s = 0.2499; // (GeV)
138  static constexpr double B_plus = 0.190; // (GeV)
139 
140  // From 1503.05762
141  static constexpr double B_c = 0.434; // (GeV)
142 
143  // From 1805.08567
144  static constexpr double eta = 0.0817; // (GeV)
145  static constexpr double eta_prime = -0.0947; // (GeV)
146  static constexpr double eta_c = 0.237; // (GeV)
147 
148  // From 0602110, using tau decays for rho
149  static constexpr double rho_plus = 0.209; // (GeV)
150  static constexpr double rho0 = 0.209; // (GeV)
151  static constexpr double phi = 0.229; // (GeV)
152 
153  // From 1708.07274, average of theoretical computations
154  static constexpr double Dstar_plus = 0.24675;// (GeV)
155  static constexpr double Dstar_s = 0.284; // (GeV)
156 
157  // From 0901.3589, thought not sure where they got it from
158  static constexpr double omega = 0.195; // (GeV)
159 
160  // From 1312.2858
161  static constexpr double Jpsi = 0.418; // (GeV)
162 
163  } meson_decay_constants;
164 
165 
169  const double mw_central_observed = 80.385;
170  const double mw_err_observed = 0.015;
172 
173 }
174 
175 #endif //#defined __numerical_constants_hpp__
const double m_planck
const double c_SI
Values from Particle Data Group 2018 (http://pdg.lbl.gov/2018/reviews/rpp2018-rev-phys-constants.pdf)
const double m_proton
const double mw_err_observed
const double alpha_EM
const double m_electron
const double GN_SI
const double s2cm
const double m_neutron
const double sigmaB_SI
const double kB_SI
const double atomic_mass_unit
const double hbar
const double m_proton_amu
const double kB_eV_over_K
const double m_planck_red
const double root2
const double pi
const double hP_SI
const double m_neutron_amu
const double hc_eVcm
const double zeta3
const double gev2cm2
const double Mpc_SI
const double GN_cgs
const double eV_to_J
constexpr double alpha_e_OS_MZ
in OS scheme from 1105.3149
const double gev2cm
double pow(const double &a)
Outputs a^i.
const double Neff_SM
const double gev2pb
const double K2eV
constexpr double delta_alpha_OS
in OS scheme.
constexpr double alpha_e_OS_thomson_limit
Thomson limit of in OS scheme from PDG 2017
const double gev2tocm3s1
const double eV2g
const double mw_central_observed
M_W (Breit-Wigner mass parameter ~ pole) = 80.385 +/- 0.015 GeV (1 sigma), Gaussian.
TODO: see if we can use this one:
Definition: Analysis.hpp:33
const double hP_eVs