StandardModelMatching.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef STANDARDMODELMATCHING_H
#define STANDARDMODELMATCHING_H
 
#include "ModelMatching.h"
#include "gslpp.h"
 
#define LEPS 1.e-5 // tolerance in the limit of S(x,y) to S(x)
 
class StandardModel;
 
class StandardModelMatching : public ModelMatching {
public:
 
    StandardModelMatching(const StandardModel & SM_i);
    
    virtual ~StandardModelMatching();
    
    void updateSMParameters();
    
    virtual std::vector<WilsonCoefficient>& CMdbd2();
    
    virtual std::vector<WilsonCoefficient>& CMdbs2() ;
 
    virtual std::vector<WilsonCoefficient>& CMdd2() ;
    
    virtual std::vector<WilsonCoefficient>& CMdk2() ;
        
    virtual std::vector<WilsonCoefficient>& CMbnlep( int a) ;
    
    virtual std::vector<WilsonCoefficient>& CMbnlepCC( int a) ;
    
    virtual std::vector<WilsonCoefficient>& CMbsg() ;
    
    virtual std::vector<WilsonCoefficient>& CMprimebsg() ;
    
    virtual std::vector<WilsonCoefficient>& CMBMll(QCD::lepton lepton) ;
    virtual std::vector<WilsonCoefficient>& CMprimeBMll(QCD::lepton lepton) ;
    
    virtual std::vector<WilsonCoefficient>& CMK() ;
    
    virtual std::vector<WilsonCoefficient>& CMKCC() ;
    
    virtual  std::vector<WilsonCoefficient>& CMd1() ;
    
    virtual  std::vector<WilsonCoefficient>& CMd1Buras() ;
    
    virtual std::vector<WilsonCoefficient>& CMkpnn() ;
    virtual std::vector<WilsonCoefficient>& CMkmm() ;
    
    virtual std::vector<WilsonCoefficient>& CMbsmm() ;
    
    virtual  std::vector<WilsonCoefficient>& CMbdmm() ;
        
    virtual std::vector<WilsonCoefficient>& CMdiujleptonknu(int i, int j, int k) ;
    
    
    virtual std::vector<WilsonCoefficient>& CMBXsnn(QCD::lepton lepton) ;
    
    virtual  std::vector<WilsonCoefficient>& CMBXdnn() ;
    
    virtual  std::vector<WilsonCoefficient>& CMDLij(int li_lj) ;
 
    virtual std::vector<WilsonCoefficient>& CMDLi3j(int li_lj) ;
 
    virtual std::vector<WilsonCoefficient>& CMmueconv() ;
 
    virtual std::vector<WilsonCoefficient>& CMgminus2mu() ;
 
    const double x_c(const double mu, const orders order = FULLNNLO) const;
    const double x_t(const double mu, const orders order = FULLNNLO) const;
    const double mt2omh2(const double mu, const orders order = FULLNNLO) const;
    
    const double A0t(double x) const;
    
    const double B0t(double x) const;
    
    const double C0t(double x) const;
    
    const double D0t(double x) const;
    
    const double F0t(double x) const;
    
    const double E0t(double x) const;
    
    const double A1t(double x, double mu) const;
    
    const double B1t(double x, double mu) const;
    
    const double C1t(double x, double mu) const;
    
    const double D1t(double x, double mu) const;
    
    const double F1t(double x, double mu) const;
    
    const double E1t(double x, double mu) const;
    
    const double G1t(double x, double mu) const;
    
    const double C7c_3L_at_mW(double x) const;
    
    const double C7t_3L_at_mt(double x) const;
    
    const double C7t_3L_func(double x, double mu) const;
    
    const double C8c_3L_at_mW(double x) const;
    
    const double C8t_3L_at_mt(double x) const;
    
    const double C8t_3L_func(double x, double mu) const;
 
    const double Tt(double x) const;
    
    const double Wt(double x) const;
    
    const double Eet(double x) const;
    
    const double Rest(double x, double mu) const;
 
