THDMW.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2017 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef THDMW_H
#define THDMW_H
 
#include "StandardModel.h"
#include "THDMWMatching.h"
#include "NPbase.h"
 
class THDMWcache; //forward reference to THDMWcache class
 
class THDMW: public NPbase {
public:
 
    static const int NTHDMWvars = 34;
    static const std::string THDMWvars[NTHDMWvars];
    
    THDMW();
    
    ~THDMW();
    
    virtual bool InitializeModel();
    
    virtual bool Init(const std::map<std::string, double>& DPars);
    
    virtual bool PreUpdate();
    
    virtual bool Update(const std::map<std::string, double>& DPars);
    
    virtual bool PostUpdate();
    
    virtual bool CheckParameters(const std::map<std::string, double>& DPars);
 
    virtual THDMWMatching& getMatching() const
    {
        return THDMWM.getObj();
    }
 
    // Flags
 
    virtual bool setFlagStr(const std::string name, const std::string value);
 
    THDMWcache* getMyTHDMWCache() const
    {
        return myTHDMWcache;
    }
 
 
    std::string getModelTypeTHDMWflag() const {
        return flag_modelTHDMW;
    }
 
    std::string getRGEorderflag() const {
        return flag_RGEorder;
    }
 
    const double getTHDMW_logtb() const {
        return THDMW_logtb;
    }
 
    const double getTHDMW_tanb() const {
        return THDMW_tanb;
    }
 
    const double getTHDMW_sinb() const {
        return THDMW_sinb;
    }
 
    const double getTHDMW_cosb() const {
        return THDMW_cosb;
    }
 
    const double getTHDMW_bma() const {
        return THDMW_bma;
    }
 
    const double getTHDMW_sin_ba() const {
        return THDMW_sin_ba;
    }
 
    const double getTHDMW_cosa() const{
        return cos(atan(pow(10.,THDMW_logtb))-THDMW_bma);
    }
 
    const double getTHDMW_sina() const{
        return sin(atan(pow(10.,THDMW_logtb))-THDMW_bma);
    }
 
    const double getTHDMW_lambda1() const {
        return THDMW_lambda1;
    }
 
    const double getTHDMW_lambda2() const {
        if (flag_modelTHDMW == "ManoharWise" || flag_modelTHDMW == "custodialMW") {
            return 0.0;
        }
        else {
            return THDMW_lambda2;
        }
    }
 
    const double getTHDMW_lambda3() const {
        if (flag_modelTHDMW == "ManoharWise" || flag_modelTHDMW == "custodialMW") {
            return 0.0;
        }
        else {
            return THDMW_lambda3;
        }
    }
 
    const double getTHDMW_lambda4() const {
        if (flag_modelTHDMW == "ManoharWise" || flag_modelTHDMW == "custodialMW") {
            return 0.0;
        }
        else {
            return THDMW_lambda4;
        }
    }
 
    const double getTHDMW_lambda5() const {
        if (flag_modelTHDMW == "custodial1") {
            return THDMW_lambda4;
        }
        else if (flag_modelTHDMW == "ManoharWise" || flag_modelTHDMW == "custodialMW") {
            return 0.0;
        }
        else {
            return THDMW_lambda5;
        }
    }
 
    const double getTHDMW_mS2() const {
        return THDMW_mS2;
    }
 
    const double getTHDMW_mu1() const {
        return THDMW_mu1;
    }
 
    const double getTHDMW_mu2() const {
        if (flag_modelTHDMW == "custodial1" || flag_modelTHDMW == "custodialMW") {
            return THDMW_mu1;
        }
        else {
            return THDMW_mu2;
        }
    }
 
    const double getTHDMW_mu3() const {
        return THDMW_mu3;
    }
    const double getTHDMW_mu4() const {
        return THDMW_mu4;
    }
 
    const double getTHDMW_mu5() const {
        if (flag_modelTHDMW == "custodial1" || flag_modelTHDMW == "custodialMW") {
            return THDMW_mu4;
        }
        else {
            return THDMW_mu5;
        }
    }
 
