#include <Bdmumu.h>

Inheritance diagram for Bdmumu:

Detailed Description

Definition at line 19 of file Bdmumu.h.

Public Member Functions
	Bdmumu (const StandardModel &SM_i, int obsFlag, QCD::lepton lep_i=QCD::MU)

double	computeAmumu (orders order)

double	computeSmumu (orders order)

double	computeThValue ()
	hep-ph/9512380v2 More...

Public Member Functions inherited from ThObservable
double	getBinMax ()
	A get method to get the maximum value of the bin. More...

double	getBinMin ()
	A get method to get the minimum value of the bin. More...

const StandardModel &	getModel ()
	A get method to get the model. More...

const std::vector< std::string >	getParametersForObservable ()
	A get method to get the parameters for the specific observable. More...

void	setBinMax (double max)
	A set method to set the maximum value of the bin. More...

void	setBinMin (double min)
	A set method to set the minimum value of the bin. More...

void	setParametersForObservable (std::vector< std::string > parametersForObservable_i)
	A set method to get the parameters for the specific observable. More...

	ThObservable (const StandardModel &SM_i)
	Constructor. More...

	ThObservable (const ThObservable &orig)
	The copy constructor. More...

virtual	~ThObservable ()
	The default destructor. More...

Protected Member Functions
void	computeAmpSq (orders order, orders_qed order_qed, double mu)

void	computeObs (orders order, orders_qed order_qed)

Private Attributes
double	absP

double	absS

gslpp::vector< gslpp::complex > **	allcoeff

double	ampSq

double	Amumu

double	argP

double	argS

double	beta

gslpp::complex	C_10

gslpp::complex	C_10p

gslpp::complex	C_P

gslpp::complex	C_Pp

gslpp::complex	C_S

gslpp::complex	C_Sp

double	chiral

double	coupling

std::unique_ptr< EvolBsmm >	evolbdmm

QCD::lepton	lep

int	leptonindex

double	mb

double	mBd

double	md

double	mlep

double	mW

int	obs

double	phiNP

double	Smumu

double	timeInt

double	yd

Additional Inherited Members
Protected Attributes inherited from ThObservable
double	max
	the bin maximum. More...

double	min
	The bin minimum. More...

std::vector< std::string >	parametersForObservable
	a vector of parameter namesfor the specific observable More...

const StandardModel &	SM
	A reference to an object of StandardMode class. More...

Constructor & Destructor Documentation

◆ Bdmumu()

Bdmumu::Bdmumu	(	const StandardModel &	SM_i,
		int	obsFlag,
		QCD::lepton	lep_i = `QCD::MU`
	)

constructor

Parameters

Flavour

Definition at line 14 of file Bdmumu.cpp.

: ThObservable(SM_i),
  evolbdmm(new EvolBsmm(8, NDR, NNLO, NLO_QED22, SM))
{  
    lep = lep_i;
    if(lep == QCD::MU) leptonindex = 1;
    else if(lep == QCD::ELECTRON) leptonindex = 0;
    else if(lep == QCD::TAU) leptonindex = 2;
    else throw std::runtime_error("Bsmumu::Bsmumu(): lepton type not defined. Can be only MU or ELECTRON");
    if (obsFlag > 0 and obsFlag < 5) obs = obsFlag;
    else throw std::runtime_error("obsFlag in Bdmumu(myFlavour, obsFlag) called from ThFactory::ThFactory() can only be 1 (BR) or 2 (BRbar) or 3 (Amumu) or 4 (Smumu)");
    SM.initializeMeson(QCD::B_D);
};

Member Function Documentation

◆ computeAmpSq()

void Bdmumu::computeAmpSq	(	orders	order,
		orders_qed	order_qed,
		double	mu
	)

protected

Parameters

order
order_qed

Returns: the short distance contribution to the | \( BR(B_s \rightarrow \mu \bar{\mu}) \)|

Definition at line 80 of file Bdmumu.cpp.

