#include <Bsmumu.h>

Inheritance diagram for Bsmumu:

Detailed Description

Definition at line 19 of file Bsmumu.h.

Public Member Functions
	Bsmumu (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 >	evolbsmm

QCD::lepton	lep

int	leptonindex

double	mb

double	mBs

double	mlep

double	ms

double	mW

int	obs

double	phiNP

double	Smumu

double	timeInt

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

◆ Bsmumu()

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

constructor

Parameters

Flavour

Definition at line 14 of file Bsmumu.cpp.

: ThObservable(SM_i),
  evolbsmm(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 Bsmumu(myFlavour, obsFlag) called from ThFactory::ThFactory() can only be 1 (BR) or 2 (BRbar) or 3 (Amumu) or 4 (Smumu)");
    SM.initializeMeson(QCD::B_S);
};

Member Function Documentation

◆ computeAmpSq()

void Bsmumu::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 Bsmumu.cpp.

{  
    if (SM.getFlavour().getHDB1().getCoeffsmumu().getOrder() < order % 3){
        std::stringstream out;
        out << order;
        throw std::runtime_error("Bsmumu::computeAmpSq(): required cofficient of "
                                 "order " + out.str() + " not computed");
    }
    allcoeff = SM.getFlavour().ComputeCoeffsmumu(mu, NDR);
 
    double alsmu = evolbsmm->alphatilde_s(mu);
    double alemu = evolbsmm->alphatilde_e(mu);
//    double alemu = SM.ale_OS(mu)/4./M_PI; // to be checked
    gslpp::matrix<gslpp::complex> Vckm = SM.getVCKM();
//    double sw2 =  (M_PI * SM.getAle() ) / ( sqrt(2.) * SM.getGF() * SM.Mw() * SM.Mw()) ;
    
    C_10p = 0.;
    C_S = 0.;
    C_Sp = 0.;
    C_P = 0.;
    C_Pp = 0.;
   
    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(1,2,leptonindex,leptonindex) - 
                    dynamic_cast<const NPSMEFTd6GeneralMatching&>(SM.getMatching()).getCedVLL(leptonindex,leptonindex,0,2)) / NPfactor; 
            }
 
            gslpp::complex CC_P = C_10 + NPfactor * ( /*C10_NP*/ - C_10p + mBs*mBs*mb / ( 2.*mlep*(mb+ms)*mW ) * (C_P - C_Pp) );
 
                
            absP = CC_P.abs(); //contains only SM contributions (P, P', S, S' not added)
            argP = CC_P.arg();
            
            gslpp::complex CC_S = NPfactor * ( beta * mBs*mBs*mb / ( 2.*mlep*(mb+ms)*mW ) * (C_S - C_Sp) );
 
            absS = CC_S.abs();
            argS = CC_S.arg();
           
            phiNP = 0.;
            
            ampSq = absP * absP ;
                  
        }
        break;
        default:
            std::stringstream out;
            out << order;
            throw std::runtime_error("Bsmumu::computeAmpSq(): order " + out.str() + " not implemented");;
        }
    }
        
}

◆ computeAmumu()

double Bsmumu::computeAmumu ( orders order )

Definition at line 68 of file Bsmumu.cpp.

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

◆ computeObs()

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

protected

Definition at line 52 of file Bsmumu.cpp.

