a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
Public Types | Public Member Functions | Private Member Functions | Private Attributes | List of all members
LEP2GIMR Class Reference
EW

A class for NP analyses of LEP-II observables with the dimension 6 operators in the GIMR basis. More...

#include <LEP2GIMR.h>

Detailed Description

A class for NP analyses of LEP-II observables with the dimension 6 operators in the GIMR basis.

Author
HEPfit Collaboration
Copyright
GNU General Public License

Definition at line 26 of file LEP2GIMR.h.

Public Types

enum  Param { C_LL =0 , C_LR , C_RL , C_RR , delta_GammaZ , delta_gLf , delta_gRf , delta_gLe , delta_gRe , delta_Mz2 }
 

Public Member Functions

 LEP2GIMR (const StandardModel &SM_i)
 LEP2GIMR constructor. More...
 
double sigma_l_LEP2_GIMR (const QCD::lepton l, const double s, const double GIMRParam_i[]) const
 
double sigma_q_LEP2_GIMR (const QCD::quark q, const double s, const double GIMRParam_i[]) const
 
double sigmaFminusB_l_LEP2_GIMR (const QCD::lepton l, const double s, const double GIMRParam_i[]) const
 
double sigmaFminusB_q_LEP2_GIMR (const QCD::quark q, const double s, const double GIMRParam_i[]) const
 

Private Member Functions

double deltaA1l (const QCD::lepton l, const double GIMRParam_i[]) const
 
double deltaA1q (const QCD::quark q, const double GIMRParam_i[]) const
 
double deltaA2l (const QCD::lepton l, const double GIMRParam_i[]) const
 
double deltaA2q (const QCD::quark q, const double GIMRParam_i[]) const
 
double deltaB1l (const QCD::lepton l, const double GIMRParam_i[]) const
 
double deltaB1q (const QCD::quark q, const double GIMRParam_i[]) const
 
double deltaB2l (const QCD::lepton l, const double GIMRParam_i[]) const
 
double deltaB2q (const QCD::quark q, const double GIMRParam_i[]) const
 
double gL_l (const QCD::lepton l) const
 
double gL_q (const QCD::quark q) const
 
double gR_l (const QCD::lepton l) const
 
double gR_q (const QCD::quark q) const
 

Private Attributes

const StandardModel & SM
 

Member Enumeration Documentation

◆ Param

enum LEP2GIMR::Param
Enumerator
C_LL 
C_LR 
C_RL 
C_RR 
delta_GammaZ 
delta_gLf 
delta_gRf 
delta_gLe 
delta_gRe 
delta_Mz2 

Definition at line 29 of file LEP2GIMR.h.

29{C_LL=0, C_LR, C_RL, C_RR, delta_GammaZ, delta_gLf, delta_gRf, delta_gLe, delta_gRe, delta_Mz2};
LEP2GIMR::delta_gLf
@ delta_gLf
Definition: LEP2GIMR.h:29
LEP2GIMR::C_RL
@ C_RL
Definition: LEP2GIMR.h:29
LEP2GIMR::C_LR
@ C_LR
Definition: LEP2GIMR.h:29
LEP2GIMR::C_RR
@ C_RR
Definition: LEP2GIMR.h:29
LEP2GIMR::delta_GammaZ
@ delta_GammaZ
Definition: LEP2GIMR.h:29
LEP2GIMR::delta_Mz2
@ delta_Mz2
Definition: LEP2GIMR.h:29
LEP2GIMR::C_LL
@ C_LL
Definition: LEP2GIMR.h:29
LEP2GIMR::delta_gLe
@ delta_gLe
Definition: LEP2GIMR.h:29
LEP2GIMR::delta_gRe
@ delta_gRe
Definition: LEP2GIMR.h:29
LEP2GIMR::delta_gRf
@ delta_gRf
Definition: LEP2GIMR.h:29

Constructor & Destructor Documentation

◆ LEP2GIMR()

LEP2GIMR::LEP2GIMR ( const StandardModel &  SM_i)

LEP2GIMR constructor.

Parameters
[in]SM_ia reference to an object of type StandardModel

Definition at line 11 of file LEP2GIMR.cpp.

