17 double m1sq=msQhat2(0,0).abs();
18 double m2sq=msQhat2(1,1).abs();
19 return (m2sq-m1sq)/
sqrt(m2sq*m1sq);
29 double m1sq=msQhat2(0,0).abs();
30 double m3sq=msQhat2(2,2).abs();
31 return (m3sq-m1sq)/
sqrt(m3sq*m1sq);
41 double m2sq=msQhat2(1,1).abs();
42 double m3sq=msQhat2(2,2).abs();
43 return (m3sq-m2sq)/
sqrt(m3sq*m2sq);
53 double m1sq=msUhat2(0,0).abs();
54 double m2sq=msUhat2(1,1).abs();
55 return (m2sq-m1sq)/
sqrt(m2sq*m1sq);
65 double m1sq=msUhat2(0,0).abs();
66 double m3sq=msUhat2(2,2).abs();
67 return (m3sq-m1sq)/
sqrt(m3sq*m1sq);
77 double m2sq=msUhat2(1,1).abs();
78 double m3sq=msUhat2(2,2).abs();
79 return (m3sq-m2sq)/
sqrt(m3sq*m2sq);
89 double m1sq=msDhat2(0,0).abs();
90 double m2sq=msDhat2(1,1).abs();
91 return (m2sq-m1sq)/
sqrt(m2sq*m1sq);
101 double m1sq=msDhat2(0,0).abs();
102 double m3sq=msDhat2(2,2).abs();
103 return (m3sq-m1sq)/
sqrt(m3sq*m1sq);
113 double m2sq=msDhat2(1,1).abs();
114 double m3sq=msDhat2(2,2).abs();
115 return (m3sq-m2sq)/
sqrt(m3sq*m2sq);
125 double m1sq=msLhat2(0,0).abs();
126 double m2sq=msLhat2(1,1).abs();
127 return (m2sq-m1sq)/
sqrt(m2sq*m1sq);
137 double m1sq=msLhat2(0,0).abs();
138 double m3sq=msLhat2(2,2).abs();
139 return (m3sq-m1sq)/
sqrt(m3sq*m1sq);
149 double m2sq=msLhat2(1,1).abs();
150 double m3sq=msLhat2(2,2).abs();
151 return (m3sq-m2sq)/
sqrt(m3sq*m2sq);
161 double m1sq=msEhat2(0,0).abs();
162 double m2sq=msEhat2(1,1).abs();
163 return (m2sq-m1sq)/
sqrt(m2sq*m1sq);
173 double m1sq=msEhat2(0,0).abs();
174 double m3sq=msEhat2(2,2).abs();
175 return (m3sq-m1sq)/
sqrt(m3sq*m1sq);
185 double m2sq=msEhat2(1,1).abs();
186 double m3sq=msEhat2(2,2).abs();
187 return (m3sq-m2sq)/
sqrt(m3sq*m2sq);
199 double Mu2average_21=
sqrt(msQhat2(1,1)*msUhat2(0,0)).
abs();
202 return TUhat(1,0).real()*v2/Mu2average_21;
214 double Mu2average_31=
sqrt(msQhat2(2,2)*msUhat2(0,0)).
abs();
217 return TUhat(2,0).real()*v2/Mu2average_31;
229 double Mu2average_32=
sqrt(msQhat2(2,2)*msUhat2(1,1)).
abs();
232 return TUhat(2,1).real()*v2/Mu2average_32;
244 double Me2average_21=
sqrt(msLhat2(1,1)*msEhat2(0,0)).
abs();
247 return TEhat(1,0).real()*v1/Me2average_21;
259 double Me2average_12=
sqrt(msLhat2(0,0)*msEhat2(1,1)).
abs();
262 return TEhat(0,1).real()*v1/Me2average_12;
274 double Me2average_31=
sqrt(msLhat2(2,2)*msEhat2(0,0)).
abs();
277 return TEhat(2,0).real()*v1/Me2average_31;
289 double Me2average_13=
sqrt(msLhat2(0,0)*msEhat2(2,2)).
abs();
292 return TEhat(0,2).real()*v1/Me2average_13;
304 double Me2average_32=
sqrt(msLhat2(2,2)*msEhat2(1,1)).
abs();
307 return TEhat(2,1).real()*v1/Me2average_32;
319 double Me2average_23=
sqrt(msLhat2(1,1)*msEhat2(2,2)).
abs();
322 return TEhat(1,2).real()*v1/Me2average_23;
334 double Me2average_31=
sqrt(msLhat2(2,2)*msEhat2(0,0)).
abs();
337 return log10(TEhat(2,0).abs()*v1/Me2average_31);
349 double Me2average_32=
sqrt(msLhat2(2,2)*msEhat2(1,1)).
abs();
352 return log10(TEhat(2,1).abs()*v1/Me2average_32);
363 double Me2average_32=
sqrt(msLhat2(2,2)*msEhat2(1,1)).
abs();
376 double Me2average_32=(
sqrt(msLhat2(2,2)).
abs()+
sqrt(msEhat2(1,1)).
abs())/2.0;
378 return Me2average_32;
388 double denominator =
sqrt(TEhat(2,1).abs()*TEhat(1,2).abs());
390 return (TEhat(2,1).abs()-TEhat(1,2).abs())/denominator;
401 double denominator =
sqrt(msLhat2(2,2).abs()*msEhat2(2,2).abs());
403 return (msLhat2(2,2).abs()-msEhat2(2,2).abs())/denominator;