10 const char* name,
const char* title,
18 _measurables, _parameters,
23 g1Val(
"g1Val",
"g1Val", this, (RooAbsReal&)*(_couplings.g1List[0][0])),
24 g2Val(
"g2Val",
"g2Val", this, (RooAbsReal&)*(_couplings.g2List[0][0])),
25 g3Val(
"g3Val",
"g3Val", this, (RooAbsReal&)*(_couplings.g3List[0][0])),
26 g4Val(
"g4Val",
"g4Val", this, (RooAbsReal&)*(_couplings.g4List[0][0])),
28 g1_primeVal(
"g1_primeVal",
"g1_primeVal", this, (RooAbsReal&)*(_couplings.g1List[1][0])),
29 g2_primeVal(
"g2_primeVal",
"g2_primeVal", this, (RooAbsReal&)*(_couplings.g2List[1][0])),
30 g3_primeVal(
"g3_primeVal",
"g3_primeVal", this, (RooAbsReal&)*(_couplings.g3List[1][0])),
31 g4_primeVal(
"g4_primeVal",
"g4_primeVal", this, (RooAbsReal&)*(_couplings.g4List[1][0])),
33 g1_prime2Val(
"g1_prime2Val",
"g1_prime2Val", this, (RooAbsReal&)*(_couplings.g1List[2][0])),
34 g2_prime2Val(
"g2_prime2Val",
"g2_prime2Val", this, (RooAbsReal&)*(_couplings.g2List[2][0])),
35 g3_prime2Val(
"g3_prime2Val",
"g3_prime2Val", this, (RooAbsReal&)*(_couplings.g3List[2][0])),
36 g4_prime2Val(
"g4_prime2Val",
"g4_prime2Val", this, (RooAbsReal&)*(_couplings.g4List[2][0])),
38 g1_prime3Val(
"g1_prime3Val",
"g1_prime3Val", this, (RooAbsReal&)*(_couplings.g1List[3][0])),
39 g2_prime3Val(
"g2_prime3Val",
"g2_prime3Val", this, (RooAbsReal&)*(_couplings.g2List[3][0])),
40 g3_prime3Val(
"g3_prime3Val",
"g3_prime3Val", this, (RooAbsReal&)*(_couplings.g3List[3][0])),
41 g4_prime3Val(
"g4_prime3Val",
"g4_prime3Val", this, (RooAbsReal&)*(_couplings.g4List[3][0])),
43 g1_prime4Val(
"g1_prime4Val",
"g1_prime4Val", this, (RooAbsReal&)*(_couplings.g1List[4][0])),
44 g2_prime4Val(
"g2_prime4Val",
"g2_prime4Val", this, (RooAbsReal&)*(_couplings.g2List[4][0])),
45 g3_prime4Val(
"g3_prime4Val",
"g3_prime4Val", this, (RooAbsReal&)*(_couplings.g3List[4][0])),
46 g4_prime4Val(
"g4_prime4Val",
"g4_prime4Val", this, (RooAbsReal&)*(_couplings.g4List[4][0])),
48 g1_prime5Val(
"g1_prime5Val",
"g1_prime5Val", this, (RooAbsReal&)*(_couplings.g1List[5][0])),
49 g2_prime5Val(
"g2_prime5Val",
"g2_prime5Val", this, (RooAbsReal&)*(_couplings.g2List[5][0])),
50 g3_prime5Val(
"g3_prime5Val",
"g3_prime5Val", this, (RooAbsReal&)*(_couplings.g3List[5][0])),
51 g4_prime5Val(
"g4_prime5Val",
"g4_prime5Val", this, (RooAbsReal&)*(_couplings.g4List[5][0])),
53 g1_prime6Val(
"g1_prime6Val",
"g1_prime6Val", this, (RooAbsReal&)*(_couplings.g1List[6][0])),
54 g2_prime6Val(
"g2_prime6Val",
"g2_prime6Val", this, (RooAbsReal&)*(_couplings.g2List[6][0])),
55 g3_prime6Val(
"g3_prime6Val",
"g3_prime6Val", this, (RooAbsReal&)*(_couplings.g3List[6][0])),
56 g4_prime6Val(
"g4_prime6Val",
"g4_prime6Val", this, (RooAbsReal&)*(_couplings.g4List[6][0])),
58 g1_prime7Val(
"g1_prime7Val",
"g1_prime7Val", this, (RooAbsReal&)*(_couplings.g1List[7][0])),
59 g2_prime7Val(
"g2_prime7Val",
"g2_prime7Val", this, (RooAbsReal&)*(_couplings.g2List[7][0])),
60 g3_prime7Val(
"g3_prime7Val",
"g3_prime7Val", this, (RooAbsReal&)*(_couplings.g3List[7][0])),
61 g4_prime7Val(
"g4_prime7Val",
"g4_prime7Val", this, (RooAbsReal&)*(_couplings.g4List[7][0])),
63 gzgs1_prime2Val(
"gzgs1_prime2Val",
"gzgs1_prime2Val", this, (RooAbsReal&)*(_couplings.gzgs1List[0][0])),
64 gzgs2Val(
"gzgs2Val",
"gzgs2Val", this, (RooAbsReal&)*(_couplings.gzgs2List[0][0])),
65 gzgs3Val(
"gzgs3Val",
"gzgs3Val", this, (RooAbsReal&)*(_couplings.gzgs3List[0][0])),
66 gzgs4Val(
"gzgs4Val",
"gzgs4Val", this, (RooAbsReal&)*(_couplings.gzgs4List[0][0])),
67 ggsgs2Val(
"ggsgs2Val",
"ggsgs2Val", this, (RooAbsReal&)*(_couplings.ggsgs2List[0][0])),
68 ggsgs3Val(
"ggsgs3Val",
"ggsgs3Val", this, (RooAbsReal&)*(_couplings.ggsgs3List[0][0])),
69 ggsgs4Val(
"ggsgs4Val",
"ggsgs4Val", this, (RooAbsReal&)*(_couplings.ggsgs4List[0][0])),
71 g1ValIm(
"g1ValIm",
"g1ValIm", this, (RooAbsReal&)*(_couplings.g1List[0][1])),
72 g2ValIm(
"g2ValIm",
"g2ValIm", this, (RooAbsReal&)*(_couplings.g2List[0][1])),
73 g3ValIm(
"g3ValIm",
"g3ValIm", this, (RooAbsReal&)*(_couplings.g3List[0][1])),
74 g4ValIm(
"g4ValIm",
"g4ValIm", this, (RooAbsReal&)*(_couplings.g4List[0][1])),
76 g1_primeValIm(
"g1_primeValIm",
"g1_primeValIm", this, (RooAbsReal&)*(_couplings.g1List[1][1])),
77 g2_primeValIm(
"g2_primeValIm",
"g2_primeValIm", this, (RooAbsReal&)*(_couplings.g2List[1][1])),
78 g3_primeValIm(
"g3_primeValIm",
"g3_primeValIm", this, (RooAbsReal&)*(_couplings.g3List[1][1])),
79 g4_primeValIm(
"g4_primeValIm",
"g4_primeValIm", this, (RooAbsReal&)*(_couplings.g4List[1][1])),
81 g1_prime2ValIm(
"g1_prime2ValIm",
"g1_prime2ValIm", this, (RooAbsReal&)*(_couplings.g1List[2][1])),
82 g2_prime2ValIm(
"g2_prime2ValIm",
"g2_prime2ValIm", this, (RooAbsReal&)*(_couplings.g2List[2][1])),
83 g3_prime2ValIm(
"g3_prime2ValIm",
"g3_prime2ValIm", this, (RooAbsReal&)*(_couplings.g3List[2][1])),
84 g4_prime2ValIm(
"g4_prime2ValIm",
"g4_prime2ValIm", this, (RooAbsReal&)*(_couplings.g4List[2][1])),
86 g1_prime3ValIm(
"g1_prime3ValIm",
"g1_prime3ValIm", this, (RooAbsReal&)*(_couplings.g1List[3][1])),
87 g2_prime3ValIm(
"g2_prime3ValIm",
"g2_prime3ValIm", this, (RooAbsReal&)*(_couplings.g2List[3][1])),
88 g3_prime3ValIm(
"g3_prime3ValIm",
"g3_prime3ValIm", this, (RooAbsReal&)*(_couplings.g3List[3][1])),
89 g4_prime3ValIm(
"g4_prime3ValIm",
"g4_prime3ValIm", this, (RooAbsReal&)*(_couplings.g4List[3][1])),
91 g1_prime4ValIm(
"g1_prime4ValIm",
"g1_prime4ValIm", this, (RooAbsReal&)*(_couplings.g1List[4][1])),
92 g2_prime4ValIm(
"g2_prime4ValIm",
"g2_prime4ValIm", this, (RooAbsReal&)*(_couplings.g2List[4][1])),
93 g3_prime4ValIm(
"g3_prime4ValIm",
"g3_prime4ValIm", this, (RooAbsReal&)*(_couplings.g3List[4][1])),
94 g4_prime4ValIm(
"g4_prime4ValIm",
"g4_prime4ValIm", this, (RooAbsReal&)*(_couplings.g4List[4][1])),
96 g1_prime5ValIm(
"g1_prime5ValIm",
"g1_prime5ValIm", this, (RooAbsReal&)*(_couplings.g1List[5][1])),
97 g2_prime5ValIm(
"g2_prime5ValIm",
"g2_prime5ValIm", this, (RooAbsReal&)*(_couplings.g2List[5][1])),
98 g3_prime5ValIm(
"g3_prime5ValIm",
"g3_prime5ValIm", this, (RooAbsReal&)*(_couplings.g3List[5][1])),
99 g4_prime5ValIm(
"g4_prime5ValIm",
"g4_prime5ValIm", this, (RooAbsReal&)*(_couplings.g4List[5][1])),
101 g1_prime6ValIm(
"g1_prime6ValIm",
"g1_prime6ValIm", this, (RooAbsReal&)*(_couplings.g1List[6][1])),
102 g2_prime6ValIm(
"g2_prime6ValIm",
"g2_prime6ValIm", this, (RooAbsReal&)*(_couplings.g2List[6][1])),
103 g3_prime6ValIm(
"g3_prime6ValIm",
"g3_prime6ValIm", this, (RooAbsReal&)*(_couplings.g3List[6][1])),
104 g4_prime6ValIm(
"g4_prime6ValIm",
"g4_prime6ValIm", this, (RooAbsReal&)*(_couplings.g4List[6][1])),
106 g1_prime7ValIm(
"g1_prime7ValIm",
"g1_prime7ValIm", this, (RooAbsReal&)*(_couplings.g1List[7][1])),
107 g2_prime7ValIm(
"g2_prime7ValIm",
"g2_prime7ValIm", this, (RooAbsReal&)*(_couplings.g2List[7][1])),
108 g3_prime7ValIm(
"g3_prime7ValIm",
"g3_prime7ValIm", this, (RooAbsReal&)*(_couplings.g3List[7][1])),
109 g4_prime7ValIm(
"g4_prime7ValIm",
"g4_prime7ValIm", this, (RooAbsReal&)*(_couplings.g4List[7][1])),
111 gzgs1_prime2ValIm(
"gzgs1_prime2ValIm",
"gzgs1_prime2ValIm", this, (RooAbsReal&)*(_couplings.gzgs1List[0][1])),
112 gzgs2ValIm(
"gzgs2ValIm",
"gzgs2ValIm", this, (RooAbsReal&)*(_couplings.gzgs2List[0][1])),
113 gzgs3ValIm(
"gzgs3ValIm",
"gzgs3ValIm", this, (RooAbsReal&)*(_couplings.gzgs3List[0][1])),
114 gzgs4ValIm(
"gzgs4ValIm",
"gzgs4ValIm", this, (RooAbsReal&)*(_couplings.gzgs4List[0][1])),
115 ggsgs2ValIm(
"ggsgs2ValIm",
"ggsgs2ValIm", this, (RooAbsReal&)*(_couplings.ggsgs2List[0][1])),
116 ggsgs3ValIm(
"ggsgs3ValIm",
"ggsgs3ValIm", this, (RooAbsReal&)*(_couplings.ggsgs3List[0][1])),
117 ggsgs4ValIm(
"ggsgs4ValIm",
"ggsgs4ValIm", this, (RooAbsReal&)*(_couplings.ggsgs4List[0][1])),
119 Lambda(
"Lambda",
"Lambda", this, (RooAbsReal&)*(_couplings.Lambda)),
125 Lambda_Q(
"Lambda_Q",
"Lambda_Q", this, (RooAbsReal&)*(_couplings.
Lambda_Q)),
127 Lambda_z11(
"Lambda_z11",
"Lambda_z11", this, (RooAbsReal&)*(_couplings.Lambda_z1qsq[0])),
128 Lambda_z21(
"Lambda_z21",
"Lambda_z21", this, (RooAbsReal&)*(_couplings.Lambda_z2qsq[0])),
129 Lambda_z31(
"Lambda_z31",
"Lambda_z31", this, (RooAbsReal&)*(_couplings.Lambda_z3qsq[0])),
130 Lambda_z41(
"Lambda_z41",
"Lambda_z41", this, (RooAbsReal&)*(_couplings.Lambda_z4qsq[0])),
132 Lambda_z12(
"Lambda_z12",
"Lambda_z12", this, (RooAbsReal&)*(_couplings.Lambda_z1qsq[1])),
133 Lambda_z22(
"Lambda_z22",
"Lambda_z22", this, (RooAbsReal&)*(_couplings.Lambda_z2qsq[1])),
134 Lambda_z32(
"Lambda_z32",
"Lambda_z32", this, (RooAbsReal&)*(_couplings.Lambda_z3qsq[1])),
135 Lambda_z42(
"Lambda_z42",
"Lambda_z42", this, (RooAbsReal&)*(_couplings.Lambda_z4qsq[1])),
137 Lambda_z10(
"Lambda_z10",
"Lambda_z10", this, (RooAbsReal&)*(_couplings.Lambda_z1qsq[2])),
138 Lambda_z20(
"Lambda_z20",
"Lambda_z20", this, (RooAbsReal&)*(_couplings.Lambda_z2qsq[2])),
139 Lambda_z30(
"Lambda_z30",
"Lambda_z30", this, (RooAbsReal&)*(_couplings.Lambda_z3qsq[2])),
140 Lambda_z40(
"Lambda_z40",
"Lambda_z40", this, (RooAbsReal&)*(_couplings.Lambda_z4qsq[2])),
142 cz_q1sq(
"cz_q1sq",
"cz_q1sq", this, (RooAbsReal&)*(_couplings.cLambda_qsq[0])),
143 cz_q2sq(
"cz_q2sq",
"cz_q2sq", this, (RooAbsReal&)*(_couplings.cLambda_qsq[1])),
144 cz_q12sq(
"cz_q12sq",
"cz_q12sq", this, (RooAbsReal&)*(_couplings.cLambda_qsq[2])),
146 gvvp1Val(
"gvvp1Val",
"gvvp1Val", this, (RooAbsReal&)*(_couplings.gvvp1List[0][0])),
147 gvpvp1Val(
"gvpvp1Val",
"gvpvp1Val", this, (RooAbsReal&)*(_couplings.gvpvp1List[0][0])),
149 gvvp1ValIm(
"gvvp1ValIm",
"gvvp1ValIm", this, (RooAbsReal&)*(_couplings.gvvp1List[0][1])),
150 gvpvp1ValIm(
"gvpvp1ValIm",
"gvpvp1ValIm", this, (RooAbsReal&)*(_couplings.gvpvp1List[0][1]))
157 g1Val(
"g1Val", this, other.g1Val),
158 g2Val(
"g2Val", this, other.g2Val),
159 g3Val(
"g3Val", this, other.g3Val),
160 g4Val(
"g4Val", this, other.g4Val),
162 g1_primeVal(
"g1_primeVal", this, other.g1_primeVal),
163 g2_primeVal(
"a2_primeVal", this, other.g2_primeVal),
164 g3_primeVal(
"g3_primeVal", this, other.g3_primeVal),
165 g4_primeVal(
"g4_primeVal", this, other.g4_primeVal),
167 g1_prime2Val(
"g1_prime2Val", this, other.g1_prime2Val),
168 g2_prime2Val(
"a2_prime2Val", this, other.g2_prime2Val),
169 g3_prime2Val(
"g3_prime2Val", this, other.g3_prime2Val),
170 g4_prime2Val(
"g4_prime2Val", this, other.g4_prime2Val),
172 g1_prime3Val(
"g1_prime3Val", this, other.g1_prime3Val),
173 g2_prime3Val(
"a2_prime3Val", this, other.g2_prime3Val),
174 g3_prime3Val(
"g3_prime3Val", this, other.g3_prime3Val),
175 g4_prime3Val(
"g4_prime3Val", this, other.g4_prime3Val),
177 g1_prime4Val(
"g1_prime4Val", this, other.g1_prime4Val),
178 g2_prime4Val(
"a2_prime4Val", this, other.g2_prime4Val),
179 g3_prime4Val(
"g3_prime4Val", this, other.g3_prime4Val),
180 g4_prime4Val(
"g4_prime4Val", this, other.g4_prime4Val),
182 g1_prime5Val(
"g1_prime5Val", this, other.g1_prime5Val),
183 g2_prime5Val(
"a2_prime5Val", this, other.g2_prime5Val),
184 g3_prime5Val(
"g3_prime5Val", this, other.g3_prime5Val),
185 g4_prime5Val(
"g4_prime5Val", this, other.g4_prime5Val),
187 g1_prime6Val(
"g1_prime6Val", this, other.g1_prime6Val),
188 g2_prime6Val(
"a2_prime6Val", this, other.g2_prime6Val),
189 g3_prime6Val(
"g3_prime6Val", this, other.g3_prime6Val),
190 g4_prime6Val(
"g4_prime6Val", this, other.g4_prime6Val),
192 g1_prime7Val(
"g1_prime7Val", this, other.g1_prime7Val),
193 g2_prime7Val(
"a2_prime7Val", this, other.g2_prime7Val),
194 g3_prime7Val(
"g3_prime7Val", this, other.g3_prime7Val),
195 g4_prime7Val(
"g4_prime7Val", this, other.g4_prime7Val),
197 gzgs1_prime2Val(
"gzgs1_prime2Val", this, other.gzgs1_prime2Val),
198 gzgs2Val(
"gzgs2Val", this, other.gzgs2Val),
199 gzgs3Val(
"gzgs3Val", this, other.gzgs3Val),
200 gzgs4Val(
"gzgs4Val", this, other.gzgs4Val),
201 ggsgs2Val(
"ggsgs2Val", this, other.ggsgs2Val),
202 ggsgs3Val(
"ggsgs3Val", this, other.ggsgs3Val),
203 ggsgs4Val(
"ggsgs4Val", this, other.ggsgs4Val),
206 g1ValIm(
"g1ValIm", this, other.g1ValIm),
207 g2ValIm(
"g2ValIm", this, other.g2ValIm),
208 g3ValIm(
"g3ValIm", this, other.g3ValIm),
209 g4ValIm(
"g4ValIm", this, other.g4ValIm),
211 g1_primeValIm(
"g1_primeValIm", this, other.g1_primeValIm),
212 g2_primeValIm(
"a2_primeValIm", this, other.g2_primeValIm),
213 g3_primeValIm(
"g3_primeValIm", this, other.g3_primeValIm),
214 g4_primeValIm(
"g4_primeValIm", this, other.g4_primeValIm),
216 g1_prime2ValIm(
"g1_prime2ValIm", this, other.g1_prime2ValIm),
217 g2_prime2ValIm(
"a2_prime2ValIm", this, other.g2_prime2ValIm),
218 g3_prime2ValIm(
"g3_prime2ValIm", this, other.g3_prime2ValIm),
219 g4_prime2ValIm(
"g4_prime2ValIm", this, other.g4_prime2ValIm),
221 g1_prime3ValIm(
"g1_prime3ValIm", this, other.g1_prime3ValIm),
222 g2_prime3ValIm(
"a2_prime3ValIm", this, other.g2_prime3ValIm),
223 g3_prime3ValIm(
"g3_prime3ValIm", this, other.g3_prime3ValIm),
224 g4_prime3ValIm(
"g4_prime3ValIm", this, other.g4_prime3ValIm),
226 g1_prime4ValIm(
"g1_prime4ValIm", this, other.g1_prime4ValIm),
227 g2_prime4ValIm(
"a2_prime4ValIm", this, other.g2_prime4ValIm),
228 g3_prime4ValIm(
"g3_prime4ValIm", this, other.g3_prime4ValIm),
229 g4_prime4ValIm(
"g4_prime4ValIm", this, other.g4_prime4ValIm),
231 g1_prime5ValIm(
"g1_prime5ValIm", this, other.g1_prime5ValIm),
232 g2_prime5ValIm(
"a2_prime5ValIm", this, other.g2_prime5ValIm),
233 g3_prime5ValIm(
"g3_prime5ValIm", this, other.g3_prime5ValIm),
234 g4_prime5ValIm(
"g4_prime5ValIm", this, other.g4_prime5ValIm),
236 g1_prime6ValIm(
"g1_prime6ValIm", this, other.g1_prime6ValIm),
237 g2_prime6ValIm(
"a2_prime6ValIm", this, other.g2_prime6ValIm),
238 g3_prime6ValIm(
"g3_prime6ValIm", this, other.g3_prime6ValIm),
239 g4_prime6ValIm(
"g4_prime6ValIm", this, other.g4_prime6ValIm),
241 g1_prime7ValIm(
"g1_prime7ValIm", this, other.g1_prime7ValIm),
242 g2_prime7ValIm(
"a2_prime7ValIm", this, other.g2_prime7ValIm),
243 g3_prime7ValIm(
"g3_prime7ValIm", this, other.g3_prime7ValIm),
244 g4_prime7ValIm(
"g4_prime7ValIm", this, other.g4_prime7ValIm),
246 gzgs1_prime2ValIm(
"gzgs1_prime2ValIm", this, other.gzgs1_prime2ValIm),
247 gzgs2ValIm(
"gzgs2ValIm", this, other.gzgs2ValIm),
248 gzgs3ValIm(
"gzgs3ValIm", this, other.gzgs3ValIm),
249 gzgs4ValIm(
"gzgs4ValIm", this, other.gzgs4ValIm),
250 ggsgs2ValIm(
"ggsgs2ValIm", this, other.ggsgs2ValIm),
251 ggsgs3ValIm(
"ggsgs3ValIm", this, other.ggsgs3ValIm),
252 ggsgs4ValIm(
"ggsgs4ValIm", this, other.ggsgs4ValIm),
254 Lambda(
"Lambda", this, other.Lambda),
281 gvvp1Val(
"gvvp1Val", this, other.gvvp1Val),
282 gvpvp1Val(
"gvpvp1Val", this, other.gvpvp1Val),
284 gvvp1ValIm(
"gvvp1ValIm", this, other.gvvp1ValIm),
285 gvpvp1ValIm(
"gvpvp1ValIm", this, other.gvpvp1ValIm)
289 Double_t& a1Re, Double_t& a1Im, Double_t& a2Re, Double_t& a2Im, Double_t& a3Re, Double_t& a3Im,
290 int VGammaVpmode1,
int VGammaVpmode2
298 Double_t s = (pow(
m12, 2) - pow(m1_, 2) - pow(m2_, 2))/2.;
299 if (m1_>m2_+
m12 || m2_>m1_+
m12) s = -s;
311 if (VGammaVpmode1==0 && VGammaVpmode2==0){
417 else if ((VGammaVpmode1==0 && VGammaVpmode2==2) || (VGammaVpmode1==2 && VGammaVpmode2==0)){
421 else if (VGammaVpmode1==2 && VGammaVpmode2==2){
457 a3Re = -2.*g4_dyn*pow(
m12, 2);
458 a2Re = -(2.*g2_dyn + g3_dyn*
kappa)*pow(
m12, 2);
459 a1Re = g1_dyn*pow(mV, 2) - a2Re*s/pow(
m12, 2);
460 a3Im = -2.*g4_dynIm*pow(
m12, 2);
461 a2Im = -(2.*g2_dynIm + g3_dynIm*
kappa)*pow(
m12, 2);
462 a1Im = g1_dynIm*pow(mV, 2) - a2Im*s/pow(
m12, 2);
465 MELAout <<
"RooSpinZero::calculateAi( "
466 << a1Re <<
" , " << a1Im <<
" , " << a2Re <<
" , " << a2Im <<
" , " << a3Re <<
" , " << a3Im <<
" , " << VGammaVpmode1 <<
" , " << VGammaVpmode2
469 <<
"\t- mV, m1, m2, m12, s = " << mV <<
" , " <<
m1 <<
" , " <<
m2 <<
" , " <<
m12 <<
" , " << s <<
'\n'
470 <<
"\t- g1_dyn = (" << g1_dyn <<
" , " << g1_dynIm <<
")" <<
'\n'
471 <<
"\t- g2_dyn = (" << g2_dyn <<
" , " << g2_dynIm <<
")" <<
'\n'
472 <<
"\t- g3_dyn = (" << g3_dyn <<
" , " << g3_dynIm <<
")" <<
'\n'
473 <<
"\t- g4_dyn = (" << g4_dyn <<
" , " << g4_dynIm <<
")"
478 Double_t& A00Re, Double_t& A00Im, Double_t& AppRe, Double_t& AppIm, Double_t& AmmRe, Double_t& AmmIm,
479 int VGammaVpmode1,
int VGammaVpmode2
481 if (
verbosity>=
TVar::DEBUG)
MELAout <<
"Begin RooSpinZero::calculateAmplitudes( " << VGammaVpmode1 <<
" , " << VGammaVpmode2 <<
" )" << endl;
486 Double_t a1Re, a2Re, a3Re, a1Im, a2Im, a3Im;
487 calculateAi(a1Re, a1Im, a2Re, a2Im, a3Re, a3Im, VGammaVpmode1, VGammaVpmode2);
489 Double_t propV1Re=1, propV2Re=1, propHRe=1;
490 Double_t propV1Im=0, propV2Im=0, propHIm=0;
493 switch (VGammaVpmode1){
510 switch (VGammaVpmode2){
529 Double_t eta1 = m1_ /
m12;
530 Double_t eta2 = m2_ /
m12;
531 Double_t eta1p2 = 1.;
536 Double_t etas = (1. - pow(eta1, 2) - pow(eta2, 2))/2.;
537 if (pow(eta1+eta2, 2)>1.) etas = -etas;
539 Double_t etasp = pow(etas, 2) - pow(eta1*eta2, 2);
540 if (etasp<0) etasp=0;
543 <<
"\t- eta1, eta2, eta1p2, etasp = "
544 << eta1 <<
" , " << eta2 <<
" , " << eta1p2 <<
" , " << etasp <<
'\n'
545 <<
"\t- ampScale, prop1, prop2 = "
546 << ampScale <<
" , " << propV1Re <<
"," << propV1Im <<
" , " << propV2Re <<
"," << propV2Im
549 Double_t A00Re_tmp=0, A00Im_tmp=0, AppRe_tmp=0, AppIm_tmp=0, AmmRe_tmp=0, AmmIm_tmp=0;
551 A00Re_tmp = -2.*(a1Re*etas + a2Re*etasp)*ampScale;
552 A00Im_tmp = -2.*(a1Im*etas + a2Im*etasp)*ampScale;
554 AppRe_tmp = (a1Re - a3Im*sqrt(etasp))*ampScale;
555 AppIm_tmp = (a1Im + a3Re*sqrt(etasp))*ampScale;
556 AmmRe_tmp = (a1Re + a3Im*sqrt(etasp))*ampScale;
557 AmmIm_tmp = (a1Im - a3Re*sqrt(etasp))*ampScale;
565 <<
"\t- A00 factor = " << A00Re_tmp <<
" , " << A00Im_tmp <<
'\n'
566 <<
"\t- A++ factor = " << AppRe_tmp <<
" , " << AppIm_tmp <<
'\n'
567 <<
"\t- A-- factor = " << AmmRe_tmp <<
" , " << AmmIm_tmp
571 std::vector<Double_t> A00_reals, A00_imags; A00_reals.push_back(A00Re_tmp); A00_imags.push_back(A00Im_tmp); A00_reals.push_back(propV1Re); A00_imags.push_back(propV1Im); A00_reals.push_back(propV2Re); A00_imags.push_back(propV2Im); A00_reals.push_back(propHRe); A00_imags.push_back(propHIm);
AnaMelaHelpers::multiplyComplexNumbers(A00_reals, A00_imags, A00Re, A00Im);
572 std::vector<Double_t> App_reals, App_imags; App_reals.push_back(AppRe_tmp); App_imags.push_back(AppIm_tmp); App_reals.push_back(propV1Re); App_imags.push_back(propV1Im); App_reals.push_back(propV2Re); App_imags.push_back(propV2Im); App_reals.push_back(propHRe); App_imags.push_back(propHIm);
AnaMelaHelpers::multiplyComplexNumbers(App_reals, App_imags, AppRe, AppIm);
573 std::vector<Double_t> Amm_reals, Amm_imags; Amm_reals.push_back(AmmRe_tmp); Amm_imags.push_back(AmmIm_tmp); Amm_reals.push_back(propV1Re); Amm_imags.push_back(propV1Im); Amm_reals.push_back(propV2Re); Amm_imags.push_back(propV2Im); Amm_reals.push_back(propHRe); Amm_imags.push_back(propHIm);
AnaMelaHelpers::multiplyComplexNumbers(Amm_reals, Amm_imags, AmmRe, AmmIm);
575 if (!(A00Re==A00Re) || !(AppRe==AppRe) || !(AmmRe==AmmRe) || !(A00Im==A00Im) || !(AppIm==AppIm) || !(AmmIm==AmmIm)){
576 MELAout <<
"A00 = " << A00Re <<
", " << A00Im << endl;
577 MELAout <<
"A++ = " << AppRe <<
", " << AppIm << endl;
578 MELAout <<
"A-- = " << AmmRe <<
", " << AmmIm << endl;
579 MELAout << eta1 <<
'\t' << eta2 <<
'\t' << etas <<
'\t' << etasp <<
'\t' << eta1p2 << endl;
586 (VGammaVpmode1==VGammaVpmode2 && VGammaVpmode1==0) && (
656 ((VGammaVpmode1==0 && VGammaVpmode2==1) || (VGammaVpmode2==0 && VGammaVpmode1==1)) && (
670 (VGammaVpmode1==VGammaVpmode2 && VGammaVpmode1==1) && (
682 ((VGammaVpmode1==0 && VGammaVpmode2==2) || (VGammaVpmode2==0 && VGammaVpmode1==2)) && (
690 (VGammaVpmode1==VGammaVpmode2 && VGammaVpmode1==2) && (