12: SM(SM_i)
13{
14}
LEP2GIMR::SM
const StandardModel & SM
Definition: LEP2GIMR.h:48

Member Function Documentation

◆ deltaA1l()

double LEP2GIMR::deltaA1l ( const QCD::lepton  l,
const double  GIMRParam_i[] 
) const
private

Definition at line 305 of file LEP2GIMR.cpp.

306{
307 double A1l;
308 double deltaGLf = GIMRParam_i[delta_gLf];
309 double deltaGRf = GIMRParam_i[delta_gRf];
310 double deltaGLe = GIMRParam_i[delta_gLe];
311 double deltaGRe = GIMRParam_i[delta_gRe];
312 A1l = deltaGLe * gL_l(l) + deltaGRe * gR_l(l)
313 + gL_l(l) * deltaGLf + gR_l(l) * deltaGRf;
314
315 return A1l;
316}
LEP2GIMR::gR_l
double gR_l(const QCD::lepton l) const
Definition: LEP2GIMR.cpp:211
LEP2GIMR::gL_l
double gL_l(const QCD::lepton l) const
Definition: LEP2GIMR.cpp:201

◆ deltaA1q()

double LEP2GIMR::deltaA1q ( const QCD::quark  q,
const double  GIMRParam_i[] 
) const
private

Definition at line 241 of file LEP2GIMR.cpp.

242{
243 double A1q;
244 double deltaGLf = GIMRParam_i[delta_gLf];
245 double deltaGRf = GIMRParam_i[delta_gRf];
246 double deltaGLe = GIMRParam_i[delta_gLe];
247 double deltaGRe = GIMRParam_i[delta_gRe];
248 A1q = deltaGLe * gL_q(q) + deltaGRe * gR_q(q)
249 + gL_l(StandardModel::ELECTRON) * deltaGLf
250 + gR_l(StandardModel::ELECTRON) * deltaGRf;
251
252 return A1q;
253}
LEP2GIMR::gL_q
double gL_q(const QCD::quark q) const
Definition: LEP2GIMR.cpp:221
LEP2GIMR::gR_q
double gR_q(const QCD::quark q) const
Definition: LEP2GIMR.cpp:231
QCD::ELECTRON
@ ELECTRON
Definition: QCD.h:312

◆ deltaA2l()

double LEP2GIMR::deltaA2l ( const QCD::lepton  l,
const double  GIMRParam_i[] 
) const
private

Definition at line 318 of file LEP2GIMR.cpp.

319{
320 double A2l;
321 double deltaGLf = GIMRParam_i[delta_gLf];
322 double deltaGRf = GIMRParam_i[delta_gRf];
323 double deltaGLe = GIMRParam_i[delta_gLe];
324 double deltaGRe = GIMRParam_i[delta_gRe];
325 A2l = deltaGLe * gR_l(l) + deltaGRe * gL_l(l)
326 + gL_l(l) * deltaGRf + gR_l(l) * deltaGLf;
327
328 return A2l;
329}

◆ deltaA2q()

double LEP2GIMR::deltaA2q ( const QCD::quark  q,
const double  GIMRParam_i[] 
) const
private

Definition at line 255 of file LEP2GIMR.cpp.

256{
257 double A2q;
258 double deltaGLf = GIMRParam_i[delta_gLf];
259 double deltaGRf = GIMRParam_i[delta_gRf];
260 double deltaGLe = GIMRParam_i[delta_gLe];
261 double deltaGRe = GIMRParam_i[delta_gRe];
262 A2q = deltaGLe * gR_q(q) + deltaGRe * gL_q(q)
263 + gL_l(StandardModel::ELECTRON) * deltaGRf
264 + gR_l(StandardModel::ELECTRON) * deltaGLf;
265
266 return A2q;
267}

◆ deltaB1l()

double LEP2GIMR::deltaB1l ( const QCD::lepton  l,
const double  GIMRParam_i[] 
) const
private

Definition at line 331 of file LEP2GIMR.cpp.

