PhysConst_constants.hpp

#pragma once
#include <cmath> // PI
 
namespace PhysConst {
 
constexpr double c_SI = 299792458.0;
 
constexpr double h_SI = (6.62607015e-34);
 
constexpr double hbar_SI = h_SI / (2.0 * M_PI);
 
constexpr double e_C = 1.602176634e-19;
 
constexpr double alpha = 1.0 / 137.035999177;
 
constexpr double g_e = -2.00231930436092;
 
constexpr double e_on_me_SI = 1.75882000838e11;
 
constexpr double m_e_kg = e_C / e_on_me_SI;
 
constexpr double m_p = 1836.152673426;
 
constexpr double m_muon = 206.7682827;
 
constexpr double m_tau = 3477.23;
 
constexpr double u_NMU_kg = 1.66053906892e-27;
 
constexpr double aB_m = hbar_SI / (m_e_kg * c_SI * alpha);
 
constexpr double Hartree_eV = 27.211386245981;
// nb: this one different from m*c^2*alpha^2 - probably because it measured better than alpha??
 
constexpr double GF_GeV2 = 1.1663787e-5;
 
//======================================================================
//======================================================================
 
constexpr double alpha2 = alpha * alpha;
 
constexpr double c = 1.0 / alpha;
constexpr double c2 = c * c;
 
constexpr double hbarc_MeVfm = (hbar_SI * c_SI / e_C) * 1.0e9;
 
constexpr double m_e_MeV = Hartree_eV * c2 / 1.0e6;
 
constexpr double u_NMU = u_NMU_kg / m_e_kg;
 
// Length:
 
constexpr double aB_cm = aB_m * (1.0e+2);
constexpr double aB_fm = aB_m * (1.0e+15);
constexpr double aB_nm = aB_m * (1.0e+9);
 
constexpr double Hartree_Hz = Hartree_eV * e_C / h_SI;
 
constexpr double Hartree_MHz = Hartree_Hz * (1.0e-6);
constexpr double Hartree_GHz = Hartree_Hz * (1.0e-9);
constexpr double Hartree_invcm = 1.0 / (2.0 * M_PI * c * aB_cm);
constexpr double HartreeWL_nm = 2.0 * M_PI * c * aB_m * (1.0e+9);
 
constexpr double hbar_on_EH = hbar_SI / e_C / Hartree_eV;
 
constexpr double muB_SI = 0.5; //
constexpr double muB_CGS = 0.5 * alpha;
constexpr double muN_SI = muB_SI / m_p;
constexpr double muN_CGS = muB_CGS / m_p;
constexpr double muN_CGS_MHz = Hartree_MHz * muB_CGS / m_p;
 
constexpr double barn_m2 = 1.0e-28;
constexpr double barn_au = barn_m2 / (aB_m * aB_m);
constexpr double barn_MHz = barn_au * Hartree_MHz;
 
 
constexpr double GF = GF_GeV2 * alpha * m_e_MeV * m_e_MeV * 1e-6;
constexpr double GFe11 = GF * 1.0e11;
 
} // namespace PhysConst