NPSMEFTd6.h

Go to the documentation of this file.

/*
 * Copyright (C) 2014 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef NPSMEFTD6_H
#define NPSMEFTD6_H
 
#include "gslpp.h"
#include "NPbase.h"
#include "NPSMEFTd6Matching.h"
#include <string.h>
#include <stdexcept>
#include <gsl/gsl_integration.h>
 
 
class NPSMEFTd6 : public NPbase {
public:
 
    static const int NNPSMEFTd6Vars = 446;
 
    static const std::string NPSMEFTd6Vars[NNPSMEFTd6Vars];
      
    static const std::string NPSMEFTd6VarsRot[NNPSMEFTd6Vars];
 
    static const int NNPSMEFTd6Vars_LFU_QFU = 250;
 
    static const std::string NPSMEFTd6Vars_LFU_QFU[NNPSMEFTd6Vars_LFU_QFU];
      
    static const std::string NPSMEFTd6VarsRot_LFU_QFU[NNPSMEFTd6Vars_LFU_QFU];
 
    NPSMEFTd6(const bool FlagLeptonUniversal_in = false, const bool FlagQuarkUniversal_in = false);
 
    virtual bool PostUpdate();
 
    virtual bool CheckParameters(const std::map<std::string, double>& DPars);
 
    virtual bool setFlag(const std::string name, const bool value);
    
    virtual NPSMEFTd6Matching& getMatching() const
    {
        return NPSMEFTd6M.getObj();
    }
    
 
    double getLambda_NP() const{
        return Lambda_NP;
    }
 
    double getCQe_2311() const{
        return CQe_2311;
    }
    
    double getCQe_2322() const{
        return CQe_2322;
    }
 
    double getCLQ1_1123() const{
        return CLQ1_1123;
    }
    
    double getCLQ1_2223() const{
        return CLQ1_2223;
    }    
 
    double getCLQ3_1123() const{
        return CLQ3_1123;
    }
    
    double getCLQ3_2223() const{
        return CLQ3_2223;
    }
    
    double getCLd_1123() const{
        return CLd_1123;
    }
    
    double getCLd_2223() const{
        return CLd_2223;
    }
    
    double getCed_1123() const{
        return Ced_1123;
    }
    
    double getCed_2223() const{
        return Ced_2223;
    }
    
    double getCHL1_11() const{
        return CHL1_11;
    }
    
    double getCHL1_22() const{
        return CHL1_22;
    }
 
    double getCHL3_11() const{
        return CHL3_11;
    }
    
    double getCHL3_22() const{
        return CHL3_22;
    }
    
    double getCHe_11() const{
        return CHe_11;
    }
        
    double getCHe_22() const{
        return CHe_22;
    }
    
    double getCeu_1133() const{
        return Ceu_1133;
    }
        
    double getCeu_2233() const{
        return Ceu_2233;
    }
    
    double getCLu_1133() const{
        return CLu_1133;
    }
        
    double getCLu_2233() const{
        return CLu_2233;
    }
    
    double getCLedQ_11() const{
        return CLedQ_11;
    }
 
    double getCLedQ_22() const{
        return CLedQ_22;
    }
    
    double getCpLedQ_11() const{
        return CpLedQ_11;
    }
    
    double getCpLedQ_22() const{
        return CpLedQ_22;
    }
 
 
    virtual double DeltaGF() const;
 
    virtual double obliqueS() const;
 
    virtual double obliqueT() const;
 
    virtual double obliqueU() const;
    
    virtual double obliqueW() const;
 
    virtual double obliqueY() const;
    
    
    
    virtual double deltaMz() const;
    
    virtual double deltaMz2() const;
        
    virtual double deltaMh() const;
    
    virtual double deltaMh2() const;
    
    virtual double deltamt() const;
    
    virtual double deltamt2() const;
    
    virtual double deltamb() const;
    
    virtual double deltamb2() const;
    
    virtual double deltamc() const;
    
    virtual double deltamc2() const;
   
    virtual double deltamtau() const;
    
    virtual double deltamtau2() const;
    
    virtual double deltaGmu() const;
    
    virtual double deltaGmu2() const;
    
    virtual double deltaaMZ() const;
    
    virtual double deltaaMZ2() const;
    
    virtual double deltaa0() const;
    
    virtual double deltaa02() const;
    
    virtual double deltaaSMZ() const;
    
    virtual double deltaaSMZ2() const;
    
    
    virtual double Mw() const;
    
    virtual double deltaMwd6() const;
    
    virtual double deltaMwd62() const;
    
    virtual double deltaGamma_Wff(const Particle fi, const Particle fj) const;
 
    virtual double GammaW(const Particle fi, const Particle fj) const;
    
    virtual double deltaGamma_W() const;
    
    virtual double GammaW() const;
    
    virtual double deltaGwd6() const;
    
    virtual double deltaGwd62() const;
    
    virtual double deltaGzd6() const;
    
    virtual double deltaGzd62() const;
 
    virtual double deltaGV_f(const Particle p) const;
 
    virtual double deltaGA_f(const Particle p) const;
 
        // no generation mixing
    double deltaGL_f(const Particle p) const;
 
        // no generation mixing
    double deltaGR_f(const Particle p) const;
 
 
 
    double CLL_mu() const;
    double CLL_tau() const;
    double CLL_up() const;
    double CLL_down() const;
    double CLL_charm() const;
    double CLL_strange() const;
    double CLL_bottom() const;
    
    double CLR_mu() const;
    double CLR_tau() const;
    double CLR_up() const;
    double CLR_down() const;
    double CLR_charm() const;
    double CLR_strange() const;
    double CLR_bottom() const;
    
    double CRL_mu() const;
    double CRL_tau() const;
    double CRL_up() const;
    double CRL_down() const;
    double CRL_charm() const;
    double CRL_strange() const;
    double CRL_bottom() const;
    
    double CRR_mu() const;
    double CRR_tau() const;
    double CRR_up() const;
    double CRR_down() const;
    double CRR_charm() const;
    double CRR_strange() const;
    double CRR_bottom() const;
 
    
 
    // no generation mixing
    virtual gslpp::complex deltaGL_Wff(const Particle pbar, const Particle p) const;
    // no generation mixing
    virtual gslpp::complex deltaGR_Wff(const Particle pbar, const Particle p) const;
 
    virtual double deltaG_hgg() const;
    virtual double deltaG_hggRatio() const;
    virtual double deltaG1_hWW() const;
    virtual double deltaG2_hWW() const;
    virtual double deltaG3_hWW() const;
    virtual double deltaG1_hZZ() const;
    virtual double deltaG2_hZZ() const;
    virtual double deltaG3_hZZ() const;
    virtual double deltaG1_hZA() const;
    virtual double deltaG1_hZARatio() const;
    virtual double deltaG2_hZA() const;
    virtual double deltaG_hAA() const;
    virtual double deltaG_hAARatio() const;
    // no generation mixing
    virtual gslpp::complex deltaG_hff(const Particle p) const;
    
    virtual double deltaG_hhhRatio() const;
 
    // no generation mixing
    gslpp::complex deltaGL_Wffh(const Particle pbar, const Particle p) const;
    
    // no generation mixing
    gslpp::complex deltaGR_Wffh(const Particle pbar, const Particle p) const;
    
    // no generation mixing
    double deltaGL_Zffh(const Particle p) const;
    
    // no generation mixing 
    double deltaGR_Zffh(const Particle p) const;
      
    // no generation mixing
    gslpp::complex deltaG_hGff(const Particle p) const;
      
    // no generation mixing
    gslpp::complex deltaG_hZff(const Particle p) const;
      
    // no generation mixing
    gslpp::complex deltaG_hAff(const Particle p) const;
      
    // no generation mixing
    gslpp::complex deltaG_Gff(const Particle p) const;
      
    // no generation mixing
    gslpp::complex deltaG_Zff(const Particle p) const;
      
    // no generation mixing
    gslpp::complex deltaG_Aff(const Particle p) const;
      
    double deltag3G() const;
 
 
 
    gslpp::complex f_triangle(const double tau) const;
    gslpp::complex g_triangle(const double tau) const;
    gslpp::complex I_triangle_1(const double tau, const double lambda) const;
    gslpp::complex I_triangle_2(const double tau, const double lambda) const;
    gslpp::complex AH_f(const double tau) const;
    
    gslpp::complex AH_W(const double tau) const;
    
    gslpp::complex AHZga_f(const double tau, const double lambda) const;
    
    gslpp::complex AHZga_W(const double tau, const double lambda) const;
 
    virtual double muggH(const double sqrt_s) const;
    virtual double muggHH(const double sqrt_s) const;
    virtual double muVBF(const double sqrt_s) const;
    virtual double muVBFgamma(const double sqrt_s) const;
    virtual double mueeWBF(const double sqrt_s) const;
    virtual double mueeWBFPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual double mueeHvv(const double sqrt_s) const;
    virtual double mueeHvvPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual double mueeZBF(const double sqrt_s) const;
    virtual double mueeZBFPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual double muepWBF(const double sqrt_s) const;
    virtual double muepZBF(const double sqrt_s) const;
    virtual double muWH(const double sqrt_s) const;
    virtual double muZH(const double sqrt_s) const;
    virtual double mueeZH(const double sqrt_s) const;
    virtual double mueeZllH(const double sqrt_s) const;
    virtual double mueeZqqH(const double sqrt_s) const;
    virtual double mueeZHPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual double mueeZllHPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const;    
    virtual double mueeZqqHPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual double aPskPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual double bPskPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual double muVH(const double sqrt_s) const;
    virtual double muVBFpVH(const double sqrt_s) const;
    virtual double muttH(const double sqrt_s) const;
    virtual double mutHq(const double sqrt_s) const;
    virtual double muggHpttH(const double sqrt_s) const;
    virtual double mueettH(const double sqrt_s) const;
    virtual double mueettHPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    virtual double mummH(const double sqrt_s) const;
    virtual double BrHggRatio() const;
    virtual double BrHWWRatio() const;
    virtual double BrHWlvRatio() const;
    virtual double BrHWW2l2vRatio() const;
    virtual double BrHWjjRatio() const;
    virtual double BrHWW4jRatio() const;
    virtual double BrHWffRatio() const;
    virtual double BrHWW4fRatio() const;
    virtual double BrHZZRatio() const;
    virtual double BrHZllRatio() const;
    virtual double BrHZZ4lRatio() const;  
    virtual double BrHZZ4eRatio() const;
    virtual double BrHZZ2e2muRatio() const;
    virtual double BrHZZ4muRatio() const;
    virtual double BrHZvvRatio() const;
    virtual double BrHZZ4vRatio() const;
    
    virtual double BrHZuuRatio() const;
    virtual double BrHZZ4uRatio() const;
    virtual double BrHZddRatio() const;
    virtual double BrHZZ4dRatio() const;
    virtual double BrHZffRatio() const;
    virtual double BrHZZ4fRatio() const;
    virtual double BrHZgaRatio() const;
    virtual double BrHZgallRatio() const;
    virtual double BrHZgaeeRatio() const;
    virtual double BrHZgamumuRatio() const;
    virtual double BrHgagaRatio() const;
    virtual double BrHmumuRatio() const;
    virtual double BrHtautauRatio() const;
    virtual double BrHccRatio() const;
    virtual double BrHbbRatio() const;
    virtual double computeGammaTotalRatio() const;
    virtual double deltaGammaTotalRatio1() const;
    virtual double deltaGammaTotalRatio1noError() const;
    virtual double deltaGammaTotalRatio2() const;
    double GammaHggRatio() const;
    double GammaHWWRatio() const;
    double GammaHWlvRatio() const;
    double GammaHWW2l2vRatio() const;
    double GammaHWjjRatio() const;
    double GammaHWW4jRatio() const;
    
    double GammaHWffRatio() const;
    double GammaHWW4fRatio() const;
    double GammaHZZRatio() const;
    double GammaHZllRatio() const;
    double GammaHZeeRatio() const;
    double GammaHZmumuRatio() const;
    double GammaHZZ4lRatio() const;  
    double GammaHZZ4eRatio() const;
    double GammaHZZ2e2muRatio() const;
    double GammaHZZ4muRatio() const;
    double GammaHZvvRatio() const;
    double GammaHZZ4vRatio() const;
    double GammaHZuuRatio() const;
    double GammaHZZ4uRatio() const
    {
        return 1.0;
    };
    double GammaHZddRatio() const;
    double GammaHZZ4dRatio() const
    {
        return 1.0;
    };
    double GammaHZffRatio() const;
    double GammaHZZ4fRatio() const;
    double GammaHZgaRatio() const;
    double GammaHgagaRatio() const;
    double GammaHmumuRatio() const;
    double GammaHtautauRatio() const;
    double GammaHccRatio() const;
    double GammaHbbRatio() const;
    double deltaGammaHggRatio1() const;
    double deltaGammaHWWRatio1() const;
    double deltaGammaHWlvRatio1() const;
    double deltaGammaHWW2l2vRatio1() const;    
    double deltaGammaHWjjRatio1() const;
    double deltaGammaHWW4jRatio1() const;
    double deltaGammaHWffRatio1() const;
    double deltaGammaHWW4fRatio1() const;
    double deltaGammaHZZRatio1() const;    
    double deltaGammaHZllRatio1() const;
    double deltaGammaHZeeRatio1() const;
    double deltaGammaHZmumuRatio1() const;
    double deltaGammaHZZ4lRatio1() const;   
    double deltaGammaHZZ4eRatio1() const;
    double deltaGammaHZZ2e2muRatio1() const;
    double deltaGammaHZZ4muRatio1() const;
    double deltaGammaHZvvRatio1() const;
    double deltaGammaHZZ4vRatio1() const;
    double deltaGammaHZuuRatio1() const;
    double deltaGammaHZZ4uRatio1() const
    {
        return 1.0;
    };
    double deltaGammaHZddRatio1() const;
    double deltaGammaHZZ4dRatio1() const
    {
        return 1.0;
    };
    double deltaGammaHZffRatio1() const;
    double deltaGammaHZZ4fRatio1() const;
    double deltaGammaHZgaRatio1() const;
    double deltaGammaHgagaRatio1() const;
    double deltaGammaHmumuRatio1() const;
    double deltaGammaHtautauRatio1() const;
    double deltaGammaHccRatio1() const;
    double deltaGammaHbbRatio1() const;
    double deltaGammaHggRatio2() const;
    double deltaGammaHWWRatio2() const;    
    double deltaGammaHWlvRatio2() const;
    double deltaGammaHWW2l2vRatio2() const;    
    double deltaGammaHWjjRatio2() const;
    double deltaGammaHWW4jRatio2() const;
    double deltaGammaHWffRatio2() const;
    double deltaGammaHWW4fRatio2() const;      
    double deltaGammaHZZRatio2() const;
    double deltaGammaHZllRatio2() const;
    double deltaGammaHZeeRatio2() const;
    double deltaGammaHZmumuRatio2() const;
    double deltaGammaHZZ4lRatio2() const;
    double deltaGammaHZZ4eRatio2() const;
    double deltaGammaHZZ2e2muRatio2() const;
    double deltaGammaHZZ4muRatio2() const;
    double deltaGammaHZvvRatio2() const;
    double deltaGammaHZZ4vRatio2() const;
    double deltaGammaHZuuRatio2() const;
    double deltaGammaHZZ4uRatio2() const;
    double deltaGammaHZddRatio2() const;
    double deltaGammaHZZ4dRatio2() const;
    double deltaGammaHZffRatio2() const;
    double deltaGammaHZZ4fRatio2() const;    
    double deltaGammaHZgaRatio2() const;
    double deltaGammaHgagaRatio2() const;
    double deltaGammaHmumuRatio2() const;
    double deltaGammaHtautauRatio2() const;
    double deltaGammaHccRatio2() const;
    double deltaGammaHbbRatio2() const;
    
    virtual double Br_H_exo() const;
    
    virtual double Br_H_inv() const;
    
    virtual double Br_H_inv_NP() const;
    
    virtual double BrHvisRatio() const;
    
    virtual double BrHtoinvRatio() const;
 
    
    virtual double muttHZbbboost(const double sqrt_s) const;
 
    virtual double muggHgaga(const double sqrt_s) const;    
    virtual double muVBFHgaga(const double sqrt_s) const;
    virtual double muZHgaga(const double sqrt_s) const;
    virtual double muWHgaga(const double sqrt_s) const;
    virtual double muVHgaga(const double sqrt_s) const;
    virtual double muttHgaga(const double sqrt_s) const;  
    virtual double muggHZga(const double sqrt_s) const;
    virtual double muVBFHZga(const double sqrt_s) const;    
    virtual double muZHZga(const double sqrt_s) const;
    virtual double muWHZga(const double sqrt_s) const;    
    virtual double muVHZga(const double sqrt_s) const;
    virtual double muttHZga(const double sqrt_s) const;    
    virtual double muggHZZ(const double sqrt_s) const;
    virtual double muVBFHZZ(const double sqrt_s) const;    
    virtual double muZHZZ(const double sqrt_s) const;
    virtual double muWHZZ(const double sqrt_s) const;    
    virtual double muVHZZ(const double sqrt_s) const;
    virtual double muttHZZ(const double sqrt_s) const;
        
    virtual double muggHZZ4l(const double sqrt_s) const;
    virtual double muVBFHZZ4l(const double sqrt_s) const;    
    virtual double muZHZZ4l(const double sqrt_s) const;
    virtual double muWHZZ4l(const double sqrt_s) const;    
    virtual double muVHZZ4l(const double sqrt_s) const;
    virtual double muttHZZ4l(const double sqrt_s) const;
       
    virtual double muggHWW(const double sqrt_s) const;
    virtual double muVBFHWW(const double sqrt_s) const;   
    virtual double muZHWW(const double sqrt_s) const;
    virtual double muWHWW(const double sqrt_s) const;    
    virtual double muVHWW(const double sqrt_s) const;
    virtual double muttHWW(const double sqrt_s) const;   
    virtual double muggHWW2l2v(const double sqrt_s) const;
    virtual double muVBFHWW2l2v(const double sqrt_s) const;   
    virtual double muZHWW2l2v(const double sqrt_s) const;
    virtual double muWHWW2l2v(const double sqrt_s) const;    
    virtual double muVHWW2l2v(const double sqrt_s) const;
    virtual double muttHWW2l2v(const double sqrt_s) const;  
    virtual double muggHmumu(const double sqrt_s) const;
    virtual double muVBFHmumu(const double sqrt_s) const;    
    virtual double muZHmumu(const double sqrt_s) const;
    virtual double muWHmumu(const double sqrt_s) const;    
    virtual double muVHmumu(const double sqrt_s) const;
    virtual double muttHmumu(const double sqrt_s) const;
    virtual double muggHtautau(const double sqrt_s) const;
    virtual double muVBFHtautau(const double sqrt_s) const;    
    virtual double muZHtautau(const double sqrt_s) const;
    virtual double muWHtautau(const double sqrt_s) const;    
    virtual double muVHtautau(const double sqrt_s) const;
    virtual double muttHtautau(const double sqrt_s) const;
    virtual double muggHbb(const double sqrt_s) const;
    virtual double muVBFHbb(const double sqrt_s) const;    
    virtual double muZHbb(const double sqrt_s) const;
    virtual double muWHbb(const double sqrt_s) const;
    virtual double muVHbb(const double sqrt_s) const;
    virtual double muttHbb(const double sqrt_s) const;
    
//-----------------------------------------------------------------------------------------
//-- Special Hadron collider signal strengths with separate full TH unc U(prod x decay) ---
//-----------------------------------------------------------------------------------------
    
    virtual double muTHUggHgaga(const double sqrt_s) const;    
    virtual double muTHUVBFHgaga(const double sqrt_s) const;
    virtual double muTHUZHgaga(const double sqrt_s) const;
    virtual double muTHUWHgaga(const double sqrt_s) const;
    virtual double muTHUVHgaga(const double sqrt_s) const;
    virtual double muTHUttHgaga(const double sqrt_s) const;
    virtual double muTHUggHZga(const double sqrt_s) const;
    virtual double muTHUVBFHZga(const double sqrt_s) const;    
    virtual double muTHUZHZga(const double sqrt_s) const;
    virtual double muTHUWHZga(const double sqrt_s) const;    
    virtual double muTHUVHZga(const double sqrt_s) const;
    virtual double muTHUttHZga(const double sqrt_s) const; 
    virtual double muTHUggHZZ(const double sqrt_s) const;
    virtual double muTHUVBFHZZ(const double sqrt_s) const;
    virtual double muTHUZHZZ(const double sqrt_s) const;
    virtual double muTHUWHZZ(const double sqrt_s) const;
    virtual double muTHUVHZZ(const double sqrt_s) const;
    virtual double muTHUttHZZ(const double sqrt_s) const;
    
    virtual double muTHUggHZZ4l(const double sqrt_s) const;
    virtual double muTHUVBFHZZ4l(const double sqrt_s) const;    
    virtual double muTHUZHZZ4l(const double sqrt_s) const;
    virtual double muTHUWHZZ4l(const double sqrt_s) const;    
    virtual double muTHUVHZZ4l(const double sqrt_s) const;
    virtual double muTHUttHZZ4l(const double sqrt_s) const;
    
    virtual double muTHUggHWW(const double sqrt_s) const;
    virtual double muTHUVBFHWW(const double sqrt_s) const;
    virtual double muTHUZHWW(const double sqrt_s) const;
    virtual double muTHUWHWW(const double sqrt_s) const; 
    virtual double muTHUVHWW(const double sqrt_s) const;
    virtual double muTHUttHWW(const double sqrt_s) const;
    virtual double muTHUggHWW2l2v(const double sqrt_s) const;
    virtual double muTHUVBFHWW2l2v(const double sqrt_s) const;   
    virtual double muTHUZHWW2l2v(const double sqrt_s) const;
    virtual double muTHUWHWW2l2v(const double sqrt_s) const;    
    virtual double muTHUVHWW2l2v(const double sqrt_s) const;
    virtual double muTHUttHWW2l2v(const double sqrt_s) const;
    virtual double muTHUggHmumu(const double sqrt_s) const;
    virtual double muTHUVBFHmumu(const double sqrt_s) const;    
    virtual double muTHUZHmumu(const double sqrt_s) const;
    virtual double muTHUWHmumu(const double sqrt_s) const;    
    virtual double muTHUVHmumu(const double sqrt_s) const;
    virtual double muTHUttHmumu(const double sqrt_s) const;
    virtual double muTHUggHtautau(const double sqrt_s) const;
    virtual double muTHUVBFHtautau(const double sqrt_s) const;
    virtual double muTHUZHtautau(const double sqrt_s) const;
    virtual double muTHUWHtautau(const double sqrt_s) const;
    virtual double muTHUVHtautau(const double sqrt_s) const;
    virtual double muTHUttHtautau(const double sqrt_s) const;
    virtual double muTHUggHbb(const double sqrt_s) const;
    virtual double muTHUVBFHbb(const double sqrt_s) const;
    virtual double muTHUZHbb(const double sqrt_s) const;
    virtual double muTHUWHbb(const double sqrt_s) const;
    virtual double muTHUVHbb(const double sqrt_s) const;
    virtual double muTHUttHbb(const double sqrt_s) const;
    
    virtual double muTHUVBFBRinv(const double sqrt_s) const;
    virtual double muTHUVBFHinv(const double sqrt_s) const;
    virtual double muTHUVHBRinv(const double sqrt_s) const;
    virtual double muTHUVHinv(const double sqrt_s) const;
    
    virtual double muTHUggHZZ4mu(const double sqrt_s) const;    
    virtual double muTHUggHZgamumu(const double sqrt_s) const;
      
      
    virtual double deltag1ZNP() const;
      
    virtual double deltaKgammaNP() const;
      
    virtual double lambdaZNP() const;
    
      
    virtual double deltag1ZNPEff() const;
      
    virtual double deltaKgammaNPEff() const;
    
 
    virtual double dxseeWWdcos(const double sqrt_s, const double cos) const;
    
    virtual double dxseeWWdcosBin(const double sqrt_s, const double cos1, const double cos2) const;
    
    virtual double xseeWW(const double sqrt_s) const;
    
    virtual double mueeWW(const double sqrt_s) const;
    
    virtual double mueeWWPol(const double sqrt_s, const double Pol_em, const double Pol_ep) const;
    
    
    //----- High Energy diboson observables at hadron colliders
 
    virtual double ppZHprobe(const double sqrt_s) const;
    
    virtual double mupTVppWZ(const double sqrt_s, const double pTV1, const double pTV2) const;
    
    
    
    //----- Simplified Template Cross Sections Bins
    
    virtual double STXS_ggH_VBFtopo_j3v(const double sqrt_s) const;
    
    
    virtual double STXS_ggH_VBFtopo_j3(const double sqrt_s) const;
    
    virtual double STXS_ggH0j(const double sqrt_s) const;
    
    virtual double STXS_ggH1j_pTH_0_60(const double sqrt_s) const;
    
    virtual double STXS_ggH1j_pTH_60_120(const double sqrt_s) const;
    
    virtual double STXS_ggH1j_pTH_120_200(const double sqrt_s) const;
    
    virtual double STXS_ggH1j_pTH_200(const double sqrt_s) const;
    
    virtual double STXS_ggH2j_pTH_0_200(const double sqrt_s) const;
    
    virtual double STXS_ggH2j_pTH_0_60(const double sqrt_s) const;
    
    
    virtual double STXS_ggH2j_pTH_60_120(const double sqrt_s) const;
    
        
    virtual double STXS_ggH2j_pTH_120_200(const double sqrt_s) const;
    
    
    virtual double STXS_ggH2j_pTH_200(const double sqrt_s) const;
    
    virtual double STXS_qqHqq_VBFtopo_Rest(const double sqrt_s) const;
    
    
    virtual double STXS_qqHqq_VBFtopo_j3v(const double sqrt_s) const;
    
    
    virtual double STXS_qqHqq_VBFtopo_j3(const double sqrt_s) const;
    
    virtual double STXS_qqHqq_VHtopo(const double sqrt_s) const;
    
    virtual double STXS_qqHqq_Rest(const double sqrt_s) const;
    
    virtual double STXS_qqHqq_pTj_200(const double sqrt_s) const;
    
    virtual double STXS_qqHlv_pTV_0_250(const double sqrt_s) const;
    
 
    virtual double STXS_qqHlv_pTV_0_150(const double sqrt_s) const;
    
    
    virtual double STXS_qqHlv_pTV_150_250_0j(const double sqrt_s) const;
    
    
    virtual double STXS_qqHlv_pTV_150_250_1j(const double sqrt_s) const;
    
    
    virtual double STXS_qqHlv_pTV_250(const double sqrt_s) const;
    
    virtual double STXS_qqHll_pTV_0_150(const double sqrt_s) const;
    
    virtual double STXS_qqHll_pTV_150_250(const double sqrt_s) const;
    
    
    virtual double STXS_qqHll_pTV_150_250_0j(const double sqrt_s) const;
    
    
    virtual double STXS_qqHll_pTV_150_250_1j(const double sqrt_s) const;
    
    virtual double STXS_qqHll_pTV_250(const double sqrt_s) const;
    
    
    virtual double STXS_ttHtH(const double sqrt_s) const;
    
    virtual double STXS_WHqqHqq_VBFtopo_j3v(const double sqrt_s) const;
    
    virtual double STXS_WHqqHqq_VBFtopo_j3(const double sqrt_s) const;
    
    virtual double STXS_WHqqHqq_VH2j(const double sqrt_s) const;
    
    virtual double STXS_WHqqHqq_Rest(const double sqrt_s) const;
    
    virtual double STXS_WHqqHqq_pTj1_200(const double sqrt_s) const;
    
    virtual double STXS_ZHqqHqq_VBFtopo_j3v(const double sqrt_s) const;
    
    virtual double STXS_ZHqqHqq_VBFtopo_j3(const double sqrt_s) const;
    
    virtual double STXS_ZHqqHqq_VH2j(const double sqrt_s) const;
    
    virtual double STXS_ZHqqHqq_Rest(const double sqrt_s) const;
    
    virtual double STXS_ZHqqHqq_pTj1_200(const double sqrt_s) const;
    
    
    
    virtual double kappamueff() const;
    
    virtual double kappataueff() const;
    
    virtual double kappaceff() const;
    
    virtual double kappabeff() const;
    
    virtual double kappaGeff() const;
    
    virtual double kappaZeff() const;
    
    virtual double kappaWeff() const;
    
    virtual double kappaAeff() const;
    
    virtual double kappaZAeff() const;
    
    
    virtual double deltayt_HB() const;
    
    virtual double deltayb_HB() const;
    
    virtual double deltaytau_HB() const;
    
    virtual double deltayc_HB() const;
    
    virtual double deltaymu_HB() const;
    
    virtual double deltacZ_HB() const;
    
    virtual double cZBox_HB() const;
    
    virtual double cZZ_HB() const;
    
    virtual double cZga_HB() const;
    
    virtual double cgaga_HB() const;
    
    virtual double cgg_HB() const;
    
    virtual double cggEff_HB() const;
    
    virtual double lambz_HB() const;
    
    
    virtual double AuxObs_NP1() const;
    
    virtual double AuxObs_NP2() const;
    
    virtual double AuxObs_NP3() const;
    
    virtual double AuxObs_NP4() const;
    
    virtual double AuxObs_NP5() const;
    
    virtual double AuxObs_NP6() const;
    
    virtual double AuxObs_NP7() const;
    
    virtual double AuxObs_NP8() const;
    
    virtual double AuxObs_NP9() const;
    
    virtual double AuxObs_NP10() const;
    
    virtual double AuxObs_NP11() const;
        
    virtual double AuxObs_NP12() const;
    
    virtual double AuxObs_NP13() const;
    
    virtual double AuxObs_NP14() const;
    
    virtual double AuxObs_NP15() const;
    
    virtual double AuxObs_NP16() const;
    
    virtual double AuxObs_NP17() const;
    
    virtual double AuxObs_NP18() const;
    
    virtual double AuxObs_NP19() const;
    
    virtual double AuxObs_NP20() const;
    
protected:
 
    virtual void setParameter(const std::string name, const double& value);
    mutable Matching<NPSMEFTd6Matching,NPSMEFTd6> NPSMEFTd6M;
 
    double CG; 
    double CW; 
    double C2B; 
    double C2W; 
    double C2BS; 
    double C2WS; 
    double CHG; 
    double CHW; 
    double CHB; 
    double CHWHB_gaga; 
    double CHWHB_gagaorth; 
    double CDHB; 
    double CDHW; 
    double CDB; 
    double CDW; 
    double CHWB; 
    double CHD; 
    double CT; 
    double CHbox; 
    double CH; 
    double CHL1_11; 
    double CHL1_12r; 
    double CHL1_13r; 
    double CHL1_22; 
    double CHL1_23r; 
    double CHL1_33; 
    double CHL1_12i; 
    double CHL1_13i; 
    double CHL1_23i; 
    double CHL3_11; 
    double CHL3_12r; 
    double CHL3_13r; 
    double CHL3_22; 
    double CHL3_23r; 
    double CHL3_33; 
    double CHL3_12i; 
    double CHL3_13i; 
    double CHL3_23i; 
    double CHe_11; 
    double CHe_12r; 
    double CHe_13r; 
    double CHe_22; 
    double CHe_23r; 
    double CHe_33; 
    double CHe_12i; 
    double CHe_13i; 
    double CHe_23i; 
    double CHQ1_11; 
    double CHQ1_12r; 
    double CHQ1_13r; 
    double CHQ1_22; 
    double CHQ1_23r; 
    double CHQ1_33; 
    double CHQ1_12i; 
    double CHQ1_13i; 
    double CHQ1_23i; 
    double CHQ3_11; 
    double CHQ3_12r; 
    double CHQ3_13r; 
    double CHQ3_22; 
    double CHQ3_23r; 
    double CHQ3_33; 
    double CHQ3_12i; 
    double CHQ3_13i; 
    double CHQ3_23i; 
    double CHu_11; 
    double CHu_12r; 
    double CHu_13r; 
    double CHu_22; 
    double CHu_23r; 
    double CHu_33; 
    double CHu_12i; 
    double CHu_13i; 
    double CHu_23i; 
    double CHd_11; 
    double CHd_12r; 
    double CHd_13r; 
    double CHd_22; 
    double CHd_23r; 
    double CHd_33; 
    double CHd_12i; 
    double CHd_13i; 
    double CHd_23i; 
    double CHud_11r; 
    double CHud_12r; 
    double CHud_13r; 
    double CHud_22r; 
    double CHud_23r; 
    double CHud_33r; 
    double CHud_11i; 
    double CHud_12i; 
    double CHud_13i; 
    double CHud_22i; 
    double CHud_23i; 
    double CHud_33i; 
    double CeH_11r; 
    double CeH_12r; 
    double CeH_13r; 
    double CeH_22r; 
    double CeH_23r; 
    double CeH_33r; 
    double CeH_11i; 
    double CeH_12i; 
    double CeH_13i; 
    double CeH_22i; 
    double CeH_23i; 
    double CeH_33i; 
    double CuH_11r; 
    double CuH_12r; 
    double CuH_13r; 
    double CuH_22r; 
    double CuH_23r; 
    double CuH_33r; 
    double CuH_11i; 
    double CuH_12i; 
    double CuH_13i; 
    double CuH_22i; 
    double CuH_23i; 
    double CuH_33i; 
    double CdH_11r; 
    double CdH_12r; 
    double CdH_13r; 
    double CdH_22r; 
    double CdH_23r; 
    double CdH_33r; 
    double CdH_11i; 
    double CdH_12i; 
    double CdH_13i; 
    double CdH_22i; 
    double CdH_23i; 
    double CdH_33i; 
    double CuG_11r; 
    double CuG_12r; 
    double CuG_13r; 
    double CuG_22r; 
    double CuG_23r; 
    double CuG_33r; 
    double CuG_11i; 
    double CuG_12i; 
    double CuG_13i; 
    double CuG_22i; 
    double CuG_23i; 
    double CuG_33i; 
    double CuW_11r; 
    double CuW_12r; 
    double CuW_13r; 
    double CuW_22r; 
    double CuW_23r; 
    double CuW_33r; 
    double CuW_11i; 
    double CuW_12i; 
    double CuW_13i; 
    double CuW_22i; 
    double CuW_23i; 
    double CuW_33i; 
    double CuB_11r; 
    double CuB_12r; 
    double CuB_13r; 
    double CuB_22r; 
    double CuB_23r; 
    double CuB_33r; 
    double CuB_11i; 
    double CuB_12i; 
    double CuB_13i; 
    double CuB_22i; 
    double CuB_23i; 
    double CuB_33i; 
    double CLL_1111;
    double CLL_1221, CLL_2112, CLL_1122, CLL_2211;
    double CLL_1331, CLL_3113, CLL_1133, CLL_3311;
    double CLQ1_1111;
    double CLQ1_1122, CLQ1_2211, CLQ1_1221, CLQ1_2112;
    double CLQ1_1133, CLQ1_3311, CLQ1_1331, CLQ1_3113;
    double CLQ1_1123, CLQ1_2223, CLQ1_3323;
    double CLQ1_1132, CLQ1_2232, CLQ1_3332;
    double CLQ3_1111;
    double CLQ3_1122, CLQ3_2211, CLQ3_1221, CLQ3_2112;
    double CLQ3_1133, CLQ3_3311, CLQ3_1331, CLQ3_3113;
    double CLQ3_1123, CLQ3_2223, CLQ3_3323;
    double CLQ3_1132, CLQ3_2232, CLQ3_3332;
    double Cee_1111;
    double Cee_1122, Cee_2211;
    double Cee_1133, Cee_3311;
    double Ceu_1111;
    double Ceu_1122, Ceu_2211;
    double Ceu_1133, Ceu_3311;
    double Ceu_2233;
    double Ced_1111;
    double Ced_1122, Ced_2211;
    double Ced_1133, Ced_3311;
    double Ced_1123, Ced_2223, Ced_3323;
    double Ced_1132, Ced_2232, Ced_3332;
    double CLe_1111;
    double CLe_1122, CLe_2211;
    double CLe_1133, CLe_3311;
    double CLu_1111;
    double CLu_1122, CLu_2211;
    double CLu_1133, CLu_3311;
    double CLu_2233;
    double CLd_1111;
    double CLd_1122, CLd_2211;
    double CLd_1133, CLd_3311;
    double CLd_1123, CLd_2223, CLd_3323;
    double CLd_1132, CLd_2232, CLd_3332;
    double CQe_1111;
    double CQe_1122, CQe_2211;
    double CQe_1133, CQe_3311;
    double CQe_2311, CQe_2322, CQe_2333;
    double CQe_3211, CQe_3222, CQe_3233;
    double CLedQ_11, CLedQ_22, CpLedQ_11, CpLedQ_22;
    double Lambda_NP; 
// 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 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; 
    
// The error in the parameters multiplying the dimension-6 operator coefficients in the production cross sections.
    double eVBF_2_Hbox;
    double eVBF_2_HQ1_11;
    double eVBF_2_Hu_11;
    double eVBF_2_Hd_11;
    double eVBF_2_HQ3_11;
    double eVBF_2_HD;
    double eVBF_2_HB;
    double eVBF_2_HW;
    double eVBF_2_HWB;
    double eVBF_2_HG;
    double eVBF_2_DHB;
    double eVBF_2_DHW;
    double eVBF_2_DeltaGF;
      
    double eVBF_78_Hbox;
    double eVBF_78_HQ1_11;
    double eVBF_78_Hu_11;
    double eVBF_78_Hd_11;
    double eVBF_78_HQ3_11;
    double eVBF_78_HD;
    double eVBF_78_HB;
    double eVBF_78_HW;
    double eVBF_78_HWB;
    double eVBF_78_HG;
    double eVBF_78_DHB;
    double eVBF_78_DHW;
    double eVBF_78_DeltaGF;
      
    double eVBF_1314_Hbox;
    double eVBF_1314_HQ1_11;
    double eVBF_1314_Hu_11;
    double eVBF_1314_Hd_11;
    double eVBF_1314_HQ3_11;
    double eVBF_1314_HD;
    double eVBF_1314_HB;
    double eVBF_1314_HW;
    double eVBF_1314_HWB;
    double eVBF_1314_HG;
    double eVBF_1314_DHB;
    double eVBF_1314_DHW;
    double eVBF_1314_DeltaGF;
 
    double eWH_2_Hbox;
    double eWH_2_HQ3_11;
    double eWH_2_HD;
    double eWH_2_HW;
    double eWH_2_HWB;
    double eWH_2_DHW;
    double eWH_2_DeltaGF;
      
    double eWH_78_Hbox;
    double eWH_78_HQ3_11;
    double eWH_78_HD;
    double eWH_78_HW;
    double eWH_78_HWB;
    double eWH_78_DHW;
    double eWH_78_DeltaGF;
      
    double eWH_1314_Hbox;
    double eWH_1314_HQ3_11;
    double eWH_1314_HD;
    double eWH_1314_HW;
    double eWH_1314_HWB;
    double eWH_1314_DHW;
    double eWH_1314_DeltaGF;
      
    double eZH_2_Hbox;
    double eZH_2_HQ1_11;
    double eZH_2_Hu_11;
    double eZH_2_Hd_11;
    double eZH_2_HQ3_11;
    double eZH_2_HD;
    double eZH_2_HB;
    double eZH_2_HW;
    double eZH_2_HWB;
    double eZH_2_DHB;
    double eZH_2_DHW;
    double eZH_2_DeltaGF;
      
    double eZH_78_Hbox;
    double eZH_78_HQ1_11;
    double eZH_78_Hu_11;
    double eZH_78_Hd_11;
    double eZH_78_HQ3_11;
    double eZH_78_HD;
    double eZH_78_HB;
    double eZH_78_HW;
    double eZH_78_HWB;
    double eZH_78_DHB;
    double eZH_78_DHW;
    double eZH_78_DeltaGF;
      
    double eZH_1314_Hbox;
    double eZH_1314_HQ1_11;
    double eZH_1314_Hu_11;
    double eZH_1314_Hd_11;
    double eZH_1314_HQ3_11;
    double eZH_1314_HD;
    double eZH_1314_HB;
    double eZH_1314_HW;
    double eZH_1314_HWB;
    double eZH_1314_DHB;
    double eZH_1314_DHW;
    double eZH_1314_DeltaGF;
 
    double ettH_2_HG;
    double ettH_2_G;
    double ettH_2_uG_33r;
    double ettH_2_DeltagHt;
      
    double ettH_78_HG;
    double ettH_78_G;
    double ettH_78_uG_33r;
    double ettH_78_DeltagHt;
      
    double ettH_1314_HG;
    double ettH_1314_G;
    double ettH_1314_uG_33r;
    double ettH_1314_DeltagHt;
        
    double BrHinv; 
    double BrHexo; 
    
    double dg1Z; 
    double dKappaga; 
    double lambZ; 
 
//  ----------------------------------------------------------------------------
//  Internal parameters    
//  ----------------------------------------------------------------------------
    
    double LambdaNP2;
    
    // Internal values for some of the dimension-6 coefficients (to allow changes of bases in Post-Update)
    double CiHL1_11;
    double CiHL1_22;
    double CiHL1_33;
    double CiHL3_11;
    double CiHL3_22;
    double CiHL3_33;
    
    double CiHQ1_11;
    double CiHQ1_22;
    double CiHQ1_33;
    double CiHQ3_11;
    double CiHQ3_22;
    double CiHQ3_33;
    
    double CiHe_11;
    double CiHe_22;
    double CiHe_33;
    
    double CiHu_11;
    double CiHu_22;
    double CiHu_33;
    
    double CiHd_11;
    double CiHd_22;
    double CiHd_33;
    
    double CiW;
    
    double CiHW;
    double CiHB; 
    double CiDHB; 
    double CiDHW;
    double CiHWB;
    
    double CiHbox;
    double CiHD;
    double CiH;
    
    double CieH_11r;
    double CieH_22r;
    double CieH_33r;
    
    double CiuH_11r;
    double CiuH_22r;
    double CiuH_33r;
    
    double CidH_11r;
    double CidH_22r;
    double CidH_33r;
    
    double CiuG_11r;
    double CiuG_22r;
    double CiuG_33r;
    
    double CiuW_11r;
    double CiuW_22r;
    double CiuW_33r;
    
    double CiuB_11r;
    double CiuB_22r;
    double CiuB_33r;
    
    double CiLL_1221;
    double CiLL_2112;    
    
    double v2;
    double v2_over_LambdaNP2;
    double aleMz;
    double eeMz;
    double eeMz2;
    double cW_tree;
    double sW_tree;
    double cW2_tree;
    double sW2_tree;
    
    double g1_tree;
    double g2_tree;
    double g3_tree;
    
    double gZvL;
    double gZlL,gZlR;
    double gZuL,gZuR;
    double gZdL,gZdR;
    
    double UevL;
    double VudL;
        
    double delta_ZZ;
    double delta_AA;
    double delta_AZ;
    double delta_h;
    
    double lambdaH_tree;
    
    double dZH;
    
    double cHSM;
    
    double cLHd6;
    
    double cLH3d62;
 
    double Yuke,Yukmu,Yuktau;
    double Yuku,Yukc,Yukt;
    double Yukd,Yuks,Yukb;
    
    // STXS ci: denote them as ai to separate from the normal conventions 
    double aiG,ai3G,ai2G;
    double aiT,aiH,aiWW,aiB,aiHW,aiHB;
    double aiA;
    double aiHQ,aipHQ,aiHL,aipHL,aiHu,aiHd,aiHe;
    double aiu,aiuG;
    
    // Cache variables
    double dGammaHTotR1,dGammaHTotR2,GammaHTotR;
    
    double CHF1_diag(const Particle F) const;
    
    double CHF3_diag(const Particle F) const;
 
    double CHf_diag(const Particle f) const;
    
    gslpp::complex CHud_diag(const Particle u) const;
    
    gslpp::complex CfH_diag(const Particle f) const;
      
    gslpp::complex CfG_diag(const Particle f) const;
      
    gslpp::complex CfW_diag(const Particle f) const;
      
    gslpp::complex CfB_diag(const Particle f) const;
    
    
private:
 
    bool FlagQuadraticTerms; 
    bool FlagRotateCHWCHB; 
    bool FlagPartialQFU; 
    bool FlagFlavU3OfX; 
    bool FlagUnivOfX; 
    bool FlagHiggsSM; 
    bool FlagLoopHd6; 
    bool FlagLoopH3d6Quad; 
 
    const bool FlagLeptonUniversal;
 
    const bool FlagQuarkUniversal;
    
    gsl_integration_cquad_workspace * w_WW;
};
 
#endif  /* NPSMEFTD6_H */