332{
333 double B1l;
334 double deltaGLf = GIMRParam_i[delta_gLf];
335 double deltaGRf = GIMRParam_i[delta_gRf];
336 double deltaGLe = GIMRParam_i[delta_gLe];
337 double deltaGRe = GIMRParam_i[delta_gRe];
338 double gL_l2 = gL_l(l) * gL_l(l);
339 double gR_l2 = gR_l(l) * gR_l(l);
340 B1l = gL_l(l) * gL_l2 * deltaGLe + gR_l(l) * gR_l2 * deltaGRe
341 + gL_l2 * gL_l(l) * deltaGLf + gR_l2 * gR_l(l) * deltaGRf;
342
343 return B1l;
344}

◆ deltaB1q()

double LEP2GIMR::deltaB1q ( const QCD::quark  q,
const double  GIMRParam_i[] 
) const
private

Definition at line 269 of file LEP2GIMR.cpp.

270{
271 double B1q;
272 double deltaGLf = GIMRParam_i[delta_gLf];
273 double deltaGRf = GIMRParam_i[delta_gRf];
274 double deltaGLe = GIMRParam_i[delta_gLe];
275 double deltaGRe = GIMRParam_i[delta_gRe];
276 double gL_q2 = gL_q(q) * gL_q(q);
277 double gR_q2 = gR_q(q) * gR_q(q);
278 double gL_l2 = gL_l(StandardModel::ELECTRON) * gL_l(StandardModel::ELECTRON);
279 double gR_l2 = gR_l(StandardModel::ELECTRON) * gR_l(StandardModel::ELECTRON);
280 B1q = gL_l(StandardModel::ELECTRON) * gL_q2 * deltaGLe
281 + gR_l(StandardModel::ELECTRON) * gR_q2 * deltaGRe
282 + gL_l2 * gL_q(q) * deltaGLf + gR_l2 * gR_q(q) * deltaGRf;
283
284 return B1q;
285}

◆ deltaB2l()

double LEP2GIMR::deltaB2l ( const QCD::lepton  l,
const double  GIMRParam_i[] 
) const
private

Definition at line 346 of file LEP2GIMR.cpp.

347{
348 double B2l;
349 double deltaGLf = GIMRParam_i[delta_gLf];
350 double deltaGRf = GIMRParam_i[delta_gRf];
351 double deltaGLe = GIMRParam_i[delta_gLe];
352 double deltaGRe = GIMRParam_i[delta_gRe];
353 double gL_l2 = gL_l(l) * gL_l(l);
354 double gR_l2 = gR_l(l) * gR_l(l);
355 B2l = gL_l(l) * gR_l2 * deltaGLe + gR_l(l) * gL_l2 * deltaGRe
356 + gR_l2 * gL_l(l) * deltaGLf + gL_l2 * gR_l(l) * deltaGRf;
357
358 return B2l;
359}

◆ deltaB2q()

double LEP2GIMR::deltaB2q ( const QCD::quark  q,
const double  GIMRParam_i[] 
) const
private

Definition at line 287 of file LEP2GIMR.cpp.

288{
289 double B2q;
290 double deltaGLf = GIMRParam_i[delta_gLf];
291 double deltaGRf = GIMRParam_i[delta_gRf];
292 double deltaGLe = GIMRParam_i[delta_gLe];
293 double deltaGRe = GIMRParam_i[delta_gRe];
294 double gL_q2 = gL_q(q) * gL_q(q);
295 double gR_q2 = gR_q(q) * gR_q(q);
296 double gL_l2 = gL_l(StandardModel::ELECTRON) * gL_l(StandardModel::ELECTRON);
297 double gR_l2 = gR_l(StandardModel::ELECTRON) * gR_l(StandardModel::ELECTRON);
298 B2q = gL_l(StandardModel::ELECTRON) * gR_q2 * deltaGLe
299 + gR_l(StandardModel::ELECTRON) * gL_q2 * deltaGRe
300 + gR_l2 * gL_q(q) * deltaGLf + gL_l2 * gR_q(q) * deltaGRf;
301
302 return B2q;
303}

◆ gL_l()

double LEP2GIMR::gL_l ( const QCD::lepton  l) const
private

Definition at line 201 of file LEP2GIMR.cpp.

