HiggsChiral.h

Go to the documentation of this file.

/*
 * Copyright (C) 2017 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef HIGGSCHIRAL_H
#define HIGGSCHIRAL_H
#include "NPbase.h"
 
class HiggsChiral : public NPbase {
public:
 
    static const int NHChiralvars = 18; 
 
    static const std::string HChiralvars[NHChiralvars];
 
    HiggsChiral();
 
    virtual ~HiggsChiral()
    {};
    
    virtual bool PostUpdate();
    
    virtual bool setFlag(const std::string name, const bool value);
 
    double getcv() const
    {
        return cv;
    }
 
    double getct() const
    {
        return ct;
    }
 
    double getcb() const
    {
        return cb;
    }
    
    double getcc() const
    {
        return cc;
    }
    
    double getcs() const
    {
        return cs;
    }
 
 
    double getctau() const
    {
        return ctau;
    }
 
    double getcmu() const
    {
        return cmu;
    }
 
    double getcg() const
    {
        return cg;
    }
 
    double getcga() const
    {
        return cga;
    }
 
    double getcZga() const
    {
        return cZga;
    }
 
    double getobsZgaLimitATLAS13() const
    {
        return obsZgaLimitATLAS13;
    }
 
    double getobsZgaLimitCMS13() const
    {
        return obsZgaLimitCMS13;
    }
 
    double getobsZgaLimitATLAS() const
    {
        return obsZgaLimitATLAS;
    }
 
    double getobsZgaLimitCMS() const
    {
        return obsZgaLimitCMS;
    }
 
    double getexpZgaLimitATLAS13() const
    {
        return expZgaLimitATLAS13;
    }
 
    double getexpZgaLimitCMS13() const
    {
        return expZgaLimitCMS13;
    }
 
    double getexpZgaLimitATLAS() const
    {
        return expZgaLimitATLAS;
    }
 
    double getexpZgaLimitCMS() const
    {
        return expZgaLimitCMS;
    }
 
    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 muWH(const double sqrt_s) const;
    virtual const double muWHpT250(const double sqrt_s) const
    {
        return muWH(sqrt_s); 
    }
    virtual const double muZH(const double sqrt_s) const;
    virtual const double muZHpT250(const double sqrt_s) const
    {
        return muZH(sqrt_s); 
    }
    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
    {
        return muVH(sqrt_s); 
    }
    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 Gammagg() const;
    virtual const double GammaWW() const;
    virtual const double GammaZZ() const;
    virtual const double GammaZga() const;
    virtual const double Gammagaga() const;
    virtual const double Gammamumu() const;
    virtual const double Gammatautau() const;
    virtual const double Gammacc() const;
    virtual const double Gammass() const;
    virtual const double Gammabb() const;
    virtual const double GammaTotal() 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 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();
    }
    
 
    virtual const double BrHtoinvRatio() 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 muppHmumu(const double sqrt_s) const;
 
    virtual const double computeGammaTotalRatio() const;
 
    virtual const double UpperLimitZgammaA13(const double sqrt_s) const;
 
    virtual const double UpperLimitZgammaC13(const double sqrt_s) const;
 
    virtual const double UpperLimitZgammaA(const double sqrt_s) const;
 
    virtual const double UpperLimitZgammaC(const double sqrt_s) const;
 
    virtual const double cgplusct() const;
 
    virtual const double cgaplusct() const;
 
    virtual const double cgminuscga() const;
 
    virtual const double cVpluscb() const;
 
    virtual const double cVplusctau() const;
 
    virtual const double cbminuscc() const;
 
    virtual const double cbminusctau() const;
 
    virtual const double ccminusctau() const;
    
    
//-----------------------------------------------------------------------------------------
//-- 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 kappaseff() 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 computecg() const;
 
    virtual const double computecV() const;
 
    virtual const double computecZga() const;
 
    virtual const double computecgaga() const;
    
    virtual const double computecmu() const;
 
    virtual const double computectau() const;
 
    virtual const double computecc() const;
    
    virtual const double computecs() const;
 
    virtual const double computect() const;
 
    virtual const double computecb() 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;
 
    double cv; 
    double ct; 
    double cb; 
    double cc; 
    double cs; 
    double ctau; 
    double cmu; 
    double cg; 
    double cga; 
    double cZga; 
    double obsZgaLimitATLAS13; 
    double obsZgaLimitCMS13; 
    double obsZgaLimitATLAS; 
    double obsZgaLimitCMS; 
    double expZgaLimitATLAS13; 
    double expZgaLimitCMS13; 
    double expZgaLimitATLAS; 
    double expZgaLimitCMS; 
    
    double cg_loop;
    double cga_loop;
    double cZga_loop;
    bool loopComputed;
    
    bool FlagUniversalcf; 
    bool FlagUniversalcvcf; 
 
};
 
#endif  /* HIGGSCHIRAL_H */