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 <TF1.h>
#include <TGraph.h>
#include <TFitResultPtr.h>
#include <gsl/gsl_spline.h>
#include "Math/GSLIntegrator.h"
#include "Math/Functor.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 BGLflag; 
    bool DMflag; 
    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; 
    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 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 afplus1, afplus2, afplus3, afplus4, afplus5, afplus6,
        afplus7, afplus8, afplus9, afplus10, afzero1, afzero2, afzero3,
        afzero4, afzero5, afzero6, afzero7, afzero8, afzero9, afzero10,
        bfplus1, bfplus2, bfplus3, bfplus4, bfplus5, bfplus6, bfplus7, 
        bfplus8, bfplus9, bfplus10, bfzero1, bfzero2, bfzero3, bfzero4, 
        bfzero5, bfzero6, bfzero7, bfzero8, bfzero9, bfzero10; 
    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;
    double afplus1_cache, afplus2_cache, afplus3_cache, afplus4_cache, afplus5_cache, afplus6_cache,
        afplus7_cache, afplus8_cache, afplus9_cache, afplus10_cache, afzero1_cache, afzero2_cache, afzero3_cache,
        afzero4_cache, afzero5_cache, afzero6_cache, afzero7_cache, afzero8_cache, afzero9_cache, afzero10_cache,
        bfplus1_cache, bfplus2_cache, bfplus3_cache, bfplus4_cache, bfplus5_cache, bfplus6_cache, bfplus7_cache, 
        bfplus8_cache, bfplus9_cache, bfplus10_cache, bfzero1_cache, bfzero2_cache, bfzero3_cache, bfzero4_cache, 
        bfzero5_cache, bfzero6_cache, bfzero7_cache, bfzero8_cache, bfzero9_cache, bfzero10_cache;
#if NBGL == 3
    double af0_3_cache,afplus_3_cache;
#endif    
    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_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);
    
    double J_res;
    ROOT::Math::GSLIntegrator ig;
    ROOT::Math::Functor1D wf;
};
 
#endif /* MPLNU_H */