202{
203 double gA = SM.getLeptons(l).getIsospin();
204 double gV = SM.getLeptons(l).getIsospin()-2.0*SM.getLeptons(l).getCharge()*SM.sW2();
205 double gL = (gA + gV)/2.;
206
207 return gL;
208 }
Particle::getIsospin
double getIsospin() const
A get method to access the particle isospin.
Definition: Particle.h:115
Particle::getCharge
double getCharge() const
A get method to access the particle charge.
Definition: Particle.h:97
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::sW2
virtual const double sW2(const double Mw_i) const
The square of the sine of the weak mixing angle in the on-shell scheme, denoted as .
Definition: StandardModel/src/StandardModel.cpp:1087

◆ gL_q()

double LEP2GIMR::gL_q ( const QCD::quark  q) const
private

Definition at line 221 of file LEP2GIMR.cpp.

222{
223 double gA = SM.getQuarks(q).getIsospin();
224 double gV = SM.getQuarks(q).getIsospin()-2.0*SM.getQuarks(q).getCharge()*SM.sW2();
225 double gL = (gA + gV)/2.;
226
227 return gL;
228 }
QCD::getQuarks
const Particle & getQuarks(const QCD::quark q) const
A get method to access a quark as an object of the type Particle.
Definition: QCD.h:536

◆ gR_l()

double LEP2GIMR::gR_l ( const QCD::lepton  l) const
private

Definition at line 211 of file LEP2GIMR.cpp.

212{
213 double gA = SM.getLeptons(l).getIsospin();
214 double gV = SM.getLeptons(l).getIsospin()-2.0*SM.getLeptons(l).getCharge()*SM.sW2();
215 double gR = (gV - gA)/2.;
216
217 return gR;
218 }

◆ gR_q()

double LEP2GIMR::gR_q ( const QCD::quark  q) const
private

Definition at line 231 of file LEP2GIMR.cpp.

232{
233 double gA = SM.getQuarks(q).getIsospin();
234 double gV = SM.getQuarks(q).getIsospin()-2.0*SM.getQuarks(q).getCharge()*SM.sW2();
235 double gR = (gV - gA)/2.;
236
237 return gR;
238 }

◆ sigma_l_LEP2_GIMR()

double LEP2GIMR::sigma_l_LEP2_GIMR ( const QCD::lepton  l,
const double  s,
const double  GIMRParam_i[] 
) const

Definition at line 18 of file LEP2GIMR.cpp.

