QCD.h

Go to the documentation of this file.

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

#ifndef QCD_H
#define QCD_H

#include "Model.h"
#include "Meson.h"
#include "OrderScheme.h"
#define MEPS 1.e-10 // mass precision

class QCD : public Model {
public:

    enum meson {
        P_0, 
        P_P, 
        K_0, 
        K_P, 
        D_0, 
        B_D, 
        B_P, 
        B_S, 
        PHI, 
        K_star, 
        MESON_END 
    };

    enum quark {
        UP, 
        DOWN, 
        CHARM, 
        STRANGE, 
        TOP, 
        BOTTOM 
    };

    static const int NQCDvars = 186; 

    static const std::string QCDvars[NQCDvars];

    QCD();

    std::string orderToString(const orders order) const;

    // Parameters

    virtual bool Init(const std::map<std::string, double>& DPars);

    virtual bool PreUpdate();

    virtual bool Update(const std::map<std::string, double>& DPars);

    virtual bool PostUpdate();

    virtual bool CheckParameters(const std::map<std::string, double>& DPars);


    // Flags

    virtual bool setFlag(const std::string name, const bool value);

    virtual bool setFlagStr(const std::string name, const std::string value);

    virtual bool CheckFlags() const;


    // get and set methods for class members

    double getNc() const
    {
        return Nc;
    }

    void setNc(double Nc)
    {
        this->Nc = Nc;
    }

    Meson getMesons(const meson m) const
    {
        return mesons[m];
    }

    Particle getQuarks(const quark q) const
    {
        return quarks[q];
    }

    double getAlsM() const
    {
        return AlsM;
    }

    double getMAls() const
    {
        return MAls;
    }

    double getMut() const
    {
        return mut;
    }

    double getMub() const
    {
        return mub;
    }

    double getMuc() const
    {
        return muc;
    }

    double getMtpole() const
    {
        return mtpole;
    }

    double getCF() const
    {
        return CF;
    }

    BParameter getBBd() const
    {
        return BBd;
    }

    BParameter getBBs() const
    {
        return BBs;
    }

    BParameter getBD() const
    {
        return BD;
    }

    BParameter getBK() const
    {
        return BK;
    }

    BParameter getBKd1() const
    {
        return BKd1;
    }

    BParameter getBKd3() const
    {
        return BKd3;
    }

    double getReA0_Kd() const
    {
        return ReA0_Kd;
    }

    double getReA2_Kd() const
    {
        return ReA2_Kd;
    }

    double getOmega_eta_etap() const
    {
        return Omega_eta_etap;
    }

    double getBr_Kp_P0enu() const
    {
        return Br_Kp_P0enu;
    }

    double getBr_Kp_munu() const
    {
        return Br_Kp_munu;
    }

    double getBr_B_Xcenu() const
    {
        return Br_B_Xcenu;
    }

    double getDeltaP_cu() const
    {
        return DeltaP_cu;
    }

    double getIB_Kl() const
    {
        return IB_Kl;
    }

    double getIB_Kp() const
    {
        return IB_Kp;
    }
        
    gslpp::complex geth_0() const
    {
        return myh_0;
    }

    gslpp::complex geth_p() const
    {
        return myh_p;
    }
    
    gslpp::complex geth_m() const
    {
        return myh_m;
    }
    
    gslpp::complex geth_0_1() const
    {
        return myh_0_1;
    }

    gslpp::complex geth_p_1() const
    {
        return myh_p_1;
    }
    
    gslpp::complex geth_m_1() const
    {
        return myh_m_1;
    }
    
    gslpp::complex geth_0_2() const
    {
        return myh_0_2;
    }

    gslpp::complex geth_p_2() const
    {
        return myh_p_2;
    }
    
    gslpp::complex geth_m_2() const
    {
        return myh_m_2;
    }
    
    gslpp::complex geth_0_MP() const
    {
        return gslpp::complex(absh_0_MP,argh_0_MP,true);
    }
    
    gslpp::complex geth_0_1_MP() const
    {
        return gslpp::complex(absh_0_1_MP,argh_0_1_MP,true);
    }
    
    double geta_0V() const
    {
        return a_0V;
    }
    
