MPlnu.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2014 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef MPLNU_H
#define MPLNU_H
 
class StandardModel;
#include <gsl/gsl_integration.h>
#include <TF1.h>
#include <TGraph.h>
#include <TFitResultPtr.h>
#include <gsl/gsl_spline.h>
#include <memory>
 
#define NBGL 3 /* ONLY 3 or 2*/
 
class MPlnu {
public:
    MPlnu(const StandardModel& SM_i, QCD::meson meson_i, QCD::meson pseudoscalar_i, QCD::lepton lep_i);
    
    virtual ~MPlnu();
    std::vector<std::string> initializeMPlnuParameters();
    
    double getDeltaGammaDeltaw(double w_min, double w_max);
 
    double get_unitarity_1min_BGL();
 
    double get_unitarity_0plus_BGL();
    
    double get_strong_unitarity_BGL();
 
    double get_fplus(double q2);
    
    double get_f0(double q2);
    
    double get_fT(double q2);
 
    double getMwidth()
    {
        updateParameters();
        return width;
    }
    
    double getaf0_0()
    {
        updateParameters();
        return af0_0;
    }
    
private:
    const StandardModel& mySM;
    QCD::lepton lep;
    QCD::meson meson;
    QCD::meson pseudoscalarM;
    std::vector<std::string> mplnuParameters;
    bool CLNflag; 
    bool btocNPpmflag; 
    bool NPanalysis; 
    double GF;            
    double Mlep;          
    double Mnu;           
    double MM;            
    double MP;            
    double w0;            
    double z0;            
    double RV;            
    double mu_b;          
    double Mb;            
    double Mc;            
    double width;         
    double ale_mub;   
    gslpp::complex Vcb;   
    double amplsq_factor;   
    double q2min, q2max; 
    double eta_EW; 
    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 fplusz0,rho1to2; 
    double N_0, alpha_0, alpha_p, beta_0, beta_p, gamma_0, gamma_p; 
    double af0_0,af0_1,af0_2,afplus_0,afplus_1,afplus_2; 
#if NBGL == 3
    double af0_3, afplus_3;
#endif    
    double mBc1m_1,mBc1m_2,mBc1m_3,mBc1m_4,mBc0p_1,mBc0p_2,chitildeT,chiL,nI; 
    double z1m_1,z1m_2,z1m_3,z0p_1,z0p_2;
    gslpp::complex z1m_4;
    double cached_intJ1_tau, cached_intJ2_tau, cached_intJ3_tau,
                cached_intJ1_mu, cached_intJ2_mu, cached_intJ3_mu,
                cached_intJ1_el, cached_intJ2_el, cached_intJ3_el; 
    double fplusz0_cache,rho1to2_cache;
    double N_0_cache, alpha_0_cache, alpha_p_cache, beta_0_cache, beta_p_cache, gamma_0_cache, gamma_p_cache;
    double af0_1_cache,af0_2_cache,afplus_0_cache,afplus_1_cache,afplus_2_cache;
#if NBGL == 3
    double af0_3_cache,afplus_3_cache;
#endif    
    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_fplus(double q2);
 
    double fplus(double q2);
    
    double phi_f0(double z);
 
    double f0(double q2);
    
    double fT(double q2);
 
    gslpp::complex HV(double q2);
    
    gslpp::complex HA(double q2);
    
    gslpp::complex HP(double q2);
    
    gslpp::complex HS(double q2);
    
    gslpp::complex HT(double q2);
    
    gslpp::complex HTt(double q2);  
 
    gslpp::complex G0(double q2);
 
    gslpp::complex G1(double q2); 
    
    gslpp::complex G2(double q2);
    
    double  J1(double q2);
    
    double  J2(double q2);
    
    double  J3(double q2);
    
    double  dGammadq2(double q2);
 
    double integrateJ(int i, double q2_min, double q2_max) ;
    
    double dGammadw(double w);
    
    gsl_error_handler_t * old_handler; 
    gsl_function FJ;
    double J_res;
    double J_err;
    gsl_integration_cquad_workspace * w_J;
};
 
#endif /* MPLNU_H */