MPll.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2014 HEPfit Collaboration
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */

#ifndef MPLL_H
#define MPLL_H

#include <math.h>
#include <StandardModel.h>
#include <ThObservable.h>
#include <gsl/gsl_integration.h>
#include <assert.h>


#define CUTOFF 10    //cutoff between LCSR and lattice values for Form Factors, in GeV^2

class MPll{
public:
    MPll(const StandardModel& SM_i, StandardModel::meson meson_i, StandardModel::meson pseudoscalar_i, StandardModel::lepton lep_i);
    
    virtual ~MPll();
    
    void updateParameters();
    
    void checkCache();
    
    double GF;            
    double ale;           
    double Mlep;          
    double MM;            
    double MP;            
    double Mb;            
    double mu_b;          
    double Mc;            
    double Ms;            
    double width;         
    double MW;            
    gslpp::complex lambda_t;     
    double b;             
    gslpp::complex h_0;          
    gslpp::complex h_0_1;        
    double q2;            
    /*LCSR fit parameters*/
    double r_1_fplus;
    double r_2_fplus;
    double m_fit2_fplus;
    double r_1_fT;
    double r_2_fT;
    double m_fit2_fT;
    double r_2_f0;
    double m_fit2_f0;
    gslpp::vector<gslpp::complex> ** allcoeff;
    gslpp::vector<gslpp::complex> ** allcoeffprime;
    gslpp::complex C_1;
    gslpp::complex C_2;
    gslpp::complex C_3;
    gslpp::complex C_4;
    gslpp::complex C_5;
    gslpp::complex C_6;
    gslpp::complex C_7;
    gslpp::complex C_9;
    gslpp::complex C_10;
    gslpp::complex C_S;
    gslpp::complex C_P;
    gslpp::complex C_7p;
    gslpp::complex C_9p;
    gslpp::complex C_10p;
    gslpp::complex C_Sp;
    gslpp::complex C_Pp;
    double LCSR_fit1(double q2, double r_1, double r_2, double m_fit2);
    
    
    double LCSR_fit2(double q2, double r_2, double m_fit2);
    
    
    double f_plus(double q2);
    
    
    double f_T(double q2);
    
    
    double f_0(double q2);
    
    gslpp::complex V_L(double q2);

    
    gslpp::complex V_R(double q2);


    gslpp::complex T_L(double q2);


    gslpp::complex T_R(double q2);


    double S_L(double q2);


    double S_R(double q2);


    gslpp::complex N();
    
    
    gslpp::complex H(double q2, double m);
    
    
    gslpp::complex Y(double q2);
    
    
    gslpp::complex H_V(double q2, int bar);


    gslpp::complex H_A(double q2, int bar);


    gslpp::complex H_S(double q2, int bar);


    gslpp::complex H_P(double q2, int bar);
    
    
    double k2 (double q2);
    
    
    double beta (double q2);
    
    
    double lambda(double q2);

    
    double F(double q2);
    
    
    double  I(int i, double q2, int bar);
    
    
    double Sigma(int i, double q2);
    
    
    double Delta(int i, double q2);
    
    double getSigma0(double q2)
    {
        return Sigma(0, q2);
    };
    
    double getSigma2(double q2)
    {
        return Sigma(2, q2);
    };
    
    double getDelta0(double q2)
    {
        return Delta(0, q2);
    };
    
    double getDelta2(double q2)
    {
        return Delta(2, q2);
    };
    
    double integrateSigma(int i, double q_min, double q_max);
    
    double integrateDelta(int i, double q_min, double q_max);
    

private:
    const StandardModel& mySM;
    StandardModel::lepton lep;
    StandardModel::meson meson;
    StandardModel::meson pseudoscalar;
    unsigned int fplus_updated;
    gslpp::vector<double> fplus_cache;
    unsigned int fT_updated;
    gslpp::vector<double> fT_cache;
    unsigned int f0_updated;
    double f0_cache;
    unsigned int k2_updated;
    gslpp::vector<double> k2_cache;
    unsigned int beta_updated;
    double beta_cache;
    unsigned int lambda_updated;
    double lambda_cache;
    unsigned int F_updated;
    unsigned int VL_updated;
    unsigned int VR_updated;
    unsigned int TL_updated;
    unsigned int TR_updated;
    unsigned int SL_updated;
    gslpp::vector<double> SL_cache;
    unsigned int SR_updated;
    unsigned int N_updated;
    gslpp::vector<double> N_cache;
    gslpp::complex Nc_cache;
    unsigned int C_1_updated;
    gslpp::complex C_1_cache;
    unsigned int C_2_updated;
    gslpp::complex C_2_cache;
    unsigned int C_3_updated;
    gslpp::complex C_3_cache;
    unsigned int C_4_updated;
    gslpp::complex C_4_cache;
    unsigned int C_5_updated;
    gslpp::complex C_5_cache;
    unsigned int C_6_updated;
    gslpp::complex C_6_cache;
    unsigned int C_7_updated;
    gslpp::complex C_7_cache;
    unsigned int C_9_updated;
    gslpp::complex C_9_cache;
    unsigned int C_10_updated;
    gslpp::complex C_10_cache;
    unsigned int C_7p_updated;
    gslpp::complex C_7p_cache;
    unsigned int C_9p_updated;
    gslpp::complex C_9p_cache;
    unsigned int C_10p_updated;
    gslpp::complex C_10p_cache;
    unsigned int C_S_updated;
    gslpp::complex C_S_cache;
    unsigned int C_P_updated;
    gslpp::complex C_P_cache;
    unsigned int C_Sp_updated;
    gslpp::complex C_Sp_cache;
    unsigned int C_Pp_updated;
    gslpp::complex C_Pp_cache;
    unsigned int Yupdated;
    gslpp::vector<double> Ycache;
    unsigned int H_V0updated;
    gslpp::vector<double> H_V0cache;
    gslpp::complex H_V0Ccache[2];
    unsigned int H_A0updated;
    unsigned int H_Supdated;
    gslpp::vector<double> H_Scache;
    unsigned int H_P_updated;
    gslpp::vector<double> H_P_cache;
    unsigned int I0_updated;
    unsigned int I2_updated;
    unsigned int I8_updated;
    std::map<std::pair<double, double>, unsigned int > sigma0Cached;
    std::map<std::pair<double, double>, unsigned int > sigma2Cached;
    std::map<std::pair<double, double>, unsigned int > delta0Cached;
    std::map<std::pair<double, double>, unsigned int > delta2Cached;
    double avaSigma0;
    double avaSigma2;
    double errSigma0;
    double errSigma2;
    double avaDelta0;
    double avaDelta2;
    double errDelta0;
    double errDelta2;
    gsl_function FS0;
    gsl_function FS2;
    gsl_function FD0;
    gsl_function FD2;
    gsl_integration_workspace * w_sigma0;
    gsl_integration_workspace * w_sigma2;
    gsl_integration_workspace * w_delta0;
    gsl_integration_workspace * w_delta2;
    std::map<std::pair<double, double>, double > cacheSigma0;
    std::map<std::pair<double, double>, double > cacheSigma2;
    std::map<std::pair<double, double>, double > cacheDelta0;
    std::map<std::pair<double, double>, double > cacheDelta2;
};

#endif  /* MPLL_H */