Bsmumu Class Reference

#include <Bsmumu.h>

Inheritance diagram for Bsmumu:

Collaboration diagram for Bsmumu:

Detailed Description

Definition at line 15 of file Bsmumu.h.

Public Member Functions
	Bsmumu (const StandardModel &SM_i, int obsFlag)
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...
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...
	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_ew order_ew, double mu)
void	computeObs (orders order, orders_ew order_ew)

Private Attributes
double	absP
double	absS
double	ampSq
double	Amumu
double	argP
double	argS
double	beta
double	chiral
EvolBsmm	evolbsmm
double	mb
double	mBs
double	mmu
double	ms
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...
const StandardModel &	SM
	A reference to an object of StandardMode class. More...

Constructor & Destructor Documentation

Bsmumu::Bsmumu	(	const StandardModel &	SM_i,
		int	obsFlag
	)

constructor

Parameters

Flavour

Definition at line 10 of file Bsmumu.cpp.

: ThObservable(SM_i),
  evolbsmm(8, NDR, NNLO, NLO_ewt4, SM)
{
    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)");
};

Member Function Documentation

void Bsmumu::computeAmpSq ( orders  order, orders_ew  order_ew, double  mu  )

protected

Parameters

order
order_ew

Returns

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

Definition at line 65 of file Bsmumu.cpp.

{  
    if (SM.getMyFlavour()->getHDB1().getCoeffsmumu().getOrder() < order % 3){
        std::stringstream out;
        out << order;
        throw std::runtime_error("Bsmumu::computeAmpSq(): required cofficient of "
                                 "order " + out.str() + " not computed");
    }
    gslpp::vector<gslpp::complex> ** allcoeff = SM.getMyFlavour()->ComputeCoeffsmumu(mu, NDR);

    double alsmu = evolbsmm.alphatilde_s(mu);
    double alemu = evolbsmm.alphatilde_e(mu);
   
    if((order == FULLNLO) && (order_ew == FULLNLO_ew)){
    
    switch(order_ew) {
        case FULLNLO_ew:
        {
            gslpp::complex CC = (*(allcoeff[LO]))(7) /alemu  + (*(allcoeff[NLO]))(7) * alsmu/alemu 
                    + (*(allcoeff[NNLO]))(7) * alsmu * alsmu/alemu + (*(allcoeff[LO_ew ]))(7) /alsmu
                    + (*(allcoeff[NLO_ew]))(7) + (*(allcoeff[NLO_ewt1]))(7) * alemu /alsmu /alsmu 
                    + (*(allcoeff[NLO_ewt2]))(7) * alsmu 
                    + (*(allcoeff[NLO_ewt3]))(7) * alemu /alsmu+ (*(allcoeff[NLO_ewt4]))(7) * alemu;
            absP = CC.abs(); //contains only SM contributions (P, P', S, S' not added)
            argP = CC.arg();
            
            absS = 0.;
            argS = 0.;
           
            phiNP = 0.;
            
            ampSq = absP * absP ;
                  
        }
        break;
        default:
            std::stringstream out;
            out << order;
            throw std::runtime_error("Bsmumu::computeAmpSq(): order " + out.str() + " not implemented");;
        }
    }
        
}

double Bsmumu::computeAmumu

(

orders

order

)

Definition at line 53 of file Bsmumu.cpp.

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

void Bsmumu::computeObs ( orders  order, orders_ew  order_ew  )

protected

Definition at line 38 of file Bsmumu.cpp.

{   
    double mu = SM.getMub();  
    
    mmu = SM.getLeptons(StandardModel::MU).getMass();
    mBs = SM.getMesons(QCD::B_S).getMass();
    mb = SM.getQuarks(QCD::BOTTOM).getMass();
    ms = SM.getQuarks(QCD::STRANGE).getMass();
    chiral = pow(mBs, 2.) / 2. / mmu * mb / (mb + ms);
    beta = sqrt(1. - pow(2. * mmu / mBs, 2.));
    computeAmpSq(order, order_ew, 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);
}

double Bsmumu::computeSmumu

(

orders

order

)

Definition at line 59 of file Bsmumu.cpp.

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

double Bsmumu::computeThValue ( )

virtual

hep-ph/9512380v2

Returns

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

Implements ThObservable.

Definition at line 18 of file Bsmumu.cpp.

{
    computeObs(FULLNLO, FULLNLO_ew);
    double FBs = SM.getMesons(QCD::B_S).getDecayconst();
    
    double coupling = SM.getGF() * SM.getGF() * SM.Mw() * SM.Mw() /M_PI /M_PI ; 
 
    double PRF = pow(coupling, 2.) / M_PI /8. / SM.getMesons(QCD::B_S).computeWidth() * pow(FBs, 2.) * pow(mmu, 2.) * mBs * beta;
    double ys = SM.getMesons(QCD::B_S).getDgamma_gamma()/2.; // For now. To be explicitly calculated.
    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

double Bsmumu::absP

private

Definition at line 53 of file Bsmumu.h.

double Bsmumu::absS

private

Definition at line 55 of file Bsmumu.h.

double Bsmumu::ampSq

private

Definition at line 57 of file Bsmumu.h.

double Bsmumu::Amumu

private

Definition at line 58 of file Bsmumu.h.

double Bsmumu::argP

private

Definition at line 54 of file Bsmumu.h.

double Bsmumu::argS

private

Definition at line 56 of file Bsmumu.h.

double Bsmumu::beta

private

Definition at line 47 of file Bsmumu.h.

double Bsmumu::chiral

private

Definition at line 52 of file Bsmumu.h.

EvolBsmm Bsmumu::evolbsmm

private

Definition at line 63 of file Bsmumu.h.

double Bsmumu::mb

private

Definition at line 50 of file Bsmumu.h.

double Bsmumu::mBs

private

Definition at line 48 of file Bsmumu.h.

double Bsmumu::mmu

private

Definition at line 49 of file Bsmumu.h.

double Bsmumu::ms

private

Definition at line 51 of file Bsmumu.h.

int Bsmumu::obs

private

Definition at line 62 of file Bsmumu.h.

double Bsmumu::phiNP

private

Definition at line 60 of file Bsmumu.h.

double Bsmumu::Smumu

private

Definition at line 59 of file Bsmumu.h.

double Bsmumu::timeInt

private

Definition at line 61 of file Bsmumu.h.

---

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

• Bsmumu.h
• Bsmumu.cpp