{   
    double mu = SM.getMub();  
    
    mlep = SM.getLeptons(lep).getMass();
    mBs = SM.getMesons(QCD::B_S).getMass();
    mW = SM.Mw();
    mb = SM.getQuarks(QCD::BOTTOM).getMass();
    ms = SM.getQuarks(QCD::STRANGE).getMass();
    chiral = pow(mBs, 2.) / 2. / mlep * mb / (mb + ms);
    beta = sqrt(1. - pow(2. * mlep / mBs, 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 Bsmumu::computeSmumu ( orders order )

Definition at line 74 of file Bsmumu.cpp.

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

◆ computeThValue()

double Bsmumu::computeThValue ( )

virtual

hep-ph/9512380v2

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

Implements ThObservable.

Definition at line 28 of file Bsmumu.cpp.

{
    double FBs = SM.getMesons(QCD::B_S).getDecayconst();
    
    double Mw = SM.Mw();
    double GF = SM.getGF();
    coupling = GF * GF * Mw * Mw /M_PI /M_PI ; 
    //convertFromSingletoDoubleGF = 2. * M_SQRT2 / GF * SM.getAle() / 4. * M_PI / Mw / Mw; /* Single GF for excluding EW corrections*/
 
    computeObs(FULLNLO, FULLNLO_QED);
 
    double PRF = pow(coupling, 2.) / M_PI /8. / SM.getMesons(QCD::B_S).computeWidth() * pow(FBs, 2.) * pow(mlep, 2.) * mBs * beta;
    double ys = SM.getMesons(QCD::B_S).getDgamma_gamma()/2.; // Using the experimental number here
    timeInt = (1. + Amumu * ys) / (1. - ys * ys); // 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("Bsmumu::computeThValue(): Observable type not defined. Can be only any of (1,2,3,4)");
    return (EXIT_FAILURE);
}

Member Data Documentation

◆ absP

double Bsmumu::absP

private

Definition at line 59 of file Bsmumu.h.

◆ absS

double Bsmumu::absS

private

Definition at line 61 of file Bsmumu.h.

◆ allcoeff

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

private

Definition at line 80 of file Bsmumu.h.

◆ ampSq

double Bsmumu::ampSq

private

Definition at line 63 of file Bsmumu.h.

◆ Amumu

double Bsmumu::Amumu

private

Definition at line 64 of file Bsmumu.h.

◆ argP

double Bsmumu::argP

private

Definition at line 60 of file Bsmumu.h.

◆ argS

double Bsmumu::argS

private

Definition at line 62 of file Bsmumu.h.

◆ beta

double Bsmumu::beta

private

Definition at line 52 of file Bsmumu.h.

◆ C_10

gslpp::complex Bsmumu::C_10

private

Definition at line 72 of file Bsmumu.h.

◆ C_10p

gslpp::complex Bsmumu::C_10p

private

Definition at line 73 of file Bsmumu.h.

◆ C_P

gslpp::complex Bsmumu::C_P

private

Definition at line 76 of file Bsmumu.h.

◆ C_Pp

gslpp::complex Bsmumu::C_Pp

private

Definition at line 77 of file Bsmumu.h.

◆ C_S

gslpp::complex Bsmumu::C_S

private

Definition at line 74 of file Bsmumu.h.

◆ C_Sp

gslpp::complex Bsmumu::C_Sp

private

Definition at line 75 of file Bsmumu.h.

◆ chiral

double Bsmumu::chiral

private

Definition at line 58 of file Bsmumu.h.

◆ coupling

double Bsmumu::coupling

private

Definition at line 69 of file Bsmumu.h.

◆ evolbsmm

std::unique_ptr<EvolBsmm> Bsmumu::evolbsmm

private

Definition at line 79 of file Bsmumu.h.

◆ lep

QCD::lepton Bsmumu::lep

private

Final leptons type.

Definition at line 51 of file Bsmumu.h.

◆ leptonindex

int Bsmumu::leptonindex

private

Definition at line 71 of file Bsmumu.h.

◆ mb

double Bsmumu::mb

private

Definition at line 56 of file Bsmumu.h.

◆ mBs

double Bsmumu::mBs

private

Definition at line 53 of file Bsmumu.h.

◆ mlep

double Bsmumu::mlep

private

Definition at line 55 of file Bsmumu.h.

◆ ms

double Bsmumu::ms

private

Definition at line 57 of file Bsmumu.h.

◆ mW

double Bsmumu::mW

private

Definition at line 54 of file Bsmumu.h.

◆ obs

int Bsmumu::obs

private

Definition at line 70 of file Bsmumu.h.

◆ phiNP

double Bsmumu::phiNP

private

Definition at line 66 of file Bsmumu.h.

◆ Smumu

double Bsmumu::Smumu

private

Definition at line 65 of file Bsmumu.h.

◆ timeInt

double Bsmumu::timeInt

private

Definition at line 67 of file Bsmumu.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

◆ Bsmumu()

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

◆ evolbsmm

◆ lep

◆ leptonindex

◆ mb

◆ mBs

◆ mlep

◆ ms

◆ mW

◆ obs

◆ phiNP

◆ Smumu

◆ timeInt