11 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
17 double m1sq=msQhat2(0,0).abs();
18 double m2sq=msQhat2(1,1).abs();
19 return (m2sq-m1sq)/
sqrt(m2sq*m1sq);
23 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
29 double m1sq=msQhat2(0,0).abs();
30 double m3sq=msQhat2(2,2).abs();
31 return (m3sq-m1sq)/
sqrt(m3sq*m1sq);
35 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
41 double m2sq=msQhat2(1,1).abs();
42 double m3sq=msQhat2(2,2).abs();
43 return (m3sq-m2sq)/
sqrt(m3sq*m2sq);
47 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
53 double m1sq=msUhat2(0,0).abs();
54 double m2sq=msUhat2(1,1).abs();
55 return (m2sq-m1sq)/
sqrt(m2sq*m1sq);
59 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
65 double m1sq=msUhat2(0,0).abs();
66 double m3sq=msUhat2(2,2).abs();
67 return (m3sq-m1sq)/
sqrt(m3sq*m1sq);
71 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
77 double m2sq=msUhat2(1,1).abs();
78 double m3sq=msUhat2(2,2).abs();
79 return (m3sq-m2sq)/
sqrt(m3sq*m2sq);
83 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
89 double m1sq=msDhat2(0,0).abs();
90 double m2sq=msDhat2(1,1).abs();
91 return (m2sq-m1sq)/
sqrt(m2sq*m1sq);
95 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
101 double m1sq=msDhat2(0,0).abs();
102 double m3sq=msDhat2(2,2).abs();
103 return (m3sq-m1sq)/
sqrt(m3sq*m1sq);
107 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
113 double m2sq=msDhat2(1,1).abs();
114 double m3sq=msDhat2(2,2).abs();
115 return (m3sq-m2sq)/
sqrt(m3sq*m2sq);
119 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
125 double m1sq=msLhat2(0,0).abs();
126 double m2sq=msLhat2(1,1).abs();
127 return (m2sq-m1sq)/
sqrt(m2sq*m1sq);
131 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
137 double m1sq=msLhat2(0,0).abs();
138 double m3sq=msLhat2(2,2).abs();
139 return (m3sq-m1sq)/
sqrt(m3sq*m1sq);
143 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
149 double m2sq=msLhat2(1,1).abs();
150 double m3sq=msLhat2(2,2).abs();
151 return (m3sq-m2sq)/
sqrt(m3sq*m2sq);
155 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
161 double m1sq=msEhat2(0,0).abs();
162 double m2sq=msEhat2(1,1).abs();
163 return (m2sq-m1sq)/
sqrt(m2sq*m1sq);
167 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
173 double m1sq=msEhat2(0,0).abs();
174 double m3sq=msEhat2(2,2).abs();
175 return (m3sq-m1sq)/
sqrt(m3sq*m1sq);
179 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
185 double m2sq=msEhat2(1,1).abs();
186 double m3sq=msEhat2(2,2).abs();
187 return (m3sq-m2sq)/
sqrt(m3sq*m2sq);
191 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
199 double Mu2average_21=
sqrt(msQhat2(1,1)*msUhat2(0,0)).
abs();
202 return TUhat(1,0).real()*v2/Mu2average_21;
206 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
214 double Mu2average_31=
sqrt(msQhat2(2,2)*msUhat2(0,0)).
abs();
217 return TUhat(2,0).real()*v2/Mu2average_31;
221 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
229 double Mu2average_32=
sqrt(msQhat2(2,2)*msUhat2(1,1)).
abs();
232 return TUhat(2,1).real()*v2/Mu2average_32;
236 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
244 double Me2average_21=
sqrt(msLhat2(1,1)*msEhat2(0,0)).
abs();
247 return TEhat(1,0).real()*v1/Me2average_21;
251 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
259 double Me2average_12=
sqrt(msLhat2(0,0)*msEhat2(1,1)).
abs();
262 return TEhat(0,1).real()*v1/Me2average_12;
266 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
274 double Me2average_31=
sqrt(msLhat2(2,2)*msEhat2(0,0)).
abs();
277 return TEhat(2,0).real()*v1/Me2average_31;
281 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
289 double Me2average_13=
sqrt(msLhat2(0,0)*msEhat2(2,2)).
abs();
292 return TEhat(0,2).real()*v1/Me2average_13;
296 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
304 double Me2average_32=
sqrt(msLhat2(2,2)*msEhat2(1,1)).
abs();
307 return TEhat(2,1).real()*v1/Me2average_32;
311 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
319 double Me2average_23=
sqrt(msLhat2(1,1)*msEhat2(2,2)).
abs();
322 return TEhat(1,2).real()*v1/Me2average_23;
326 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
334 double Me2average_31=
sqrt(msLhat2(2,2)*msEhat2(0,0)).
abs();
337 return log10(TEhat(2,0).abs()*v1/Me2average_31);
341 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
349 double Me2average_32=
sqrt(msLhat2(2,2)*msEhat2(1,1)).
abs();
352 return log10(TEhat(2,1).abs()*v1/Me2average_32);
356 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
363 double Me2average_32=
sqrt(msLhat2(2,2)*msEhat2(1,1)).
abs();
369 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
376 double Me2average_32=(
sqrt(msLhat2(2,2)).
abs()+
sqrt(msEhat2(1,1)).
abs())/2.0;
378 return Me2average_32;
382 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
388 double denominator =
sqrt(TEhat(2,1).abs()*TEhat(1,2).abs());
390 return (TEhat(2,1).abs()-TEhat(1,2).abs())/denominator;
394 :
ThObservable(SM_i), mySUSY(static_cast<const SUSY&> (SM_i))
401 double denominator =
sqrt(msLhat2(2,2).abs()*msEhat2(2,2).abs());
403 return (msLhat2(2,2).abs()-msEhat2(2,2).abs())/denominator;