a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
SFermiondeltas.cpp
Go to the documentation of this file.
1/*
2 * Copyright (C) 2016 HEPfit Collaboration
3 *
4 *
5 * For the licensing terms see doc/COPYING.
6 */
7
8#include "SFermiondeltas.h"
9
11: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
12{}
13
15{
16 gslpp::matrix<gslpp::complex> msQhat2( mySUSY.getMsQhat2() );
17 double m1sq=msQhat2(0,0).abs();
18 double m2sq=msQhat2(1,1).abs();
19 return (m2sq-m1sq)/sqrt(m2sq*m1sq);
20}
21
23: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
24{}
25
27{
28 gslpp::matrix<gslpp::complex> msQhat2( mySUSY.getMsQhat2() );
29 double m1sq=msQhat2(0,0).abs();
30 double m3sq=msQhat2(2,2).abs();
31 return (m3sq-m1sq)/sqrt(m3sq*m1sq);
32}
33
35: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
36{}
37
39{
40 gslpp::matrix<gslpp::complex> msQhat2( mySUSY.getMsQhat2() );
41 double m2sq=msQhat2(1,1).abs();
42 double m3sq=msQhat2(2,2).abs();
43 return (m3sq-m2sq)/sqrt(m3sq*m2sq);
44}
45
47: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
48{}
49
51{
52 gslpp::matrix<gslpp::complex> msUhat2( mySUSY.getMsUhat2() );
53 double m1sq=msUhat2(0,0).abs();
54 double m2sq=msUhat2(1,1).abs();
55 return (m2sq-m1sq)/sqrt(m2sq*m1sq);
56}
57
59: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
60{}
61
63{
64 gslpp::matrix<gslpp::complex> msUhat2( mySUSY.getMsUhat2() );
65 double m1sq=msUhat2(0,0).abs();
66 double m3sq=msUhat2(2,2).abs();
67 return (m3sq-m1sq)/sqrt(m3sq*m1sq);
68}
69
71: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
72{}
73
75{
76 gslpp::matrix<gslpp::complex> msUhat2( mySUSY.getMsUhat2() );
77 double m2sq=msUhat2(1,1).abs();
78 double m3sq=msUhat2(2,2).abs();
79 return (m3sq-m2sq)/sqrt(m3sq*m2sq);
80}
81
83: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
84{}
85
87{
88 gslpp::matrix<gslpp::complex> msDhat2( mySUSY.getMsDhat2() );
89 double m1sq=msDhat2(0,0).abs();
90 double m2sq=msDhat2(1,1).abs();
91 return (m2sq-m1sq)/sqrt(m2sq*m1sq);
92}
93
95: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
96{}
97
99{
100 gslpp::matrix<gslpp::complex> msDhat2( mySUSY.getMsDhat2() );
101 double m1sq=msDhat2(0,0).abs();
102 double m3sq=msDhat2(2,2).abs();
103 return (m3sq-m1sq)/sqrt(m3sq*m1sq);
104}
105
107: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
108{}
109
111{
112 gslpp::matrix<gslpp::complex> msDhat2( mySUSY.getMsDhat2() );
113 double m2sq=msDhat2(1,1).abs();
114 double m3sq=msDhat2(2,2).abs();
115 return (m3sq-m2sq)/sqrt(m3sq*m2sq);
116}
117
119: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
120{}
121
123{
124 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
125 double m1sq=msLhat2(0,0).abs();
126 double m2sq=msLhat2(1,1).abs();
127 return (m2sq-m1sq)/sqrt(m2sq*m1sq);
128}
129
131: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
132{}
133
135{
136 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
137 double m1sq=msLhat2(0,0).abs();
138 double m3sq=msLhat2(2,2).abs();
139 return (m3sq-m1sq)/sqrt(m3sq*m1sq);
140}
141
143: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
144{}
145
147{
148 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
149 double m2sq=msLhat2(1,1).abs();
150 double m3sq=msLhat2(2,2).abs();
151 return (m3sq-m2sq)/sqrt(m3sq*m2sq);
152}
153
155: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
156{}
157
159{
160 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
161 double m1sq=msEhat2(0,0).abs();
162 double m2sq=msEhat2(1,1).abs();
163 return (m2sq-m1sq)/sqrt(m2sq*m1sq);
164}
165
167: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
168{}
169
171{
172 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
173 double m1sq=msEhat2(0,0).abs();
174 double m3sq=msEhat2(2,2).abs();
175 return (m3sq-m1sq)/sqrt(m3sq*m1sq);
176}
177
179: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
180{}
181
183{
184 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
185 double m2sq=msEhat2(1,1).abs();
186 double m3sq=msEhat2(2,2).abs();
187 return (m3sq-m2sq)/sqrt(m3sq*m2sq);
188}
189
191: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
192{}
193
195{
196 gslpp::matrix<gslpp::complex> TUhat = mySUSY.getTUhat();
197 gslpp::matrix<gslpp::complex> msQhat2( mySUSY.getMsQhat2() );
198 gslpp::matrix<gslpp::complex> msUhat2( mySUSY.getMsUhat2() );
199 double Mu2average_21=sqrt(msQhat2(1,1)*msUhat2(0,0)).abs();
200 double v2 = mySUSY.v2();
201
202 return TUhat(1,0).real()*v2/Mu2average_21; //deltaRL_12_u
203}
204
206: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
207{}
208
210{
211 gslpp::matrix<gslpp::complex> TUhat = mySUSY.getTUhat();
212 gslpp::matrix<gslpp::complex> msQhat2( mySUSY.getMsQhat2() );
213 gslpp::matrix<gslpp::complex> msUhat2( mySUSY.getMsUhat2() );
214 double Mu2average_31=sqrt(msQhat2(2,2)*msUhat2(0,0)).abs();
215 double v2 = mySUSY.v2();
216
217 return TUhat(2,0).real()*v2/Mu2average_31; //deltaRL_13_u
218}
219
221: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
222{}
223
225{
226 gslpp::matrix<gslpp::complex> TUhat = mySUSY.getTUhat();
227 gslpp::matrix<gslpp::complex> msQhat2( mySUSY.getMsQhat2() );
228 gslpp::matrix<gslpp::complex> msUhat2( mySUSY.getMsUhat2() );
229 double Mu2average_32=sqrt(msQhat2(2,2)*msUhat2(1,1)).abs();
230 double v2 = mySUSY.v2();
231
232 return TUhat(2,1).real()*v2/Mu2average_32; //deltaRL_23_u
233}
234
236: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
237{}
238
240{
241 gslpp::matrix<gslpp::complex> TEhat = mySUSY.getTEhat();
242 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
243 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
244 double Me2average_21=sqrt(msLhat2(1,1)*msEhat2(0,0)).abs();
245 double v1 = mySUSY.v1();
246
247 return TEhat(1,0).real()*v1/Me2average_21; //deltaRL_12_e
248}
249
251: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
252{}
253
255{
256 gslpp::matrix<gslpp::complex> TEhat = mySUSY.getTEhat();
257 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
258 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
259 double Me2average_12=sqrt(msLhat2(0,0)*msEhat2(1,1)).abs();
260 double v1 = mySUSY.v1();
261
262 return TEhat(0,1).real()*v1/Me2average_12; //deltaRL_12_e
263}
264
266: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
267{}
268
270{
271 gslpp::matrix<gslpp::complex> TEhat = mySUSY.getTEhat();
272 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
273 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
274 double Me2average_31=sqrt(msLhat2(2,2)*msEhat2(0,0)).abs();
275 double v1 = mySUSY.v1();
276
277 return TEhat(2,0).real()*v1/Me2average_31; //deltaRL_13_e
278}
279
281: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
282{}
283
285{
286 gslpp::matrix<gslpp::complex> TEhat = mySUSY.getTEhat();
287 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
288 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
289 double Me2average_13=sqrt(msLhat2(0,0)*msEhat2(2,2)).abs();
290 double v1 = mySUSY.v1();
291
292 return TEhat(0,2).real()*v1/Me2average_13; //deltaRL_13_e
293}
294
296: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
297{}
298
300{
301 gslpp::matrix<gslpp::complex> TEhat = mySUSY.getTEhat();
302 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
303 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
304 double Me2average_32=sqrt(msLhat2(2,2)*msEhat2(1,1)).abs();
305 double v1 = mySUSY.v1();
306
307 return TEhat(2,1).real()*v1/Me2average_32; //deltaRL_23_e
308}
309
311: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
312{}
313
315{
316 gslpp::matrix<gslpp::complex> TEhat = mySUSY.getTEhat();
317 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
318 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
319 double Me2average_23=sqrt(msLhat2(1,1)*msEhat2(2,2)).abs();
320 double v1 = mySUSY.v1();
321
322 return TEhat(1,2).real()*v1/Me2average_23; //deltaRL_32_e
323}
324
326: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
327{}
328
330{
331 gslpp::matrix<gslpp::complex> TEhat = mySUSY.getTEhat();
332 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
333 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
334 double Me2average_31=sqrt(msLhat2(2,2)*msEhat2(0,0)).abs();
335 double v1 = mySUSY.v1();
336
337 return log10(TEhat(2,0).abs()*v1/Me2average_31); //logdeltaRL_13_e
338}
339
341: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
342{}
343
345{
346 gslpp::matrix<gslpp::complex> TEhat = mySUSY.getTEhat();
347 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
348 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
349 double Me2average_32=sqrt(msLhat2(2,2)*msEhat2(1,1)).abs();
350 double v1 = mySUSY.v1();
351
352 return log10(TEhat(2,1).abs()*v1/Me2average_32); //logdeltaRL_23_e
353}
354
356: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
357{}
358
360{
361 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
362 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
363 double Me2average_32=sqrt(msLhat2(2,2)*msEhat2(1,1)).abs();
364
365 return log10(sqrt(Me2average_32)); //logmslepton
366}
367
369: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
370{}
371
373{
374 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
375 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
376 double Me2average_32=(sqrt(msLhat2(2,2)).abs()+sqrt(msEhat2(1,1)).abs())/2.0;
377
378 return Me2average_32; //mslepton
379}
380
382: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
383{}
384
386{
387 gslpp::matrix<gslpp::complex> TEhat = mySUSY.getTEhat();
388 double denominator = sqrt(TEhat(2,1).abs()*TEhat(1,2).abs());
389
390 return (TEhat(2,1).abs()-TEhat(1,2).abs())/denominator;
391}
392
394: ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
395{}
396
398{
399 gslpp::matrix<gslpp::complex> msLhat2( mySUSY.getMsLhat2() );
400 gslpp::matrix<gslpp::complex> msEhat2( mySUSY.getMsEhat2() );
401 double denominator = sqrt(msLhat2(2,2).abs()*msEhat2(2,2).abs());
402
403 return (msLhat2(2,2).abs()-msEhat2(2,2).abs())/denominator;
404}
A base class for SUSY models.
Definition: SUSY.h:33
gslpp::matrix< gslpp::complex > getMsLhat2() const
Definition: SUSY.h:426
double v2() const
Definition: SUSY.cpp:384
gslpp::matrix< gslpp::complex > getMsDhat2() const
Definition: SUSY.h:344
gslpp::matrix< gslpp::complex > getTEhat() const
Gets the trilinear-coupling matrix for charged sleptons.
Definition: SUSY.h:454
gslpp::matrix< gslpp::complex > getMsUhat2() const
Definition: SUSY.h:339
gslpp::matrix< gslpp::complex > getMsQhat2() const
Definition: SUSY.h:334
gslpp::matrix< gslpp::complex > getMsEhat2() const
Definition: SUSY.h:436
gslpp::matrix< gslpp::complex > getTUhat() const
Gets the trilinear-coupling matrix for up-type squarks.
Definition: SUSY.h:353
double v1() const
Definition: SUSY.cpp:379
A model class for the Standard Model.
A class for a model prediction of an observable.
Definition: ThObservable.h:25
virtual double computeThValue()=0
A member to be overloaded by the respective theory observable. class that calculates the value of the...
const SUSY & mySUSY
deltaLL1_l(const StandardModel &SM_i)
deltaLL1_l constructor.
double computeThValue()
deltaLL1_q(const StandardModel &SM_i)
deltaLL1_q constructor.
double computeThValue()
const SUSY & mySUSY
double computeThValue()
const SUSY & mySUSY
deltaLL2_l(const StandardModel &SM_i)
deltaLL2_l constructor.
double computeThValue()
const SUSY & mySUSY
deltaLL2_q(const StandardModel &SM_i)
deltaLL2_q constructor.
const SUSY & mySUSY
deltaLL3_l(const StandardModel &SM_i)
deltaLL3_l constructor.
double computeThValue()
double computeThValue()
deltaLL3_q(const StandardModel &SM_i)
deltaLL3_q constructor.
const SUSY & mySUSY
deltaLLRR_l(const StandardModel &SM_i)
deltaLLRR_l constructor.
double computeThValue()
deltaRL_12_e(const StandardModel &SM_i)
deltaRL_12_e constructor.
const SUSY & mySUSY
const SUSY & mySUSY
double computeThValue()
deltaRL_12_u(const StandardModel &SM_i)
deltaRL_12_u constructor.
double computeThValue()
const SUSY & mySUSY
deltaRL_13_e(const StandardModel &SM_i)
deltaRL_13_e constructor.
deltaRL_13_u(const StandardModel &SM_i)
deltaRL_13_u constructor.
double computeThValue()
const SUSY & mySUSY
deltaRL_21_e(const StandardModel &SM_i)
deltaRL_21_e constructor.
double computeThValue()
const SUSY & mySUSY
const SUSY & mySUSY
double computeThValue()
deltaRL_23_e(const StandardModel &SM_i)
deltaRL_23_e constructor.
const SUSY & mySUSY
double computeThValue()
deltaRL_23_u(const StandardModel &SM_i)
deltaRL_23_u constructor.
const SUSY & mySUSY
deltaRL_31_e(const StandardModel &SM_i)
deltaRL_31_e constructor.
double computeThValue()
double computeThValue()
deltaRL_32_e(const StandardModel &SM_i)
deltaRL_32_e constructor.
const SUSY & mySUSY
double computeThValue()
deltaRR1_d(const StandardModel &SM_i)
deltaRR1_d constructor.
const SUSY & mySUSY
deltaRR1_e(const StandardModel &SM_i)
deltaRR1_e constructor.
const SUSY & mySUSY
double computeThValue()
double computeThValue()
deltaRR1_u(const StandardModel &SM_i)
deltaRR1_u constructor.
const SUSY & mySUSY
const SUSY & mySUSY
deltaRR2_d(const StandardModel &SM_i)
deltaRR2_d constructor.
double computeThValue()
deltaRR2_e(const StandardModel &SM_i)
deltaRR2_e constructor.
const SUSY & mySUSY
double computeThValue()
const SUSY & mySUSY
double computeThValue()
deltaRR2_u(const StandardModel &SM_i)
deltaRR2_u constructor.
double computeThValue()
deltaRR3_d(const StandardModel &SM_i)
deltaRR3_d constructor.
const SUSY & mySUSY
deltaRR3_e(const StandardModel &SM_i)
deltaRR3_e constructor.
const SUSY & mySUSY
double computeThValue()
double computeThValue()
const SUSY & mySUSY
deltaRR3_u(const StandardModel &SM_i)
deltaRR3_u constructor.
deltaTEhat23(const StandardModel &SM_i)
deltaTEhat23 constructor.
logdeltaRL_13_e(const StandardModel &SM_i)
logdeltaRL_13_e constructor.
const SUSY & mySUSY
const SUSY & mySUSY
logdeltaRL_23_e(const StandardModel &SM_i)
logdeltaRL_23_e constructor.
logmslepton(const StandardModel &SM_i)
logmslepton constructor.
double computeThValue()
const SUSY & mySUSY
const SUSY & mySUSY
mslepton(const StandardModel &SM_i)
mslepton constructor.
double computeThValue()