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 <gsl/gsl_integration.h>
#include <gsl/gsl_spline.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_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 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; 
    double CS; 
    double CP; 
    double CSp; 
    double CPp; 
    double CV; 
    double CA; 
    double CVp; 
    double CAp; 
    double C7; 
    double C7p; 
    double CT; 
    double 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,aF21,aF22; 
    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 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,aF21_cache,aF22_cache;
    double 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) ;
    
    gsl_error_handler_t * old_handler; 
    gsl_function FJ;
    double J_res;
    double J_err;
    gsl_integration_cquad_workspace * w_J;
};
 
#endif /* MVLNU_H */