ew_input.hpp

Go to the documentation of this file.

// ====================================================================
// This file is part of FlexibleSUSY.
//
// FlexibleSUSY 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.
//
// FlexibleSUSY is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with FlexibleSUSY.  If not, see
// <http://www.gnu.org/licenses/>.
// ====================================================================
 
#ifndef ELECTROWEAK_INPUT_H
#define ELECTROWEAK_INPUT_H
 
#include "wrappers.hpp"
 
#include <cmath>
 
namespace flexiblesusy {
 
namespace Electroweak_constants {
   namespace {
      constexpr double vev = 246.22;
      constexpr double root2 = 1.4142135623730950;
      constexpr double mtoprun = 165;
      constexpr double mbrun = 2.9;
      constexpr double mtau = 1.77699;
      constexpr double yt = mtoprun * root2 / vev;
      constexpr double yb = mbrun * root2 / vev;
      constexpr double ytau = mtau * root2 / vev;
      constexpr double MZ = 91.1876;
      constexpr double Error_MZ = 0.0021; 
      constexpr double MW = 80.385;
      constexpr double MWSM = 80.357; 
      constexpr double MH = 125.25; 
      constexpr double Error_MH = 0.24; 
      constexpr double MUP = 2.4e-3; 
      constexpr double MDOWN = 4.75e-3; 
      constexpr double MSTRANGE = 0.104; 
      constexpr double MCHARM = 1.27; 
      constexpr double MBOTTOM = 4.18; 
      constexpr double MTOP = 165.0; 
      constexpr double MELECTRON = 5.10998902e-4; 
      constexpr double MMUON = 1.056583715e-1; 
      constexpr double MTAU = 1.77699; 
      constexpr double PMTOP = 173.34; 
      constexpr double PMBOTTOM = 4.9; 
      constexpr double aem = 1.0 / 127.916; // at MZ
      constexpr double sinThetaW2 = 0.23122;
      const double sinThetaW = std::sqrt(sinThetaW2);
      constexpr double cosThetaW2 = 1 - sinThetaW2;
      const double cosThetaW = std::sqrt(cosThetaW2);
      constexpr double alpha1 = 5.0 * aem / (3.0 * (1.0 - sinThetaW2));
      constexpr double alpha2 = aem / sinThetaW2;
      constexpr double alpha3 = 0.1184; // at MZ from PDG
      const double e  = std::sqrt(4.0 * Pi * aem);
      const double g1 = std::sqrt(4.0 * Pi * alpha1);
      const double g2 = std::sqrt(4.0 * Pi * alpha2);
      const double g3 = std::sqrt(4.0 * Pi * alpha3);
      constexpr double gYSM = 3.57232027E-01;     
      const double g1SM = std::sqrt(5./3.) * gYSM; 
      constexpr double g2SM = 6.51103848E-01;     
      constexpr double g3SM = 1.21087245E+00;     
      constexpr double CKM_THETA12 = 0.229206; 
      constexpr double CKM_THETA13 = 0.003960; 
      constexpr double CKM_THETA23 = 0.042223; 
      constexpr double CKM_DELTA   = 0.;
      const double PMNS_THETA12 = 0.5 * std::asin(std::sqrt(0.846));
      const double PMNS_THETA13 = 0.5 * std::asin(std::sqrt(0.093));
      const double PMNS_THETA23 = 0.5 * std::asin(std::sqrt(0.999));
      constexpr double PMNS_DELTA   = 0.;
      constexpr double PMNS_ALPHA1  = 0.;
      constexpr double PMNS_ALPHA2  = 0.;
      constexpr double gfermi = 1.1663787e-5; 
      constexpr double yeSM = 2.85784368E-06; 
      constexpr double ymSM = 5.90911374E-04; 
      constexpr double ylSM = 9.95869693E-03; 
      constexpr double yuSM = 7.89527379E-06; 
      constexpr double ycSM = 3.60854291E-03; 
      constexpr double ytSM = 9.76017610E-01; 
      constexpr double ydSM = 1.56989573E-05; 
      constexpr double ysSM = 3.43724539E-04; 
      constexpr double ybSM = 1.64406299E-02; 
      constexpr double mu2SM = 7.67488232E+03; 
      constexpr double lamSM = 2.79613357E-01; 
      constexpr double vSM = 2.48997424E+02;   
      constexpr double delta_alpha_s_5_had = 0.02766; 
   } // namespace
} // namespace Electroweak_constants
 
} // namespace flexiblesusy
 
#endif