BXqll.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2016 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef BXQLL_H
#define BXQLL_H
 
class StandardModel;
 
#include "QCD.h"
#include "ThObservable.h"
//#include "Particle.h"
#include "gslpp.h"
#include "Expanded.h"
#include <gsl/gsl_integration.h>
#include "HeffDF1.h"
#include "F_1.h"
#include "F_2.h"
 
class BXqll {
public:
    
    BXqll(const StandardModel& SM_i, QCD::quark quark_i, QCD::lepton lep_i);
    
    virtual ~BXqll();  
 
    std::vector<std::string> initializeBXqllParameters();
    
    double getR_LOWQ2(double sh);
    
    double getH(std::string obs, double sh);
    
    double integrateH(std::string obs, double q_min, double q_max);
 
private:
    const StandardModel& mySM;
    F_1 myF_1;
    F_2 myF_2;
    HeffDF1 myHeff;
    QCD::lepton lep;
    QCD::quark quark;
    double CF, GF, ale, alsmu, alsmuc, alstilde, aletilde, kappa;
    double Mlep, mu_b, mu_c, Mb, Mc, Mtau, Mb_pole, Mc_pole, Ms, MW;
    double abslambdat_over_Vcb, Vts_over_Vcb, z, muh, lambda_1, lambda_2, Lbl;
    double phi1, phi2, phi00, phi00_2, phi20, phi01;
    double phinv00, phinv10, phinv20, phinv01, phinv11, phinv21, phinv_00_01;
    double BR_BXcenu, C_ratio, pre;
    unsigned int QCD_max, QED_max;
   
    std::vector< gslpp::vector<gslpp::complex> > M_7;
    std::vector< gslpp::vector<gslpp::complex> > M_9;
    std::vector< gslpp::vector<gslpp::complex> > M_10;
    
    std::vector< gslpp::matrix<gslpp::complex> > Hij_T;
    std::vector< gslpp::matrix<gslpp::complex> > Hij_L;
    std::vector< gslpp::matrix<gslpp::complex> > Hij_A;
 
    std::vector<std::string> BXqllParameters;
//    gslpp::vector<gslpp::complex> ** allcoeff_smm;/**<Vector that contains the Wilson coeffients */
    Expanded<gslpp::vector<gslpp::complex> > allcoeff;
    gslpp::matrix<gslpp::complex> WC;
    double aveH;
    double errH;
    gsl_function FH;
    gsl_integration_cquad_workspace * w_H;
    gsl_error_handler_t * old_handler; 
    void updateParameters();
 
    double F_17re(double muh, double z, double sh, int maxpow=20);
    double F_17im(double muh, double z, double sh, int maxpow=20);
    double F_19re(double muh, double z, double sh, int maxpow=20);
    double F_19im(double muh, double z, double sh, int maxpow=20);
    double F_27re(double muh, double z, double sh, int maxpow=20);
    double F_27im(double muh, double z, double sh, int maxpow=20);
    double F_29re(double muh, double z, double sh, int maxpow=20);
    double F_29im(double muh, double z, double sh, int maxpow=20);
    
    double DeltaF_19re(double muh, double z, double sh, int maxpow=20);
    double DeltaF_19im(double muh, double z, double sh, int maxpow=20);
    double DeltaF_29re(double muh, double z, double sh, int maxpow=20);
    double DeltaF_29im(double muh, double z, double sh, int maxpow=20);
    gslpp::complex F17(double sh);
 
    gslpp::complex F27(double sh);
 
    gslpp::complex F19(double sh);
 
    gslpp::complex F29(double sh);
 
    gslpp::complex F87(double sh);
 
    double F89(double sh);
    
    double S77_T(double sh, orders order);
    double S79_T(double sh, orders order);
    double S99_T(double sh, orders order);
    double S1010_T(double sh, orders order);
    
    double S77_L(double sh, orders order);
    double S79_L(double sh, orders order);
    double S99_L(double sh, orders order);
    double S1010_L(double sh, orders order);
    
    double S710_A(double sh, orders order);
    double S910_A(double sh, orders order);
    
    gslpp::complex cij_T(unsigned int i, unsigned int j, double sh, unsigned int ord);
    gslpp::complex cij_L(unsigned int i, unsigned int j, double sh, unsigned int ord);
    gslpp::complex cij_A(unsigned int i, unsigned int j, double sh);
    
    gslpp::complex eij_T(unsigned int i, unsigned int j, double sh);
    gslpp::complex eij_L(unsigned int i, unsigned int j, double sh);
    gslpp::complex eij_A(unsigned int i, unsigned int j, double sh);
    
    double omega77_T(double sh);
    double omega79_T(double sh);
    double omega99_T(double sh);
    
    double omega77_L(double sh);
    double omega79_L(double sh);
    double omega99_L(double sh);
    
    double omega710_A(double sh);
    double omega910_A(double sh);
    
    double omega77em_T(double sh);
    double omega79em_T(double sh);
    double omega99em_T(double sh);
    double omega22em_T(double sh);
    gslpp::complex omega27em_T(double sh);
    gslpp::complex omega29em_T(double sh);
    
    double omega77em_L(double sh);
    double omega79em_L(double sh);
    double omega99em_L(double sh);
    double omega22em_L(double sh);
    gslpp::complex omega27em_L(double sh);
    gslpp::complex omega29em_L(double sh);
    
    double omega710em_A(double sh);
    double omega910em_A(double sh);
    gslpp::complex omega210em_A(double sh);
    
    gslpp::complex f_Huber(double sh, double gamma_9, double rho_c, double rho_b, double rho_0, double rho_num);
    
    gslpp::complex f9pen_Huber(double sh);
    
    gslpp::complex g_Huber(double y);
    
    gslpp::complex KS_cc(double sh);
    
    gslpp::complex KS_aux(double sh, double m, double Gamma, double Br_ll, double Br_had);
    
    gslpp::complex F_BIR(double r);
    
    void computeMi(double sh);
    
    double H_T (double sh);
    
    double H_L (double sh);
    
    double H_A (double sh);
 
    double Phi_u(orders ord);
    double Phi_u(orders_qed ord_qed);
    
    double Phi_u_inv(unsigned int ord_qcd, unsigned int ord_qed);
    
    double PhiTL_brems(double sh);
    
    double PhiA_brems(double sh);
    
    gslpp::complex C9mod(double sh);
    
    gslpp::complex h_z(double zed, double sh);
    
    double tau78(double sh);
    
    double tau89(double sh);
    
    double tau88(double sh);
    
    double tau22fit(double sh);
    
    double tau27fit_Re(double sh);
    
    double tau27fit_Im(double sh);
    
    double tau28fit_Re(double sh);
    
    double tau28fit_Im(double sh);
    
    double tau29fit_Re(double sh);
    
    double tau29fit_Im(double sh);
    
    double t810(double sh);
    
    double t210fit(double sh);
    
    unsigned int int_qed(orders_qed order_qed);
    
    double CCH_multiplication(std::vector< gslpp::matrix<gslpp::complex> >& Hij);
    
    void Test_WC_DF1();
    
    friend double gslpp_special_functions::dilog(double x);
};
#endif  /* BXqLL_H */