GeneralTHDMMatching.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2017 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef GENERALTHDMMATCHING_H
#define GENERALTHDMMATCHING_H
 
 #include <Polylogarithms.h>
#include "gslpp.h"
#include "StandardModelMatching.h"
 
class GeneralTHDM;
 
class GeneralTHDMMatching : public StandardModelMatching {
public:
    GeneralTHDMMatching(const GeneralTHDM & GeneralTHDM_i);
 
    
    
    
    
    
    virtual std::vector<WilsonCoefficient>& CMgminus2mu();
 
    virtual double gminus2muLO();
    //FOR THE MOMENT WE LEAVE THIS DEFINITION OF THE LO CONTRIBUTION TO THE g-2, MUST BE CHECKED
    
    
    
    double F1oneloopgm2(const double ratio_sq);
    
    
    
    double F2oneloopgm2(const double ratio_sq);
    
    
    
    double F3oneloopgm2(const double ratio_sq);
    
    
    
    double F1twoloopgm2(const double ratio_sq);
    
    
    
    double F1tildetwoloopgm2(const double ratio_sq);
    
    
    double F2twoloopgm2(const double ratio_sq);
    
    
    
    double F3twoloopgm2(const double ratio_sq);
    
    
    
    double F4twoloopgm2(const double ratio_sq);
    
    
    
    
    double F5twoloopgm2(const double ratio_sq);
    
    
    
    virtual double gminus2muNLO();
    
    
    
    
    
    //virtual double gminus2muNLOF();
    //WE'RE REMOVING THIS DEFINITIONS AND SUBSTITUTING THEM FROM THOSE OF 1502.04199
    
    //virtual double gminus2muNLOB();
    //WE'RE REMOVING THIS DEFINITIONS AND SUBSTITUTING THEM FROM THOSE OF 1502.04199
    
    
    
    
 
    virtual gslpp::complex negsquareroot(double x);
    
    virtual gslpp::complex negpow(double basis, double exp);
    
    virtual gslpp::complex neglog(gslpp::complex argument);
 
    
    
    
    
    
    
    //virtual gslpp::complex TF(double m1, double m2, double m3);
    //WE'RE REMOVING THIS DEFINITIONS AND SUBSTITUTING THEM FROM THOSE OF 1502.04199
    
    
    
    
    
    
    
    
    virtual std::vector<WilsonCoefficient>& CMdbs2();
//    virtual std::vector<WilsonCoefficient>& CMdbsp2();
 
    virtual std::vector<WilsonCoefficient>& CMdiujleptonknu(int i, int j, int k) ;
    
    
    virtual std::vector<WilsonCoefficient>& CMBMll(QCD::lepton lepton);
    
    
    virtual double C10Bll(double xt, double xHp, gslpp::complex su);
    
    virtual  gslpp::complex CSboxBll(double xt, double xHp, gslpp::complex su, gslpp::complex sd, gslpp::complex sl);
    
    virtual  gslpp::complex CPboxBll(double xt, double xHp, gslpp::complex su, gslpp::complex sd, gslpp::complex sl);
    
    virtual  gslpp::complex CPZUBll(double xt, double xHp, double sW2, gslpp::complex su, gslpp::complex sd);
    
   
    virtual  double f1(double xHp, double xt);
    
        
    virtual  double f2(double xHp, double xt);
    
        
    virtual  double f3(double xHp, double xt);
    
        
    virtual  double f4(double xHp, double xt);
    
        
    virtual  double f5(double xHp, double xt);
    
        
    virtual  double f6(double xHp, double xt);
    
        
    virtual  double f7(double xHp, double xt);
    
        
    virtual  double f8(double xHp, double xt);
   
    
        
    virtual  double f9(double xHp, double xt);
    
    virtual  double f10(double xHp, double xt);
    
    
    virtual gslpp::complex  g0(double xHp, double xt, gslpp::complex su, gslpp::complex sd);
    
    virtual gslpp::complex  g1a(double xHp, double xt, gslpp::complex su, gslpp::complex sd);
    
    virtual gslpp::complex  g2a(double xHp, double xt, gslpp::complex su, gslpp::complex sd);
    
    virtual gslpp::complex  g3a(double xHp, double xt, gslpp::complex su, gslpp::complex sd);
    
    
    virtual gslpp::complex lambdaHHphi(double lambda3, double Relambda7,double Imlambda7, double Ri1, double Ri2, double Ri3 );
    
    virtual gslpp::complex CphiU(double xHp, double xt, double vev, double xphi, double mu, double Ri1, double Ri2, double Ri3, double mHi_2, double lambda3, double Relambda7,double Imlambda7, gslpp::complex su, gslpp::complex sd);
    
    virtual std::vector<WilsonCoefficient>& CMbsg() ;
    
    gslpp::complex setWCbsg (int i, gslpp::complex sigmau, gslpp::complex sigmad, double mt, double mhp, double mu, orders order);
    
     const Polylogarithms getPolyLog() const
       {
           return PolyLog;
       }
 
    void updateGTHDMParameters();
 
private:
    const GeneralTHDM & myGTHDM;
 
    gslpp::matrix<gslpp::complex> myCKM;
    WilsonCoefficient mcdbs2, mculeptonnu, mcBMll, mcbsg, mcgminus2mu, mcbsmm;
 
    double GF, mMU;
    gslpp::complex CWbsgArrayLO[8], CWbsgArrayNLO[8], CWbsgArrayNNLO[8];
    double mtbsg, mhpbsg, mubsg; // caching
    gslpp::complex su, sd, sl; // caching
 
    const Polylogarithms PolyLog;
};
 
#endif  /* GENERALTHDMMATCHING_H */