Physics.h

Go to the documentation of this file.

//
// Class Physics
//   A namespace defining various mathematical and physical constants.
//
// Copyright (c) 2015-2021, Christof Metzger-Kraus, Helmholtz-Zentrum Berlin
//                          Pedro Calvo, CIEMAT, Spain
// All rights reserved
//
// Implemented as part of the PhD thesis
// "Optimizing the radioisotope production of the novel AMIT
// superconducting weak focusing cyclotron"
//
// This file is part of OPAL.
//
// OPAL is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// You should have received a copy of the GNU General Public License
// along with OPAL. If not, see <https://www.gnu.org/licenses/>.
//
#ifndef OPALX_Physics_HH
#define OPALX_Physics_HH
 
#ifdef OPALX_USE_KOKKOS_MATH_CONSTANTS
#include <Kokkos_MathematicalConstants.hpp>
#endif
 
namespace Physics {
 
#ifdef OPALX_USE_KOKKOS_MATH_CONSTANTS
    constexpr double pi = Kokkos::numbers::pi_v<double>;
#else
    constexpr double pi = 3.14159265358979323846;
#endif
 
    constexpr double two_pi = 2 * pi;
 
    constexpr double u_two_pi = 1.0 / two_pi;
 
#ifdef OPALX_USE_KOKKOS_MATH_CONSTANTS
    constexpr double e = Kokkos::numbers::e_v<double>;
#else
    constexpr double e = 2.7182818284590452354;
#endif
 
#ifdef OPALX_USE_KOKKOS_MATH_CONSTANTS
    constexpr double log10e = Kokkos::numbers::log10e_v<double>;
#else
    constexpr double log10e = 0.43429448190325182765;
#endif
 
    constexpr double c = 299792458.0;
 
    constexpr double mu_0 = 1.25663706212e-06;
 
    constexpr double epsilon_0 = 8.8541878128e-12;
 
    constexpr double h_bar = 6.582119569e-25;
 
    constexpr double Avo = 6.02214076e23;
 
    constexpr double kB = 8.617333262e-5;
 
    constexpr double E_ryd = 13.605693122994e-9;
 
    constexpr double a0 = 5.29177210903e-11;
 
    constexpr double q_e = 1.602176634e-19;
 
    constexpr double alpha = 7.2973525693e-03;
 
    constexpr double amu = 0.93149410242;
 
    constexpr double m_e = 0.51099895000e-03;
 
    constexpr double r_e = 2.8179403262e-15;
 
    constexpr double lamda_e = 3.8615926796e-13;
 
    constexpr double a_e = 1.15965218128e-03;
 
    constexpr double m_p = 0.93827208816e+00;
 
    constexpr double r_p = 1.53469857e-18;
 
    constexpr double lamda_p = 2.10308910336e-16;
 
    constexpr double a_p = 1.792847386e+00;
 
    constexpr double z_p = 1;
 
    constexpr double m_c = 11.9967074146787 * amu;
 
    constexpr double m_hm = 1.00837 * amu;
 
    constexpr double m_u = 237.999501 * amu;
 
    constexpr double m_mu = 0.1056583755;
 
    constexpr double tau_mu = 2.1969811e-6;
 
    constexpr double a_mu = 1.16592061e-03;
 
    constexpr double m_pi = 0.13957039;
 
    constexpr double tau_pi = 2.6033e-8;
 
    constexpr double m_d = 2.013553212745 * amu;
 
    constexpr double m_xe = 128.87494026 * amu;
 
    constexpr double m_alpha = 4.001506179127 * amu;
 
    constexpr double m_h = 1.00782503224 * amu;
 
    constexpr double m_h2p = 2.01510 * amu;
 
    constexpr double m_h3p = 3.02293 * amu;
 
    constexpr double PMASS = 1.67262192369e-27;  // kg
 
    constexpr double EMASS = 9.1093837015e-31;  // kg
 
    constexpr double PCHARGE = 1.602176634e-19;  // C
 
    // alfven current
    constexpr double Ia = 17.045148e+03;
    // e/m - electron charge to mass quotient
    constexpr double e0m = 1.75882001076e+11;
    // e/mc
    constexpr double e0mc = e0m / c;
};  // namespace Physics
 
#endif  // OPALX_Physics_HH