    const double getTHDMW_mu6() const {
        if (flag_modelTHDMW == "custodial1" || flag_modelTHDMW == "custodialMW") {
            return 2.0*THDMW_mu1;
        }
        else {
            return THDMW_mu6;
        }
    }
 
    const double getTHDMW_nu1() const {
        return THDMW_nu1;
    }
 
    const double getTHDMW_nu2() const {
        return THDMW_nu2;
    }
 
    const double getTHDMW_nu3() const {
        if (flag_modelTHDMW == "custodial1" || flag_modelTHDMW == "custodialMW") {
            return 0.5*THDMW_nu2;
        }
        else {
            return THDMW_nu3;
        }
    }
    const double getTHDMW_nu4() const {
        return THDMW_nu4;
    }
 
    const double getTHDMW_nu5() const {
        if (flag_modelTHDMW == "custodial1" || flag_modelTHDMW == "custodialMW") {
            return THDMW_nu4;
        }
        else {
            return THDMW_nu5;
        }
    }
 
    const double getTHDMW_omega1() const {
        if (flag_modelTHDMW == "ManoharWise" || flag_modelTHDMW == "custodialMW") {
            return 0.0;
        }
        else {
            return THDMW_omega1;
        }
    }
 
    const double getTHDMW_omega2() const {
        if (flag_modelTHDMW == "ManoharWise" || flag_modelTHDMW == "custodialMW") {
            return 0.0;
        }
        else {
            return THDMW_omega2;
        }
    }
 
    const double getTHDMW_omega3() const {
        if (flag_modelTHDMW == "custodial1") {
            return 0.5*THDMW_omega2;
        }
        else if (flag_modelTHDMW == "ManoharWise" || flag_modelTHDMW == "custodialMW") {
            return 0.0;
        }
        else {
            return THDMW_omega3;
        }
    }
 
    const double getTHDMW_omega4() const {
        if (flag_modelTHDMW == "ManoharWise" || flag_modelTHDMW == "custodialMW") {
            return 0.0;
        }
        else {
            return THDMW_omega4;
        }
    }
 
    const double getTHDMW_omega5() const {
        if (flag_modelTHDMW == "custodial1") {
            return THDMW_omega4;
        }
        else if (flag_modelTHDMW == "ManoharWise" || flag_modelTHDMW == "custodialMW") {
            return 0.0;
        }
        else {
            return THDMW_omega5;
        }
    }
 
    const double getTHDMW_kappa1() const {
        if (flag_modelTHDMW == "ManoharWise" || flag_modelTHDMW == "custodialMW") {
            return 0.0;
        }
        else {
            return THDMW_kappa1;
        }
    }
 
    const double getTHDMW_kappa2() const {
        if (flag_modelTHDMW == "ManoharWise" || flag_modelTHDMW == "custodialMW") {
            return 0.0;
        }
        else {
            return THDMW_kappa2;
        }
    }
 
    const double getTHDMW_kappa3() const {
        if (flag_modelTHDMW == "custodial1") {
            return THDMW_kappa2;
        }
        else if (flag_modelTHDMW == "ManoharWise" || flag_modelTHDMW == "custodialMW") {
            return 0.0;
        }
        else {
            return THDMW_kappa3;
        }
    }
 
    const double getTHDMW_etaU() const {
        return THDMW_etaU;
    }
 
    const double getTHDMW_etaD() const {
        return THDMW_etaD;
    }
    
    
    const double getTHDMW_rho_b() const {
        return THDMW_rho_b;
    }
 
    const double getTHDMW_S_b() const {
        return THDMW_S_b;
    }
    
    
 
    const double getQ_THDMW() const {
        return Q_THDMW;
    }
 
    const double getRpepsTHDMW() const {
        return RpepsTHDMW;
    }
 
    const double getNLOuniscaleTHDMW() const {
        return NLOuniscaleTHDMW;
    }
 