    double C10_OS1(double x, double mu);
    
    const double Y0(double x) const;
    
    const double Y1(double x, double mu) const;
    
    const double C7LOeff(double x) const;
    
    const double C7NLOeff(double x) const;
    
    const double C8LOeff(double x) const;
    
    const double C8NLOeff(double x) const;
 
    const double B0b(double x) const;
    
    const double C0b(double x) const;
    
    const double D0b(double x) const;
    
    const double D0b_tilde(double x) const;
    
    const double E0b(double x) const;
    
    const double B1d(double x, double mu) const;
    
    const double B1d_tilde(double x, double mu) const;
    
    const double B1u(double x, double mu) const;
    
    const double B1u_tilde(double x, double mu) const;
    
    const double C1ew(double x) const;
    
    const double Zew(double xt, double xz) const;
        
    const double Gew(double xt, double xz, double mu) const;
    
    const double Hew(double xt, double xz, double mu) const;
    
    gslpp::complex S0c() const;
    
    gslpp::complex S0ct() const;
    
    gslpp::complex S0tt() const;
    
    const double X0t(double x) const;
    
    const double X1t(double x) const;
    
    const double Xewt(double x, double a, double mu) const;
 
    const double S1(double x) const;
 
    virtual const double S0(double, double) const;
 
    /*
     * Wilson coefficients Misiak basis
     * Operator block: L (2)
     * Normalization: 4 G_F / sqrt(2) x CKM
     */
    virtual std::vector<WilsonCoefficientNew>& CMDF1(std::string blocks, unsigned int nops);
 
    const double getMt_mut() const {
        return Mt_mut;
    }
     
protected:
    std::vector<WilsonCoefficient> vmcdb, vmcds, vmcd2, vmck2, vmck, vmckcc;
    std::vector<WilsonCoefficient> vmcbsg, vmcprimebsg, vmcBMll, vmcprimeBMll, vmcbnlep, vmcbnlepCC, vmcd1, vmcd1Buras;
    std::vector<WilsonCoefficient> vmckpnn, vmckmm, vmcbsnn, vmcbdnn, vmcbsmm, vmcbdmm, vmculeptonnu;
    std::vector<WilsonCoefficient> vmcDLij, vmcDLi3j, vmcmueconv, vmcgminus2mu;
    std::vector<WilsonCoefficientNew> vmcDF1;
 
    
private:    
    const StandardModel & SM; 
    const double S0(double) const;
    const double S0p(double x) const;
    const double S11(double x) const;
    const double S18(double x) const;
    WilsonCoefficient mcdbd2, mcdbs2, mcdd2, mcdk2, mck, mckcc;
    WilsonCoefficient mcbsg, mcprimebsg, mcBMll, mcprimeBMll, mcbnlep, mcbnlepCC, mcd1, mcd1Buras;
    WilsonCoefficient mckpnn, mckmm, mcbsnn, mcbdnn, mcbsmm, mcbdmm, mculeptonnu;
    WilsonCoefficient mcDLij, mcDLi3j, mcmueconv, mcgminus2mu;
    WilsonCoefficientNew mcC, mcP, mcM, mcL, mcQ, mcB;
    
    double Mut, Muw, Ale, GF, Mw_tree, Nc, CF, Mt_muw, Mt_mut;
    double gamma0, J5, BtNDR, Mw, sW2, mu_b, sw;
    double L, Lz;
    double alstilde, aletilde;
    
    gslpp::complex lam_t;
    gslpp::matrix<gslpp::complex> Vckm;
 
    double setWCBsmm(int i, double x, orders order);
    
    double setWCBsmmEW(int i, double x, orders_qed order_qed);
    
    double setWCBdmm(int i, double x, orders order);
    
    double setWCBdmmEW(int i, double x, orders_qed order_qed);
    
    double setWCbsg (int i, double x, orders order);
    
    double setWCBMll (int i, double x, orders order);
    
    double setWCbnlep (int i, double x, orders order);
    
