Flavour/src/bsgamma.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2015 HEPfit Collaboration
 * All rights reserved.
 *
 * 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>

#define FOUR_BODY false

class Bsgamma : public ThObservable {
public: 
    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 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_re_t_1mt(double t)
    {
        return kappa(Mc,t).real() * t * (1. - t);
    };
    
    
    double getKc_re_t_1mt2(double t)
    {
        return kappa(Mc,t).real() * 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_re_Kb_1mt2(double t)
    {
        return kappa(Mc,t).real() * 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_re_Kb_t_1mt2(double t)
    {
        return kappa(Mc,t).real() * 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);
    
    
    double Int_bc1(double E0);
    
    
    double Int_bc2(double E0);
    
    
    double Int_c1(double E0);
    
    
    double Int_c1_im(double E0);
    
    
    double Int_c2(double E0);
    
    
    double 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);
    
    
    double Phi13_1(double E0);
    
    
    double Phi14_1(double E0);
    
    
    double Phi15_1(double E0);
    
    
    double Phi16_1(double E0);
    
    
    double Phi17_1(double E0, double z);
    
    
    double Phi18_1(double E0, double z);
    
    
    double Phi22_1(double E0);
    
    
    double Phi23_1_4body(double E0);
    
    
    double Phi23_1(double E0);
    
    
    double Phi24_1_4body(double E0);
    
    
    double Phi24_1(double E0);
    
    
    double Phi25_1_4body(double E0);
    
    
    double Phi25_1(double E0);
    
    
    double Phi26_1_4body(double E0);
    
    
    double Phi26_1(double E0);
    
    
    double Phi27_1(double E0, double z);
    
    
    double Phi27_1_im(double E0, double z);
    
    
    double Phi28_1(double E0, double z);
    
    double Phi33_1(double E0);
    
    double Phi34_1(double E0);
    
    double Phi35_1(double E0);
    
    double Phi36_1(double E0);
    
    double Phi37_1(double E0);
    
    double Phi38_1(double E0);
    
    double Phi44_1(double E0);
    
    double Phi45_1(double E0);
    
    double Phi46_1(double E0);
    
    
    double Phi47_1(double E0);
    
    double Phi48_1(double E0);
    
    double Phi55_1(double E0);
    
    double Phi56_1(double E0);
    
    
    double Phi57_1(double E0);
    
    double Phi58_1(double E0);
    
    double Phi66_1(double E0);
    
    
    double Phi67_1(double E0);
    
    double Phi68_1(double E0);
    
    
    double Phi77_1(double E0);
    
    
    double Phi78_1(double E0);
    
    
    double Phi88_1(double E0);
    
    
    double Kij_1(int i, int j, double E0, double mu);
    
    
    void computeCoeff(double mu);
    
    
    double P0(double E0);
    
    
    double P11();
    
    
    double P21(double E0, double mu);
    
    
    double P12();
    
    
    double P22(double E0, double mu_b, double mu_c);
    
    
    double P32(double E0, double mu);
    
    
    double Vub_NLO_2body(bool CPodd);
    
    
    double Vub_NLO_3body(double E0, bool CPodd);
    
    
    double Vub_NLO_4body(double E0, bool CPodd);
    
    
    double Vub_NLO(double E0, bool CPodd);
    
    
    double P(double E0, double mu_b, double mu_c, orders order, bool CPodd);
    
    
    double N_27();
    
    
    double N_77(double E0, double mu);
    
    
    double N(double E0, double mu);
    
    
    double C_sem();
    
    
    void updateParameters();
    
    
    double computeThValue();
    
    
private:
    double ale; 
    double alsUps; 
    double Alstilde; 
    double E0; 
    double mu_b; 
    double mu_c; 
    double mu_kin; 
    double Mb_kin; 
    double Mc; 
    double Ms;
    double BRsl; 
    double C; 
    gslpp::complex lambda_t; 
    gslpp::complex CKMu; 
    gslpp::complex V_cb; 
    double overall; 
    double mu_pi2; 
    double mu_G2; 
    double rho_D3; 
    double rho_LS3; 
    double BLNPcorr; 
    int obs; 
    double BR; 
    double BR_CPodd; 
    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_2;
    gslpp::complex C7p_0;
    gslpp::complex C7p_1;
    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;
    double CacheIntbc1;
    double CacheIntbc2;
    double CacheIntc1;
    double CacheIntc1im;
    double CacheIntc2;
    double 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 */