a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
Robs_tau_3mu Class Reference

A class for calculating the upper limit ratio for the process \( \tau \to \mu \mu \mu \). More...

#include <li_3lj.h>

+ Inheritance diagram for Robs_tau_3mu:

Detailed Description

A class for calculating the upper limit ratio for the process \( \tau \to \mu \mu \mu \).

Author
HEPfit Collaboration

Definition at line 121 of file li_3lj.h.

Public Member Functions

double computeThValue ()
 
 Robs_tau_3mu (const StandardModel &SM_i)
 Calculates the upper limit ratio of the process \( \tau \to \mu \mu \mu \). 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 StandardModelgetModel ()
 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...
 

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 StandardModelSM
 A reference to an object of StandardMode class. More...
 

Constructor & Destructor Documentation

◆ Robs_tau_3mu()

Robs_tau_3mu::Robs_tau_3mu ( const StandardModel SM_i)

Calculates the upper limit ratio of the process \( \tau \to \mu \mu \mu \).

Definition at line 214 of file li_3lj.cpp.

215 : ThObservable(SM_i)
216 {}

Member Function Documentation

◆ computeThValue()

double Robs_tau_3mu::computeThValue ( )
virtual
Returns
\(R=1+(BR_{\text{theo}}-BR_{\text{obs,90\%}})/BR_{\text{exp,90\%}}\cdot \sqrt{2}\text{InvErf}(0.90)\)

Implements ThObservable.

Definition at line 218 of file li_3lj.cpp.

