Flavour.h

Go to the documentation of this file.

/*
 * Copyright (C) 2012 HEPfit Collaboration
 *
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef FLAVOUR_H
#define FLAVOUR_H
 
class StandardModel;
class HeffDF2;
class HeffDF1_diujlknu;
class HeffDB1;
class HeffDC1;
class HeffDS1;
class MVll;
class MPll;
class MVgamma;
class MVlnu;
class MPlnu;
class AmpDB2;
#include "QCD.h"
#include <boost/tuple/tuple.hpp>
#include <memory>
 
class Flavour {
public:
 
    Flavour(const StandardModel& SM_i);
 
    HeffDF2& getHDF2() const;
 
 
    HeffDF1_diujlknu& getHDF1_diujlknu() const;
 
    HeffDS1& getHDS1() const;
 
 
    HeffDC1& getHDC1() const;
 
    HeffDB1& getHDB1() const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffBd(double mu, schemes scheme = NDR, bool SM = false) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffDS1pnunu() const;
    gslpp::vector<gslpp::complex>** ComputeCoeffDS1pnunuC() const;
    gslpp::vector<gslpp::complex>** ComputeCoeffDS1mumu() const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffBs(double mu, schemes scheme = NDR, bool SM = false) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffdd(double mu, schemes scheme = NDR) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffK(double mu, schemes scheme = NDR) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffmK(double mu, schemes scheme = NDR) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffDS1PPv(double mu, schemes scheme = NDR) const;
    gslpp::vector<gslpp::complex>** ComputeCoeffDS1PPz(double muc, schemes scheme = NDR) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffsmumu(double mu, schemes scheme = NDR) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffdmumu(double mu, schemes scheme = NDR) const;
 
    
    gslpp::vector<gslpp::complex>** ComputeCoeffdiujlknu(int i, int j, int k, double mu) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffsnunu(QCD::lepton lepton = QCD::NOLEPTON, bool noSM = false) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffdnunu() const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffsgamma(double mu, bool noSM = false, schemes scheme = NDR) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffprimesgamma(double mu, schemes scheme = NDR) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffsgamma_Buras(double mu, bool noSM = false, schemes scheme = NDR) const;
    
    gslpp::vector<gslpp::complex>** ComputeCoeffBMll(double mu, QCD::lepton lepton, bool noSM = false, schemes scheme = NDR) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffBMll_Buras(double mu, QCD::lepton lepton, bool noSM = false, schemes scheme = NDR) const;
 
    gslpp::vector<gslpp::complex>** ComputeCoeffprimeBMll(double mu, QCD::lepton lepton, schemes scheme = NDR) const;
 
    AmpDB2& getDB2(int BMeson_i, bool flag_fixmub = false, bool flag_RI = false) const;
 
    MVll& getMVll(QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i) const;
 
    MPll& getMPll(QCD::meson meson_i, QCD::meson pseudoscalar_i, QCD::lepton lep_i) const;
 
    MVgamma& getMVgamma(QCD::meson meson_i, QCD::meson vector_i) const;
 
    MVlnu& getMVlnu(QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i) const;
 
    MPlnu& getMPlnu(QCD::meson meson_i, QCD::meson pseudoscalar_i, QCD::lepton lep_i) const;
 
    void setUpdateFlag(QCD::meson meson_i, QCD::meson meson_j, QCD::lepton lep_i, bool updated_i) const;
 
    bool getUpdateFlag(QCD::meson meson_i, QCD::meson meson_j, QCD::lepton lep_i) const;
 
    void setSMupdated() const;
 
    bool setFlag(const std::string name, const bool value);
 
    bool setFlagUseDispersionRelation(bool dispersion) {
        return (this->dispersion = dispersion);
    }
 
    bool setFlagUsezExpansion(bool zExpansion) {
        return (this->zExpansion = zExpansion);
    }
 
    bool setFlagCLN(bool CLNflag) {
        return (this->CLNflag = CLNflag);
    }
    
    bool setFlagBGL(bool BGLflag) {
        return (this->BGLflag = BGLflag);
    }
    
    bool setFlagDM(bool DMflag) {
        return (this->DMflag = DMflag);
    }
 
    bool setFlagFixedWCbtos(bool FixedWCbtosflag) {
        return (this->FixedWCbtosflag = FixedWCbtosflag);
    }
 
    bool setFlagMPll_FNALMILC(bool MPll_FNALMILC_flag) {
        return (this->MPll_FNALMILC_flag = MPll_FNALMILC_flag);
    }
 
    bool setFlagMPll_GRvDV(bool MPll_GRvDV_flag) {
        return (this->MPll_GRvDV_flag = MPll_GRvDV_flag);
    }
 
    bool setFlagMPll_DM(bool MPll_DM_flag) {
        return (this->MPll_DM_flag = MPll_DM_flag);
    }
 
    bool setFlagMVll_DM(bool MVll_DM_flag) {
        return (this->MVll_DM_flag = MVll_DM_flag);
    }
 
    bool setFlagNeutrinoTree(bool NeutrinoTree_flag) const {
        return (this->NeutrinoTree_flag = NeutrinoTree_flag);
    }
 
    bool setFlagBXsnunu_LFUNP(bool BXsnunu_LFUNP_flag) const {
        return (this->BXsnunu_LFUNP_flag = BXsnunu_LFUNP_flag);
    }
 
    bool getFlagUseDispersionRelation() const {
        return dispersion;
    }
 
    bool getFlagUsezExpansion() const {
        return zExpansion;
    }
 
    bool getFlagCLN() const {
        return CLNflag;
    }
    
    bool getFlagBGL() const {
        return BGLflag;
    }
    
    bool getFlagDM() const {
        return DMflag;
    }
 
    bool getFlagFixedWCbtos() const {
        return FixedWCbtosflag;
    }
 
    bool getFlagMPll_FNALMILC() const {
        return MPll_FNALMILC_flag;
    }
 
    bool getFlagMPll_GRvDV() const {
        return MPll_GRvDV_flag;
    }
 
    bool getFlagMPll_DM() const {
        return MPll_DM_flag;
    }
 
    bool getFlagMVll_DM() const {
        return MVll_DM_flag;
    }
 
    bool getFlagNeutrinoTree() const {
        return NeutrinoTree_flag;
    }
 
    bool getFlagBXsnunu_LFUNP() const {
        return BXsnunu_LFUNP_flag;
    }
 
private:
    template<typename T, typename... Args> std::shared_ptr<T>& getPtr(std::shared_ptr<T>& x, Args... args) const;
    template <typename T, typename... Args> T& getM(std::map<std::vector<int>, std::shared_ptr<T> >& map, Args ... args) const;
    const StandardModel & mySM;
    mutable std::shared_ptr<HeffDF2> HDF2; 
    mutable std::shared_ptr<HeffDF1_diujlknu> HDF1_diujlknu; 
    mutable std::shared_ptr<HeffDB1> HDB1; 
    mutable std::shared_ptr<HeffDS1> HDS1; 
    mutable std::map<std::vector<int>, std::shared_ptr<MVll> > MVllMap;
    mutable std::map<std::vector<int>, std::shared_ptr<MVlnu> > MVlnuMap;
    mutable std::map<std::vector<int>, std::shared_ptr<MVgamma> > MVgammaMap;
    mutable std::map<std::vector<int>, std::shared_ptr<MPll> > MPllMap;
    mutable std::map<std::vector<int>, std::shared_ptr<MPlnu> > MPlnuMap;
    mutable std::map<std::vector<int>, std::shared_ptr<AmpDB2> > AmpDB2Map;
    mutable std::map<std::vector<int>, bool> flagUpdateMap;
 
    mutable bool dispersion;
    mutable bool zExpansion;
    mutable bool CLNflag;
    mutable bool BGLflag;
    mutable bool DMflag;
    mutable bool FixedWCbtosflag;
    mutable bool MPll_FNALMILC_flag;
    mutable bool MPll_GRvDV_flag;
    mutable bool MPll_DM_flag;
    mutable bool MVll_DM_flag;
    mutable bool NeutrinoTree_flag;
    mutable bool BXsnunu_LFUNP_flag;
};
 
#endif /* FLAVOUR_H */