Robs_mu_3e Class Reference

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

#include <li_3lj.h>

Inheritance diagram for Robs_mu_3e:

Detailed Description

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

Author

HEPfit Collaboration

Copyright

GNU General Public License

Definition at line 69 of file li_3lj.h.

Public Member Functions
double	computeThValue ()
	Robs_mu_3e (const StandardModel &SM_i)
	Calculates the upper limit ratio of the process \( \mu \to eee \). 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...

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

Robs_mu_3e()

Robs_mu_3e::Robs_mu_3e

(

const StandardModel &

SM_i

)

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

Definition at line 127 of file li_3lj.cpp.

: ThObservable(SM_i)
{}

Member Function Documentation

computeThValue()

double Robs_mu_3e::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 131 of file li_3lj.cpp.

{
    double alph = SM.getAle();
    double mE = SM.getLeptons(StandardModel::ELECTRON).getMass();
    double mMU = SM.getLeptons(StandardModel::MU).getMass();
    double GammaMU = 2.99598e-19;
    gslpp::vector<gslpp::complex> ** allcoeff_m3e = SM.getMyLeptonFlavour()->ComputeCoeffli_3lj(1);
 
    double Brmu3e = alph*alph/(32.0*M_PI) * pow(mMU,5.0)
                    * ((*(allcoeff_m3e[LO]))(1).abs2()+(*(allcoeff_m3e[LO]))(0).abs2()
                       -4.0*((*(allcoeff_m3e[LO]))(1)*(*(allcoeff_m3e[LO]))(2).conjugate()+(*(allcoeff_m3e[LO]))(3)*(*(allcoeff_m3e[LO]))(0).conjugate()).real()
                       +((*(allcoeff_m3e[LO]))(3).abs2()+(*(allcoeff_m3e[LO]))(2).abs2())*(16.0*log(mMU/mE)-22.0)/3.0
                       +((*(allcoeff_m3e[LO]))(5).abs2()+(*(allcoeff_m3e[LO]))(4).abs2())/6.0
                       +(((*(allcoeff_m3e[LO]))(7)+(*(allcoeff_m3e[LO]))(13)).abs2()+((*(allcoeff_m3e[LO]))(6)+(*(allcoeff_m3e[LO]))(12)).abs2())/3.0
                       +(((*(allcoeff_m3e[LO]))(9)+(*(allcoeff_m3e[LO]))(15)).abs2()+((*(allcoeff_m3e[LO]))(8)+(*(allcoeff_m3e[LO]))(14)).abs2())/24.0
                       +((*(allcoeff_m3e[LO]))(11).abs2()+(*(allcoeff_m3e[LO]))(10).abs2())*6.0
                       -(((*(allcoeff_m3e[LO]))(9)+(*(allcoeff_m3e[LO]))(15))*(*(allcoeff_m3e[LO]))(11).conjugate()+((*(allcoeff_m3e[LO]))(8)+(*(allcoeff_m3e[LO]))(14))*(*(allcoeff_m3e[LO]))(10).conjugate()).real()
                       +((*(allcoeff_m3e[LO]))(1)*(*(allcoeff_m3e[LO]))(5).conjugate()+(*(allcoeff_m3e[LO]))(0)*(*(allcoeff_m3e[LO]))(4).conjugate()
                         +(*(allcoeff_m3e[LO]))(1)*((*(allcoeff_m3e[LO]))(7)+(*(allcoeff_m3e[LO]))(13)).conjugate()+(*(allcoeff_m3e[LO]))(0)*((*(allcoeff_m3e[LO]))(6)+(*(allcoeff_m3e[LO]))(12)).conjugate()).real()*2.0/3.0
                       -4.0*((*(allcoeff_m3e[LO]))(2)*(*(allcoeff_m3e[LO]))(5).conjugate()+(*(allcoeff_m3e[LO]))(3)*(*(allcoeff_m3e[LO]))(4).conjugate()
                             +(*(allcoeff_m3e[LO]))(3)*((*(allcoeff_m3e[LO]))(6)+(*(allcoeff_m3e[LO]))(12)).conjugate()+(*(allcoeff_m3e[LO]))(2)*((*(allcoeff_m3e[LO]))(7)+(*(allcoeff_m3e[LO]))(13)).conjugate()).real()/3.0
                       +2.0*((*(allcoeff_m3e[LO]))(19).abs2()+(*(allcoeff_m3e[LO]))(16).abs2())/3.0
                       +((*(allcoeff_m3e[LO]))(18).abs2()+(*(allcoeff_m3e[LO]))(17).abs2())/3.0
                       +((*(allcoeff_m3e[LO]))(5)*(*(allcoeff_m3e[LO]))(19).conjugate()+(*(allcoeff_m3e[LO]))(4)*(*(allcoeff_m3e[LO]))(16).conjugate()
                         +((*(allcoeff_m3e[LO]))(7)+(*(allcoeff_m3e[LO]))(13))*(*(allcoeff_m3e[LO]))(18).conjugate()+((*(allcoeff_m3e[LO]))(6)+(*(allcoeff_m3e[LO]))(12))*(*(allcoeff_m3e[LO]))(17).conjugate()).real()*2.0/3.0
                       +((*(allcoeff_m3e[LO]))(1)*(*(allcoeff_m3e[LO]))(19).conjugate()+(*(allcoeff_m3e[LO]))(0)*(*(allcoeff_m3e[LO]))(16).conjugate()).real()*4.0/3.0
                       +((*(allcoeff_m3e[LO]))(1)*(*(allcoeff_m3e[LO]))(18).conjugate()+(*(allcoeff_m3e[LO]))(0)*(*(allcoeff_m3e[LO]))(17).conjugate()).real()*2.0/3.0
                       -((*(allcoeff_m3e[LO]))(2)*(*(allcoeff_m3e[LO]))(19).conjugate()+(*(allcoeff_m3e[LO]))(3)*(*(allcoeff_m3e[LO]))(16).conjugate()).real()*8.0/3.0
                       -((*(allcoeff_m3e[LO]))(3)*(*(allcoeff_m3e[LO]))(17).conjugate()+(*(allcoeff_m3e[LO]))(2)*(*(allcoeff_m3e[LO]))(18).conjugate()).real()*4.0/3.0) / GammaMU;
    double Brmu3e_obs=1.;  //from ?
    double Brmu3e_exp=1.;  //from ?
    double ntos=1.64485362695147;
    
    return (1+(Brmu3e-Brmu3e_obs)/Brmu3e_exp)*ntos;
}

---

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

• li_3lj.h
• li_3lj.cpp