20{
21 double Mz = SM.getMz();
22 double Nf = 1.0;
23 double Qe = SM.getLeptons(StandardModel::ELECTRON).getCharge();
24 double Ql = SM.getLeptons(l).getCharge();
25 double GammaZ = SM.Gamma_Z();
26 double deltaGammaZ = GIMRParam_i[delta_GammaZ];
27 double alpha = SM.ale_OS(sqrt(s), FULLNLO);
28 double alpha2 = alpha*alpha;
29 double deltaMzsq = GIMRParam_i[delta_Mz2];
30 double GRl = gR_l(l);
31 double GLl = gL_l(l);
32 double GRe = gR_l(StandardModel::ELECTRON);
33 double GLe = gL_l(StandardModel::ELECTRON);
34 double sW2 = SM.sW2();
35 double cW = sqrt(1. - SM.cW2());
36 double cW2 = cW*cW;
37 double CLL = GIMRParam_i[C_LL];
38 double CLR = GIMRParam_i[C_LR];
39 double CRL = GIMRParam_i[C_RL];
40 double CRR = GIMRParam_i[C_RR];
41 complex denom = complex(s - Mz*Mz, Mz*GammaZ, false);
42 complex chiZ = s/denom;
43 complex deltachiZ = chiZ/denom*(complex(deltaMzsq,-Mz*deltaGammaZ-GammaZ*deltaMzsq/(2.0*Mz)));
44 complex chideltachi = (chiZ*deltachiZ.conjugate()+chiZ.conjugate()*deltachiZ);
45 double dA1l=deltaA1l(l,GIMRParam_i);
46 double dA2l=deltaA2l(l,GIMRParam_i);
47 double dB1l=deltaB1l(l,GIMRParam_i);
48 double dB2l=deltaB2l(l,GIMRParam_i);
49
50
51 double ds = alpha/(6.0)*Nf*(CLL*(GLl*GLe/(cW2*sW2)*chiZ.real()+Ql*Qe)
52 + CLR*(GLe*GRl/(cW2*sW2)*chiZ.real()+Ql*Qe)
53 + CRL*(GLl*GRe/(cW2*sW2)*chiZ.real()+Qe*Ql)
54 + CRR*(GRl*GRe/(cW2*sW2)*chiZ.real()+Qe*Ql))
55 + 2.*alpha2*Nf*M_PI*Ql*Qe/(3.*cW2*sW2*s)*(chiZ.real()*(dA1l+dA2l)
56 + deltachiZ.real()*(GLl+GRl)*(GLe+GRe))
57 + alpha2*Nf*M_PI/(3.*cW2*cW2*sW2*sW2*s)*(2.*chiZ.abs2()*(dB1l+dB2l)
58 + chideltachi.real()*(GLl*GLl+GRl*GRl)*(GLe*GLe+GRe*GRe));
59
60 return ds;
61 }
FULLNLO
@ FULLNLO
Definition: OrderScheme.h:38
GammaZ
An observable class for the total decay width of the boson.
Definition: GammaZ.h:32
LEP2GIMR::deltaB1l
double deltaB1l(const QCD::lepton l, const double GIMRParam_i[]) const
Definition: LEP2GIMR.cpp:331
LEP2GIMR::deltaB2l
double deltaB2l(const QCD::lepton l, const double GIMRParam_i[]) const
Definition: LEP2GIMR.cpp:346
LEP2GIMR::deltaA1l
double deltaA1l(const QCD::lepton l, const double GIMRParam_i[]) const
Definition: LEP2GIMR.cpp:305
LEP2GIMR::deltaA2l
double deltaA2l(const QCD::lepton l, const double GIMRParam_i[]) const
Definition: LEP2GIMR.cpp:318
StandardModel::getMz
const double getMz() const
A get method to access the mass of the boson .
Definition: StandardModel.h:753
StandardModel::ale_OS
const double ale_OS(const double mu, orders order=FULLNLO) const
The running electromagnetic coupling in the on-shell scheme.
Definition: StandardModel/src/StandardModel.cpp:605
StandardModel::Gamma_Z
virtual const double Gamma_Z() const
The total decay width of the boson, .
Definition: StandardModel/src/StandardModel.cpp:1433
StandardModel::cW2
virtual const double cW2(const double Mw_i) const
The square of the cosine of the weak mixing angle in the on-shell scheme, denoted as .
Definition: StandardModel/src/StandardModel.cpp:1076
s
Test Observable.

◆ sigma_q_LEP2_GIMR()

double LEP2GIMR::sigma_q_LEP2_GIMR ( const QCD::quark  q,
const double  s,
const double  GIMRParam_i[] 
) const

Definition at line 111 of file LEP2GIMR.cpp.

