bsgamma.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2015 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef BSGAMMA_H
#define BSGAMMA_H
 
#include "ThObservable.h"
#include <gsl/gsl_integration.h>
#include <Polylogarithms.h>
#include <ClausenFunctions.h>
 
class Bsgamma : public ThObservable {
public: 
    Bsgamma(const StandardModel& SM_i, QCD::quark quark_i, int obsFlag);
    
    
    Bsgamma(const StandardModel& SM_i, int obsFlag);
    
    
    double delta(double E0);
    
    
    double rho(double E0);
    
    
    double omega(double E0);
    
    
    double T1(double E0, double t);
    
    
    double T2(double E0, double t);
    
    
    double T3(double E0, double t);
    
    
    double P0_4body(double E0, double t);
    
    
    double zeta();
    
    
    gslpp::complex a(double z);
    
    
    gslpp::complex b(double z);
    
    
    gslpp::complex r1(int i, double z);
    
    
    gslpp::complex r1_ew(int i, double z);
    
    
    gslpp::complex Gamma_t(double t);
    
    
    gslpp::complex kappa(double Mq, double t);
    
    
    double getKc_abs2_t(double t)
    {
        return kappa(Mc,t).abs2() * t;
    };
    
    
    double getKc_abs2_1mt(double t)
    {
        return kappa(Mc,t).abs2() * (1. - t);
    };
    
    
    double getKc_abs2_t_1mt(double t)
    {
        return kappa(Mc,t).abs2() * t * (1. - t);
    };
    
    
    double getKc_abs2_1mt2(double t)
    {
        return kappa(Mc,t).abs2() * (1. - t) * (1. - t);
    };
    
    
    double getKc_re_t(double t)
    {
        return kappa(Mc,t).real() * t ;
    };
    
    
    double getKc_im_t(double t)
    {
        return kappa(Mc,t).imag() * t ;
    };
    
    
    double getKc_re_t_1mt(double t)
    {
        return kappa(Mc,t).real() * t * (1. - t);
    };
    
    
    double getKc_im_t_1mt(double t)
    {
        return kappa(Mc,t).imag() * t * (1. - t);
    };
    
    
    double getKc_re_t_1mt2(double t)
    {
        return kappa(Mc,t).real() * t * (1. - t) * (1. - t);
    };
    
    
    double getKc_im_t_1mt2(double t)
    {
        return kappa(Mc,t).imag() * t * (1. - t) * (1. - t);
    };
    
    
    double getKc_re_1mt(double t)
    {
        return kappa(Mc,t).real() * (1. - t);
    };
    
    
    double getKc_im_1mt(double t)
    {
        return kappa(Mc,t).imag() * (1. - t);
    };
    
    
    double getKc_re_1mt2(double t)
    {
        return kappa(Mc,t).real() * (1. - t) * (1. - t);
    };
    
    
    double getKc_im_1mt2(double t)
    {
        return kappa(Mc,t).imag() * (1. - t) * (1. - t);
    };
    
    
    double getKb_abs2_1mt(double t)
    {
        return kappa(Mb_kin,t).abs2() * (1. - t);
    };
    
    
    double getKb_abs2_1mt2(double t)
    {
        return kappa(Mb_kin,t).abs2() * (1. - t) * (1. - t);
    };
    
    
    double getKb_abs2_t_1mt(double t)
    {
        return kappa(Mb_kin,t).abs2() * t * (1. - t);
    };
    
    
    double getKb_abs2_t_1mt2(double t)
    {
        return kappa(Mb_kin,t).abs2() * t * (1. - t) * (1. - t);
    };
    
    
    double getKb_abs2_t2_1mt(double t)
    {
        return kappa(Mb_kin,t).abs2() * t * t * (1. - t);
    };
    
    
    double getKb_abs2_t2_1mt2(double t)
    {
        return kappa(Mb_kin,t).abs2() * t * t * (1. - t) * (1. - t);
    };
    
    
    double getKb_re_t(double t)
    {
        return kappa(Mb_kin,t).real() * t ;
    };
    
    
    double getKb_re_t_1mt(double t)
    {
        return kappa(Mb_kin,t).real() * t * (1. - t);
    };
    
    
    double getKb_re_t2_1mt(double t)
    {
        return kappa(Mb_kin,t).real() * t * t * (1. - t);
    };
    
    
    double getKb_re_t2_1mt2(double t)
    {
        return kappa(Mb_kin,t).real() * t * t * (1. - t) * (1. - t);
    };
    
    
    double getKb_re_t_1mt2(double t)
    {
        return kappa(Mb_kin,t).real() * t * (1. - t) * (1. - t);
    };
    
    
    double getKb_re_1mt(double t)
    {
        return kappa(Mb_kin,t).real() * (1. - t);
    };
    
    
    double getKb_re_1mt2(double t)
    {
        return kappa(Mb_kin,t).real() * (1. - t) * (1. - t);
    };
    
    
    double getKc_re_Kb_1mt(double t)
    {
        return kappa(Mc,t).real() * kappa(Mb_kin,t).real() * (1. - t);
    };
    
    
    double getKc_im_Kb_1mt(double t)
    {
        return kappa(Mc,t).imag() * kappa(Mb_kin,t).real() * (1. - t);
    };
    
    
    double getKc_re_Kb_1mt2(double t)
    {
        return kappa(Mc,t).real() * kappa(Mb_kin,t).real() * (1. - t) * (1. - t);
    };
    
    
    double getKc_im_Kb_1mt2(double t)
    {
        return kappa(Mc,t).imag() * kappa(Mb_kin,t).real() * (1. - t) * (1. - t);
    };
    
    
    double getKc_re_Kb_t_1mt(double t)
    {
        return kappa(Mc,t).real() * kappa(Mb_kin,t).real() * t * (1. - t);
    };
    
    
    double getKc_im_Kb_t_1mt(double t)
    {
        return kappa(Mc,t).imag() * kappa(Mb_kin,t).real() * t * (1. - t);
    };
    
    
    double getKc_re_Kb_t_1mt2(double t)
    {
        return kappa(Mc,t).real() * kappa(Mb_kin,t).real() * t * (1. - t) * (1. - t);
    };
    
    
    double getKc_im_Kb_t_1mt2(double t)
    {
        return kappa(Mc,t).imag() * kappa(Mb_kin,t).real() * t * (1. - t) * (1. - t);
    };
    
    
    double Int_b1(double E0);
    
    
    double Int_b2(double E0);
    
    
    double Int_b3(double E0);
    
    
    double Int_b4(double E0);
    
    
    double Int_bb1(double E0);
    
    
    double Int_bb2(double E0);
    
    
    double Int_bb4(double E0);
    
    
    gslpp::complex Int_bc1(double E0);
    
    
    gslpp::complex Int_bc2(double E0);
    
    
    gslpp::complex Int_c1(double E0);
    
    
    gslpp::complex Int_c2(double E0);
    
    
    gslpp::complex Int_c3(double E0);
    
    
    double Int_cc(double E0);
    
    
    double Int_cc1(double E0);
    
    
    double Int_cc1_part1(double E0);
    
    
    double ff7_dMP(double E0);
    
    
    double ff7_sMP(double E0);
    
    
    double ff8_dMP(double E0);
    
    
    double ff8_sMP(double E0);
    
    
    double Phi11_1(double E0);
    
    
    double Phi12_1(double E0);
    
    
    gslpp::complex Phi13_1(double E0);
    
    
    gslpp::complex Phi14_1(double E0);
    
    
    gslpp::complex Phi15_1(double E0);
    
    
    gslpp::complex Phi16_1(double E0);
    
    
    gslpp::complex Phi17_1(double E0, double z);
    
    
    gslpp::complex Phi18_1(double E0, double z);
    
    
    double Phi22_1(double E0);
    
    
    double Phi23_1_4body(double E0);
    
    
    gslpp::complex Phi23_1(double E0);
    
    
    double Phi24_1_4body(double E0);
    
    
    gslpp::complex Phi24_1(double E0);
    
    
    double Phi25_1_4body(double E0);
    
    
    gslpp::complex Phi25_1(double E0);
    
    
    double Phi26_1_4body(double E0);
    
    
    gslpp::complex Phi26_1(double E0);
    
    
    gslpp::complex Phi27_1(double E0, double z);
    
    
    gslpp::complex Phi28_1(double E0, double z);
    