    virtual const double muggH(const double sqrt_s) const;
    virtual const double muVBF(const double sqrt_s) const;
    virtual const double mueeWBF(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual const double muWH(const double sqrt_s) const;
    virtual const double muZH(const double sqrt_s) const;
    virtual const double mueeZH(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual const double muVH(const double sqrt_s) const;
    virtual const double muVBFpVH(const double sqrt_s) const;
    virtual const double muttH(const double sqrt_s) const;
    virtual const double GammaTotal() const;
    virtual const double BrHggRatio() const;
    virtual const double BrHWWRatio() const;
    virtual const double BrHZZRatio() const;
    virtual const double BrHZgaRatio() const;
    virtual const double BrHgagaRatio() const;
    virtual const double BrHmumuRatio() const;
    virtual const double BrHtautauRatio() const;
    virtual const double BrHccRatio() const;
    virtual const double BrHbbRatio() const;
    virtual const double muggHgaga(const double sqrt_s) const;
    virtual const double muVBFHgaga(const double sqrt_s) const;
    virtual const double muVHgaga(const double sqrt_s) const;
    virtual const double muttHgaga(const double sqrt_s) const;
    virtual const double muggHZZ(const double sqrt_s) const;
    virtual const double muVBFHZZ(const double sqrt_s) const;
    virtual const double muVHZZ(const double sqrt_s) const;
    virtual const double muttHZZ(const double sqrt_s) const;
    virtual const double muggHWW(const double sqrt_s) const;
    virtual const double muVBFHWW(const double sqrt_s) const;
    virtual const double muVHWW(const double sqrt_s) const;
    virtual const double muttHWW(const double sqrt_s) const;
    virtual const double muggHtautau(const double sqrt_s) const;
    virtual const double muVBFHtautau(const double sqrt_s) const;
    virtual const double muVHtautau(const double sqrt_s) const;
    virtual const double muttHtautau(const double sqrt_s) const;
    virtual const double muggHbb(const double sqrt_s) const;
    virtual const double muVBFHbb(const double sqrt_s) const;
    virtual const double muVHbb(const double sqrt_s) const;
    virtual const double muttHbb(const double sqrt_s) const;
    virtual const double muppHmumu(const double sqrt_s) const;
    virtual const double muppHZga(const double sqrt_s) const;
    virtual const double computeGammaTotalRatio() const;
 
    /*
     * @brief Method overriding the NPbase Mw with the tree-level W mass
     */
    virtual const double Mw() const;
 
protected:
 
    virtual void setParameter(const std::string, const double&);
 
    bool checkmodelTypeTHDMW(const std::string modeltype) const
    {
        if (modeltype.compare("custodial1") == 0)
            return true;
        else if (modeltype.compare("ManoharWise") == 0)
            return true;
        else if (modeltype.compare("custodialMW") == 0)
            return true;
        else if (modeltype.compare("custodial2") == 0)
            return true;
        else
            return false;
    }
 
    mutable Matching<THDMWMatching,THDMW> THDMWM; 
 
    bool checkRGEorder(const std::string RGEorder) const
    {
        if (RGEorder.compare("LO") == 0
                || RGEorder.compare("approxNLO") == 0)
            return true;
        else
            return false;
    }
 
private:
 
    THDMWcache* myTHDMWcache;
 
    double THDMW_logtb, THDMW_tanb, THDMW_sinb, THDMW_cosb, THDMW_bma, THDMW_sin_ba;
    double THDMW_lambda1, THDMW_lambda2, THDMW_lambda3, THDMW_lambda4, THDMW_lambda5;
    double THDMW_mS2, THDMW_mu1, THDMW_mu2, THDMW_mu3, THDMW_mu4, THDMW_mu5, THDMW_mu6;
    double THDMW_nu1, THDMW_nu2, THDMW_nu3, THDMW_nu4, THDMW_nu5;
    double THDMW_omega1, THDMW_omega2, THDMW_omega3, THDMW_omega4, THDMW_omega5;
    double THDMW_kappa1, THDMW_kappa2, THDMW_kappa3;
    double THDMW_etaU, THDMW_etaD;
    double THDMW_rho_b, THDMW_S_b;
    double Q_THDMW, RpepsTHDMW, NLOuniscaleTHDMW;
    std::string flag_modelTHDMW, flag_RGEorder;
};
 
#endif  /* THDMW_H */