192 double GammaTAU = 2.26735e-12;
196 double Brtau3e = alph*alph/(32.0*M_PI) *
pow(mTAU,5.0)
197 * ((*(allcoeff_t3e[
LO]))(1).abs2()+(*(allcoeff_t3e[
LO]))(0).abs2()
198 -4.0*((*(allcoeff_t3e[
LO]))(1)*(*(allcoeff_t3e[
LO]))(2).conjugate()+(*(allcoeff_t3e[
LO]))(3)*(*(allcoeff_t3e[LO]))(0).conjugate()).real()
199 +((*(allcoeff_t3e[
LO]))(3).abs2()+(*(allcoeff_t3e[
LO]))(2).abs2())*(16.0*
log(mTAU/mE)-22.0)/3.0
200 +((*(allcoeff_t3e[
LO]))(5).abs2()+(*(allcoeff_t3e[
LO]))(4).abs2())/6.0
201 +(((*(allcoeff_t3e[LO]))(7)+(*(allcoeff_t3e[
LO]))(13)).abs2()+((*(allcoeff_t3e[
LO]))(6)+(*(allcoeff_t3e[LO]))(12)).abs2())/3.0
202 +(((*(allcoeff_t3e[LO]))(9)+(*(allcoeff_t3e[
LO]))(15)).abs2()+((*(allcoeff_t3e[
LO]))(8)+(*(allcoeff_t3e[LO]))(14)).abs2())/24.0
203 +((*(allcoeff_t3e[LO]))(11).abs2()+(*(allcoeff_t3e[
LO]))(10).abs2())*6.0
204 -(((*(allcoeff_t3e[LO]))(9)+(*(allcoeff_t3e[
LO]))(15))*(*(allcoeff_t3e[
LO]))(11).conjugate()+((*(allcoeff_t3e[
LO]))(8)+(*(allcoeff_t3e[LO]))(14))*(*(allcoeff_t3e[LO]))(10).conjugate()).real()
205 +((*(allcoeff_t3e[
LO]))(1)*(*(allcoeff_t3e[LO]))(5).conjugate()+(*(allcoeff_t3e[
LO]))(0)*(*(allcoeff_t3e[LO]))(4).conjugate()
206 +(*(allcoeff_t3e[
LO]))(1)*((*(allcoeff_t3e[LO]))(7)+(*(allcoeff_t3e[
LO]))(13)).conjugate()+(*(allcoeff_t3e[
LO]))(0)*((*(allcoeff_t3e[LO]))(6)+(*(allcoeff_t3e[
LO]))(12)).conjugate()).real()*2.0/3.0
207 -4.0*((*(allcoeff_t3e[
LO]))(2)*(*(allcoeff_t3e[LO]))(5).conjugate()+(*(allcoeff_t3e[
LO]))(3)*(*(allcoeff_t3e[LO]))(4).conjugate()
208 +(*(allcoeff_t3e[
LO]))(3)*((*(allcoeff_t3e[LO]))(6)+(*(allcoeff_t3e[
LO]))(12)).conjugate()+(*(allcoeff_t3e[
LO]))(2)*((*(allcoeff_t3e[LO]))(7)+(*(allcoeff_t3e[
LO]))(13)).conjugate()).real()/3.0
209 +2.0*((*(allcoeff_t3e[
LO]))(19).abs2()+(*(allcoeff_t3e[
LO]))(16).abs2())/3.0
210 +((*(allcoeff_t3e[LO]))(18).abs2()+(*(allcoeff_t3e[
LO]))(17).abs2())/3.0
211 +((*(allcoeff_t3e[LO]))(5)*(*(allcoeff_t3e[
LO]))(19).conjugate()+(*(allcoeff_t3e[
LO]))(4)*(*(allcoeff_t3e[LO]))(16).conjugate()
212 +((*(allcoeff_t3e[
LO]))(7)+(*(allcoeff_t3e[LO]))(13))*(*(allcoeff_t3e[LO]))(18).conjugate()+((*(allcoeff_t3e[
LO]))(6)+(*(allcoeff_t3e[LO]))(12))*(*(allcoeff_t3e[LO]))(17).conjugate()).real()*2.0/3.0
213 +((*(allcoeff_t3e[
LO]))(1)*(*(allcoeff_t3e[LO]))(19).conjugate()+(*(allcoeff_t3e[
LO]))(0)*(*(allcoeff_t3e[LO]))(16).conjugate()).real()*4.0/3.0
214 +((*(allcoeff_t3e[
LO]))(1)*(*(allcoeff_t3e[LO]))(18).conjugate()+(*(allcoeff_t3e[
LO]))(0)*(*(allcoeff_t3e[LO]))(17).conjugate()).real()*2.0/3.0
215 -((*(allcoeff_t3e[
LO]))(2)*(*(allcoeff_t3e[LO]))(19).conjugate()+(*(allcoeff_t3e[
LO]))(3)*(*(allcoeff_t3e[LO]))(16).conjugate()).real()*8.0/3.0
216 -((*(allcoeff_t3e[
LO]))(3)*(*(allcoeff_t3e[LO]))(17).conjugate()+(*(allcoeff_t3e[
LO]))(2)*(*(allcoeff_t3e[LO]))(18).conjugate()).real()*4.0/3.0) / GammaTAU;
Particle getLeptons(const StandardModel::lepton p) const
A get method to retrieve the member object of a lepton.
complex pow(const complex &z1, const complex &z2)
gslpp::vector< gslpp::complex > ** ComputeCoeffli_3lj(int li_lj)
Computes the Wilson coefficient for the process .
LeptonFlavour * getMyLeptonFlavour() const
const StandardModel & mySM
Constructor containing the Wilson coefficient.
double getAle() const
A get method to retrieve the fine-structure constant .
const double & getMass() const
A get method to access the particle mass.
complex log(const complex &z)