    double Phi33_1(double E0);
    
    double Phi34_1(double E0);
    
    double Phi35_1(double E0);
    
    gslpp::complex Phi36_1(double E0);
    
    double Phi37_1(double E0);
    
    double Phi38_1(double E0);
    
    double Phi44_1(double E0);
    
    double Phi45_1(double E0);
    
    gslpp::complex Phi46_1(double E0);
    
    
    double Phi47_1(double E0);
    
    double Phi48_1(double E0);
    
    double Phi55_1(double E0);
    
    gslpp::complex Phi56_1(double E0);
    
    
    double Phi57_1(double E0);
    
    double Phi58_1(double E0);
    
    gslpp::complex Phi66_1(double E0);
    
    
    gslpp::complex Phi67_1(double E0);
    
    gslpp::complex Phi68_1(double E0);
    
    
    double Phi77_1(double E0);
    
    
    double Phi78_1(double E0);
    
    
    double Phi88_1(double E0);
    
    
    gslpp::complex Kij_1(int i, int j, double E0, double mu);
    
    
    double Rer22(double z);
    
    
    double Phi22_2beta0(double E0, double mu);
    
    
    double Phi28_2beta0(double E0, double mu);
    
    
    double Phi77_2beta0(double E0, double mu);
    
    
    double Phi88_2beta0(double E0, double mu);
    
    
    double dY1(double E0);
    
    
    double Y1(double E0, double mu);
    
    
    double Y2CF(double E0, double mu);
    
    
    double Y2CA(double E0, double mu);
    
    
    double Y2NL(double E0, double mu);
    
    
    double Y2NV_PHI1(double rho);
    
    
    double Y2NV_PHI2(double rho);
    
    
    double Y2NV_PHI3(double rho);
    
    
    double Y2NV_PHI4(double rho);
    
    
    double Y2NV(double E0, double mu);
    
    
    double Y2NH(double E0, double mu);
    
    
    double Y2(double E0, double mu);
    
    
    double f_NLO_1(double z);
    
    
    double zdz_f_NLO(double z, double E0);
    
    
    double mddel_f_NLO(double z, double E0);
    
    
    double h27_2(double z, double E0);
    
    
    double f_q(double z, double E0);
    
    
    double f_b(double z);
    
    
    double f_c(double z);
    
    
    double F_1(double z);
    
    
    double F_2(double z);
    
    
    double delddel_Phi22_1(double E0);
    
    
    double zdz_Phi22_1(double E0);
    
    
    double delddel_Phi28_1(double z, double E0);
    
    
    double zdz_Phi28_1(double z, double E0);
    
    
    double delddel_Phi88_1(double E0);
    
    
    double f(double r);
    
    
    double Delta(double r);
    
    
    double f_u(double r);
    
    
    double omega77(double z);
    
    
    double Int_Phi77_2rem(double E0);
    
    
    double Phi77_2rem(double E0);
    
    
    double K77_2_z1(double E0, double mu);
    
    
    double Kij_2(int i, int j, double E0, double mu_b, double mu_c);
    
    
    void computeCoeff(double mu);
    
    
    double P0(double E0);
    
    
    double P11();
    
    
    double P21(double E0, double mu);
    
    
    double P21_CPodd(double E0, double mu);
    
    
    double P12();
    
    
    double P22(double E0, double mu_b, double mu_c);
    
    
    double P32(double E0, double mu);
    
    
    double EW_NLO(double mu);
    
    
    double Vub_NLO_2body();
    
    
    double Vub_NLO_2body_CPodd();
    
    
    double Vub_NLO_3body_A(double E0);
    
    
    double Vub_NLO_3body_A_CPodd(double E0);
    
    
    double Vub_NLO_3body_B(double E0);
    
    
    double Vub_NLO_3body_B_CPodd(double E0);
    
    
    double Vub_NLO_4body(double E0);
    
    
    double Vub_NLO_4body_CPodd(double E0);
    
    
    double Vub_NLO(double E0);
    
    
    double Vub_NLO_CPodd(double E0);
    
    
    double Vub_NNLO(double E0);
    
    
    double P(double E0, double mu_b, double mu_c, orders order);
    
    
    double N_27();
    
    
    double N_77(double E0, double mu);
    
    
    double N(double E0, double mu);
    
    
    double C_sem();
    
    
    void updateParameters();
    
    
    double computeThValue();
    
    
private:
    QCD::quark quark;
    bool SUM;
    bool EWflag;
    bool FOUR_BODY;
    bool WET_NP_btos = false, SMEFT_NP_btos = false;
    