    double geta_1V() const
    {
        return a_1V;
    }
    
    double geta_2V() const
    {
        return a_2V;
    }
    
    double getMRV() const
    {
        return MRV;
    }
    
    double geta_0A0() const
    {
        return a_0A0;
    }
    
    double geta_1A0() const
    {
        return a_1A0;
    }
    
    double geta_2A0() const
    {
        return a_2A0;
    }
    
    double getMRA0() const
    {
        return MRA0;
    }
    
    double geta_0A1() const
    {
        return a_0A1;
    }
    
    double geta_1A1() const
    {
        return a_1A1;
    }
    
    double geta_2A1() const
    {
        return a_2A1;
    }
    
    double getMRA1() const
    {
        return MRA1;
    }
    
    double geta_0A12() const
    {
        return a_0A12;
    }
    
    double geta_1A12() const
    {
        return a_1A12;
    }
    
    double geta_2A12() const
    {
        return a_2A12;
    }
    
    double getMRA12() const
    {
        return MRA12;
    }
    
    double geta_0T1() const
    {
        return a_0T1;
    }
    
    double geta_1T1() const
    {
        return a_1T1;
    }
    
    double geta_2T1() const
    {
        return a_2T1;
    }
    
    double getMRT1() const
    {
        return MRT1;
    }
    
    double geta_0T2() const
    {
        return a_0T2;
    }
    
    double geta_1T2() const
    {
        return a_1T2;
    }
    
    double geta_2T2() const
    {
        return a_2T2;
    }
    
    double getMRT2() const
    {
        return MRT2;
    }
    
    double geta_0T23() const
    {
        return a_0T23;
    }
    
    double geta_1T23() const
    {
        return a_1T23;
    }
    
    double geta_2T23() const
    {
        return a_2T23;
    }
    
    double getMRT23() const
    {
        return MRT23;
    }

    double geta_0Vphi() const
    {
        return a_0Vphi;
    }
    
    double geta_1Vphi() const
    {
        return a_1Vphi;
    }
    
    double geta_2Vphi() const
    {
        return a_2Vphi;
    }
    
    double getMRVphi() const
    {
        return MRVphi;
    }
    
    double geta_0A0phi() const
    {
        return a_0A0phi;
    }
    
    double geta_1A0phi() const
    {
        return a_1A0phi;
    }
    
    double geta_2A0phi() const
    {
        return a_2A0phi;
    }
    
    double getMRA0phi() const
    {
        return MRA0phi;
    }
    
    double geta_0A1phi() const
    {
        return a_0A1phi;
    }
    
    double geta_1A1phi() const
    {
        return a_1A1phi;
    }
    
    double geta_2A1phi() const
    {
        return a_2A1phi;
    }
    
    double getMRA1phi() const
    {
        return MRA1phi;
    }
    
    double geta_0A12phi() const
    {
        return a_0A12phi;
    }
    
    double geta_1A12phi() const
    {
        return a_1A12phi;
    }
    
    double geta_2A12phi() const
    {
        return a_2A12phi;
    }
    
    double getMRA12phi() const
    {
        return MRA12phi;
    }
    
    double geta_0T1phi() const
    {
        return a_0T1phi;
    }
    
    double geta_1T1phi() const
    {
        return a_1T1phi;
    }
    
    double geta_2T1phi() const
    {
        return a_2T1phi;
    }
    
    double getMRT1phi() const
    {
        return MRT1phi;
    }
    
    double geta_0T2phi() const
    {
        return a_0T2phi;
    }
    
    double geta_1T2phi() const
    {
        return a_1T2phi;
    }
    
    double geta_2T2phi() const
    {
        return a_2T2phi;
    }
    
    double getMRT2phi() const
    {
        return MRT2phi;
    }
    
    double geta_0T23phi() const
    {
        return a_0T23phi;
    }
    
    double geta_1T23phi() const
    {
        return a_1T23phi;
    }
    
    double geta_2T23phi() const
    {
        return a_2T23phi;
    }
    
    double getMRT23phi() const
    {
        return MRT23phi;
    }
    
    double getr_1_fplus() const
    {
        return r_1_fplus;
    }