219 {
220  double alph = SM.getAle();
221  double mMU = SM.getLeptons(StandardModel::MU).getMass();
222  double mTAU = SM.getLeptons(StandardModel::TAU).getMass();
223  double GammaTAU = 2.26735e-12;
224 
226 
227  double Brtau3mu = alph*alph/(32.0*M_PI) * pow(mTAU,5.0)
228  * ((*(allcoeff_t3m[LO]))(1).abs2()+(*(allcoeff_t3m[LO]))(0).abs2()
229  -4.0*((*(allcoeff_t3m[LO]))(1)*(*(allcoeff_t3m[LO]))(2).conjugate()+(*(allcoeff_t3m[LO]))(3)*(*(allcoeff_t3m[LO]))(0).conjugate()).real()
230  +((*(allcoeff_t3m[LO]))(3).abs2()+(*(allcoeff_t3m[LO]))(2).abs2())*(16.0*log(mTAU/mMU)-22.0)/3.0
231  +((*(allcoeff_t3m[LO]))(5).abs2()+(*(allcoeff_t3m[LO]))(4).abs2())/6.0
232  +(((*(allcoeff_t3m[LO]))(7)+(*(allcoeff_t3m[LO]))(13)).abs2()+((*(allcoeff_t3m[LO]))(6)+(*(allcoeff_t3m[LO]))(12)).abs2())/3.0
233  +(((*(allcoeff_t3m[LO]))(9)+(*(allcoeff_t3m[LO]))(15)).abs2()+((*(allcoeff_t3m[LO]))(8)+(*(allcoeff_t3m[LO]))(14)).abs2())/24.0
234  +((*(allcoeff_t3m[LO]))(11).abs2()+(*(allcoeff_t3m[LO]))(10).abs2())*6.0
235  -(((*(allcoeff_t3m[LO]))(9)+(*(allcoeff_t3m[LO]))(15))*(*(allcoeff_t3m[LO]))(11).conjugate()+((*(allcoeff_t3m[LO]))(8)+(*(allcoeff_t3m[LO]))(14))*(*(allcoeff_t3m[LO]))(10).conjugate()).real()
236  +((*(allcoeff_t3m[LO]))(1)*(*(allcoeff_t3m[LO]))(5).conjugate()+(*(allcoeff_t3m[LO]))(0)*(*(allcoeff_t3m[LO]))(4).conjugate()
237  +(*(allcoeff_t3m[LO]))(1)*((*(allcoeff_t3m[LO]))(7)+(*(allcoeff_t3m[LO]))(13)).conjugate()+(*(allcoeff_t3m[LO]))(0)*((*(allcoeff_t3m[LO]))(6)+(*(allcoeff_t3m[LO]))(12)).conjugate()).real()*2.0/3.0
238  -4.0*((*(allcoeff_t3m[LO]))(2)*(*(allcoeff_t3m[LO]))(5).conjugate()+(*(allcoeff_t3m[LO]))(3)*(*(allcoeff_t3m[LO]))(4).conjugate()
239  +(*(allcoeff_t3m[LO]))(3)*((*(allcoeff_t3m[LO]))(6)+(*(allcoeff_t3m[LO]))(12)).conjugate()+(*(allcoeff_t3m[LO]))(2)*((*(allcoeff_t3m[LO]))(7)+(*(allcoeff_t3m[LO]))(13)).conjugate()).real()/3.0
240  +2.0*((*(allcoeff_t3m[LO]))(19).abs2()+(*(allcoeff_t3m[LO]))(16).abs2())/3.0
241  +((*(allcoeff_t3m[LO]))(18).abs2()+(*(allcoeff_t3m[LO]))(17).abs2())/3.0
242  +((*(allcoeff_t3m[LO]))(5)*(*(allcoeff_t3m[LO]))(19).conjugate()+(*(allcoeff_t3m[LO]))(4)*(*(allcoeff_t3m[LO]))(16).conjugate()
243  +((*(allcoeff_t3m[LO]))(7)+(*(allcoeff_t3m[LO]))(13))*(*(allcoeff_t3m[LO]))(18).conjugate()+((*(allcoeff_t3m[LO]))(6)+(*(allcoeff_t3m[LO]))(12))*(*(allcoeff_t3m[LO]))(17).conjugate()).real()*2.0/3.0
244  +((*(allcoeff_t3m[LO]))(1)*(*(allcoeff_t3m[LO]))(19).conjugate()+(*(allcoeff_t3m[LO]))(0)*(*(allcoeff_t3m[LO]))(16).conjugate()).real()*4.0/3.0
245  +((*(allcoeff_t3m[LO]))(1)*(*(allcoeff_t3m[LO]))(18).conjugate()+(*(allcoeff_t3m[LO]))(0)*(*(allcoeff_t3m[LO]))(17).conjugate()).real()*2.0/3.0
246  -((*(allcoeff_t3m[LO]))(2)*(*(allcoeff_t3m[LO]))(19).conjugate()+(*(allcoeff_t3m[LO]))(3)*(*(allcoeff_t3m[LO]))(16).conjugate()).real()*8.0/3.0
247  -((*(allcoeff_t3m[LO]))(3)*(*(allcoeff_t3m[LO]))(17).conjugate()+(*(allcoeff_t3m[LO]))(2)*(*(allcoeff_t3m[LO]))(18).conjugate()).real()*4.0/3.0) / GammaTAU;
248  double Brtau3mu_obs=4.6e-8; //from JHEP 02 (2015) 121
249  double Brtau3mu_exp=5.0e-8; //from JHEP 02 (2015) 121
250  double ntos=1.64485362695147;
251 
252  return (1+(Brtau3mu-Brtau3mu_obs)/Brtau3mu_exp)*ntos;
253 }

The documentation for this class was generated from the following files:
QCD::TAU
Definition: QCD.h:316
StandardModel::getMyLeptonFlavour
LeptonFlavour * getMyLeptonFlavour() const
Definition: StandardModel.h:1011
LO
Definition: OrderScheme.h:33
gslpp::log
complex log(const complex &z)
Definition: gslpp_complex.cpp:342
ThObservable::ThObservable
ThObservable(const StandardModel &SM_i)
Constructor.
Definition: ThObservable.h:32
Particle::getMass
const double & getMass() const
A get method to access the particle mass.
Definition: Particle.h:61
LeptonFlavour::ComputeCoeffli_3lj
gslpp::vector< gslpp::complex > ** ComputeCoeffli_3lj(int li_lj)
Computes the Wilson coefficient for the process .
Definition: LeptonFlavour.h:71
gslpp::pow
complex pow(const complex &z1, const complex &z2)
Definition: gslpp_complex.cpp:395
ThObservable::SM
const StandardModel & SM
A reference to an object of StandardMode class.
Definition: ThObservable.h:121
StandardModel::getAle
double getAle() const
A get method to retrieve the fine-structure constant .
Definition: StandardModel.h:745
gslpp::vector< gslpp::complex >
QCD::MU
Definition: QCD.h:314
StandardModel::getLeptons
Particle getLeptons(const QCD::lepton p) const
A get method to retrieve the member object of a lepton.
Definition: StandardModel.h:709