    double ale; 
    double AleatMztilde; 
    double alsUps; 
    double Alstilde; 
    double mu_b; 
    double mu_c; 
    double mu_kin; 
    double Mb_kin; 
    double Mc; 
    double Ms;
    double Mb;
    double Mz;
    double BRsl; 
    double C; 
    double CKMratio; 
    gslpp::complex V_ub; 
    gslpp::complex V_cb; 
    gslpp::complex V_tb; 
    gslpp::complex CKMu; 
    double CKMusq; 
    double overall; 
    double mu_pi2; 
    double mu_G2; 
    double rho_D3; 
    double rho_LS3; 
    double BLNPcorr; 
    int obs; 
//    double BR; /**<BR of the decay */
//    double BR_CPodd; /**<BR of the decay */
    
    gslpp::vector<gslpp::complex> ** allcoeff;
    gslpp::vector<gslpp::complex> ** allcoeffprime;
    gslpp::complex C1_0;
    gslpp::complex C2_0;
    gslpp::complex C3_0;
    gslpp::complex C4_0;
    gslpp::complex C5_0;
    gslpp::complex C6_0;
    gslpp::complex C7_0;
    gslpp::complex C8_0;
    gslpp::complex C1_1;
    gslpp::complex C2_1;
    gslpp::complex C3_1;
    gslpp::complex C4_1;
    gslpp::complex C5_1;
    gslpp::complex C6_1;
    gslpp::complex C7_1;
    gslpp::complex C8_1;
    gslpp::complex C7_1ew;
    gslpp::complex C7_2;
    gslpp::complex C7p_0;
    gslpp::complex C7p_1;
    gslpp::complex C_7_NP;
    gslpp::complex C_7p_NP;
    
