7constexpr double PI = 3.1415926535897932;
8constexpr double PI2 = 9.869604401089358;
9constexpr double PI3 = 31.006276680299816;
10constexpr double INV_PI = 0.318309886183791;
11constexpr double INV_PI2 = 0.101321183642338;
12constexpr double INV_PI3 = 0.032251534433199;
13constexpr double E = 2.7182818284590452;
14constexpr double ZETA3 = 1.2020569031595942;
15constexpr double RT2 = 1.414213562373095;
16constexpr double INV_RT2 = 0.707106781186547;
17constexpr double GAMMA = 0.5772156649015328;
20constexpr std::complex<double>
I = std::complex<double>(0, 1);
23constexpr double HBAR = 6.58211889e-25;
25constexpr double M_PLANCK = 1.2209102930946623e+19;
27constexpr double M_P = 0.9382720;
28constexpr double M_N = 0.9395654;
constexpr double GEV_TO_INV_S
constexpr double ERF_INV_RT2
constexpr std::complex< double > I
constexpr double GEV_TO_INV_PS
constexpr double M_PLANCK