NPEpsilons_pureNP.cpp

Go to the documentation of this file.

/*
 * Copyright (C) 2013 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#include "NPEpsilons_pureNP.h"
#include <stdexcept>
 
const std::string NPEpsilons_pureNP::EPSILONpureNPvars[NEPSILONpureNPvars]
        = {"delEps_1", "delEps_2", "delEps_3", "delEps_b"};
 
NPEpsilons_pureNP::NPEpsilons_pureNP()
: NPbase()
{
    ModelParamMap.insert(std::make_pair("delEps_1", std::cref(deltaEps_1)));
    ModelParamMap.insert(std::make_pair("delEps_2", std::cref(deltaEps_2)));
    ModelParamMap.insert(std::make_pair("delEps_3", std::cref(deltaEps_3)));
    ModelParamMap.insert(std::make_pair("delEps_b", std::cref(deltaEps_b)));
}
 
void NPEpsilons_pureNP::setParameter(const std::string name, const double& value)
{
    if (name.compare("delEps_1") == 0)
        deltaEps_1 = value;
    else if (name.compare("delEps_2") == 0)
        deltaEps_2 = value;
    else if (name.compare("delEps_3") == 0)
        deltaEps_3 = value;
    else if (name.compare("delEps_b") == 0)
        deltaEps_b = value;
    else
        NPbase::setParameter(name, value);
}
 
bool NPEpsilons_pureNP::CheckParameters(const std::map<std::string, double>& DPars)
{
    for (int i = 0; i < NEPSILONpureNPvars; i++) {
        if (DPars.find(EPSILONpureNPvars[i]) == DPars.end()) {
            std::cout << "ERROR: Missing mandatory NPEpsilons_pureNP parameter "
                    << EPSILONpureNPvars[i] << std::endl;
            raiseMissingModelParameterCount();
            addMissingModelParameter(EPSILONpureNPvars[i]);
        }
    }
    return (NPbase::CheckParameters(DPars));
}
 
 
double NPEpsilons_pureNP::epsilon1() const
{
    return (trueSM.epsilon1() + deltaEps_1);
}
 
double NPEpsilons_pureNP::epsilon2() const
{
    return (trueSM.epsilon2() + deltaEps_2);
}
 
double NPEpsilons_pureNP::epsilon3() const
{
    return (trueSM.epsilon3() + deltaEps_3);
}
 
double NPEpsilons_pureNP::epsilonb() const
{
    return (trueSM.epsilonb() + deltaEps_b);
}
 
 
 
double NPEpsilons_pureNP::Mw() const
{
    double Mw_SM = trueSM.Mw();
    double sW2_SM = trueSM.sW2();
    double cW2_SM = trueSM.cW2();
    return ( Mw_SM * (1.0 - (-c02() * deltaEps_1
            + (c02() - s02()) * deltaEps_2
            + 2.0 * s02() * deltaEps_3) / (cW2_SM - sW2_SM) / 2.0));
}
 
double NPEpsilons_pureNP::GammaW() const
{
    throw std::runtime_error("NPEpsilons_pureNP::GammaW() is not available.");
}
 
 
 
double NPEpsilons_pureNP::deltaGV_f(const Particle f) const
{
    if (f.is("TOP")) return 0.0;
 
    /* SM values */
    double gV_SM = trueSM.gV_f(f).real();
    double gA_SM = trueSM.gA_f(f).real();
 
    if (f.is("BOTTOM"))
        return ( gV_SM * deltaEps_1 / 2.0 + gV_SM * deltaEps_b
            + (gV_SM - gA_SM)*(deltaEps_3 - c02() * deltaEps_1)
            / (c02() - s02())
            - (gV_SM - gA_SM) * deltaEps_b);
    else
        return ( gV_SM * deltaEps_1 / 2.0
            + (gV_SM - gA_SM)*(deltaEps_3 - c02() * deltaEps_1)
            / (c02() - s02()));
}
 
double NPEpsilons_pureNP::deltaGA_f(const Particle f) const
{
    if (f.is("TOP")) return 0.0;
 
    if (f.is("BOTTOM"))
        return ( f.getIsospin()*(deltaEps_1 / 2.0 + deltaEps_b));
 
    return ( f.getIsospin() * deltaEps_1 / 2.0);
}