    gsl_function INT;
    gsl_integration_cquad_workspace * w_INT;
    double avaINT;
    double errINT;
    unsigned int Intb1Cached;
    unsigned int Intb2Cached;
    unsigned int Intb3Cached;
    unsigned int Intb4Cached;
    unsigned int Intbb1Cached;
    unsigned int Intbb2Cached;
    unsigned int Intbb4Cached;
    unsigned int Intbc1Cached;
    unsigned int Intbc2Cached;
    unsigned int Intc1Cached;
    unsigned int Intc1imCached;
    unsigned int Intc2Cached;
    unsigned int Intc3Cached;
    unsigned int IntccCached;
    unsigned int Intcc1Cached;
    unsigned int Intcc1p1Cached;
    unsigned int IntPhi772rCached;
    double CacheIntb1;
    double CacheIntb2;
    double CacheIntb3;
    double CacheIntb4;
    double CacheIntbb1;
    double CacheIntbb2;
    double CacheIntbb4;
    gslpp::complex CacheIntbc1;
    gslpp::complex CacheIntbc2;
    gslpp::complex CacheIntc1;
    gslpp::complex CacheIntc2;
    gslpp::complex CacheIntc3;
    double CacheIntcc;
    double CacheIntcc1;
    double CacheIntcc1p1;
    double CacheIntPhi772r;
    unsigned int Intb_updated;
    unsigned int Intbc_updated;
    double Intb_cache;
    gslpp::vector<double> Intbc_cache;
    void checkCache();
    
};
 
#endif  /* BSGAMMA_H */