113{
114 double Mz = SM.getMz();
115 double Nf = 3.0;
116 double Qe = SM.getLeptons(StandardModel::ELECTRON).getCharge();
117 double Qq = SM.getQuarks(q).getCharge();
118 double GammaZ = SM.Gamma_Z();
119 double deltaGammaZ = GIMRParam_i[delta_GammaZ];
120 double alpha = SM.ale_OS(sqrt(s), FULLNLO);
121 double alpha2 = alpha*alpha;
122 double deltaMzsq = GIMRParam_i[delta_Mz2];
123 double GRq = gR_q(q);
124 double GLq = gL_q(q);
125 double GRe = gR_l(StandardModel::ELECTRON);
126 double GLe = gL_l(StandardModel::ELECTRON);
127 double sW2 = SM.sW2();
128 double cW = sqrt(1. - SM.cW2());
129 double cW2 = cW*cW;
130 double CLL = GIMRParam_i[C_LL];
131 double CLR = GIMRParam_i[C_LR];
132 double CRL = GIMRParam_i[C_RL];
133 double CRR = GIMRParam_i[C_RR];
134 complex denom = complex(s - Mz*Mz, Mz*GammaZ, false);
135 complex chiZ = s/denom;
136 complex deltachiZ = chiZ/denom*(complex(deltaMzsq,-Mz*deltaGammaZ-GammaZ*deltaMzsq/(2.0*Mz)));
137 complex chideltachi = (chiZ*deltachiZ.conjugate()+chiZ.conjugate()*deltachiZ);
138 double dA1q=deltaA1q(q,GIMRParam_i);
139 double dA2q=deltaA2q(q,GIMRParam_i);
140 double dB1q=deltaB1q(q,GIMRParam_i);
141 double dB2q=deltaB2q(q,GIMRParam_i);
142
143 double ds = alpha/(6.0)*Nf*(CLL*(GLq*GLe/(cW2*sW2)*chiZ.real()+Qq*Qe)
144 + CLR*(GLe*GRq/(cW2*sW2)*chiZ.real()+Qq*Qe)
145 + CRL*(GLq*GRe/(cW2*sW2)*chiZ.real()+Qe*Qq)
146 + CRR*(GRq*GRe/(cW2*sW2)*chiZ.real()+Qe*Qq))
147 + 2.*alpha2*Nf*M_PI*Qq*Qe/(3.*cW2*sW2*s)*(chiZ.real()*(dA1q+dA2q)
148 + deltachiZ.real()*(GLq+GRq)*(GLe+GRe))
149 + alpha2*Nf*M_PI/(3.*cW2*cW2*sW2*sW2*s)*(2.*chiZ.abs2()*(dB1q+dB2q)
150 + chideltachi.real()*(GLq*GLq+GRq*GRq)*(GLe*GLe+GRe*GRe));
151
152 return ds;
153 }
LEP2GIMR::deltaA2q
double deltaA2q(const QCD::quark q, const double GIMRParam_i[]) const
Definition: LEP2GIMR.cpp:255
LEP2GIMR::deltaA1q
double deltaA1q(const QCD::quark q, const double GIMRParam_i[]) const
Definition: LEP2GIMR.cpp:241
LEP2GIMR::deltaB1q
double deltaB1q(const QCD::quark q, const double GIMRParam_i[]) const
Definition: LEP2GIMR.cpp:269
LEP2GIMR::deltaB2q
double deltaB2q(const QCD::quark q, const double GIMRParam_i[]) const
Definition: LEP2GIMR.cpp:287

◆ sigmaFminusB_l_LEP2_GIMR()

double LEP2GIMR::sigmaFminusB_l_LEP2_GIMR ( const QCD::lepton  l,
const double  s,
const double  GIMRParam_i[] 
) const

Definition at line 64 of file LEP2GIMR.cpp.

66{
67 double Mz = SM.getMz();
68 double Nf = 1.0;
69 double Qe = SM.getLeptons(StandardModel::ELECTRON).getCharge();
70 double Ql = SM.getLeptons(l).getCharge();
71 double GammaZ = SM.Gamma_Z();
72 double deltaGammaZ = GIMRParam_i[delta_GammaZ];
73 double alpha = SM.ale_OS(sqrt(s), FULLNLO);
74 double alpha2 = alpha*alpha;
75 double deltaMzsq = GIMRParam_i[delta_Mz2];
76 double GRl = gR_l(l);
77 double GLl = gL_l(l);
78 double GRe = gR_l(StandardModel::ELECTRON);
79 double GLe = gL_l(StandardModel::ELECTRON);
80 double sW2 = SM.sW2();
81 double cW = sqrt(1. - SM.cW2());
82 double cW2 = cW*cW;
83 double CLL = GIMRParam_i[C_LL];
84 double CLR = GIMRParam_i[C_LR];
85 double CRL = GIMRParam_i[C_RL];
86 double CRR = GIMRParam_i[C_RR];
87 complex denom = complex(s - Mz*Mz, Mz*GammaZ, false);
88 complex chiZ = s/denom;
89 complex deltachiZ = chiZ/denom*(complex(deltaMzsq,-Mz*deltaGammaZ-GammaZ*deltaMzsq/(2.0*Mz)));
90 complex chideltachi = (chiZ*deltachiZ.conjugate()+chiZ.conjugate()*deltachiZ);
91 double dA1l=deltaA1l(l,GIMRParam_i);
92 double dA2l=deltaA2l(l,GIMRParam_i);
93 double dB1l=deltaB1l(l,GIMRParam_i);
94 double dB2l=deltaB2l(l,GIMRParam_i);
95
96
97 double ds = alpha/(8.0)*Nf*(CLL*(GLl*GLe/(cW2*sW2)*chiZ.real()+Ql*Qe)
98 - CLR*(GLe*GRl/(cW2*sW2)*chiZ.real()+Ql*Qe)
99 - CRL*(GLl*GRe/(cW2*sW2)*chiZ.real()+Qe*Ql)
100 + CRR*(GRl*GRe/(cW2*sW2)*chiZ.real()+Qe*Ql))
101 + alpha2*Nf*M_PI*Ql*Qe/(2.*cW2*sW2*s)*(chiZ.real()*(dA1l-dA2l)
102 + deltachiZ.real()*(GLl-GRl)*(GLe-GRe))
103 + alpha2*Nf*M_PI/(4.*cW2*cW2*sW2*sW2*s)*(2.*chiZ.abs2()*(dB1l-dB2l)
104 + chideltachi.real()*(GLl*GLl-GRl*GRl)*(GLe*GLe-GRe*GRe));
105
106 return ds;
107 }