{
    if (SM.getFlavour().getHDB1().getCoeffdmumu().getOrder() < order % 3){
        std::stringstream out;
        out << order;
        throw std::runtime_error("Bdmumu::computeAmpSq(): required cofficient of "
                                 "order " + out.str() + " not computed");
    }
    allcoeff = SM.getFlavour().ComputeCoeffdmumu(mu, NDR);
    
    double alsmu = evolbdmm->alphatilde_s(mu);
    double alemu = evolbdmm->alphatilde_e(mu);
//    double alemu = SM.ale_OS(mu)/4./M_PI; // to be checked
    gslpp::matrix<gslpp::complex> Vckm = SM.getVCKM();
    
    if((order == FULLNLO) && (order_qed == FULLNLO_QED)){
    
    switch(order_qed) {
        case FULLNLO_QED:
        {
            C_10 = (*(allcoeff[LO]))(7) /alemu  + (*(allcoeff[NLO]))(7) * alsmu/alemu 
                    + (*(allcoeff[NNLO]))(7) * alsmu * alsmu/alemu + (*(allcoeff[LO_QED ]))(7) /alsmu
                    + (*(allcoeff[NLO_QED11]))(7) + (*(allcoeff[NLO_QED02]))(7) * alemu /alsmu /alsmu 
                    + (*(allcoeff[NLO_QED21]))(7) * alsmu 
                    + (*(allcoeff[NLO_QED12]))(7) * alemu /alsmu+ (*(allcoeff[NLO_QED22]))(7) * alemu;
 
            gslpp::complex NPfactor = coupling;
 
            if(SM.getModelName().compare("NPSMEFTd6U2") == 0 || SM.getModelName().compare("NPSMEFTd6U3") == 0)
            {
                C_10 = C_10 + (dynamic_cast<const NPSMEFTd6GeneralMatching&>(SM.getMatching()).getCdeVLR(0,2,leptonindex,leptonindex) - 
                    dynamic_cast<const NPSMEFTd6GeneralMatching&>(SM.getMatching()).getCedVLL(leptonindex,leptonindex,0,2)) / NPfactor; 
            }
 
 
 
            absP = C_10.abs();
            argP = C_10.arg();
           
            phiNP = 0.;
            
            ampSq = absP * absP ;
                  
        }
        break;
        default:
            std::stringstream out;
            out << order;
            throw std::runtime_error("Bdmumu::computeAmpSq(): order " + out.str() + " not implemented");;
    }
    }
    
    
}

◆ computeAmumu()

double Bdmumu::computeAmumu ( orders order )

Definition at line 68 of file Bdmumu.cpp.

{
    computeObs(FULLNLO, FULLNLO_QED);
    return(Amumu);
}

◆ computeObs()

void Bdmumu::computeObs	(	orders	order,
		orders_qed	order_qed
	)

protected

Definition at line 52 of file Bdmumu.cpp.

{
    double mu = SM.getMub();  
        
    mlep = SM.getLeptons(lep).getMass();
    mBd = SM.getMesons(QCD::B_D).getMass();
    mW = SM.Mw();
    mb = SM.getQuarks(QCD::BOTTOM).getMass();
    md = SM.getQuarks(QCD::DOWN).getMass();
    chiral = pow(mBd, 2.) / 2. / mlep * mb / (mb + md);
    beta = sqrt(1. - pow(2. * mlep / mBd, 2.));
    computeAmpSq(order, order_qed, mu);
    Amumu = (absP * absP * cos(2. * argP - phiNP) -  absS * absS * cos(2. * argS - phiNP)) / (absP * absP + absS * absS);
    Smumu = (absP * absP * sin(2. * argP - phiNP) -  absS * absS * sin(2. * argS - phiNP)) / (absP * absP + absS * absS);
}

◆ computeSmumu()

double Bdmumu::computeSmumu ( orders order )

Definition at line 74 of file Bdmumu.cpp.

{
    computeObs(FULLNLO, FULLNLO_QED);
    return(Smumu);
}

◆ computeThValue()

double Bdmumu::computeThValue ( )

virtual

hep-ph/9512380v2

Returns: theoretical value of | \( BR(B_d \rightarrow \mu \bar{\mu}) \)|

Implements ThObservable.

Definition at line 28 of file Bdmumu.cpp.

{   
    double FBd = SM.getMesons(QCD::B_D).getDecayconst();
    
    double Mw = SM.Mw();
    double GF = SM.getGF();
    coupling = GF * GF * Mw * Mw /M_PI /M_PI ; 
 
    computeObs(FULLNLO, FULLNLO_QED);
    
 
    double PRF = pow(coupling, 2.) / M_PI /8. / SM.getMesons(QCD::B_D).computeWidth() * pow(FBd, 2.) * pow(mlep, 2.) * mBd * beta;
    yd = 0.; // Use the experimental number here
    timeInt = (1. + Amumu * yd) / (1. - yd * yd); // Note modification in form due to algorithm
    
    if (obs == 1) return( PRF * ampSq);
    if (obs == 2) return( PRF * ampSq * timeInt);
    if (obs == 3) return( Amumu );
    if (obs == 4) return( Smumu );
 
    throw std::runtime_error("Bdmumu::computeThValue(): Observable type not defined. Can be only any of (1,2,3,4)");
    return (EXIT_FAILURE);
}