    double setWCbnlepEW (int i, double x);
 
    const double phi1(double z) const;
    
    const double phi2 (double x, double y) const;
    
    double CWbsgArrayLO[8], CWbsgArrayNLO[8], CWbsgArrayNNLO[8];
    double CWprimebsgArrayLO[8], CWprimebsgArrayNLO[8];
    double CWBMllArrayLO[19], CWBMllArrayNLO[19];
    double CWD1ArrayLO[10], CWD1ArrayNLO[10];
    double CWbnlepArrayLOqcd[10], CWbnlepArrayNLOqcd[10];
    double CWbnlepArrayLOew[10], CWbnlepArrayNLOew[10];
    
    double CWBsmmArrayNNLOqcd[8], CWBsmmArrayNLOqcd[8], CWBsmmArrayLOqcd[8];
    double CWBsmmArrayNLOewt4[8], CWBsmmArrayNLOewt2[8], CWBsmmArrayNLOew[8];
    
    double CWBdmmArrayNNLOqcd[8], CWBdmmArrayNLOqcd[8], CWBdmmArrayLOqcd[8];
    double CWBdmmArrayNLOewt4[8], CWBdmmArrayNLOewt2[8], CWBdmmArrayNLOew[8];
    
    double sw2, swa, sw2a, swb, sw2b, swc, sw2c, swd, sw2d, swe, sw2e, swf, sw2f; //sen(theta_W) tree level
    double xcachea, xcacheb, xcachec, xcached, xcachee, xcachef; // caching
 
    
    const double phi_z(double z) const;
    
    const double phi_xy (double x, double y) const;
    
    double C3funNNLO(double x);
    
    double C4fun(double x, orders ord);
 
    double C5funNNLO(double x);
 
    double C6funNNLO(double x);
    
    double C7funLO(double x);
    
    double C8funLO(double x);
 
    unsigned int setCMDF1(WilsonCoefficientNew& CMDF1, WilsonCoefficientNew& DF1block, unsigned int tot, schemes scheme, qcd_orders order_qcd, qed_orders order_qed);
 
    double fbb(double x);
    
    double gbb(double x);
 
    double taub2(double x);
 
    double Delta_t(double mu, double x);
 
    /*
     * Wilson coefficients Misiak basis
     * Operator block: C (2)
     * Normalization: 4 G_F / sqrt(2) x CKM
     */
    WilsonCoefficientNew& mc_C();
 
    /*
     * Wilson coefficients Misiak basis
     * Operator block: P (4)
     * Normalization: 4 G_F / sqrt(2) x CKM
     */
    WilsonCoefficientNew& mc_P();
 
    /*
     * Wilson coefficients Misiak basis
     * Operator block: M (2)
     * Normalization: 4 G_F / sqrt(2) x CKM
     * QED only available at NLO and in approximate formulas
     * QED ref.: Gambino, Haisch, JHEP 0110, 020, hep-ph/0109058
     */
    WilsonCoefficientNew& mc_M();
 
    /*
     * Wilson coefficients Misiak basis
     * Operator block: L (2)
     * Normalization: 4 G_F / sqrt(2) x CKM
     */
    WilsonCoefficientNew& mc_L();
 
    /*
     * Wilson coefficients Misiak basis
     * Operator block: Q (4)
     * Normalization: 4 G_F / sqrt(2) x CKM
     * QED_NLO ref.: Gambino, Haisch, JHEP 0110, 020, hep-ph/0109058 - COULD BE CHANGED TO X,Y,W
     */
    WilsonCoefficientNew& mc_Q();
 
    /*
     * Wilson coefficients Misiak basis
     * Operator block: B (1)
     * Normalization: 4 G_F / sqrt(2) x CKM
     */
    WilsonCoefficientNew& mc_B();
 
    friend double gslpp_special_functions::dilog(double x);
    friend double gslpp_special_functions::clausen(double x);
    friend double gslpp_special_functions::zeta(int i);
    
};
 
#endif  /* STANDARDMODELMATCHING_H */