HiggsKigen.h

Go to the documentation of this file.

/*
 * Copyright (C) 2014 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef HIGGSKIGEN_H
#define HIGGSKIGEN_H
#include "NPbase.h"
 
class HiggsKigen : public NPbase {
public:
 
    static const int NHKvKfgenvars = 94; 
 
    static const std::string HKvKfgenvars[NHKvKfgenvars];
 
    HiggsKigen();
 
    virtual ~HiggsKigen()
    {
    };
    
    virtual bool PostUpdate();
    
    virtual bool setFlag(const std::string name, const bool value);
 
    double getKw() const
    {
        return Kw;
    }
 
    void setKw(double Kw)
    {
        this->Kw = Kw;
    }
       
    double getKz() const
    {
        return Kz;
    }
 
    void setKz(double Kz)
    {
        this->Kz = Kz;
    }
      
    double getKg() const
    {
        return Kg;
    }
 
    void setKg(double Kg)
    {
        this->Kg = Kg;
    }
      
    double getKga() const
    {
        return Kga;
    }
 
    void setKga(double Kga)
    {
        this->Kga = Kga;
    }
 
    double getKzga() const
    {
        return Kzga;
    }
 
    void setKzga(double Kzga)
    {
        this->Kzga = Kzga;
    }
 
    double getKd() const
    {
        return Kd;
    }
 
    void setKd(double Kd)
    {
        this->Kd = Kd;
    }
        
    double getKs() const
    {
        return Ks;
    }
 
    void setKs(double Ks)
    {
        this->Ks = Ks;
    }    
    
    double getKb() const
    {
        return Kb;
    }
 
    void setKb(double Kb)
    {
        this->Kb = Kb;
    }
 
    double getKe() const
    {
        return Ke;
    }
 
    void setKe(double Ke)
    {
        this->Ke = Ke;
    }
    
    double getKmu() const
    {
        return Kmu;
    }
 
    void setKmu(double Kmu)
    {
        this->Kmu = Kmu;
    }
    
    double getKtau() const
    {
        return Ktau;
    }
 
    void setKtau(double Ktau)
    {
        this->Ktau = Ktau;
    }
 
    double getKu() const
    {
        return Ku;
    }
 
    void setKu(double Ku)
    {
        this->Ku = Ku;
    }
    
    double getKc() const
    {
        return Kc;
    }
 
    void setKc(double Kc)
    {
        this->Kc = Kc;
    }    
    
    double getKt() const
    {
        return Kt;
    }
 
    void setKt(double Kt)
    {
        this->Kt = Kt;
    }   
    
    double getKH() const
    {
        return KH;
    }
 
    void setKH(double KH)
    {
        this->KH = KH;
    }
    
    
    double getBrHinv() const
    {
        return BrHinv;
    }
 
    void setBrHinv(double BrHinv)
    {
        this->BrHinv = BrHinv;
    } 
    
    double getBrHexo() const
    {
        return BrHexo;
    }
    
    void setBrHexo(double BrHexo)
    {
        this->BrHexo = BrHexo;
    }
 
    virtual bool CheckParameters(const std::map<std::string, double>& DPars);
 
 
    virtual const double obliqueS() const;
 
    virtual const double obliqueT() const;
 
    virtual const double obliqueU() const;
 
    virtual const double muggH(const double sqrt_s) const;
    virtual const double muVBF(const double sqrt_s) const;
    virtual const double muVBFgamma(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 mueeHvv(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual const double mueeZBF(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual const double muepWBF(const double sqrt_s) const;
    virtual const double muepZBF(const double sqrt_s) const;
    virtual const double muWH(const double sqrt_s) const;
    virtual const double muWHpT250(const double sqrt_s) const;
    virtual const double muZH(const double sqrt_s) const;
    virtual const double muZHpT250(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 mueeZllH(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual const double mueeZqqH(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 muVHpT250(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 mutHq(const double sqrt_s) const;
    virtual const double mummH(const double sqrt_s) const;
    virtual const double mummHNWA(const double sqrt_s) const; 
    virtual const double mummZH(const double sqrt_s) const;
    virtual const double mummHvv(const double sqrt_s) const;
    virtual const double mummHmm(const double sqrt_s) const;
    virtual const double mummttH(const double sqrt_s) const;
    virtual const double muggHpttH(const double sqrt_s) const;
    virtual const double mueettH(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual const double muttHZbbboost(const double sqrt_s) const;
    virtual const double muttHgagaZeeboost(const double sqrt_s) const;
 
    
    
    virtual const double muggHgaga(const double sqrt_s) const;    
    virtual const double muggHgagaInt(const double sqrt_s) const;
    virtual const double muVBFHgaga(const double sqrt_s) const;
    virtual const double muZHgaga(const double sqrt_s) const;
    virtual const double muWHgaga(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 muggHZga(const double sqrt_s) const;
    virtual const double muVBFHZga(const double sqrt_s) const;    
    virtual const double muZHZga(const double sqrt_s) const;
    virtual const double muWHZga(const double sqrt_s) const;    
    virtual const double muVHZga(const double sqrt_s) const;
    virtual const double muttHZga(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 muZHZZ(const double sqrt_s) const;
    virtual const double muWHZZ(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 muggHZZ4l(const double sqrt_s) const;
    virtual const double muVBFHZZ4l(const double sqrt_s) const;    
    virtual const double muZHZZ4l(const double sqrt_s) const;
    virtual const double muWHZZ4l(const double sqrt_s) const;    
    virtual const double muVHZZ4l(const double sqrt_s) const;
    virtual const double muttHZZ4l(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 muZHWW(const double sqrt_s) const;
    virtual const double muWHWW(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 muggHWW2l2v(const double sqrt_s) const;
    virtual const double muVBFHWW2l2v(const double sqrt_s) const;   
    virtual const double muZHWW2l2v(const double sqrt_s) const;
    virtual const double muWHWW2l2v(const double sqrt_s) const;    
    virtual const double muVHWW2l2v(const double sqrt_s) const;
    virtual const double muttHWW2l2v(const double sqrt_s) const;
    virtual const double muggHmumu(const double sqrt_s) const;
    virtual const double muVBFHmumu(const double sqrt_s) const;    
    virtual const double muZHmumu(const double sqrt_s) const;
    virtual const double muWHmumu(const double sqrt_s) const;    
    virtual const double muVHmumu(const double sqrt_s) const;
    virtual const double muttHmumu(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 muZHtautau(const double sqrt_s) const;
    virtual const double muWHtautau(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 muZHbb(const double sqrt_s) const;
    virtual const double muWHbb(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 computeGammaTotalRatio() const;
    
    virtual const double computeGammaTotalnoErrorRatio() const;
    
    virtual const double computeGammaTotalSMRatio() const;
    
    virtual const double computeGammaTotalSMnoErrorRatio() const;
    
    
    virtual const double BrHggRatio() const;
    virtual const double BrHWWRatio() const;
    virtual const double BrHZZRatio() const;
    virtual const double BrHVVRatio() const;
    virtual const double BrHZgaRatio() const;
    virtual const double BrHZgallRatio() const;
    virtual const double BrHZgaeeRatio() const;
    virtual const double BrHZgamumuRatio() const;
    virtual const double BrHgagaRatio() const;
    virtual const double BrHmumuRatio() const;
    virtual const double BrHtautauRatio() const;
    virtual const double BrHccRatio() const;
    virtual const double BrHssRatio() const;
    virtual const double BrHbbRatio() const;
    
    virtual const double Br_H_exo() const;
    
    virtual const double Br_H_inv() const;
    
    virtual const double Br_H_inv_NP() const;
    
    virtual const double BrHvisRatio() const;
    
    virtual const double BrHtoinvRatio() const;
    
    
    virtual const double BrH2L2LRatio() const;
    
    virtual const double BrH2e2muRatio() const;    
 
    virtual const double BrH2v2vRatio() const;
 
    virtual const double BrH2L2vRatio() const;
    
    virtual const double BrH2L2v2Ratio() const;
    
    virtual const double BrH2e2vRatio() const;
    
    virtual const double BrH2mu2vRatio() const;
    
    virtual const double BrH2u2uRatio() const;
    
    virtual const double BrH2d2dRatio() const;
    
    virtual const double BrH2u2dRatio() const;
    
    virtual const double BrH2L2uRatio() const;
    
    virtual const double BrH2L2dRatio() const;
    
    virtual const double BrH2v2uRatio() const;
    
    virtual const double BrH2v2dRatio() const;
    
    virtual const double BrH4LRatio() const;
    
    virtual const double BrH4L2Ratio() const;
    
    virtual const double BrH4eRatio() const;
    
    virtual const double BrH4muRatio() const;
    
    virtual const double BrH4vRatio() const;
    
    virtual const double BrH4uRatio() const;
    
    virtual const double BrH4dRatio() const;
    
    virtual const double BrHLvvLRatio() const;
    
    virtual const double BrHevmuvRatio() const;
    
    virtual const double BrHudduRatio() const;
    
    virtual const double BrHLvudRatio() const;
    
    virtual const double BrH2udRatio() const;
    
    virtual const double BrH2LvRatio() const;   
    
    virtual const double BrH2Lv2Ratio() const;
    
    virtual const double BrH2evRatio() const;    
    
    virtual const double BrH2muvRatio() const;    
    
    virtual const double BrH4fRatio() const;
    
    // DECAYS INVOLVING ONLY ELECTRONS, MUONS OR NEUTRINOS IN THE FINAL STATES 
 
    virtual const double BrH4lRatio() const;
    
    virtual const double BrH2l2vRatio() const;
    
    
    virtual const double BrHlljjRatio() const
    {
        return BrHZZRatio(); 
    }
    
    virtual const double BrHlvjjRatio() const
    {
        return BrHWWRatio(); 
    }
    
    virtual const double BrHlv_lvorjjRatio() const
    {
        return BrHWWRatio();
    }   
    
    virtual const double BrHll_vvorjjRatio() const
    {
        return BrHZZRatio();
    }
    
//-----------------------------------------------------------------------------------------
//-- Special Hadron collider signal strengths with separate full TH unc U(prod x decay) ---
//-----------------------------------------------------------------------------------------
    
    virtual const double muTHUggHgaga(const double sqrt_s) const;    
    virtual const double muTHUVBFHgaga(const double sqrt_s) const;
    virtual const double muTHUZHgaga(const double sqrt_s) const;
    virtual const double muTHUWHgaga(const double sqrt_s) const;
    virtual const double muTHUVHgaga(const double sqrt_s) const;
    virtual const double muTHUttHgaga(const double sqrt_s) const;
    virtual const double muTHUggHZga(const double sqrt_s) const;
    virtual const double muTHUVBFHZga(const double sqrt_s) const;    
    virtual const double muTHUZHZga(const double sqrt_s) const;
    virtual const double muTHUWHZga(const double sqrt_s) const;    
    virtual const double muTHUVHZga(const double sqrt_s) const;
    virtual const double muTHUttHZga(const double sqrt_s) const; 
    virtual const double muTHUggHZZ(const double sqrt_s) const;
    virtual const double muTHUVBFHZZ(const double sqrt_s) const;
    virtual const double muTHUZHZZ(const double sqrt_s) const;
    virtual const double muTHUWHZZ(const double sqrt_s) const;
    virtual const double muTHUVHZZ(const double sqrt_s) const;
    virtual const double muTHUttHZZ(const double sqrt_s) const;
    
    virtual const double muTHUggHZZ4l(const double sqrt_s) const;
    virtual const double muTHUVBFHZZ4l(const double sqrt_s) const;    
    virtual const double muTHUZHZZ4l(const double sqrt_s) const;
    virtual const double muTHUWHZZ4l(const double sqrt_s) const;    
    virtual const double muTHUVHZZ4l(const double sqrt_s) const;
    virtual const double muTHUttHZZ4l(const double sqrt_s) const;
    
    virtual const double muTHUggHWW(const double sqrt_s) const;
    virtual const double muTHUVBFHWW(const double sqrt_s) const;
    virtual const double muTHUZHWW(const double sqrt_s) const;
    virtual const double muTHUWHWW(const double sqrt_s) const; 
    virtual const double muTHUVHWW(const double sqrt_s) const;
    virtual const double muTHUttHWW(const double sqrt_s) const;
    virtual const double muTHUggHWW2l2v(const double sqrt_s) const;
    virtual const double muTHUVBFHWW2l2v(const double sqrt_s) const;   
    virtual const double muTHUZHWW2l2v(const double sqrt_s) const;
    virtual const double muTHUWHWW2l2v(const double sqrt_s) const;    
    virtual const double muTHUVHWW2l2v(const double sqrt_s) const;
    virtual const double muTHUttHWW2l2v(const double sqrt_s) const;
    virtual const double muTHUggHmumu(const double sqrt_s) const;
    virtual const double muTHUVBFHmumu(const double sqrt_s) const;    
    virtual const double muTHUZHmumu(const double sqrt_s) const;
    virtual const double muTHUWHmumu(const double sqrt_s) const;    
    virtual const double muTHUVHmumu(const double sqrt_s) const;
    virtual const double muTHUttHmumu(const double sqrt_s) const;
    virtual const double muTHUggHtautau(const double sqrt_s) const;
    virtual const double muTHUVBFHtautau(const double sqrt_s) const;
    virtual const double muTHUZHtautau(const double sqrt_s) const;
    virtual const double muTHUWHtautau(const double sqrt_s) const;
    virtual const double muTHUVHtautau(const double sqrt_s) const;
    virtual const double muTHUttHtautau(const double sqrt_s) const;
    virtual const double muTHUggHbb(const double sqrt_s) const;
    virtual const double muTHUVBFHbb(const double sqrt_s) const;
    virtual const double muTHUZHbb(const double sqrt_s) const;
    virtual const double muTHUWHbb(const double sqrt_s) const;
    virtual const double muTHUVHbb(const double sqrt_s) const;
    virtual const double muTHUttHbb(const double sqrt_s) const;
    
    virtual const double muTHUVBFBRinv(const double sqrt_s) const;
    virtual const double muTHUVBFHinv(const double sqrt_s) const;
    virtual const double muTHUVHBRinv(const double sqrt_s) const;
    virtual const double muTHUVHinv(const double sqrt_s) const;
    
    virtual const double muTHUggHZZ4mu(const double sqrt_s) const;    
    virtual const double muTHUggHZgamumu(const double sqrt_s) const;    
 
 
    virtual const double deltaG_hgg() const;
    virtual const double deltaG_hggRatio() const;
    virtual const double deltaG1_hWW() const;
    virtual const double deltaG2_hWW() const;
    virtual const double deltaG3_hWW() const;
    virtual const double deltaG1_hZZ() const;
    virtual const double deltaG2_hZZ() const;
    virtual const double deltaG3_hZZ() const;
    virtual const double deltaG1_hZA() const;
    virtual const double deltaG1_hZARatio() const;
    virtual const double deltaG2_hZA() const;
    virtual const double deltaG_hAA() const;
    virtual const double deltaG_hAARatio() const;
    // no generation mixing
    virtual gslpp::complex deltaG_hff(const Particle p) const;   
   
    
    virtual const double kappamueff() const;
    
    virtual const double kappataueff() const;
    
    virtual const double kappaceff() const;
    
    virtual const double kappabeff() const;
    
    virtual const double kappaGeff() const;
    
    virtual const double kappaZeff() const;
    
    virtual const double kappaWeff() const;
    
    virtual const double kappaAeff() const;
    
    virtual const double kappaZAeff() const;
 
protected:
 
    virtual void setParameter(const std::string name, const double& value);
 
    virtual const double computeKg() const;
 
    virtual const double computeKW() const;
 
    virtual const double computeKZ() const;
 
    virtual const double computeKZga() const;
 
    virtual const double computeKgaga() const;
    
    virtual const double computeKe() const;
    
    virtual const double computeKmu() const;
 
    virtual const double computeKtau() const;
    
    virtual const double computeKu() const;
 
    virtual const double computeKc() const;
 
    virtual const double computeKt() const;
    
    virtual const double computeKd() const;
    
    virtual const double computeKs() const;
 
    virtual const double computeKb() const;
    
    virtual const double computeKgLoop() const;
    
    virtual const double computeKZgaLoop() const;
 
    virtual const double computeKgagaLoop() const;
 
private:
    gslpp::complex f_func(const double x) const;
    gslpp::complex g_func(const double x) const;
    gslpp::complex Int1(const double tau, const double lambda) const;
    gslpp::complex Int2(const double tau, const double lambda) const;
 
    // Model parameters    
    double Kw; 
    double Kz; 
    double Kg; 
    double Kga; 
    double Kzga; 
    double Ku; 
    double Kc; 
    double Kt; 
    double Kd; 
    double Ks; 
    double Kb; 
    double Ke; 
    double Kmu; 
    double Ktau; 
    double KH; 
    double BrHinv; 
    double BrHexo; 
// The intrinsic and parametric theory relative errors in the Higgs observables.
    double eggFint; 
    double eggFpar; 
    double ettHint; 
    double ettHpar; 
    double eVBFint; 
    double eVBFpar; 
    double eWHint; 
    double eWHpar; 
    double eZHint; 
    double eZHpar; 
    double eeeWBFint; 
    double eeeWBFpar; 
    double eeeZHint; 
    double eeeZHpar; 
    double eeettHint; 
    double eeettHpar; 
    double eepWBFint; 
    double eepWBFpar; 
    double eepZBFint; 
    double eepZBFpar; 
    double eHggint; 
    double eHggpar; 
    double eHWWint; 
    double eHWWpar; 
    double eHZZint; 
    double eHZZpar; 
    double eHZgaint; 
    double eHZgapar; 
    double eHgagaint; 
    double eHgagapar; 
    double eHmumuint; 
    double eHmumupar; 
    double eHtautauint; 
    double eHtautaupar; 
    double eHccint; 
    double eHccpar; 
    double eHssint; 
    double eHsspar; 
    double eHbbint; 
    double eHbbpar; 
    
    double eggFHgaga,eggFHZga,eggFHZZ,eggFHWW,eggFHtautau,eggFHbb,eggFHmumu; 
    double eVBFHgaga,eVBFHZga,eVBFHZZ,eVBFHWW,eVBFHtautau,eVBFHbb,eVBFHmumu; 
    double eWHgaga,eWHZga,eWHZZ,eWHWW,eWHtautau,eWHbb,eWHmumu; 
    double eZHgaga,eZHZga,eZHZZ,eZHWW,eZHtautau,eZHbb,eZHmumu; 
    double ettHgaga,ettHZga,ettHZZ,ettHWW,ettHtautau,ettHbb,ettHmumu; 
    double eVBFHinv,eVHinv; 
    
    double eHwidth; 
    
    // Cache variables
    double GammaHTotR;
    double GammaHTotSMR;
    
    double GammaHTotnoER;
    double GammaHTotSMnoER;
    
    // Flags    
    bool FlagKiLoop; 
    bool FlagCustodial; 
    bool FlagUniversalKf; 
    bool FlagUniversalK; 
    bool FlagUseKH; 
    
};
 
#endif  /* HIGGSKIGEN_H */