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>& CMbtaunu(QCD::meson meson_i) ;
    
    virtual  std::vector<WilsonCoefficient>& CMBXsnn() ;
    
    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() ;
 
    double x_c(const double mu, const orders order = FULLNNLO) const;
    double x_t(const double mu, const orders order = FULLNNLO) const;
    double mt2omh2(const double mu, const orders order = FULLNNLO) const;
    
    double A0t(double x)const;
    
    double B0t(double x)const;
    
    double C0t(double x)const;
    
    double D0t(double x)const;
    
    double F0t(double x)const;
    
    double E0t(double x)const;
    
    double A1t(double x, double mu) const;
    
    double B1t(double x, double mu)const;
    
    double C1t(double x, double mu)const;
    
    double D1t(double x, double mu)const;
    
    double F1t(double x, double mu) const;
    
    double E1t(double x, double mu) const;
    
    double G1t(double x, double mu) const;
    
    double C7c_3L_at_mW(double x) const;
    
    double C7t_3L_at_mt(double x) const;
    
    double C7t_3L_func(double x, double mu) const;
    
    double C8c_3L_at_mW(double x) const;
    
    double C8t_3L_at_mt(double x) const;
    
    double C8t_3L_func(double x, double mu) const;
 
    double Tt(double x) const;
    
    double Wt(double x) const;
    
    double Eet(double x) const;
    
    double Rest(double x, double mu) const;
 
    double C10_OS1(double x, double mu);
    
    double Y0(double x)const;
    
    double Y1(double x, double mu)const;
    
    double C7LOeff(double x)const;
    
    double C7NLOeff(double x)const;
    
    double C8LOeff(double x)const;
    
    double C8NLOeff(double x)const;
 
    double B0b(double x)const;
    
    double C0b(double x)const;
    
    double D0b(double x)const;
    
    double D0b_tilde(double x)const;
    
    double E0b(double x)const;
    
    double B1d(double x, double mu)const;
    
    double B1d_tilde(double x, double mu)const;
    
    double B1u(double x, double mu)const;
    
    double B1u_tilde(double x, double mu)const;
    
    double C1ew(double x)const;
    
    double Zew(double xt, double xz)const;
        
    double Gew(double xt, double xz, double mu)const;
    
    double Hew(double xt, double xz, double mu)const;
    
    gslpp::complex S0c() const;
    
    gslpp::complex S0ct() const;
    
    gslpp::complex S0tt() const;
    
    double X0t(double x)const;
    
    double X1t(double x)const;
    
    double Xewt(double x, double a, double mu)const;
 
    double S1(double x) const;
 
    virtual double S0(double, double) 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, vmcbtaunu;
    std::vector<WilsonCoefficient> vmcDLij, vmcDLi3j, vmcmueconv, vmcgminus2mu;
 
    
private:    
    const StandardModel & SM; 
    double S0(double) const;
    double S0p(double x) const;
    double S11(double x) 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, mcbtaunu;
    WilsonCoefficient mcDLij, mcDLi3j, mcmueconv, mcgminus2mu;
    
    double Mut, Muw, Ale, GF, Mw_tree, Nc, CF, Mt_muw, Mt_mut;
    double gamma0, J5, BtNDR, Mw, sW2, mu_b;
    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);
 
    double phi1(double z) 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 sw, swa, swb, swc, swd, swe, swf; //sen(theta_W) tree level
    double xcachea, xcacheb, xcachec, xcached, xcachee, xcachef; // caching
 
    
    double phi_z(double z) 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);
 
    double fbb(double x);
    
    double gbb(double x);
 
    double taub2(double x);
 
    double Delta_t(double mu, double x);
 
    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 */