NPHiggs.cpp

Go to the documentation of this file.

/*
 * Copyright (C) 2013 HEPfit Collaboration
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */

#include "NPHiggs.h"
#include <stdexcept>


const std::string NPHiggs::NPHIGGSvars[NNPHIGGSvars]
        = {"a", "b", "c_u", "c_d", "c_e", "d_3", "d_4", "LambdaNP"};

NPHiggs::NPHiggs()
: NPbase(), LambdaNP_in(0.0)
{
    ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("a", boost::cref(a)));
    ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("b", boost::cref(b)));
    ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("c_u", boost::cref(c_u)));
    ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("c_d", boost::cref(c_d)));
    ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("c_e", boost::cref(c_e)));
    ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("d_3", boost::cref(d_3)));
    ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("d_4", boost::cref(d_4)));
    ModelParamMap.insert(std::pair<std::string, boost::reference_wrapper<const double> >("LambdaNP", boost::cref(LambdaNP_in)));

}

void NPHiggs::setParameter(const std::string name, const double& value)
{
    if (name.compare("a") == 0)
        a = value;
    else if (name.compare("b") == 0)
        b = value;
    else if (name.compare("c_u") == 0)
        c_u = value;
    else if (name.compare("c_d") == 0)
        c_d = value;
    else if (name.compare("c_e") == 0)
        c_e = value;
    else if (name.compare("d_3") == 0)
        d_3 = value;
    else if (name.compare("d_4") == 0)
        d_4 = value;
    else if (name.compare("LambdaNP") == 0)
        LambdaNP_in = value;
    else
        NPbase::setParameter(name, value);
}

bool NPHiggs::CheckParameters(const std::map<std::string, double>& DPars)
{
    for (int i = 0; i < NNPHIGGSvars; i++) {
        if (DPars.find(NPHIGGSvars[i]) == DPars.end()) {
            std::cout << "ERROR: Missing mandatory NPHiggs parameter "
                    << NPHIGGSvars[i] << std::endl;
            return false;
        }
    }
    return (NPbase::CheckParameters(DPars));
}


double NPHiggs::obliqueS() const
{
    double Lambda;
    if (LambdaNP_in != 0.0)
        Lambda = LambdaNP_in;
    else if (fabs(1.0 - a * a) < pow(10.0, -32.0))
        Lambda = pow(10.0, 19.0);
    else
        Lambda = 4.0 * M_PI * v() / sqrt(fabs(1.0 - a * a));

    return ( 1.0 / 12.0 / M_PI * (1.0 - a * a) * log(Lambda * Lambda / mHl / mHl));
}

double NPHiggs::obliqueT() const
{
    double Lambda;
    double cW2_SM = trueSM.cW2(); /* This has to be the SM value. */
    if (LambdaNP_in != 0.0)
        Lambda = LambdaNP_in;
    else if (fabs(1.0 - a * a) < pow(10.0, -32.0))
        Lambda = pow(10.0, 19.0);
    else
        Lambda = 4.0 * M_PI * v() / sqrt(fabs(1.0 - a * a));

    return ( -3.0 / 16.0 / M_PI / cW2_SM * (1.0 - a * a) * log(Lambda * Lambda / mHl / mHl));
}

double NPHiggs::obliqueU() const
{
    return 0.0;
}