THDMWMatching.h

Go to the documentation of this file.

/*
 * Copyright (C) 2018 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
 
#ifndef THDMWMATCHING_H
#define THDMWMATCHING_H
 
#include <Polylogarithms.h>
#include "gslpp.h"
#include "StandardModelMatching.h"
 
class THDMW;
 
class THDMWMatching : public StandardModelMatching {
public:
    THDMWMatching(const THDMW & THDMW_i);
 
 
    virtual gslpp::complex C10NP(double xt, double xS, gslpp::complex etaU);
    
    virtual gslpp::complex CSNP(double nu1, double xh, double xt, double xS, gslpp::complex etaU,  gslpp::complex etaD);
    
    
    virtual gslpp::complex CPNP(double xt, double xS, gslpp::complex etaU,  gslpp::complex etaD);
    
    
    
    
    virtual  std::vector<WilsonCoefficient>& CMdbs2();
 
    
 
    virtual std::vector<WilsonCoefficient>& CMBMll(QCD::lepton lepton);
    
    
    
    
    virtual double CNPVLLctEtadEtasu(double xc, double xb, double xt, double xS);
    
    virtual double CNPVLLctEtau4(double xc, double xb, double xt, double xS);
    
    virtual double CNPVLLctEtau2(double xc, double xb, double xt, double xS);
    
    virtual gslpp::complex CNPVLLct(double xc, double xb, double xt, double xS,  gslpp::complex etaU, gslpp::complex etaD);
    
    virtual double CNPVLLttEtadEtasu(double xc, double xb, double xt, double xS);
    
    virtual double CNPVLLttEtau4(double xc, double xb, double xt, double xS);
    
    virtual double CNPVLLttEtau2(double xc, double xb, double xt, double xS);
    
    virtual gslpp::complex CNPVLLtt(double xc, double xb, double xt, double xS,  gslpp::complex etaU, gslpp::complex etaD);
 
    virtual double CNPVLLccEtadEtasu(double xc, double xb, double xS);
    
    virtual double CNPVLLccEtau2(double xc, double xb, double xS);
    
    virtual gslpp::complex CNPVLLcc(double xc, double xb, double xS,  gslpp::complex etaU, gslpp::complex etaD);
 
    virtual gslpp::complex CNPVLL(double xc, double xb, double xt, double xS,  gslpp::complex etaU, gslpp::complex etaD);
    
    
    
    
    
    
    virtual double CNPSRR1ctEtadEtasuEtau2(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR1ctEtad2Etasu2(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR1ctEtadEtasu(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR1ctEtau4(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR1ctEtau2(double xc, double xb, double xt, double xS);
    
    virtual gslpp::complex CNPSRR1ct(double xc, double xb, double xt, double xS,  gslpp::complex etaU, gslpp::complex etaD);
 
    
    
    
    
    
    
    virtual double CNPSRR1ttEtadEtasuEtau2(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR1ttEtad2Etasu2(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR1ttEtadEtasu(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR1ttEtau4(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR1ttEtau2(double xc, double xb, double xt, double xS);
    
    virtual gslpp::complex CNPSRR1tt(double xc, double xb, double xt, double xS,  gslpp::complex etaU, gslpp::complex etaD);
 
    virtual gslpp::complex CNPSRR1(double xc, double xb, double xt, double xS,  gslpp::complex etaU, gslpp::complex etaD);
 
    
 
    
    
    
    
    
    virtual double CNPSRR2ctEtadEtasuEtau2(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR2ctEtad2Etasu2(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR2ctEtadEtasu(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR2ctEtau4(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR2ctEtau2(double xc, double xb, double xt, double xS);
    
    virtual gslpp::complex CNPSRR2ct(double xc, double xb, double xt, double xS,  gslpp::complex etaU, gslpp::complex etaD);
 
    
    
    
    
    
    virtual double CNPSRR2ttEtadEtasuEtau2(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR2ttEtad2Etasu2(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR2ttEtadEtasu(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR2ttEtau4(double xc, double xb, double xt, double xS);
    
    virtual double CNPSRR2ttEtau2(double xc, double xb, double xt, double xS);
    
    virtual gslpp::complex CNPSRR2tt(double xc, double xb, double xt, double xS,  gslpp::complex etaU, gslpp::complex etaD);
 
    virtual gslpp::complex CNPSRR2(double xc, double xb, double xt, double xS,  gslpp::complex etaU, gslpp::complex etaD);
 
    
    
    
    
    
    
    
    
    
    virtual  double f1(double xc, double xt, double xS);
  
    virtual  double f2(double xb, double xt, double xS);
    
    virtual  double f3(double xb, double xt, double xS);
    
    virtual  double f4(double xb, double xt, double xS);
    
    virtual  double f5(double xc, double xt, double xS);
    
    virtual  double f6(double xb, double xt, double xS);
    
    virtual  double f7(double xb, double xt, double xS);
    
    virtual  double f8(double xb, double xt);
    
    virtual  double f9(double xc ,double xb, double xt, double xS);
    
    virtual  double f10( double xt, double xS);
    
    virtual  double f11( double xt, double xS);
    
    virtual  double f12( double xt, double xS);
    
    virtual  double f13( double xt, double xS);
    
    virtual  double f14(double xb, double xt, double xS);
    
    virtual  double f15(double xb, double xt, double xS);
    
    virtual  double f16(double xb, double xt, double xS);
    
    virtual  double f17(double xt, double xS);
    
    virtual  double f18(double xb, double xt, double xS);
    
    virtual  double f19(double xb, double xt, double xS);
    
    virtual  double f20(double xb, double xt, double xS);
    
    virtual  double f21(double xb, double xt);
    
    virtual  double f22(double xb, double xt, double xS);
    
    virtual  double f23(double xb, double xt, double xS);
    
    virtual  double f24(double xt, double xS);
    
    virtual  double f25(double xb, double xt, double xS);
    
    virtual  double f26(double xb, double xt, double xS);
    
    virtual  double f27(double xb, double xt, double xS);
    
    virtual  double f28(double xb, double xt);
    
    virtual  double f29(double xt, double xS);
    
    virtual  double f30(double xb, double xt, double xS);
    
    virtual  double f31(double xb, double xt, double xS);
    
    virtual  double f32(double xc, double xb, double xt, double xS);
    
    virtual  double f33(double xc, double xb, double xt, double xS);
    
    virtual  double f34(double xt, double xS);
    
    virtual  double f35(double xc, double xb, double xt, double xS);
    
    virtual  double f36(double xt, double xS);
    
    virtual  double f37(double xt, double xS);
    
    virtual  double f38(double xb, double xt, double xS);
    
    virtual  double f39(double xb, double xt, double xS);
    
    virtual  double f40(double xb, double xt, double xS);
    
    virtual  double f41(double xt, double xS);
    
    virtual  double f42(double xt, double xS);
    
    virtual  double f43(double xt, double xS);
    
    
    
    
     const Polylogarithms getPolyLog() const
       {
           return PolyLog;
       }
 
    void updateTHDMWParameters();
 
private:
    const THDMW & myTHDMW;
 
    gslpp::matrix<gslpp::complex> myCKM;
    WilsonCoefficient mcBMll,mcbsg , mcdbs2;
 
    //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 /* THDMWMATCHING_H */