The parent class in Flavour for calculating all the Wilson coefficients for various Flavor Violating processes. More...

#include <Flavour.h>

Detailed Description

The parent class in Flavour for calculating all the Wilson coefficients for various Flavor Violating processes.

Author: HEPfit Collaboration

Copyright: GNU General Public License

The Flavour class aggregates the Wilson coefficients for each of the processes generated by the model by calling the Hamiltonians.

Definition at line 33 of file Flavour.h.

Public Member Functions
gslpp::vector< gslpp::complex > **	ComputeCoeffBd (double mu, schemes scheme=NDR) const
	Computes the Wilson coefficient for the process \( B_d \to \mu \mu \). More...

gslpp::vector< gslpp::complex > **	ComputeCoeffBMll (double mu, QCD::lepton lepton, bool noSM=false, schemes scheme=NDR) const
	Computes the Wilson coefficient for the process \( B \to V/P \ell^+ \ell^- \). More...

gslpp::vector< gslpp::complex > **	ComputeCoeffBs (double mu, schemes scheme=NDR, bool SM=false) const
	Computes the Wilson coefficient for the process \( B_s \to \mu \mu \). More...

gslpp::vector< gslpp::complex > **	ComputeCoeffbtaunu (QCD::meson meson_i) const
	Computes the Wilson coefficient for the process \( B_d \to \tau \nu \). More...

gslpp::vector< gslpp::complex > **	ComputeCoeffdd (double mu, schemes scheme=NDR) const

gslpp::vector< gslpp::complex > **	ComputeCoeffdmumu (double mu, schemes scheme=NDR) const
	Computes the Wilson coefficient for the process \( B_d \to \mu \mu \). More...

gslpp::vector< gslpp::complex > **	ComputeCoeffdnunu () const

gslpp::vector< gslpp::complex > **	ComputeCoeffDS1mumu () const

gslpp::vector< gslpp::complex > **	ComputeCoeffDS1pnunu () const

