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 32 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...

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::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 19 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 53 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 113 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 58 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 88 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 63 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 83 of file Flavour.cpp.

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

◆ ComputeCoeffdnunu()

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

Definition at line 98 of file Flavour.cpp.

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

◆ ComputeCoeffDS1mumu()

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

◆ ComputeCoeffDS1pnunu()

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

◆ ComputeCoeffDS1PP()

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

◆ ComputeCoeffK()

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

Definition at line 68 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 73 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 118 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 108 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 103 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 78 of file Flavour.cpp.

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

◆ ComputeCoeffsnunu()

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

Definition at line 93 of file Flavour.cpp.

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

◆ getFlagCLN()

bool Flavour::getFlagCLN ( ) const

inline

Definition at line 247 of file Flavour.h.

     {
         return CLNflag;
     }

◆ getFlagFixedWCbtos()

bool Flavour::getFlagFixedWCbtos ( ) const

inline

Definition at line 252 of file Flavour.h.

     {
         return FixedWCbtosflag;
     }

◆ getFlagUseDispersionRelation()

bool Flavour::getFlagUseDispersionRelation ( ) const

inline

Definition at line 242 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 48 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 43 of file Flavour.cpp.

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

◆ 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 209 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 191 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 175 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 142 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 123 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 158 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 274 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 243 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 26 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 232 of file Flavour.h.

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

◆ setFlagFixedWCbtos()

bool Flavour::setFlagFixedWCbtos ( bool FixedWCbtosflag )

inline

Definition at line 237 of file Flavour.h.

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

◆ setFlagUseDispersionRelation()

bool Flavour::setFlagUseDispersionRelation ( bool dispersion )

inline

Definition at line 227 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 281 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 218 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 271 of file Flavour.h.

◆ dispersion

bool Flavour::dispersion

mutableprivate

Definition at line 270 of file Flavour.h.

◆ FixedWCbtosflag

bool Flavour::FixedWCbtosflag

mutableprivate

Definition at line 272 of file Flavour.h.

◆ flagUpdateMap

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

mutableprivate

Definition at line 268 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 262 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 261 of file Flavour.h.

◆ MPllMap

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

mutableprivate

Definition at line 266 of file Flavour.h.

◆ MPlnuMap

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

mutableprivate

Definition at line 267 of file Flavour.h.

◆ MVgammaMap

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

mutableprivate

Definition at line 265 of file Flavour.h.

◆ MVllMap

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

mutableprivate

Definition at line 263 of file Flavour.h.

◆ MVlnuMap

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

mutableprivate

Definition at line 264 of file Flavour.h.

◆ mySM

const StandardModel& Flavour::mySM

private

Definition at line 260 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()

◆ 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

◆ MPllMap

◆ MPlnuMap

◆ MVgammaMap

◆ MVllMap

◆ MVlnuMap

◆ mySM