◆ sigmaFminusB_q_LEP2_GIMR()

double LEP2GIMR::sigmaFminusB_q_LEP2_GIMR ( const QCD::quark  q,
const double  s,
const double  GIMRParam_i[] 
) const

Definition at line 155 of file LEP2GIMR.cpp.

157{
158 double Mz = SM.getMz();
159 double Nf = 3.0;
160 double Qe = SM.getLeptons(StandardModel::ELECTRON).getCharge();
161 double Qq = SM.getQuarks(q).getCharge();
162 double GammaZ = SM.Gamma_Z();
163 double deltaGammaZ = GIMRParam_i[delta_GammaZ];
164 double alpha = SM.ale_OS(sqrt(s), FULLNLO);
165 double alpha2 = alpha*alpha;
166 double deltaMzsq = GIMRParam_i[delta_Mz2];
167 double GRq = gR_q(q);
168 double GLq = gL_q(q);
169 double GRe = gR_l(StandardModel::ELECTRON);
170 double GLe = gL_l(StandardModel::ELECTRON);
171 double sW2 = SM.sW2();
172 double cW = sqrt(1. - SM.cW2());
173 double cW2 = cW*cW;
174 double CLL = GIMRParam_i[C_LL];
175 double CLR = GIMRParam_i[C_LR];
176 double CRL = GIMRParam_i[C_RL];
177 double CRR = GIMRParam_i[C_RR];
178 complex denom = complex(s - Mz*Mz, Mz*GammaZ, false);
179 complex chiZ = s/denom;
180 complex deltachiZ = chiZ/denom*(complex(deltaMzsq,-Mz*deltaGammaZ-GammaZ*deltaMzsq/(2.0*Mz)));
181 complex chideltachi = (chiZ*deltachiZ.conjugate()+chiZ.conjugate()*deltachiZ);
182 double dA1q=deltaA1q(q,GIMRParam_i);
183 double dA2q=deltaA2q(q,GIMRParam_i);
184 double dB1q=deltaB1q(q,GIMRParam_i);
185 double dB2q=deltaB2q(q,GIMRParam_i);
186
187 double ds = alpha/(8.0)*Nf*(CLL*(GLq*GLe/(cW2*sW2)*chiZ.real()+Qq*Qe)
188 - CLR*(GLe*GRq/(cW2*sW2)*chiZ.real()+Qq*Qe)
189 - CRL*(GLq*GRe/(cW2*sW2)*chiZ.real()+Qe*Qq)
190 + CRR*(GRq*GRe/(cW2*sW2)*chiZ.real()+Qe*Qq))
191 + alpha2*Nf*M_PI*Qq*Qe/(4.*cW2*sW2*s)*(chiZ.real()*(dA1q-dA2q)
192 + deltachiZ.real()*(GLq-GRq)*(GLe-GRe))
193 + alpha2*Nf*M_PI/(4.*cW2*cW2*sW2*sW2*s)*(2.*chiZ.abs2()*(dB1q-dB2q)
194 + chideltachi.real()*(GLq*GLq-GRq*GRq)*(GLe*GLe-GRe*GRe));
195
196 return ds;
197 }

Member Data Documentation

◆ SM

const StandardModel& LEP2GIMR::SM
private

Definition at line 48 of file LEP2GIMR.h.

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