a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
li_lj_gamma.cpp
Go to the documentation of this file.
1/*
2 * Copyright (C) 2012 HEPfit Collaboration
3 *
4 *
5 * For the licensing terms see doc/COPYING.
6 */
7
8#include "li_lj_gamma.h"
9#include "LeptonFlavour.h"
10
11mu_e_gamma::mu_e_gamma(const StandardModel& SM_i)
12: ThObservable(SM_i)
13{}
14
15double mu_e_gamma::computeThValue()
16{
17 double alph = SM.getAle();
18 double mMU = SM.getLeptons(StandardModel::MU).getMass();
19 gslpp::vector<gslpp::complex> ** allcoeff_me = SM.getMyLeptonFlavour()->ComputeCoeffli_lj_gamma(1);
20 double BR_mu_e_gamma = (alph*pow(mMU,5) * ((*(allcoeff_me[LO])) * (*(allcoeff_me[LO])).conjugate()).abs() / (2.99598e-19) );
21 return BR_mu_e_gamma;
22}
23
24log_meg::log_meg(const StandardModel& SM_i)
25: ThObservable(SM_i)
26{}
27
28double log_meg::computeThValue()
29{
30 double alph = SM.getAle();
31 double mMU = SM.getLeptons(StandardModel::MU).getMass();
32 gslpp::vector<gslpp::complex> ** allcoeff_me = SM.getMyLeptonFlavour()->ComputeCoeffli_lj_gamma(1);
33 double BR_mu_e_gamma = (alph*pow(mMU,5) * ((*(allcoeff_me[LO])) * (*(allcoeff_me[LO])).conjugate()).abs() / (2.99598e-19) );
34 return log10(BR_mu_e_gamma);
35}
36
37Robs_mu_e_gamma::Robs_mu_e_gamma(const StandardModel& SM_i)
38: ThObservable(SM_i)
39{}
40
41double Robs_mu_e_gamma::computeThValue()
42{
43 double alph = SM.getAle();
44 double mMU = SM.getLeptons(StandardModel::MU).getMass();
45 gslpp::vector<gslpp::complex> ** allcoeff_me = SM.getMyLeptonFlavour()->ComputeCoeffli_lj_gamma(1);
46 double BR_mu_e_gamma = (alph*pow(mMU,5) * ((*(allcoeff_me[LO])) * (*(allcoeff_me[LO])).conjugate()).abs() / (2.99598e-19) );
47 double BR_mu_e_gamma_obs=4.2e-13; //from 1605.05081v1 (TheMEG:2016wtm)
48 double BR_mu_e_gamma_exp=5.3e-13; //from 1605.05081v1 (TheMEG:2016wtm)
49 double ntos=1.64485362695147; //90% CL
50
51 return (1+(BR_mu_e_gamma-BR_mu_e_gamma_obs)/BR_mu_e_gamma_exp)*ntos;
52}
53
54tau_mu_gamma::tau_mu_gamma(const StandardModel& SM_i)
55: ThObservable(SM_i)
56{}
57
58double tau_mu_gamma::computeThValue()
59{
60 double alph = SM.getAle();
61 double mTAU = SM.getLeptons(StandardModel::TAU).getMass();
62 gslpp::vector<gslpp::complex> ** allcoeff_tm = SM.getMyLeptonFlavour()->ComputeCoeffli_lj_gamma(2);
63 double BR_tau_mu_gamma = (alph*pow(mTAU,5) * ((*(allcoeff_tm[LO])) * (*(allcoeff_tm[LO])).conjugate()).abs() / (2.26735e-12) );
64 return BR_tau_mu_gamma;
65}
66
67log_tmg::log_tmg(const StandardModel& SM_i)
68: ThObservable(SM_i)
69{}
70
71double log_tmg::computeThValue()
72{
73 double alph = SM.getAle();
74 double mTAU = SM.getLeptons(StandardModel::TAU).getMass();
75 gslpp::vector<gslpp::complex> ** allcoeff_tm = SM.getMyLeptonFlavour()->ComputeCoeffli_lj_gamma(2);
76 double BR_tau_mu_gamma = (alph*pow(mTAU,5) * ((*(allcoeff_tm[LO])) * (*(allcoeff_tm[LO])).conjugate()).abs() / (2.26735e-12) );
77 return log10(BR_tau_mu_gamma);
78}
79
80Robs_tau_mu_gamma::Robs_tau_mu_gamma(const StandardModel& SM_i)
81: ThObservable(SM_i)
82{}
83
84double Robs_tau_mu_gamma::computeThValue()
85{
86 double alph = SM.getAle();
87 double mTAU = SM.getLeptons(StandardModel::TAU).getMass();
88 gslpp::vector<gslpp::complex> ** allcoeff_tm = SM.getMyLeptonFlavour()->ComputeCoeffli_lj_gamma(2);
89 double BR_tau_mu_gamma = (alph*pow(mTAU,5) * ((*(allcoeff_tm[LO])) * (*(allcoeff_tm[LO])).conjugate()).abs() / (2.26735e-12) );
90 double BR_tau_mu_gamma_obs=4.4e-8; //from Phys.Rev.Lett. 104 (2010) 021802 (Aubert:2009ag)
91 double BR_tau_mu_gamma_exp=8.2e-8; //from Phys.Rev.Lett. 104 (2010) 021802 (Aubert:2009ag)
92 double ntos=1.64485362695147; //90% CL
93
94 return (1+(BR_tau_mu_gamma-BR_tau_mu_gamma_obs)/BR_tau_mu_gamma_exp)*ntos;
95}
96
97Robs_tau_mu_gamma_BelleII::Robs_tau_mu_gamma_BelleII(const StandardModel& SM_i)
98: ThObservable(SM_i)
99{}
100
101double Robs_tau_mu_gamma_BelleII::computeThValue()
102{
103 double alph = SM.getAle();
104 double mTAU = SM.getLeptons(StandardModel::TAU).getMass();
105 gslpp::vector<gslpp::complex> ** allcoeff_tm = SM.getMyLeptonFlavour()->ComputeCoeffli_lj_gamma(2);
106 double BR_tau_mu_gamma = (alph*pow(mTAU,5) * ((*(allcoeff_tm[LO])) * (*(allcoeff_tm[LO])).conjugate()).abs() / (2.26735e-12) );
107 double BR_tau_mu_gamma_obs=1.0e-9; //1002.5012
108 double BR_tau_mu_gamma_exp=1.0e-9; //1002.5012
109 double ntos=1.64485362695147; //90% CL
110
111 return (1+(BR_tau_mu_gamma-BR_tau_mu_gamma_obs)/BR_tau_mu_gamma_exp)*ntos;
112}
113
114tau_e_gamma::tau_e_gamma(const StandardModel& SM_i)
115: ThObservable(SM_i)
116{}
117
118double tau_e_gamma::computeThValue()
119{
120 double alph = SM.getAle();
121 double mTAU = SM.getLeptons(StandardModel::TAU).getMass();
122 gslpp::vector<gslpp::complex> ** allcoeff_te = SM.getMyLeptonFlavour()->ComputeCoeffli_lj_gamma(3);
123 double BR_tau_e_gamma = (alph*pow(mTAU,5) * ((*(allcoeff_te[LO])) * (*(allcoeff_te[LO])).conjugate()).abs() / (2.26735e-12) );
124 return BR_tau_e_gamma;
125}
126
127log_teg::log_teg(const StandardModel& SM_i)
128: ThObservable(SM_i)
129{}
130
131double log_teg::computeThValue()
132{
133 double alph = SM.getAle();
134 double mTAU = SM.getLeptons(StandardModel::TAU).getMass();
135 gslpp::vector<gslpp::complex> ** allcoeff_te = SM.getMyLeptonFlavour()->ComputeCoeffli_lj_gamma(3);
136 double BR_tau_e_gamma = (alph*pow(mTAU,5) * ((*(allcoeff_te[LO])) * (*(allcoeff_te[LO])).conjugate()).abs() / (2.26735e-12) );
137 return log10(BR_tau_e_gamma);
138}
139
140Robs_tau_e_gamma::Robs_tau_e_gamma(const StandardModel& SM_i)
141: ThObservable(SM_i)
142{}
143
144double Robs_tau_e_gamma::computeThValue()
145{
146 double alph = SM.getAle();
147 double mTAU = SM.getLeptons(StandardModel::TAU).getMass();
148 gslpp::vector<gslpp::complex> ** allcoeff_te = SM.getMyLeptonFlavour()->ComputeCoeffli_lj_gamma(3);
149 double BR_tau_e_gamma = (alph*pow(mTAU,5) * ((*(allcoeff_te[LO])) * (*(allcoeff_te[LO])).conjugate()).abs() / (2.26735e-12) );
150 double BR_tau_e_gamma_obs=3.3e-8; //from Phys.Rev.Lett. 104 (2010) 021802 (Aubert:2009ag)
151 double BR_tau_e_gamma_exp=9.8e-8; //from Phys.Rev.Lett. 104 (2010) 021802 (Aubert:2009ag)
152 double ntos=1.64485362695147; //90% CL
153
154 return (1+(BR_tau_e_gamma-BR_tau_e_gamma_obs)/BR_tau_e_gamma_exp)*ntos;
155}
LO
@ LO
Definition: OrderScheme.h:34
LeptonFlavour::ComputeCoeffli_lj_gamma
gslpp::vector< gslpp::complex > ** ComputeCoeffli_lj_gamma(int li_lj)
Computes the Wilson coefficient for the process .
Definition: LeptonFlavour.h:52
Particle::getMass
const double & getMass() const
A get method to access the particle mass.
Definition: Particle.h:61
QCD::MU
@ MU
Definition: QCD.h:314
QCD::TAU
@ TAU
Definition: QCD.h:316
Robs_mu_e_gamma::Robs_mu_e_gamma
Robs_mu_e_gamma(const StandardModel &SM_i)
Calculates the upper limit ratio of the process .
Definition: li_lj_gamma.cpp:37
Robs_mu_e_gamma::computeThValue
double computeThValue()
Definition: li_lj_gamma.cpp:41
Robs_tau_e_gamma::computeThValue
double computeThValue()
Definition: li_lj_gamma.cpp:144
Robs_tau_e_gamma::Robs_tau_e_gamma
Robs_tau_e_gamma(const StandardModel &SM_i)
Calculates the upper limit ratio of the process .
Definition: li_lj_gamma.cpp:140
Robs_tau_mu_gamma_BelleII::computeThValue
double computeThValue()
Definition: li_lj_gamma.cpp:101
Robs_tau_mu_gamma_BelleII::Robs_tau_mu_gamma_BelleII
Robs_tau_mu_gamma_BelleII(const StandardModel &SM_i)
Calculates the upper limit ratio of the process .
Definition: li_lj_gamma.cpp:97
Robs_tau_mu_gamma::computeThValue
double computeThValue()
Definition: li_lj_gamma.cpp:84
Robs_tau_mu_gamma::Robs_tau_mu_gamma
Robs_tau_mu_gamma(const StandardModel &SM_i)
Calculates the upper limit ratio of the process .
Definition: li_lj_gamma.cpp:80
StandardModel
A model class for the Standard Model.
Definition: StandardModel.h:509
StandardModel::getLeptons
const Particle & getLeptons(const QCD::lepton p) const
A get method to retrieve the member object of a lepton.
Definition: StandardModel.h:744
StandardModel::getMyLeptonFlavour
LeptonFlavour * getMyLeptonFlavour() const
Definition: StandardModel.h:1038
StandardModel::getAle
const double getAle() const
A get method to retrieve the fine-structure constant .
Definition: StandardModel.h:789
ThObservable
A class for a model prediction of an observable.
Definition: ThObservable.h:25
ThObservable::SM
const StandardModel & SM
A reference to an object of StandardMode class.
Definition: ThObservable.h:121
log_meg::computeThValue
double computeThValue()
Definition: li_lj_gamma.cpp:28
log_meg::log_meg
log_meg(const StandardModel &SM_i)
Calculates the value of the process .
Definition: li_lj_gamma.cpp:24
log_teg::computeThValue
double computeThValue()
Definition: li_lj_gamma.cpp:131
log_teg::log_teg
log_teg(const StandardModel &SM_i)
Calculates the value of the process .
Definition: li_lj_gamma.cpp:127
log_tmg::log_tmg
log_tmg(const StandardModel &SM_i)
Calculates the value of the process .
Definition: li_lj_gamma.cpp:67
log_tmg::computeThValue
double computeThValue()
Definition: li_lj_gamma.cpp:71
mu_e_gamma::mu_e_gamma
mu_e_gamma(const StandardModel &SM_i)
Calculates the value of the process .
Definition: li_lj_gamma.cpp:11
mu_e_gamma::computeThValue
double computeThValue()
Definition: li_lj_gamma.cpp:15
tau_e_gamma::computeThValue
double computeThValue()
Definition: li_lj_gamma.cpp:118
tau_e_gamma::tau_e_gamma
tau_e_gamma(const StandardModel &SM_i)
Calculates the value of the process .
Definition: li_lj_gamma.cpp:114
tau_mu_gamma::tau_mu_gamma
tau_mu_gamma(const StandardModel &SM_i)
Calculates the value of the process .
Definition: li_lj_gamma.cpp:54
tau_mu_gamma::computeThValue
double computeThValue()
Definition: li_lj_gamma.cpp:58
li_lj_gamma.h