StandardModelMatching.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2012 HEPfit Collaboration
 * All rights reserved.
 *
 * 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);
    
    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();
    virtual  std::vector<WilsonCoefficient>& CMprimeBMll();
    
    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();
    
    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 Tt(double x) const;
    
    double Wt(double x) const;
    
    double Eet(double x) const;
    
    double Rest(double x, double mu) const;
    
    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 E0b(double x)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;

    double S0(double, double) const;
    
    
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;
    double Muw;
    double Ale;
    double GF;
    double Mw_tree;
    double Nc;
    double CF;
    double gamma0;
    double J5;
    double BtNDR;
    double Mw;
    double sW2;
    double mu_b;
    //double MM;
    gslpp::matrix<gslpp::complex> Vckm;
    gslpp::complex lam_t;
    double L;
    
    double setWCBsmm(int i, double x, orders order);
    
    double setWCBsmmEW(int i, double x, orders_ew order_ew);
    
    double setWCBdmm(int i, double x, orders order);
    
    double setWCBdmmEW(int i, double x, orders_ew order_ew);
    
    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];
    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
    
    
};

#endif  /* STANDARDMODELMATCHING_H */