Member Data Documentation

◆ absP

double Bdmumu::absP

private

Definition at line 59 of file Bdmumu.h.

◆ absS

double Bdmumu::absS

private

Definition at line 61 of file Bdmumu.h.

◆ allcoeff

gslpp::vector<gslpp::complex>** Bdmumu::allcoeff

private

Definition at line 80 of file Bdmumu.h.

◆ ampSq

double Bdmumu::ampSq

private

Definition at line 63 of file Bdmumu.h.

◆ Amumu

double Bdmumu::Amumu

private

Definition at line 64 of file Bdmumu.h.

◆ argP

double Bdmumu::argP

private

Definition at line 60 of file Bdmumu.h.

◆ argS

double Bdmumu::argS

private

Definition at line 62 of file Bdmumu.h.

◆ beta

double Bdmumu::beta

private

Definition at line 52 of file Bdmumu.h.

◆ C_10

gslpp::complex Bdmumu::C_10

private

Definition at line 72 of file Bdmumu.h.

◆ C_10p

gslpp::complex Bdmumu::C_10p

private

Definition at line 73 of file Bdmumu.h.

◆ C_P

gslpp::complex Bdmumu::C_P

private

Definition at line 76 of file Bdmumu.h.

◆ C_Pp

gslpp::complex Bdmumu::C_Pp

private

Definition at line 77 of file Bdmumu.h.

◆ C_S

gslpp::complex Bdmumu::C_S

private

Definition at line 74 of file Bdmumu.h.

◆ C_Sp

gslpp::complex Bdmumu::C_Sp

private

Definition at line 75 of file Bdmumu.h.

◆ chiral

double Bdmumu::chiral

private

Definition at line 58 of file Bdmumu.h.

◆ coupling

double Bdmumu::coupling

private

Definition at line 69 of file Bdmumu.h.

◆ evolbdmm

std::unique_ptr<EvolBsmm> Bdmumu::evolbdmm

private

Definition at line 79 of file Bdmumu.h.

◆ lep

QCD::lepton Bdmumu::lep

private

Final leptons type.

Definition at line 51 of file Bdmumu.h.

◆ leptonindex

int Bdmumu::leptonindex

private

Definition at line 71 of file Bdmumu.h.

◆ mb

double Bdmumu::mb

private

Definition at line 56 of file Bdmumu.h.

◆ mBd

double Bdmumu::mBd

private

Definition at line 53 of file Bdmumu.h.

◆ md

double Bdmumu::md

private

Definition at line 57 of file Bdmumu.h.

◆ mlep

double Bdmumu::mlep

private

Definition at line 55 of file Bdmumu.h.

◆ mW

double Bdmumu::mW

private

Definition at line 54 of file Bdmumu.h.

◆ obs

int Bdmumu::obs

private

Definition at line 70 of file Bdmumu.h.

◆ phiNP

double Bdmumu::phiNP

private

Definition at line 66 of file Bdmumu.h.

◆ Smumu

double Bdmumu::Smumu

private

Definition at line 65 of file Bdmumu.h.

◆ timeInt

double Bdmumu::timeInt

private

Definition at line 67 of file Bdmumu.h.

◆ yd

double Bdmumu::yd

private

Definition at line 68 of file Bdmumu.h.

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

Detailed Description

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ Bdmumu()

Member Function Documentation

◆ computeAmpSq()

◆ computeAmumu()

◆ computeObs()

◆ computeSmumu()

◆ computeThValue()

Member Data Documentation

◆ absP

◆ absS

◆ allcoeff

◆ ampSq

◆ Amumu

◆ argP

◆ argS

◆ beta

◆ C_10

◆ C_10p

◆ C_P

◆ C_Pp

◆ C_S

◆ C_Sp

◆ chiral

◆ coupling

◆ evolbdmm

◆ lep

◆ leptonindex

◆ mb

◆ mBd

◆ md

◆ mlep

◆ mW

◆ obs

◆ phiNP

◆ Smumu

◆ timeInt

◆ yd