MVlnu.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2014 HEPfit Collaboration
 * 
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef MVLNU_H
#define MVLNU_H
 
class StandardModel;
 
#include <TF1.h>
#include <TGraph.h>
#include <TFitResultPtr.h>
#include <gsl/gsl_spline.h>
#include "Math/GSLIntegrator.h"
#include "Math/Functor.h"
#include <memory>
 
class MVlnu {
public:
    MVlnu(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i);
    
    virtual ~MVlnu();
    std::vector<std::string> initializeMVlnuParameters();
    
    double getDeltaGammaDeltaw(double w_min, double w_max);
    
    double getDeltaGammaDeltacl(double cl_min, double cl_max);
    
    double getDeltaGammaDeltacV(double cV_min, double cV_max);
    
    double getDeltaGammaDeltachi(double chi_min, double chi_max);
    
    double getFL();
 
    double get_unitarity_V_BGL();
 
    double get_unitarity_A_BGL();
 
    double get_unitarity_P_BGL();
 
    double get_hA1w1();
    
    double get_hA1(double w);
    
    double get_hA2(double w);
    
    double get_hA3(double w);
    
    double get_hV(double w);
    
    double get_R1(double w);
    
    double get_R2(double w);
    
    double get_R0(double w);
 
    double getPlep();
 
    double getMwidth(){
        updateParameters();
        return width;
    }
    
private:
    const StandardModel& mySM;
    QCD::lepton lep;
    QCD::meson meson;
    QCD::meson vectorM;
    std::vector<std::string> mvlnuParameters;
    bool CLNflag; 
    bool BGLflag; 
    bool DMflag; 
    bool btocNPpmflag; 
    bool NPanalysis; 
    double MV_o_MM;
    double sqrtMV_o_MM;
    
    double GF;            
    double Mlep;          
    double Mnu;           
    double MM;            
    double MV;            
    double w0;            
    double RV;            
    double mu_b;          
    double Mb;            
    double Mc;            
    double width;         
    double ale_mub;   
    gslpp::complex Vcb;   
    double amplsq_factor;   
    double q2min, q2max; 
    double CV_SM; 
    gslpp::complex CS; 
    gslpp::complex CP; 
    gslpp::complex CSp; 
    gslpp::complex CPp; 
    gslpp::complex CV; 
    gslpp::complex CA; 
    gslpp::complex CVp; 
    gslpp::complex CAp; 
    gslpp::complex C7; 
    gslpp::complex C7p; 
    gslpp::complex CT; 
    gslpp::complex CTp; 
    double hA1w1,rho2,R1w1,R2w1; 
    double N_A, N_1, N_2, j_A, j_0, j_1, j_2, k_A, k_0, k_1, k_2, l_A; 
    double af0,af1,af2,ag0,ag1,ag2,aF11,aF12,aF13,aF21,aF22,aF23; 
    double mBcstV1,mBcstV2,mBcstV3,mBcstV4,mBcstA1,mBcstA2,mBcstA3,mBcstA4,
        mBcstP1,mBcstP2,mBcstP3,chiTV,chiTA,chiTP,nI; 
    double zV1,zV2,zV3,zV4,zA1,zA2,zA3,zA4,zP1,zP2,zP3;
    double af_1, ag_1, aF1_1, aF2_1, af_2, ag_2, aF1_2, aF2_2, af_3, ag_3, aF1_3, aF2_3, af_4, ag_4, 
        aF1_4, aF2_4, af_5, ag_5, aF1_5, aF2_5, af_6, ag_6, aF1_6, aF2_6, af_7, ag_7, aF1_7, aF2_7,
        af_8, ag_8, aF1_8, aF2_8, af_9, ag_9, aF1_9, aF2_9, af_10, ag_10, aF1_10, aF2_10,
        bf_1, bg_1, bF1_1, bF2_1, bf_2, bg_2, bF1_2, bF2_2, bf_3, bg_3, bF1_3, bF2_3, bf_4, bg_4, 
        bF1_4, bF2_4, bf_5, bg_5, bF1_5, bF2_5, bf_6, bg_6, bF1_6, bF2_6, bf_7, bg_7, bF1_7, bF2_7,
        bf_8, bg_8, bF1_8, bF2_8, bf_9, bg_9, bF1_9, bF2_9, bf_10, bg_10, bF1_10, bF2_10; 
    double cached_intJ1s_tau, cached_intJ1c_tau, cached_intJ2s_tau, cached_intJ2c_tau, cached_intJ3_tau, cached_intJ4_tau,
                cached_intJ5_tau, cached_intJ6s_tau, cached_intJ6c_tau, cached_intJ7_tau, cached_intJ8_tau, cached_intJ9_tau,
                cached_intJ1s_mu, cached_intJ1c_mu, cached_intJ2s_mu, cached_intJ2c_mu, cached_intJ3_mu,
                cached_intJ4_mu, cached_intJ5_mu, cached_intJ6s_mu, cached_intJ6c_mu, cached_intJ7_mu, cached_intJ8_mu, cached_intJ9_mu,
                cached_intJ1s_el, cached_intJ1c_el, cached_intJ2s_el, cached_intJ2c_el, cached_intJ3_el, cached_intJ4_el,
                cached_intJ5_el, cached_intJ6s_el, cached_intJ6c_el, cached_intJ7_el, cached_intJ8_el, cached_intJ9_el; 
    double hA1w1_cache,rho2_cache,R1w1_cache,R2w1_cache;
    double N_A_cache, N_1_cache, N_2_cache, j_A_cache, j_0_cache, j_1_cache, j_2_cache, k_A_cache, k_0_cache, k_1_cache, k_2_cache, l_A_cache;
    double af0_cache,af1_cache,af2_cache,ag0_cache,ag1_cache,ag2_cache,aF11_cache,aF12_cache,aF13_cache,aF21_cache,aF22_cache,aF23_cache;
    double af_1_cache, ag_1_cache, aF1_1_cache, aF2_1_cache, af_2_cache, ag_2_cache, aF1_2_cache, aF2_2_cache, 
        af_3_cache, ag_3_cache, aF1_3_cache, aF2_3_cache, af_4_cache, ag_4_cache, aF1_4_cache, aF2_4_cache, 
        af_5_cache, ag_5_cache, aF1_5_cache, aF2_5_cache, af_6_cache, ag_6_cache, aF1_6_cache, aF2_6_cache, 
        af_7_cache, ag_7_cache, aF1_7_cache, aF2_7_cache, af_8_cache, ag_8_cache, aF1_8_cache, aF2_8_cache, 
        af_9_cache, ag_9_cache, aF1_9_cache, aF2_9_cache, af_10_cache, ag_10_cache, aF1_10_cache, aF2_10_cache,
        bf_1_cache, bg_1_cache, bF1_1_cache, bF2_1_cache, bf_2_cache, bg_2_cache, bF1_2_cache, bF2_2_cache, 
        bf_3_cache, bg_3_cache, bF1_3_cache, bF2_3_cache, bf_4_cache, bg_4_cache, bF1_4_cache, bF2_4_cache, 
        bf_5_cache, bg_5_cache, bF1_5_cache, bF2_5_cache, bf_6_cache, bg_6_cache, bF1_6_cache, bF2_6_cache, 
        bf_7_cache, bg_7_cache, bF1_7_cache, bF2_7_cache, bf_8_cache, bg_8_cache, bF1_8_cache, bF2_8_cache, 
        bf_9_cache, bg_9_cache, bF1_9_cache, bF2_9_cache, bf_10_cache, bg_10_cache, bF1_10_cache, bF2_10_cache;
    gslpp::complex CS_cache,CSp_cache,CP_cache,CPp_cache,CV_cache,CVp_cache,CA_cache,CAp_cache,CT_cache,CTp_cache;
    bool checkcache_int_tau, checkcache_int_mu, checkcache_int_el;
    