gslpp::vector< gslpp::complex > **	ComputeCoeffDS1PP (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 > **	ComputeCoeffprimeBMll (double mu, QCD::lepton lepton, schemes scheme=NDR) const
	Computes the chirality flipped Wilson coefficient for the process \( B \to V/P \ell^+ \ell^- \). More...

gslpp::vector< gslpp::complex > **	ComputeCoeffprimesgamma (double mu, schemes scheme=NDR) const
	Computes the chirality flipped Wilson coefficient for the process \( b \to s \gamma \). More...

gslpp::vector< gslpp::complex > **	ComputeCoeffsgamma (double mu, bool noSM=false, schemes scheme=NDR) const
	Computes the Wilson coefficient for the process \( b \to s \gamma \). More...

gslpp::vector< gslpp::complex > **	ComputeCoeffsmumu (double mu, schemes scheme=NDR) const
	Computes the Wilson coefficient for the process \( B_s \to \mu \mu \). More...

gslpp::vector< gslpp::complex > **	ComputeCoeffsnunu () const

	Flavour (const StandardModel &SM_i)
	The constructor. More...

bool	getFlagCLN () const

bool	getFlagFixedWCbtos () const

bool	getFlagUseDispersionRelation () const

HeffDB1 &	getHDB1 () const
	The member that returns an object of the class HeffDB1. More...

HeffDF2 &	getHDF2 () const
	The member that returns an object of the class HeffDF2. More...

HeffDS1 &	getHDS1 () const
	The member that returns an object of the class HeffDS1. More...

MPll &	getMPll (QCD::meson meson_i, QCD::meson pseudoscalar_i, QCD::lepton lep_i) const
	Returns the initial and final state dependent object for \( B \to P \ell^+ \ell^- \). More...

MPlnu &	getMPlnu (QCD::meson meson_i, QCD::meson pseudoscalar_i, QCD::lepton lep_i) const
	Returns the initial and final state dependent object for \( B \to P \ell \nu \). More...

MVgamma &	getMVgamma (QCD::meson meson_i, QCD::meson vector_i) const
	Returns the initial and final state dependent object for \( B \to V \gamma \). More...

MVll &	getMVll (QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i) const
	Returns the initial and final state dependent object for \( B \to V \ell^+ \ell^- \). More...

MVlnu &	getMVlnu (QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i) const
	Returns the initial and final state dependent object for \( B \to V \ell \nu \). More...

bool	getUpdateFlag (QCD::meson meson_i, QCD::meson meson_j, QCD::lepton lep_i) const
	gets the update flag for the initial and final state dependent object for \( B \to V \ell^+ \ell^- \). More...

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

bool	setFlagCLN (bool CLNflag)

bool	setFlagFixedWCbtos (bool FixedWCbtosflag)

bool	setFlagUseDispersionRelation (bool dispersion)

void	setSMupdated () const
	a member used for the caching for \( B \to V \ell^+ \ell^- \). More...

void	setUpdateFlag (QCD::meson meson_i, QCD::meson meson_j, QCD::lepton lep_i, bool updated_i) const
	sets the update flag for the initial and final state dependent object for \( B \to V \ell^+ \ell^- \). More...

Private Member Functions
template<typename T , typename... Args>
T &	getM (std::map< std::vector< int >, std::shared_ptr< T > > &map, Args ... args) const

template<typename T , typename... Args>
std::shared_ptr< T > &	getPtr (std::shared_ptr< T > &x, Args ... args) const

template<typename T , typename... Args>
std::shared_ptr< T > &	getPtr (std::shared_ptr< T > &x, Args... args) const

Private Attributes
bool	CLNflag

bool	dispersion

bool	FixedWCbtosflag

std::map< std::vector< int >, bool >	flagUpdateMap

std::shared_ptr< HeffDB1 >	HDB1
	An Object for the Hamiltonian of the \( \Delta B = 1 \) processes. More...

std::shared_ptr< HeffDF2 >	HDF2
	An Object for the Hamiltonian of the \( \Delta F = 2 \) processes. More...

std::shared_ptr< HeffDS1 >	HDS1
	An Object for the Hamiltonian of the \( \Delta S = 1 \) processes. More...

std::map< std::vector< int >, std::shared_ptr< MPll > >	MPllMap

std::map< std::vector< int >, std::shared_ptr< MPlnu > >	MPlnuMap

std::map< std::vector< int >, std::shared_ptr< MVgamma > >	MVgammaMap

std::map< std::vector< int >, std::shared_ptr< MVll > >	MVllMap

std::map< std::vector< int >, std::shared_ptr< MVlnu > >	MVlnuMap

const StandardModel &	mySM

Constructor & Destructor Documentation

◆ Flavour()

Flavour::Flavour ( const StandardModel & SM_i )

The constructor.

Parameters

[in] SM_i a reference to an object of the class StandardModel

Definition at line 20 of file Flavour.cpp.

 : mySM(SM_i)
 {
     dispersion = false;
     CLNflag = false;
 };

Member Function Documentation

◆ ComputeCoeffBd()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffBd	(	double	mu,
		schemes	scheme = `NDR`
	)		const

Computes the Wilson coefficient for the process \( B_d \to \mu \mu \).

Parameters

[in]	mu	the lower matching scale for the process
[in]	scheme	the scheme in which the Wilson Coefficients need to be calculated

Returns: returns the Wilson coefficients for the process \( B_d \to \mu \mu \)

Definition at line 59 of file Flavour.cpp.

 {
     return getPtr<HeffDF2>(HDF2)->ComputeCoeffBd(mu, scheme);
 }

◆ ComputeCoeffBMll()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffBMll	(	double	mu,
		QCD::lepton	lepton,
		bool	noSM = `false`,
		schemes	scheme = `NDR`
	)		const

Computes the Wilson coefficient for the process \( B \to V/P \ell^+ \ell^- \).

Parameters

[in]	mu	the lower matching scale for the process
[in]	scheme	the scheme in which the Wilson Coefficients need to be calculated

Returns: returns the Wilson coefficients for the process \( B \to V/P \ell^+ \ell^- \)

Definition at line 134 of file Flavour.cpp.

 {
     return getPtr<HeffDB1>(HDB1)->ComputeCoeffBMll(mu, lepton, noSM, scheme);
 }

◆ ComputeCoeffBs()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffBs	(	double	mu,
		schemes	scheme = `NDR`,
		bool	SM = `false`
	)		const

Computes the Wilson coefficient for the process \( B_s \to \mu \mu \).

Parameters

[in]	mu	the lower matching scale for the process
[in]	scheme	the scheme in which the Wilson Coefficients need to be calculated

Returns: returns the Wilson coefficients for the process \( B_s \to \mu \mu \)

Definition at line 64 of file Flavour.cpp.

 {
     return getPtr<HeffDF2>(HDF2)->ComputeCoeffBs(mu, scheme, SM);
 }

◆ ComputeCoeffbtaunu()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffbtaunu ( QCD::meson meson_i ) const

Computes the Wilson coefficient for the process \( B_d \to \tau \nu \).

Parameters

[in]	mu	the lower matching scale for the process
[in]	scheme	the scheme in which the Wilson Coefficients need to be calculated

Returns: returns the Wilson coefficients for the process \( B_d \to \tau \nu \)

Definition at line 109 of file Flavour.cpp.

 {
     return getPtr<HeffDB1>(HDB1)->ComputeCoeffbtaunu(meson_i);
 }

◆ ComputeCoeffdd()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffdd	(	double	mu,
		schemes	scheme = `NDR`
	)		const

Definition at line 69 of file Flavour.cpp.

 {
     return getPtr<HeffDF2>(HDF2)->ComputeCoeffdd(mu, scheme);
 }

◆ ComputeCoeffdmumu()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffdmumu	(	double	mu,
		schemes	scheme = `NDR`
	)		const

Computes the Wilson coefficient for the process \( B_d \to \mu \mu \).

Parameters

[in]	mu	the lower matching scale for the process
[in]	scheme	the scheme in which the Wilson Coefficients need to be calculated

Returns: returns the Wilson coefficients for the process \( B_d \to \mu \mu \)

Definition at line 104 of file Flavour.cpp.

 {
     return getPtr<HeffDB1>(HDB1)->ComputeCoeffdmumu(mu, scheme);
 }

◆ ComputeCoeffdnunu()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffdnunu ( ) const

Definition at line 119 of file Flavour.cpp.

 {
     return getPtr<HeffDB1>(HDB1)->ComputeCoeffdnunu();
 }

◆ ComputeCoeffDS1mumu()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffDS1mumu ( ) const

Definition at line 94 of file Flavour.cpp.

 {
     return getPtr<HeffDS1>(HDS1)->ComputeCoeffDS1mumu();
 }

◆ ComputeCoeffDS1pnunu()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffDS1pnunu ( ) const

Definition at line 89 of file Flavour.cpp.

 {
     return getPtr<HeffDS1>(HDS1)->ComputeCoeffDS1pnunu();
 }

◆ ComputeCoeffDS1PP()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffDS1PP	(	double	mu,
		schemes	scheme = `NDR`
	)		const

Definition at line 84 of file Flavour.cpp.

 {
     return getPtr<HeffDS1>(HDS1)->ComputeCoeffDS1PP(mu, scheme);
 }

◆ ComputeCoeffK()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffK	(	double	mu,
		schemes	scheme = `NDR`
	)		const

Definition at line 74 of file Flavour.cpp.

 {
     return getPtr<HeffDF2>(HDF2)->ComputeCoeffK(mu, scheme);
 }

◆ ComputeCoeffmK()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffmK	(	double	mu,
		schemes	scheme = `NDR`
	)		const

Definition at line 79 of file Flavour.cpp.

 {
     return getPtr<HeffDF2>(HDF2)->ComputeCoeffmK(mu, scheme);
 }

◆ ComputeCoeffprimeBMll()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffprimeBMll	(	double	mu,
		QCD::lepton	lepton,
		schemes	scheme = `NDR`
	)		const

Computes the chirality flipped Wilson coefficient for the process \( B \to V/P \ell^+ \ell^- \).

Parameters

[in]	mu	the lower matching scale for the process
[in]	scheme	the scheme in which the Wilson Coefficients need to be calculated

Returns: returns the chirality flipped Wilson coefficients for the process \( B \to V/P \ell^+ \ell^- \)

Definition at line 139 of file Flavour.cpp.

 {
     return getPtr<HeffDB1>(HDB1)->ComputeCoeffprimeBMll(mu, lepton, scheme);
 }

◆ ComputeCoeffprimesgamma()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffprimesgamma	(	double	mu,
		schemes	scheme = `NDR`
	)		const

Computes the chirality flipped Wilson coefficient for the process \( b \to s \gamma \).

Parameters

[in]	mu	the lower matching scale for the process
[in]	scheme	the scheme in which the Wilson Coefficients need to be calculated

Returns: returns the chirality flipped Wilson coefficients for the process \( b \to s \gamma \)

Definition at line 129 of file Flavour.cpp.

 {
     return getPtr<HeffDB1>(HDB1)->ComputeCoeffprimesgamma(mu, scheme);
 }

◆ ComputeCoeffsgamma()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffsgamma	(	double	mu,
		bool	noSM = `false`,
		schemes	scheme = `NDR`
	)		const

Computes the Wilson coefficient for the process \( b \to s \gamma \).

Parameters

[in]	mu	the lower matching scale for the process
[in]	scheme	the scheme in which the Wilson Coefficients need to be calculated

Returns: returns the Wilson coefficients for the process \( b \to s \gamma \)

Definition at line 124 of file Flavour.cpp.

 {
     return getPtr<HeffDB1>(HDB1)->ComputeCoeffsgamma(mu, noSM, scheme);
 }

◆ ComputeCoeffsmumu()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffsmumu	(	double	mu,
		schemes	scheme = `NDR`
	)		const

Computes the Wilson coefficient for the process \( B_s \to \mu \mu \).

Parameters

[in]	mu	the lower matching scale for the process
[in]	scheme	the scheme in which the Wilson Coefficients need to be calculated

Returns: returns the Wilson coefficients for the process \( B_s \to \mu \mu \)

Definition at line 99 of file Flavour.cpp.

 {
     return getPtr<HeffDB1>(HDB1)->ComputeCoeffsmumu(mu, scheme);
 }

◆ ComputeCoeffsnunu()

gslpp::vector< gslpp::complex > ** Flavour::ComputeCoeffsnunu ( ) const

Definition at line 114 of file Flavour.cpp.

 {
     return getPtr<HeffDB1>(HDB1)->ComputeCoeffsnunu();
 }

◆ getFlagCLN()

bool Flavour::getFlagCLN ( ) const

inline

Definition at line 255 of file Flavour.h.

     {
         return CLNflag;
     }

◆ getFlagFixedWCbtos()

bool Flavour::getFlagFixedWCbtos ( ) const

inline

Definition at line 260 of file Flavour.h.

     {
         return FixedWCbtosflag;
     }

◆ getFlagUseDispersionRelation()

bool Flavour::getFlagUseDispersionRelation ( ) const

inline

Definition at line 250 of file Flavour.h.

     {
         return dispersion;
     }

◆ getHDB1()

HeffDB1 & Flavour::getHDB1 ( ) const

The member that returns an object of the class HeffDB1.

Returns: returns the Hamiltonian for the \( \Delta B = 1 \) processes.

Definition at line 54 of file Flavour.cpp.

 {
     return *getPtr<HeffDB1>(HDB1);
 }

◆ getHDF2()

HeffDF2 & Flavour::getHDF2 ( ) const

The member that returns an object of the class HeffDF2.

Returns: returns the Hamiltonian for the \( \Delta F = 2 \) processes

Definition at line 44 of file Flavour.cpp.

 {
     return *getPtr<HeffDF2>(HDF2);
 }

◆ getHDS1()

HeffDS1 & Flavour::getHDS1 ( ) const

The member that returns an object of the class HeffDS1.

Returns: returns the Hamiltonian for the \( \Delta S = 1 \) processes.

Definition at line 49 of file Flavour.cpp.

 {
     return *getPtr<HeffDS1>(HDS1);
 }

◆ getM()

template<typename T , typename... Args>

T & Flavour::getM	(	std::map< std::vector< int >, std::shared_ptr< T > > &	map,
		Args ...	args
	)		const

private

Definition at line 230 of file Flavour.cpp.

 {
     std::vector<int> key({args...});
     if(map.find(key)==map.end()) {
         map.insert(std::make_pair(key,std::shared_ptr<T>(new T(mySM,args...))));
     }
     return *map.at(key);
 }

◆ getMPll()

MPll & Flavour::getMPll	(	QCD::meson	meson_i,
		QCD::meson	pseudoscalar_i,
		QCD::lepton	lep_i
	)		const

Returns the initial and final state dependent object for \( B \to P \ell^+ \ell^- \).

Parameters

[in]	meson_i	specifies the meson in the initial state
[in]	pseudoscalar_i	specifies the vector in the final state
[in]	lepton_i	specifies the lepton in the final state

Returns: returns a pointer to the initial and final state dependent object for the process \( B \to P \ell^+ \ell^- \)

Definition at line 212 of file Flavour.cpp.

 {
 //    std::reference_wrapper<std::shared_ptr<MPll> > x(MPll_BpKmu);
 //    if (meson_i == QCD::B_P && pseudoscalar_i == QCD::K_P && lep_i == QCD::MU) x = std::ref(MPll_BpKmu);
 //    else if (meson_i == QCD::B_P && pseudoscalar_i == QCD::K_P && lep_i == QCD::ELECTRON) x = std::ref(MPll_BpKel);
 //    else if (meson_i == QCD::B_D && pseudoscalar_i == QCD::K_0 && lep_i == QCD::MU) x = std::ref(MPll_B0Kmu);
 //    else if (meson_i == QCD::B_D && pseudoscalar_i == QCD::K_0 && lep_i == QCD::ELECTRON) x = std::ref(MPll_B0Kel);
 //    else throw std::runtime_error("Flavour: Decay channel not implemented.");
 //    return *getPtr<MPll>(x.get(), meson_i, pseudoscalar_i, lep_i);
 //    std::vector<int> key({meson_i, vector_i, lep_i});
 //    if(MPllMap.find(key)==MPllMap.end()) {
 //        MPllMap.insert(std::make_pair(key,std::shared_ptr<MPll>(new MPll(mySM,meson_i, vector_i, lep_i))));
 //    }
 //    return *MPllMap.at(key);
     return getM<MPll>(MPllMap, meson_i, vector_i, lep_i);
 }

◆ getMPlnu()

MPlnu & Flavour::getMPlnu	(	QCD::meson	meson_i,
		QCD::meson	pseudoscalar_i,
		QCD::lepton	lep_i
	)		const

Returns the initial and final state dependent object for \( B \to P \ell \nu \).

Parameters

[in]	meson_i	specifies the meson in the initial state
[in]	pseudoscalar_i	specifies the vector in the final state
[in]	lepton_i	specifies the lepton in the final state

Returns: returns a pointer to the initial and final state dependent object for the process \( B \to V \ell \nu \)

Definition at line 196 of file Flavour.cpp.

 {
 //    std::reference_wrapper<std::shared_ptr<MPlnu> > x(MPlnu_BdbarDtaunu);
 //    if (meson_i == StandardModel::B_D && vector_i == StandardModel::D_P && lep_i == StandardModel::TAU) x = std::ref(MPlnu_BdbarDtaunu);
 //    else if (meson_i == StandardModel::B_D && vector_i == StandardModel::D_P && lep_i == StandardModel::MU) x = std::ref(MPlnu_BdbarDmunu);
 //    else if (meson_i == StandardModel::B_D && vector_i == StandardModel::D_P && lep_i == StandardModel::ELECTRON) x = std::ref(MPlnu_BdbarDelnu);
 //    else throw std::runtime_error("Flavour: Decay channel not implemented.");
 //    return *getPtr<MPlnu>(x.get(), meson_i, vector_i, lep_i);
 //    std::vector<int> key({meson_i, vector_i, lep_i});
 //    if(MPlnuMap.find(key)==MPlnuMap.end()) {
 //        MPlnuMap.insert(std::make_pair(key,std::shared_ptr<MPlnu>(new MPlnu(mySM,meson_i, vector_i, lep_i))));
 //    }
 //    return *MPlnuMap.at(key);
     return getM<MPlnu>(MPlnuMap, meson_i, vector_i, lep_i);
 }

◆ getMVgamma()

MVgamma & Flavour::getMVgamma	(	QCD::meson	meson_i,
		QCD::meson	vector_i
	)		const

Returns the initial and final state dependent object for \( B \to V \gamma \).

Parameters

[in]	meson_i	specifies the meson in the initial state
[in]	vector_i	specifies the vector in the final state

Returns: returns a pointer to the initial and final state dependent object for the process \( B \to V \gamma \)

Definition at line 163 of file Flavour.cpp.

 {
 //    std::reference_wrapper<std::shared_ptr<MVgamma> > x(MVgamma_BdKstgamma);
 //    if (meson_i == StandardModel::B_D && vector_i == StandardModel::K_star) x = std::ref(MVgamma_BdKstgamma);
 //    else if (meson_i == StandardModel::B_P && vector_i == StandardModel::K_star_P) x = std::ref(MVgamma_BpKstgamma);
 //    else if (meson_i == StandardModel::B_S && vector_i == StandardModel::PHI) x = std::ref(MVgamma_Bsphigamma);
 //    else throw std::runtime_error("Flavour: Decay channel not implemented.");
 //    return *getPtr<MVgamma>(x.get(), meson_i, vector_i);
 //    std::vector<int> key({meson_i, vector_i});
 //    if(MVgammaMap.find(key)==MVgammaMap.end()) {
 //        MVgammaMap.insert(std::make_pair(key,std::shared_ptr<MVgamma>(new MVgamma(mySM,meson_i, vector_i))));
 //    }
 //    return *MVgammaMap.at(key);
     return getM<MVgamma>(MVgammaMap, meson_i, vector_i);
 }

◆ getMVll()

MVll & Flavour::getMVll	(	QCD::meson	meson_i,
		QCD::meson	vector_i,
		QCD::lepton	lep_i
	)		const

Returns the initial and final state dependent object for \( B \to V \ell^+ \ell^- \).

Parameters

[in]	meson_i	specifies the meson in the initial state
[in]	vector_i	specifies the vector in the final state
[in]	lepton_i	specifies the lepton in the final state

Returns: returns a pointer to the initial and final state dependent object for the process \( B \to V \ell^+ \ell^- \)

Definition at line 144 of file Flavour.cpp.

 {
 //    std::reference_wrapper<std::shared_ptr<MVll> > x(MVll_BdKstarmu);
 //    if (meson_i == StandardModel::B_D && vector_i == StandardModel::K_star && lep_i == StandardModel::MU) x = std::ref(MVll_BdKstarmu);
 //    else if (meson_i == StandardModel::B_D && vector_i == StandardModel::K_star && lep_i == StandardModel::ELECTRON) x = std::ref(MVll_BdKstarel);
 //    else if (meson_i == StandardModel::B_P && vector_i == StandardModel::K_star_P && lep_i == StandardModel::MU) x = std::ref(MVll_BpKstarmu);
 //    else if (meson_i == StandardModel::B_P && vector_i == StandardModel::K_star_P && lep_i == StandardModel::ELECTRON) x = std::ref(MVll_BpKstarel);
 //    else if (meson_i == StandardModel::B_S && vector_i == StandardModel::PHI && lep_i == StandardModel::MU) x = std::ref(MVll_Bsphimu);
 //    else if (meson_i == StandardModel::B_S && vector_i == StandardModel::PHI && lep_i == StandardModel::ELECTRON) x = std::ref(MVll_Bsphiel);
 //    else throw std::runtime_error("Flavour: Decay channel not implemented.");
 //    return *getPtr<MVll>(x.get(), meson_i, vector_i, lep_i);
 //    std::vector<int> key({meson_i, vector_i, lep_i});
 //    if(MVllMap.find(key)==MVllMap.end()) {
 //        MVllMap.insert(std::make_pair(key,std::shared_ptr<MVll>(new MVll(mySM,meson_i, vector_i, lep_i))));
 //    }
 //    return *MVllMap.at(key);
     return getM<MVll>(MVllMap, meson_i, vector_i, lep_i);
 }

◆ getMVlnu()

MVlnu & Flavour::getMVlnu	(	QCD::meson	meson_i,
		QCD::meson	vector_i,
		QCD::lepton	lep_i
	)		const

Returns the initial and final state dependent object for \( B \to V \ell \nu \).

Parameters

[in]	meson_i	specifies the meson in the initial state
[in]	vector_i	specifies the vector in the final state
[in]	lepton_i	specifies the lepton in the final state

Returns: returns a pointer to the initial and final state dependent object for the process \( B \to V \ell^+ \ell^- \)

Definition at line 179 of file Flavour.cpp.

 {
 //    std::reference_wrapper<std::shared_ptr<MVlnu> > x(MVlnu_BdbarDstartaunu);
 //    if (meson_i == StandardModel::B_D && vector_i == StandardModel::D_star_P && lep_i == StandardModel::TAU) x = std::ref(MVlnu_BdbarDstartaunu);
 //    else if (meson_i == StandardModel::B_D && vector_i == StandardModel::D_star_P && lep_i == StandardModel::MU) x = std::ref(MVlnu_BdbarDstarmunu);
 //    else if (meson_i == StandardModel::B_D && vector_i == StandardModel::D_star_P && lep_i == StandardModel::ELECTRON) x = std::ref(MVlnu_BdbarDstarelnu);
 //    else throw std::runtime_error("Flavour: Decay channel not implemented.");
 //    return *getPtr<MVlnu>(x.get(), meson_i, vector_i, lep_i);
 //    std::vector<int> key({meson_i, vector_i, lep_i});
 //    if(MVlnuMap.find(key)==MVlnuMap.end()) {
 //        MVlnuMap.insert(std::make_pair(key,std::shared_ptr<MVlnu>(new MVlnu(mySM,meson_i, vector_i, lep_i))));
 //    }
 //    return *MVlnuMap.at(key);
     return getM<MVlnu>(MVlnuMap, meson_i, vector_i, lep_i);
  
 }

◆ getPtr() [1/2]

template<typename T , typename... Args>

std::shared_ptr<T>& Flavour::getPtr	(	std::shared_ptr< T > &	x,
		Args ...	args
	)		const

private

Definition at line 295 of file Flavour.cpp.

 {
     if (x.get() == nullptr)
         x.reset(new T(mySM, args...));
     return x;
 }

◆ getPtr() [2/2]

template<typename T , typename... Args>

std::shared_ptr<T>& Flavour::getPtr	(	std::shared_ptr< T > &	x,
		Args...	args
	)		const

private

◆ getUpdateFlag()

bool Flavour::getUpdateFlag	(	QCD::meson	meson_i,
		QCD::meson	meson_j,
		QCD::lepton	lep_i
	)		const

gets the update flag for the initial and final state dependent object for \( B \to V \ell^+ \ell^- \).

Parameters

[in]	meson_i	specifies the meson in the initial state
[in]	vector_i	specifies the vector in the final state
[in]	lepton_i	specifies the lepton in the final state

Definition at line 264 of file Flavour.cpp.

 {
 //    if (meson_i == StandardModel::B_D && meson_j == StandardModel::K_star && lep_i == StandardModel::MU) return update_BdKstarmu;
 //    if (meson_i == StandardModel::B_D && meson_j == StandardModel::K_star && lep_i == StandardModel::ELECTRON) return update_BdKstarel;
 //    if (meson_i == StandardModel::B_P && meson_j == StandardModel::K_star_P && lep_i == StandardModel::MU) return update_BpKstarmu;
 //    if (meson_i == StandardModel::B_P && meson_j == StandardModel::K_star_P && lep_i == StandardModel::ELECTRON) return update_BpKstarel;
 //    if (meson_i == StandardModel::B_S && meson_j == StandardModel::PHI && lep_i == StandardModel::MU) return update_Bsphimu;
 //    if (meson_i == StandardModel::B_S && meson_j == StandardModel::PHI && lep_i == StandardModel::ELECTRON) return update_Bsphiel;
 //    if (meson_i == StandardModel::B_P && meson_j == StandardModel::K_P && lep_i == StandardModel::MU) return update_BpKmu;
 //    if (meson_i == StandardModel::B_P && meson_j == StandardModel::K_P && lep_i == StandardModel::ELECTRON) return update_BpKel;
 //    if (meson_i == StandardModel::B_D && meson_j == StandardModel::K_0 && lep_i == StandardModel::MU) return update_B0Kmu;
 //    if (meson_i == StandardModel::B_D && meson_j == StandardModel::K_0 && lep_i == StandardModel::ELECTRON) return update_B0Kel;
 //    if (meson_i == StandardModel::B_D && meson_j == StandardModel::K_star && lep_i == StandardModel::NOLEPTON) return update_BdKstgamma;
 //    if (meson_i == StandardModel::B_P && meson_j == StandardModel::K_star_P && lep_i == StandardModel::NOLEPTON) return update_BpKstgamma;
 //    if (meson_i == StandardModel::B_S && meson_j == StandardModel::PHI && lep_i == StandardModel::NOLEPTON) return update_Bsphigamma;
 //    if (meson_i == StandardModel::B_D && meson_j == StandardModel::D_star_P && lep_i == StandardModel::TAU) return update_BdDstartaunu;
 //    if (meson_i == StandardModel::B_D && meson_j == StandardModel::D_star_P && lep_i == StandardModel::MU) return update_BdDstarmunu;
 //    if (meson_i == StandardModel::B_D && meson_j == StandardModel::D_star_P && lep_i == StandardModel::ELECTRON) return update_BdDstarelnu;
 //    if (meson_i == StandardModel::B_D && meson_j == StandardModel::D_P && lep_i == StandardModel::TAU) return update_BdDtaunu;
 //    if (meson_i == StandardModel::B_D && meson_j == StandardModel::D_P && lep_i == StandardModel::MU) return update_BdDmunu;
 //    if (meson_i == StandardModel::B_D && meson_j == StandardModel::D_P && lep_i == StandardModel::ELECTRON) return update_BdDelnu;
 //    else throw std::runtime_error("Flavour: Wrong update flags requested.");
 //
     std::vector<int> key({meson_i, meson_j, lep_i});
     if(flagUpdateMap.find(key)==flagUpdateMap.end()) {
         // Default value
         flagUpdateMap[key]=true;
     }
     return flagUpdateMap.at(key);
 }

◆ setFlag()

bool Flavour::setFlag	(	const std::string	name,
		const bool	value
	)

Definition at line 27 of file Flavour.cpp.

 {
     
     if (name.compare("UseDispersionRelation") == 0) {
         setFlagUseDispersionRelation(value);
         return true;
     } else if (name.compare("CLNflag") == 0) {
         setFlagCLN(value);
         return true;
     } else if (name.compare("FixedWCbtos") == 0) {
         setFlagFixedWCbtos(value);
         return true;
     } else
         return false;
 }

◆ setFlagCLN()

bool Flavour::setFlagCLN ( bool CLNflag )

inline

Definition at line 240 of file Flavour.h.

     {
         return (this->CLNflag = CLNflag);
     }

◆ setFlagFixedWCbtos()

bool Flavour::setFlagFixedWCbtos ( bool FixedWCbtosflag )

inline

Definition at line 245 of file Flavour.h.

     {
         return (this->FixedWCbtosflag = FixedWCbtosflag);
     }

◆ setFlagUseDispersionRelation()

bool Flavour::setFlagUseDispersionRelation ( bool dispersion )

inline

Definition at line 235 of file Flavour.h.

     {
         return (this->dispersion = dispersion);
     }

◆ setSMupdated()

void Flavour::setSMupdated ( ) const

a member used for the caching for \( B \to V \ell^+ \ell^- \).

Definition at line 302 of file Flavour.cpp.

 {
     for(auto &x : flagUpdateMap) x.second = true;
 }

◆ setUpdateFlag()

void Flavour::setUpdateFlag	(	QCD::meson	meson_i,
		QCD::meson	meson_j,
		QCD::lepton	lep_i,
		bool	updated_i
	)		const

sets the update flag for the initial and final state dependent object for \( B \to V \ell^+ \ell^- \).

Parameters

[in]	meson_i	specifies the meson in the initial state
[in]	vector_i	specifies the vector in the final state
[in]	lepton_i	specifies the lepton in the final state

Definition at line 239 of file Flavour.cpp.

 {
 //    if (meson_i == StandardModel::B_D && meson_j == StandardModel::K_star && lep_i == StandardModel::MU) update_BdKstarmu = updated_i;
 //    else if (meson_i == StandardModel::B_D && meson_j == StandardModel::K_star && lep_i == StandardModel::ELECTRON) update_BdKstarel = updated_i;
 //    else if (meson_i == StandardModel::B_P && meson_j == StandardModel::K_star_P && lep_i == StandardModel::MU) update_BpKstarmu = updated_i;
 //    else if (meson_i == StandardModel::B_P && meson_j == StandardModel::K_star_P && lep_i == StandardModel::ELECTRON) update_BpKstarel = updated_i;
 //    else if (meson_i == StandardModel::B_S && meson_j == StandardModel::PHI && lep_i == StandardModel::MU) update_Bsphimu = updated_i;
 //    else if (meson_i == StandardModel::B_S && meson_j == StandardModel::PHI && lep_i == StandardModel::ELECTRON) update_Bsphiel = updated_i;
 //    else if (meson_i == StandardModel::B_P && meson_j == StandardModel::K_P && lep_i == StandardModel::MU) update_BpKmu = updated_i;
 //    else if (meson_i == StandardModel::B_P && meson_j == StandardModel::K_P && lep_i == StandardModel::ELECTRON) update_BpKel = updated_i;
 //    else if (meson_i == StandardModel::B_D && meson_j == StandardModel::K_0 && lep_i == StandardModel::MU) update_B0Kmu = updated_i;
 //    else if (meson_i == StandardModel::B_D && meson_j == StandardModel::K_0 && lep_i == StandardModel::ELECTRON) update_B0Kel = updated_i;
 //    else if (meson_i == StandardModel::B_D && meson_j == StandardModel::K_star && lep_i == StandardModel::NOLEPTON) update_BdKstgamma = updated_i;
 //    else if (meson_i == StandardModel::B_P && meson_j == StandardModel::K_star_P && lep_i == StandardModel::NOLEPTON) update_BpKstgamma = updated_i;
 //    else if (meson_i == StandardModel::B_S && meson_j == StandardModel::PHI && lep_i == StandardModel::NOLEPTON) update_Bsphigamma = updated_i;
 //    else if (meson_i == StandardModel::B_D && meson_j == StandardModel::D_star_P && lep_i == StandardModel::TAU) update_BdDstartaunu = updated_i;
 //    else if (meson_i == StandardModel::B_D && meson_j == StandardModel::D_star_P && lep_i == StandardModel::MU) update_BdDstarmunu = updated_i;
 //    else if (meson_i == StandardModel::B_D && meson_j == StandardModel::D_star_P && lep_i == StandardModel::ELECTRON) update_BdDstarelnu = updated_i;
 //    else if (meson_i == StandardModel::B_D && meson_j == StandardModel::D_P && lep_i == StandardModel::TAU) update_BdDtaunu = updated_i;
 //    else if (meson_i == StandardModel::B_D && meson_j == StandardModel::D_P && lep_i == StandardModel::MU) update_BdDmunu = updated_i;
 //    else if (meson_i == StandardModel::B_D && meson_j == StandardModel::D_P && lep_i == StandardModel::ELECTRON) update_BdDelnu = updated_i;
 //    else throw std::runtime_error("Flavour: Wrong update flag requested.");
     flagUpdateMap[std::vector<int>({meson_i, meson_j, lep_i})] = updated_i;
 }

Member Data Documentation

◆ CLNflag

bool Flavour::CLNflag

mutableprivate

Definition at line 280 of file Flavour.h.

◆ dispersion

bool Flavour::dispersion

mutableprivate

Definition at line 279 of file Flavour.h.

◆ FixedWCbtosflag

bool Flavour::FixedWCbtosflag

mutableprivate

Definition at line 281 of file Flavour.h.

◆ flagUpdateMap

std::map<std::vector<int>, bool> Flavour::flagUpdateMap

mutableprivate

Definition at line 277 of file Flavour.h.

◆ HDB1

std::shared_ptr<HeffDB1> Flavour::HDB1

mutableprivate

An Object for the Hamiltonian of the \( \Delta B = 1 \) processes.

Definition at line 270 of file Flavour.h.

◆ HDF2

std::shared_ptr<HeffDF2> Flavour::HDF2

mutableprivate

An Object for the Hamiltonian of the \( \Delta F = 2 \) processes.

Definition at line 269 of file Flavour.h.

◆ HDS1

std::shared_ptr<HeffDS1> Flavour::HDS1

mutableprivate

An Object for the Hamiltonian of the \( \Delta S = 1 \) processes.

Definition at line 271 of file Flavour.h.

◆ MPllMap

std::map<std::vector<int>, std::shared_ptr<MPll> > Flavour::MPllMap

mutableprivate

Definition at line 275 of file Flavour.h.

◆ MPlnuMap

std::map<std::vector<int>, std::shared_ptr<MPlnu> > Flavour::MPlnuMap

mutableprivate

Definition at line 276 of file Flavour.h.

◆ MVgammaMap

std::map<std::vector<int>, std::shared_ptr<MVgamma> > Flavour::MVgammaMap

mutableprivate

Definition at line 274 of file Flavour.h.

◆ MVllMap

std::map<std::vector<int>, std::shared_ptr<MVll> > Flavour::MVllMap

mutableprivate

Definition at line 272 of file Flavour.h.

◆ MVlnuMap

std::map<std::vector<int>, std::shared_ptr<MVlnu> > Flavour::MVlnuMap

mutableprivate

Definition at line 273 of file Flavour.h.

◆ mySM

const StandardModel& Flavour::mySM

private

Definition at line 268 of file Flavour.h.

The documentation for this class was generated from the following files:

Detailed Description

Public Member Functions

Private Member Functions

Private Attributes

Constructor & Destructor Documentation

◆ Flavour()

Member Function Documentation

◆ ComputeCoeffBd()

◆ ComputeCoeffBMll()

◆ ComputeCoeffBs()

◆ ComputeCoeffbtaunu()

◆ ComputeCoeffdd()

◆ ComputeCoeffdmumu()

◆ ComputeCoeffdnunu()

◆ ComputeCoeffDS1mumu()

◆ ComputeCoeffDS1pnunu()

◆ ComputeCoeffDS1PP()

◆ ComputeCoeffK()

◆ ComputeCoeffmK()

◆ ComputeCoeffprimeBMll()

◆ ComputeCoeffprimesgamma()

◆ ComputeCoeffsgamma()

◆ ComputeCoeffsmumu()

◆ ComputeCoeffsnunu()

◆ getFlagCLN()

◆ getFlagFixedWCbtos()

◆ getFlagUseDispersionRelation()

◆ getHDB1()

◆ getHDF2()

◆ getHDS1()

◆ getM()

◆ getMPll()

◆ getMPlnu()

◆ getMVgamma()

◆ getMVll()

◆ getMVlnu()

◆ getPtr() [1/2]

◆ getPtr() [2/2]

◆ getUpdateFlag()

◆ setFlag()

◆ setFlagCLN()

◆ setFlagFixedWCbtos()

◆ setFlagUseDispersionRelation()

◆ setSMupdated()

◆ setUpdateFlag()

Member Data Documentation

◆ CLNflag

◆ dispersion

◆ FixedWCbtosflag

◆ flagUpdateMap

◆ HDB1

◆ HDF2

◆ HDS1

◆ MPllMap

◆ MPlnuMap

◆ MVgammaMap

◆ MVllMap

◆ MVlnuMap

◆ mySM