d3/d50/_b_xqll_8h_source.html

/*

 * 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, phi1, phi2;

    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);


    void computeHij_T(double sh);


    void computeHij_L(double sh);


    void computeHij_A(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 Phi_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);


    unsigned int int_qed(orders_qed order_qed);


    double CCH_multiplication(std::vector< gslpp::matrix<gslpp::complex> >& Hij);


    double FULLCCH_multiplication(std::vector< gslpp::matrix<gslpp::complex> >& Hij);


    void Test_WC_DF1();


    friend double gslpp_special_functions::dilog(double x);

};

#endif  /* BXqLL_H */