    double getr_2_fplus() const
    {
        return r_2_fplus;
    }
    
    double getm_fit2_fplus() const
    {
        return m_fit2_fplus;
    }
    
    double getr_1_fT() const
    {
        return r_1_fT;
    }

    double getr_2_fT() const
    {
        return r_2_fT;
    }
    
    double getm_fit2_fT() const
    {
        return m_fit2_fT;
    }
    
    double getr_2_f0() const
    {
        return r_2_f0;
    }
    
    double getm_fit2_f0() const
    {
        return m_fit2_f0;
    }
    
    double getbsgamma_E0() const
    {
        return bsgamma_E0;
    }
    
    double getBLNPcorr() const
    {
        return BLNPcorr;
    }
    
    double getGambino_mukin() const
    {
        return Gambino_mukin;
    }
    
    double getGambino_BRsem() const
    {
        return Gambino_BRsem;
    }
    
    double getGambino_Mbkin() const
    {
        return Gambino_Mbkin;
    }
    
    double getGambino_Mcatmuc() const
    {
        return Gambino_Mcatmuc;
    }
    
    double getGambino_mupi2() const
    {
        return Gambino_mupi2;
    }
    
    double getGambino_rhoD3() const
    {
        return Gambino_rhoD3;
    }
    
    double getGambino_muG2() const
    {
        return Gambino_muG2;
    }
    
    double getGambino_rhoLS3() const
    {
        return Gambino_rhoLS3;
    }
    
    double getFKstarp() const
    {
        return FKstarp;
    }

    

    double Thresholds(const int i) const;

    double AboveTh(const double mu) const;

    double BelowTh(const double mu) const;

    double Nf(const double mu) const;


    double Beta0(const double nf) const;

    double Beta1(const double nf) const;

    double Beta2(const double nf) const;

    double AlsWithInit(const double mu, const double alsi, const double mu_i,
            const orders order) const;

    double AlsWithLambda(const double mu, const orders order) const;

    double Als(const double mu, const orders order = FULLNLO) const;

    double logLambda(const double nf, orders order) const;

    double Als4(const double mu) const;
    
    virtual double alphaMz() const = 0;

    double Alstilde5(const double mu) const;

    double Mrun4(const double mu_f, const double mu_i, const double m) const;


    double Gamma0(const double nf) const;

    double Gamma1(const double nf) const;

    double Gamma2(const double nf) const;

    double Mrun(const double mu, const double m, const orders order = FULLNLO) const;

    double Mrun(const double mu_f, const double mu_i, const double m,
            const orders order = FULLNLO) const;


    double Mbar2Mp(const double mbar, const orders order = FULLNLO) const;

    double Mp2Mbar(const double mp, const orders order = FULLNLO) const;

    double MS2DRqmass(const double MSscale, const double MSbar) const;

    double MS2DRqmass(const double MSbar) const;


protected:

    virtual void setParameter(const std::string name, const double& value);

    bool computemt; 
    bool requireYu; 
    bool requireYd; 

