file frontends/MicrOmegas_ScalarSingletDM_Z2_3_6_9_2.cpp
[No description available] More…
Detailed Description
Author: Jonathan Cornell
Date: April 2017
Frontend for MicrOmegas ScalarSingletDM_Z2 3.6.9.2 backend
Authors (add name and date if you modify):
Source code
// GAMBIT: Global and Modular BSM Inference Tool
// *********************************************
/// \file
///
/// Frontend for MicrOmegas ScalarSingletDM_Z2 3.6.9.2 backend
///
/// *********************************************
///
/// Authors (add name and date if you modify):
///
/// \author Jonathan Cornell
/// \date April 2017
///
/// *********************************************
#include "gambit/Backends/frontend_macros.hpp"
#include "gambit/Backends/frontends/MicrOmegas_ScalarSingletDM_Z2_3_6_9_2.hpp"
#include <unistd.h>
// Convenience functions (definitions)
BE_NAMESPACE
{
double dNdE(double Ecm, double E, int inP, int outN)
{
// outN 0-5: gamma, e+, p-, nu_e, nu_mu, nu_tau
// inP: 0 - 6: glu, d, u, s, c, b, t
// 7 - 9: e, m, l
// 10 - 15: Z, ZT, ZL, W, WT, WL
double tab[250]; // NZ = 250
// readSpectra() moved to initialization function.
// Must be inside critical block if used here!
// readSpectra();
mInterp(Ecm/2, inP, outN, tab);
return zInterp(log(E/Ecm*2), tab);
}
}
END_BE_NAMESPACE
// Initialisation function (definition)
BE_INI_FUNCTION
{
int error;
char cdmName[10];
const Spectrum& spec = *Dep::ScalarSingletDM_Z2_spectrum;
const SubSpectrum& he = spec.get_HE();
const SubSpectrum& le = spec.get_LE();
const SMInputs& sminputs = spec.get_SMInputs();
int VZdecayOpt, VWdecayOpt; // 0=no 3 body final states
// 1=3 body final states in annihlations
// 2=3 body final states in co-annihilations
VZdecayOpt = runOptions->getValueOrDef<int>(1, "VZdecay");
VWdecayOpt = runOptions->getValueOrDef<int>(1, "VWdecay");
*VZdecay = VZdecayOpt;
*VWdecay = VWdecayOpt;
logger() << LogTags::debug << "Initializing MicrOmegas ScalarSingletDM_Z2 with ";
logger() << "VWdecay: " << VWdecay << " VZdecay: " << VZdecay << EOM;
// Uncomment below to force MicrOmegas to do calculations in unitary gauge
// *ForceUG=1;
double mS = spec.get(Par::Pole_Mass,"S");
double mH = spec.get(Par::Pole_Mass,"h0_1");
double mW = le.get(Par::Pole_Mass, "W+");
double lambda = he.get(Par::dimensionless,"lambda_hS");
//double lambdaS = he.get(Par::dimensionless,"lambda_S");
error = assignVal((char*)"MS", mS);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set DM mass in"
"MicrOmegas. MicrOmegas error code: " + std::to_string(error));
error = assignVal((char*)"lambda", lambda);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set lambda in"
"MicrOmegas. MicrOmegas error code: " + std::to_string(error));
//error = assignVal((char*)"lamS", lambdaS);
//if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set lamS in"
// "MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// Set SM particle masses in micrOmegas_3.6.9.2
// EE = sqrt(4*pi*(1/alphainv))
error = assignVal((char*)"EE", sqrt(4*M_PI*(1/sminputs.alphainv)));
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set electric charge in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// mu(2 GeV) in MSbar scheme
error = assignVal((char*)"Mu", sminputs.mU);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set mU in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// md(2 GeV) in MSbar scheme
error = assignVal((char*)"Md", sminputs.mD);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set mD in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// ms(2 GeV) in MSbar scheme
error = assignVal((char*)"Ms", sminputs.mS);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set mS in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// mc(mc) in MSbar scheme
error = assignVal((char*)"McMc", sminputs.mCmC);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set Mc(Mc) in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// mb(mb) in MSbar scheme
error = assignVal((char*)"MbMb", sminputs.mBmB);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set Mb(Mb) in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// mtop(pole)
error = assignVal((char*)"Mtp", sminputs.mT);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set mT in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// mE(pole)
error = assignVal((char*)"Me", sminputs.mE);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set mE in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// mMu(pole)
error = assignVal((char*)"Mm", sminputs.mMu);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set mMu in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// mTau(pole)
error = assignVal((char*)"Ml", sminputs.mTau);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set mTau in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// mZ(pole)
error = assignVal((char*)"MZ", sminputs.mZ);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set mZ in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// Set Higgs boson mass in micrOmegas_3.6.9.2
error = assignVal((char*)"Mh", mH);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set mH in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// Set sin(theta_W) in micrOmegas_3.6.9.2
error = assignVal((char*)"CW", sqrt(1-mW/(sminputs.mZ)*mW/(sminputs.mZ)));
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set sin(theta_W) in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
// Set widths. Higgs width is needed from the decay table, otherwise error is
// thrown. If widths of the other SM particle are not present, the MicrOmegas
// default values are used.
const DecayTable* tbl = &(*Dep::decay_rates);
double width;
bool present = true;
try {width = tbl->at("h0_1").width_in_GeV;}
catch (std::exception& e) {backend_error().raise(LOCAL_INFO, "Higgs width missing from"
"decay table.");}
error = assignVal((char*)"wh", width);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set Higgs width in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
present = true;
try {width = tbl->at("Z0").width_in_GeV;}
catch (std::exception& e) {present = false;}
if (present)
{
error = assignVal((char*)"wZ", width);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set Z width in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
}
present = true;
try {width = tbl->at("W+").width_in_GeV;}
catch (std::exception& e) {present = false;}
if (present)
{
error = assignVal((char*)"wW", width);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set W width in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
}
present = true;
try {width = tbl->at("mu-").width_in_GeV;}
catch (std::exception& e) {present = false;}
if (present)
{
error = assignVal((char*)"wm", width);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set muon width in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
}
present = true;
try {width = tbl->at("tau-").width_in_GeV;}
catch (std::exception& e) {present = false;}
if (present)
{
error = assignVal((char*)"wl", width);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set tau width in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
}
present = true;
try {width = tbl->at("t").width_in_GeV;}
catch (std::exception& e) {present = false;}
if (present)
{
error = assignVal((char*)"wt", width);
if (error != 0) backend_error().raise(LOCAL_INFO, "Unable to set top width in"
" MicrOmegas. MicrOmegas error code: " + std::to_string(error));
}
// Initialise micrOMEGAs mass spectrum
error = sortOddParticles(byVal(cdmName));
if (error != 0) backend_error().raise(LOCAL_INFO, "MicrOmegas function "
"sortOddParticles returned error code: " + std::to_string(error));
}
END_BE_INI_FUNCTION
Updated on 2023-06-26 at 21:36:57 +0000