    void updateParameters();
 
    double lambda_half(double a, double b, double c);
    
    double phi_f(double z);
 
    double f_BGL(double q2);
    
    double phi_g(double q2);
 
    double g_BGL(double q2);
 
    double phi_F1(double q2);
    
    double F1_BGL(double q2);
    
    double phi_F2(double q2);
    
    double F2_BGL(double q2);
    
    double hA1(double q2);
    
    double R1(double q2);
    
    double R2(double q2);
    
    double R0(double q2);
    
    double V(double q2);
 
    double A0(double q2);
 
    
    double A1(double q2);
    
    double A2(double q2);
    
    double T1(double q2);
 
    
    double T2(double q2);
    
    double A12(double q2); 
    
    double T23(double q2);
 
    gslpp::complex HV0(double q2);
 
    gslpp::complex HVp(double q2);
    
    gslpp::complex HVm(double q2);
    
    gslpp::complex HA0(double q2);
 
    gslpp::complex HAp(double q2);
    
    gslpp::complex HAm(double q2);
    
    gslpp::complex HP(double q2);
    
    gslpp::complex HS(double q2);
    
    gslpp::complex HT0(double q2);
    
    gslpp::complex HT0t(double q2);
 
    gslpp::complex HTp(double q2);
    
    gslpp::complex HTpt(double q2);
    
    gslpp::complex HTm(double q2);
    
    gslpp::complex HTmt(double q2);  
 
    gslpp::complex G000(double q2);
 
    gslpp::complex G010(double q2); 
    
    gslpp::complex G020(double q2);
    
    gslpp::complex G200(double q2);
    
    gslpp::complex G210(double q2);
    
    gslpp::complex G220(double q2);
    
    gslpp::complex G211(double q2);
    
    gslpp::complex G221(double q2);
    
    gslpp::complex G222(double q2);
    
    double  J1s(double q2);
    
    double  J1c(double q2);
    
    double  J2s(double q2);
    
    double  J2c(double q2);
    
    double  J3(double q2);
    
    double  J4(double q2);
    
    double  J5(double q2);
    
    double  J6s(double q2);
    
    double  J6c(double q2);
    
    double  J7(double q2);
    
    double  J8(double q2);
    
    double  J9(double q2);
 
    double integrateJ(int i, double q2_min, double q2_max) ;
    
    double dGammadw(double w);
 
    double dGammadcldq2(double q2, double cl);
 
    double dGammadcl(double cl);
 
    double dGammadcVdq2(double q2, double cV);
    
    double dGammadcV(double cV);
    
    double dGammadchidq2(double q2, double chi);
    
    double dGammadchi(double chi);
    
    double  Hminus(double q2);
    
    double  Hplus(double q2);
    
    double  H0(double q2);
    
    double  H0t(double q2);
    
    double dGpdq2(double q2);
    
    double dGmdq2(double q2);
 
    double integrateGpm(int i, double q2_min, double q2_max) ;
    
    double J_res;
    ROOT::Math::GSLIntegrator ig;
    ROOT::Math::Functor1D wf;
};
 
#endif /* MVLNU_H */