    // model parameters
    double AlsM; 
    double MAls; 
    double mtpole; 
    double mut; 
    double mub; 
    double muc; 
    double ReA0_Kd;
    double ReA2_Kd;
    double Omega_eta_etap;
    double Br_Kp_P0enu;
    double IB_Kl;
    double IB_Kp;
    double DeltaP_cu;
    double Br_Kp_munu;
    double Br_B_Xcenu;
    double BBsoBBd; 
    double FBsoFBd; 
    double absh_0, absh_p, absh_m, argh_0, argh_p, argh_m;
    double absh_0_1, absh_p_1, absh_m_1, argh_0_1, argh_p_1, argh_m_1;
    double absh_0_2, absh_p_2, absh_m_2, argh_0_2, argh_p_2, argh_m_2;
    double absh_0_MP, argh_0_MP, absh_0_1_MP, argh_0_1_MP;
    double a_0V, a_1V, a_2V, MRV;
    double a_0A0, a_1A0, a_2A0, MRA0;
    double a_0A1, a_1A1, a_2A1, MRA1;
    double a_0A12, a_1A12, a_2A12, MRA12;
    double a_0T1, a_1T1, a_2T1, MRT1;
    double a_0T2, a_1T2, a_2T2, MRT2;
    double a_0T23, a_1T23, a_2T23, MRT23;
    double a_0Vphi, a_1Vphi, a_2Vphi, MRVphi;
    double a_0A0phi, a_1A0phi, a_2A0phi, MRA0phi;
    double a_0A1phi, a_1A1phi, a_2A1phi, MRA1phi;
    double a_0A12phi, a_1A12phi, a_2A12phi, MRA12phi;
    double a_0T1phi, a_1T1phi, a_2T1phi, MRT1phi;
    double a_0T2phi, a_1T2phi, a_2T2phi, MRT2phi;
    double a_0T23phi, a_1T23phi, a_2T23phi, MRT23phi;
    double r_1_fplus, r_2_fplus, m_fit2_fplus;
    double r_1_fT, r_2_fT, m_R_fT, m_fit2_fT;
    double r_2_f0, m_fit2_f0;
    double bsgamma_E0, BLNPcorr;
    double Gambino_mukin, Gambino_BRsem, Gambino_Mbkin, Gambino_Mcatmuc, Gambino_mupi2, Gambino_rhoD3, Gambino_muG2, Gambino_rhoLS3;
    double FKstarp; //matrix element of tensor current for transverse polarization at 1 GeV
    
    //double r_2A0, r_2T1, r_2T2, r_2A0phi, r_2T1phi, r_2T2phi removed because they are fixed by form factors relations
    double Nc; 
    Particle quarks[6]; 
    Meson mesons[MESON_END]; 
private:
    
    double CF; 
    BParameter BBs; 
    BParameter BBd; 
    BParameter BD; 
    BParameter BK; 
    BParameter BKd1;
    BParameter BKd3;
    double BBsB0, BBsB1, BBsB2, BBsB3, BBsB4, BBsscale;
    //double BBdB0, BBdB1, BBdB2, BBdB3, BBdB4, BBdscale;
    double BDB0, BDB1, BDB2, BDB3, BDB4, BDscale;
    double BKB0, BKB1, BKB2, BKB3, BKB4, BKscale;
    double BKd1B0, BKd1B1, BKd1B2, BKd1B3, BKd1B4, BKd1B5, BKd1B6, BKd1B7, BKd1B8, BKd1B9, BKd_scale;
    double BKd3B0, BKd3B1, BKd3B2, BKd3B3, BKd3B4, BKd3B5, BKd3B6, BKd3B7, BKd3B8, BKd3B9;

    double zeta2; 
    double zeta3; 
    bool computeFBd; 
    bool computeBd; 
    static const int CacheSize = 5; 
    mutable double als_cache[8][CacheSize]; 
    mutable double logLambda5_cache[4][CacheSize];
    mutable double logLambdaNLO_cache[9][CacheSize];
    mutable double mrun_cache[10][CacheSize]; 
    mutable double mp2mbar_cache[5][CacheSize]; 
    bool unknownParameterWarning; 

    double AlsWithLambda(const double mu, const double logLambda, const orders order) const;

    double ZeroNf6NLO(double *logLambda6, double *logLambda5_in) const;

    double ZeroNf5(double *logLambda5, double *order) const;

    double ZeroNf4NLO(double *logLambda4, double *logLambda5_in) const;

    double ZeroNf3NLO(double *logLambda3, double *logLambda4_in) const;

    double logLambda5(orders order) const;

    double logLambdaNLO(const double nfNEW, const double nfORG, const double logLambdaORG) const;

    double logLambda(const double muMatching, const double mf,
            const double nfNEW, const double nfORG,
            const double logLambdaORG, orders order) const;

    double threCorrForMass(const double nf_f, const double nf_i) const;

    double MrunTMP(const double mu_f, const double mu_i, const double m, const orders order) const;

    double Mp2MbarTMP(double *mu, double *params) const;

    void CacheShift(double cache[][5], int n) const;

    gslpp::complex myh_0,myh_p,myh_m,myh_0_1,myh_p_1,myh_m_1,myh_0_2,myh_p_2,myh_m_2;
};

#endif  /* QCD_H */