JHUGen MELA  v2.4.1
Matrix element calculations as used in JHUGen. MELA is an important tool that was used for the Higgs boson discovery and for precise measurements of its structure and interactions. Please see the website https://spin.pha.jhu.edu/ and papers cited there for more details, and kindly cite those papers when using this code.
RooSpinZero.cc
Go to the documentation of this file.
1 #include "RooSpinZero.h"
2 
3 
4 using namespace std;
5 using namespace MELAStreamHelpers;
6 
7 
8 RooSpinZero::RooSpinZero() : RooSpin(){}
9 RooSpinZero::RooSpinZero(
10  const char* name, const char* title,
11  modelMeasurables const& _measurables,
12  modelParameters const& _parameters,
13  modelCouplings const& _couplings,
14  RooSpin::VdecayType _Vdecay1, RooSpin::VdecayType _Vdecay2,
15  TVar::VerbosityLevel verbosity_
16 ) : RooSpin(
17  name, title,
18  _measurables, _parameters,
19  _Vdecay1, _Vdecay2,
20  verbosity_
21  ),
22 
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])),
27 
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])),
32 
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])),
37 
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])),
42 
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])),
47 
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])),
52 
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])),
57 
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])),
62 
63  gzgs1_prime2Val("gzgs1_prime2Val", "gzgs1_prime2Val", this, (RooAbsReal&)*(_couplings.gzgs1List[0][0])), // Special case!
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])),
70 
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])),
75 
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])),
80 
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])),
85 
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])),
90 
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])),
95 
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])),
100 
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])),
105 
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])),
110 
111  gzgs1_prime2ValIm("gzgs1_prime2ValIm", "gzgs1_prime2ValIm", this, (RooAbsReal&)*(_couplings.gzgs1List[0][1])), // Special case!
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])),
118 
119  Lambda("Lambda", "Lambda", this, (RooAbsReal&)*(_couplings.Lambda)),
120  Lambda_zgs1("Lambda_zgs1", "Lambda_zgs1", this, (RooAbsReal&)*(_couplings.Lambda_zgs1)),
121  Lambda_z1("Lambda_z1", "Lambda_z1", this, (RooAbsReal&)*(_couplings.Lambda_z1)),
122  Lambda_z2("Lambda_z2", "Lambda_z2", this, (RooAbsReal&)*(_couplings.Lambda_z2)),
123  Lambda_z3("Lambda_z3", "Lambda_z3", this, (RooAbsReal&)*(_couplings.Lambda_z3)),
124  Lambda_z4("Lambda_z4", "Lambda_z4", this, (RooAbsReal&)*(_couplings.Lambda_z4)),
125  Lambda_Q("Lambda_Q", "Lambda_Q", this, (RooAbsReal&)*(_couplings.Lambda_Q)),
126 
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])),
131 
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])),
136 
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])),
141 
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])),
145 
146  gvvp1Val("gvvp1Val", "gvvp1Val", this, (RooAbsReal&)*(_couplings.gvvp1List[0][0])),
147  gvpvp1Val("gvpvp1Val", "gvpvp1Val", this, (RooAbsReal&)*(_couplings.gvpvp1List[0][0])),
148 
149  gvvp1ValIm("gvvp1ValIm", "gvvp1ValIm", this, (RooAbsReal&)*(_couplings.gvvp1List[0][1])),
150  gvpvp1ValIm("gvpvp1ValIm", "gvpvp1ValIm", this, (RooAbsReal&)*(_couplings.gvpvp1List[0][1]))
151 {}
152 
153 
154 RooSpinZero::RooSpinZero(const RooSpinZero& other, const char* name) :
155 RooSpin(other, name),
156 
157 g1Val("g1Val", this, other.g1Val),
158 g2Val("g2Val", this, other.g2Val),
159 g3Val("g3Val", this, other.g3Val),
160 g4Val("g4Val", this, other.g4Val),
161 
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),
166 
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),
171 
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),
176 
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),
181 
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),
186 
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),
191 
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),
196 
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),
204 
205 
206 g1ValIm("g1ValIm", this, other.g1ValIm),
207 g2ValIm("g2ValIm", this, other.g2ValIm),
208 g3ValIm("g3ValIm", this, other.g3ValIm),
209 g4ValIm("g4ValIm", this, other.g4ValIm),
210 
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),
215 
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),
220 
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),
225 
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),
230 
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),
235 
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),
240 
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),
245 
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),
253 
254 Lambda("Lambda", this, other.Lambda),
255 Lambda_zgs1("Lambda_zgs1", this, other.Lambda_zgs1),
256 Lambda_z1("Lambda_z1", this, other.Lambda_z1),
257 Lambda_z2("Lambda_z2", this, other.Lambda_z2),
258 Lambda_z3("Lambda_z3", this, other.Lambda_z3),
259 Lambda_z4("Lambda_z4", this, other.Lambda_z4),
260 Lambda_Q("Lambda_Q", this, other.Lambda_Q),
261 
262 Lambda_z11("Lambda_z11", this, other.Lambda_z11),
263 Lambda_z21("Lambda_z21", this, other.Lambda_z21),
264 Lambda_z31("Lambda_z31", this, other.Lambda_z31),
265 Lambda_z41("Lambda_z41", this, other.Lambda_z41),
266 
267 Lambda_z12("Lambda_z12", this, other.Lambda_z12),
268 Lambda_z22("Lambda_z22", this, other.Lambda_z22),
269 Lambda_z32("Lambda_z32", this, other.Lambda_z32),
270 Lambda_z42("Lambda_z42", this, other.Lambda_z42),
271 
272 Lambda_z10("Lambda_z10", this, other.Lambda_z10),
273 Lambda_z20("Lambda_z20", this, other.Lambda_z20),
274 Lambda_z30("Lambda_z30", this, other.Lambda_z30),
275 Lambda_z40("Lambda_z40", this, other.Lambda_z40),
276 
277 cz_q1sq("cz_q1sq", this, other.cz_q1sq),
278 cz_q2sq("cz_q2sq", this, other.cz_q2sq),
279 cz_q12sq("cz_q12sq", this, other.cz_q12sq),
280 
281 gvvp1Val("gvvp1Val", this, other.gvvp1Val),
282 gvpvp1Val("gvpvp1Val", this, other.gvpvp1Val),
283 
284 gvvp1ValIm("gvvp1ValIm", this, other.gvvp1ValIm),
285 gvpvp1ValIm("gvpvp1ValIm", this, other.gvpvp1ValIm)
286 {}
287 
288 void RooSpinZero::calculateAi(
289  Double_t& a1Re, Double_t& a1Im, Double_t& a2Re, Double_t& a2Im, Double_t& a3Re, Double_t& a3Im,
290  int VGammaVpmode1, int VGammaVpmode2
291 )const{
292  Double_t mV;
293  getMVGamV(&mV);
294 
295  Double_t m1_=m1; if (Vdecay1==RooSpin::kVdecayType_GammaOnshell) m1_=0;
296  Double_t m2_=m2; if (Vdecay2==RooSpin::kVdecayType_GammaOnshell) m2_=0;
297 
298  Double_t s = (pow(m12, 2) - pow(m1_, 2) - pow(m2_, 2))/2.;
299  if (m1_>m2_+m12 || m2_>m1_+m12) s = -s;
300 
301  Double_t g1_dyn=0;
302  Double_t g2_dyn=0;
303  Double_t g3_dyn=0;
304  Double_t g4_dyn=0;
305  Double_t g1_dynIm=0;
306  Double_t g2_dynIm=0;
307  Double_t g3_dynIm=0;
308  Double_t g4_dynIm=0;
309 
310  if (Vdecay1!=RooSpin::kVdecayType_GammaOnshell && Vdecay2!=RooSpin::kVdecayType_GammaOnshell){ // VV/VVp/VpVp
311  if (VGammaVpmode1==0 && VGammaVpmode2==0){ // VV
312  if (g1_primeVal!=0) g1_dyn += g1_primeVal * pow(Lambda_z1, 4)/(pow(Lambda_z1, 2) + pow(m1_, 2))/(pow(Lambda_z1, 2) + pow(m2_, 2));
313  if (g1_prime2Val!=0) g1_dyn += g1_prime2Val* (pow(m1_, 2) + pow(m2_, 2))/pow(Lambda_z1, 2);
314  if (g1_prime3Val!=0) g1_dyn += g1_prime3Val* (pow(m1_, 2) - pow(m2_, 2))/pow(Lambda_z1, 2);
315  if (g1_prime4Val!=0) g1_dyn += g1_prime4Val* (m12*m12)/pow(Lambda_Q, 2);
316  if (g1_prime5Val!=0) g1_dyn += g1_prime5Val* (pow(m1_, 4) + pow(m2_, 4))/pow(Lambda_z1, 4);
317  if (g1_prime6Val!=0) g1_dyn += g1_prime6Val* (pow(m1_, 4) - pow(m2_, 4))/pow(Lambda_z1, 4);
318  if (g1_prime7Val!=0) g1_dyn += g1_prime7Val* (pow(m1_, 2) * pow(m2_, 2))/pow(Lambda_z1, 4);
319 
320  if (cz_q1sq!=0.) g1_dyn *= 1./(1.+ cz_q1sq*pow(m1_/Lambda_z11, 2));
321  if (cz_q2sq!=0.) g1_dyn *= 1./(1.+ cz_q2sq*pow(m2_/Lambda_z12, 2));
322  if (cz_q12sq!=0.) g1_dyn *= 1./(1.+ cz_q12sq*pow(m12/Lambda_z10, 2));
323  g1_dyn += g1Val;
324 
325  g2_dyn = g2Val;
326  if (g2_primeVal!=0) g2_dyn += g2_primeVal * pow(Lambda_z2, 4)/(pow(Lambda_z2, 2) + pow(m1_, 2))/(pow(Lambda_z2, 2) + pow(m2_, 2));
327  if (g2_prime2Val!=0) g2_dyn += g2_prime2Val* (pow(m1_, 2) + pow(m2_, 2))/pow(Lambda_z2, 2);
328  if (g2_prime3Val!=0) g2_dyn += g2_prime3Val* (pow(m1_, 2) - pow(m2_, 2))/pow(Lambda_z2, 2);
329  if (g2_prime4Val!=0) g2_dyn += g2_prime4Val* (m12*m12)/pow(Lambda_Q, 2);
330  if (g2_prime5Val!=0) g2_dyn += g2_prime5Val* (pow(m1_, 4) + pow(m2_, 4))/pow(Lambda_z2, 4);
331  if (g2_prime6Val!=0) g2_dyn += g2_prime6Val* (pow(m1_, 4) - pow(m2_, 4))/pow(Lambda_z2, 4);
332  if (g2_prime7Val!=0) g2_dyn += g2_prime7Val* (pow(m1_, 2) * pow(m2_, 2))/pow(Lambda_z2, 4);
333 
334  if (cz_q1sq!=0.) g2_dyn *= 1./(1.+ cz_q1sq*pow(m1_/Lambda_z21, 2));
335  if (cz_q2sq!=0.) g2_dyn *= 1./(1.+ cz_q2sq*pow(m2_/Lambda_z22, 2));
336  if (cz_q12sq!=0.) g2_dyn *= 1./(1.+ cz_q12sq*pow(m12/Lambda_z20, 2));
337 
338  g3_dyn = g3Val;
339  if (g3_primeVal!=0) g3_dyn += g3_primeVal * pow(Lambda_z3, 4)/(pow(Lambda_z3, 2) + pow(m1_, 2))/(pow(Lambda_z3, 2) + pow(m2_, 2));
340  if (g3_prime2Val!=0) g3_dyn += g3_prime2Val* (pow(m1_, 2) + pow(m2_, 2))/pow(Lambda_z3, 2);
341  if (g3_prime3Val!=0) g3_dyn += g3_prime3Val* (pow(m1_, 2) - pow(m2_, 2))/pow(Lambda_z3, 2);
342  if (g3_prime4Val!=0) g3_dyn += g3_prime4Val* (m12*m12)/pow(Lambda_Q, 2);
343  if (g3_prime5Val!=0) g3_dyn += g3_prime5Val* (pow(m1_, 4) + pow(m2_, 4))/pow(Lambda_z3, 4);
344  if (g3_prime6Val!=0) g3_dyn += g3_prime6Val* (pow(m1_, 4) - pow(m2_, 4))/pow(Lambda_z3, 4);
345  if (g3_prime7Val!=0) g3_dyn += g3_prime7Val* (pow(m1_, 2) * pow(m2_, 2))/pow(Lambda_z3, 4);
346 
347  if (cz_q1sq!=0.) g3_dyn *= 1./(1.+ cz_q1sq*pow(m1_/Lambda_z31, 2));
348  if (cz_q2sq!=0.) g3_dyn *= 1./(1.+ cz_q2sq*pow(m2_/Lambda_z32, 2));
349  if (cz_q12sq!=0.) g3_dyn *= 1./(1.+ cz_q12sq*pow(m12/Lambda_z30, 2));
350 
351  g4_dyn = g4Val;
352  if (g4_primeVal!=0) g4_dyn += g4_primeVal * pow(Lambda_z4, 4)/(pow(Lambda_z4, 2) + pow(m1_, 2))/(pow(Lambda_z4, 2) + pow(m2_, 2));
353  if (g4_prime2Val!=0) g4_dyn += g4_prime2Val* (pow(m1_, 2) + pow(m2_, 2))/pow(Lambda_z4, 2);
354  if (g4_prime3Val!=0) g4_dyn += g4_prime3Val* (pow(m1_, 2) - pow(m2_, 2))/pow(Lambda_z4, 2);
355  if (g4_prime4Val!=0) g4_dyn += g4_prime4Val* (m12*m12)/pow(Lambda_Q, 2);
356  if (g4_prime5Val!=0) g4_dyn += g4_prime5Val* (pow(m1_, 4) + pow(m2_, 4))/pow(Lambda_z4, 4);
357  if (g4_prime6Val!=0) g4_dyn += g4_prime6Val* (pow(m1_, 4) - pow(m2_, 4))/pow(Lambda_z4, 4);
358  if (g4_prime7Val!=0) g4_dyn += g4_prime7Val* (pow(m1_, 2) * pow(m2_, 2))/pow(Lambda_z4, 4);
359 
360  if (cz_q1sq!=0.) g4_dyn *= 1./(1.+ cz_q1sq*pow(m1_/Lambda_z41, 2));
361  if (cz_q2sq!=0.) g4_dyn *= 1./(1.+ cz_q2sq*pow(m2_/Lambda_z42, 2));
362  if (cz_q12sq!=0.) g4_dyn *= 1./(1.+ cz_q12sq*pow(m12/Lambda_z40, 2));
363 
364  g1_dynIm = 0;
365  if (g1_primeValIm!=0) g1_dynIm += g1_primeValIm * pow(Lambda_z1, 4)/(pow(Lambda_z1, 2) + pow(m1_, 2))/(pow(Lambda_z1, 2) + pow(m2_, 2));
366  if (g1_prime2ValIm!=0) g1_dynIm += g1_prime2ValIm* (pow(m1_, 2) + pow(m2_, 2))/pow(Lambda_z1, 2);
367  if (g1_prime3ValIm!=0) g1_dynIm += g1_prime3ValIm* (pow(m1_, 2) - pow(m2_, 2))/pow(Lambda_z1, 2);
368  if (g1_prime4ValIm!=0) g1_dynIm += g1_prime4ValIm* (m12*m12)/pow(Lambda_Q, 2);
369  if (g1_prime5ValIm!=0) g1_dynIm += g1_prime5ValIm* (pow(m1_, 4) + pow(m2_, 4))/pow(Lambda_z1, 4);
370  if (g1_prime6ValIm!=0) g1_dynIm += g1_prime6ValIm* (pow(m1_, 4) - pow(m2_, 4))/pow(Lambda_z1, 4);
371  if (g1_prime7ValIm!=0) g1_dynIm += g1_prime7ValIm* (pow(m1_, 2) * pow(m2_, 2))/pow(Lambda_z1, 4);
372 
373  if (cz_q1sq!=0.) g1_dynIm *= 1./(1.+ cz_q1sq*pow(m1_/Lambda_z11, 2));
374  if (cz_q2sq!=0.) g1_dynIm *= 1./(1.+ cz_q2sq*pow(m2_/Lambda_z12, 2));
375  if (cz_q12sq!=0.) g1_dynIm *= 1./(1.+ cz_q12sq*pow(m12/Lambda_z10, 2));
376  g1_dynIm += g1ValIm;
377 
378  g2_dynIm = g2ValIm;
379  if (g2_primeValIm!=0) g2_dynIm += g2_primeValIm * pow(Lambda_z2, 4)/(pow(Lambda_z2, 2) + pow(m1_, 2))/(pow(Lambda_z2, 2) + pow(m2_, 2));
380  if (g2_prime2ValIm!=0) g2_dynIm += g2_prime2ValIm* (pow(m1_, 2) + pow(m2_, 2))/pow(Lambda_z2, 2);
381  if (g2_prime3ValIm!=0) g2_dynIm += g2_prime3ValIm* (pow(m1_, 2) - pow(m2_, 2))/pow(Lambda_z2, 2);
382  if (g2_prime4ValIm!=0) g2_dynIm += g2_prime4ValIm* (m12*m12)/pow(Lambda_Q, 2);
383  if (g2_prime5ValIm!=0) g2_dynIm += g2_prime5ValIm* (pow(m1_, 4) + pow(m2_, 4))/pow(Lambda_z2, 4);
384  if (g2_prime6ValIm!=0) g2_dynIm += g2_prime6ValIm* (pow(m1_, 4) - pow(m2_, 4))/pow(Lambda_z2, 4);
385  if (g2_prime7ValIm!=0) g2_dynIm += g2_prime7ValIm* (pow(m1_, 2) * pow(m2_, 2))/pow(Lambda_z2, 4);
386 
387  if (cz_q1sq!=0.) g2_dynIm *= 1./(1.+ cz_q1sq*pow(m1_/Lambda_z21, 2));
388  if (cz_q2sq!=0.) g2_dynIm *= 1./(1.+ cz_q2sq*pow(m2_/Lambda_z22, 2));
389  if (cz_q12sq!=0.) g2_dynIm *= 1./(1.+ cz_q12sq*pow(m12/Lambda_z20, 2));
390 
391  g3_dynIm = g3ValIm;
392  if (g3_primeValIm!=0) g3_dynIm += g3_primeValIm * pow(Lambda_z3, 4)/(pow(Lambda_z3, 2) + pow(m1_, 2))/(pow(Lambda_z3, 2) + pow(m2_, 2));
393  if (g3_prime2ValIm!=0) g3_dynIm += g3_prime2ValIm* (pow(m1_, 2) + pow(m2_, 2))/pow(Lambda_z3, 2);
394  if (g3_prime3ValIm!=0) g3_dynIm += g3_prime3ValIm* (pow(m1_, 2) - pow(m2_, 2))/pow(Lambda_z3, 2);
395  if (g3_prime4ValIm!=0) g3_dynIm += g3_prime4ValIm* (m12*m12)/pow(Lambda_Q, 2);
396  if (g3_prime5ValIm!=0) g3_dynIm += g3_prime5ValIm* (pow(m1_, 4) + pow(m2_, 4))/pow(Lambda_z3, 4);
397  if (g3_prime6ValIm!=0) g3_dynIm += g3_prime6ValIm* (pow(m1_, 4) - pow(m2_, 4))/pow(Lambda_z3, 4);
398  if (g3_prime7ValIm!=0) g3_dynIm += g3_prime7ValIm* (pow(m1_, 2) * pow(m2_, 2))/pow(Lambda_z3, 4);
399 
400  if (cz_q1sq!=0.) g3_dynIm *= 1./(1.+ cz_q1sq*pow(m1_/Lambda_z31, 2));
401  if (cz_q2sq!=0.) g3_dynIm *= 1./(1.+ cz_q2sq*pow(m2_/Lambda_z32, 2));
402  if (cz_q12sq!=0.) g3_dynIm *= 1./(1.+ cz_q12sq*pow(m12/Lambda_z30, 2));
403 
404  g4_dynIm = g4ValIm;
405  if (g4_primeValIm!=0) g4_dynIm += g4_primeValIm * pow(Lambda_z4, 4)/(pow(Lambda_z4, 2) + pow(m1_, 2))/(pow(Lambda_z4, 2) + pow(m2_, 2));
406  if (g4_prime2ValIm!=0) g4_dynIm += g4_prime2ValIm* (pow(m1_, 2) + pow(m2_, 2))/pow(Lambda_z4, 2);
407  if (g4_prime3ValIm!=0) g4_dynIm += g4_prime3ValIm* (pow(m1_, 2) - pow(m2_, 2))/pow(Lambda_z4, 2);
408  if (g4_prime4ValIm!=0) g4_dynIm += g4_prime4ValIm* (m12*m12)/pow(Lambda_Q, 2);
409  if (g4_prime5ValIm!=0) g4_dynIm += g4_prime5ValIm* (pow(m1_, 4) + pow(m2_, 4))/pow(Lambda_z4, 4);
410  if (g4_prime6ValIm!=0) g4_dynIm += g4_prime6ValIm* (pow(m1_, 4) - pow(m2_, 4))/pow(Lambda_z4, 4);
411  if (g4_prime7ValIm!=0) g4_dynIm += g4_prime7ValIm* (pow(m1_, 2) * pow(m2_, 2))/pow(Lambda_z4, 4);
412 
413  if (cz_q1sq!=0.) g4_dynIm *= 1./(1.+ cz_q1sq*pow(m1_/Lambda_z41, 2));
414  if (cz_q2sq!=0.) g4_dynIm *= 1./(1.+ cz_q2sq*pow(m2_/Lambda_z42, 2));
415  if (cz_q12sq!=0.) g4_dynIm *= 1./(1.+ cz_q12sq*pow(m12/Lambda_z40, 2));
416  }
417  else if ((VGammaVpmode1==0 && VGammaVpmode2==2) || (VGammaVpmode1==2 && VGammaVpmode2==0)){ // VVp/VpV
418  g1_dyn = gvvp1Val;
419  g1_dynIm = gvvp1ValIm;
420  }
421  else if (VGammaVpmode1==2 && VGammaVpmode2==2){ // VpVp
422  g1_dyn = gvpvp1Val;
423  g1_dynIm = gvpvp1ValIm;
424  }
425  }
426  else if (!(Vdecay1==RooSpin::kVdecayType_GammaOnshell && Vdecay2==RooSpin::kVdecayType_GammaOnshell) && (VGammaVpmode1+VGammaVpmode2)==1){ // Zgs/Zg
427  if (gzgs1_prime2Val!=0 && VGammaVpmode1==0) g1_dyn += gzgs1_prime2Val* pow(m1_, 2)/pow(Lambda_zgs1, 2);
428  if (gzgs1_prime2Val!=0 && VGammaVpmode2==0) g1_dyn += gzgs1_prime2Val* pow(m2_, 2)/pow(Lambda_zgs1, 2);
429  g2_dyn = gzgs2Val;
430  g3_dyn = gzgs3Val;
431  g4_dyn = gzgs4Val;
432 
433  if (gzgs1_prime2ValIm!=0 && VGammaVpmode1==0) g1_dynIm += gzgs1_prime2ValIm* pow(m1_, 2)/pow(Lambda_zgs1, 2);
434  if (gzgs1_prime2ValIm!=0 && VGammaVpmode2==0) g1_dynIm += gzgs1_prime2ValIm* pow(m2_, 2)/pow(Lambda_zgs1, 2);
435  g2_dynIm = gzgs2ValIm;
436  g3_dynIm = gzgs3ValIm;
437  g4_dynIm = gzgs4ValIm;
438  }
439  else{ // gsgs/gsg/gg
440  g2_dyn = ggsgs2Val;
441  g3_dyn = ggsgs3Val;
442  g4_dyn = ggsgs4Val;
443 
444  g2_dynIm = ggsgs2ValIm;
445  g3_dynIm = ggsgs3ValIm;
446  g4_dynIm = ggsgs4ValIm;
447  }
448  /*
449  MELAout << "g1: " << g1_dyn << " " << g1_dynIm << '\t';
450  MELAout << "g2: " << g2_dyn << " " << g2_dynIm << '\t';
451  MELAout << "g3: " << g3_dyn << " " << g3_dynIm << '\t';
452  MELAout << "g4: " << g4_dyn << " " << g4_dynIm << '\t';
453  MELAout << endl;
454  */
455 
456  Double_t kappa = s/pow(Lambda, 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);
463 
464  if (verbosity>=TVar::DEBUG){
465  MELAout << "RooSpinZero::calculateAi( "
466  << a1Re << " , " << a1Im << " , " << a2Re << " , " << a2Im << " , " << a3Re << " , " << a3Im << " , " << VGammaVpmode1 << " , " << VGammaVpmode2
467  << " ):\n"
468  << "\t- Input:\n"
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 << ")"
474  << endl;
475  }
476 }
477 void RooSpinZero::calculateAmplitudes(
478  Double_t& A00Re, Double_t& A00Im, Double_t& AppRe, Double_t& AppIm, Double_t& AmmRe, Double_t& AmmIm,
479  int VGammaVpmode1, int VGammaVpmode2
480 )const{
481  if (verbosity>=TVar::DEBUG) MELAout << "Begin RooSpinZero::calculateAmplitudes( " << VGammaVpmode1 << " , " << VGammaVpmode2 << " )" << endl;
482 
483  Double_t m1_=m1; if (Vdecay1==RooSpin::kVdecayType_GammaOnshell) m1_=0;
484  Double_t m2_=m2; if (Vdecay2==RooSpin::kVdecayType_GammaOnshell) m2_=0;
485 
486  Double_t a1Re, a2Re, a3Re, a1Im, a2Im, a3Im;
487  calculateAi(a1Re, a1Im, a2Re, a2Im, a3Re, a3Im, VGammaVpmode1, VGammaVpmode2);
488 
489  Double_t propV1Re=1, propV2Re=1, propHRe=1;
490  Double_t propV1Im=0, propV2Im=0, propHIm=0;
491  if (Vdecay1!=RooSpin::kVdecayType_GammaOnshell){
492  int proptype=-1;
493  switch (VGammaVpmode1){
494  case 0:
495  proptype=1;
496  break;
497  case 1:
498  proptype=0;
499  break;
500  case 2:
501  proptype=3;
502  break;
503  default:
504  break;
505  }
506  calculatePropagator(propV1Re, propV1Im, m1_, proptype);
507  }
508  if (Vdecay2!=RooSpin::kVdecayType_GammaOnshell){
509  int proptype=-1;
510  switch (VGammaVpmode2){
511  case 0:
512  proptype=1;
513  break;
514  case 1:
515  proptype=0;
516  break;
517  case 2:
518  proptype=3;
519  break;
520  default:
521  break;
522  }
523  calculatePropagator(propV2Re, propV2Im, m2_, proptype);
524  }
525  calculatePropagator(propHRe, propHIm, m12, 2);
526 
527  Double_t ampScale = calculateAmplitudeScale(VGammaVpmode2, VGammaVpmode2)*pow(GeVunit, 2);
528 
529  Double_t eta1 = m1_ / m12;
530  Double_t eta2 = m2_ / m12;
531  Double_t eta1p2 = 1.;
532  if (Vdecay1!=RooSpin::kVdecayType_GammaOnshell) eta1p2 *= eta1;
533  if (Vdecay2!=RooSpin::kVdecayType_GammaOnshell) eta1p2 *= eta2;
534 
535  // etas = s/m12**2
536  Double_t etas = (1. - pow(eta1, 2) - pow(eta2, 2))/2.;
537  if (pow(eta1+eta2, 2)>1.) etas = -etas;
538  // x = (s/(m1*m2))**2 - 1 = (etas**2 - (eta1*eta2)**2)/(eta1*eta2)**2
539  Double_t etasp = pow(etas, 2) - pow(eta1*eta2, 2); // Notice how eta1p2 is not used. The second set of multiplications below is the reason to it!
540  if (etasp<0) etasp=0;
541 
542  if (verbosity>=TVar::DEBUG) MELAout << "RooSpinZero::calculateAmplitudes:\n"
543  << "\t- eta1, eta2, eta1p2, etasp = "
544  << eta1 << " , " << eta2 << " , " << eta1p2 << " , " << etasp << '\n'
545  << "\t- ampScale, prop1, prop2 = "
546  << ampScale << " , " << propV1Re << "," << propV1Im << " , " << propV2Re << "," << propV2Im
547  << endl;
548 
549  Double_t A00Re_tmp=0, A00Im_tmp=0, AppRe_tmp=0, AppIm_tmp=0, AmmRe_tmp=0, AmmIm_tmp=0;
550  if (Vdecay1!=RooSpin::kVdecayType_GammaOnshell && Vdecay2!=RooSpin::kVdecayType_GammaOnshell){
551  A00Re_tmp = -2.*(a1Re*etas + a2Re*etasp)*ampScale;
552  A00Im_tmp = -2.*(a1Im*etas + a2Im*etasp)*ampScale;
553  }
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;
558 
559  AppRe_tmp *= eta1p2;
560  AppIm_tmp *= eta1p2;
561  AmmRe_tmp *= eta1p2;
562  AmmIm_tmp *= eta1p2;
563 
564  if (verbosity>=TVar::DEBUG) MELAout << "RooSpinZero::calculateAmplitudes:\n"
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
568  << endl;
569 
570  // A_old = ARe_old+i*AIm_old => A_new = ARe_new + i*AIm_new = A_old*propV1*propV2
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);
574 
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;
580  }
581 }
582 
583 Bool_t RooSpinZero::computeNeededAmplitude(int VGammaVpmode1, int VGammaVpmode2) const{
584  return (
585  (
586  (VGammaVpmode1==VGammaVpmode2 && VGammaVpmode1==0) && (
587  g1Val!=0.
588  || g1_primeVal!=0.
589  || g1_prime2Val!=0.
590  || g1_prime3Val!=0.
591  || g1_prime4Val!=0.
592  || g1_prime5Val!=0.
593  || g1_prime6Val!=0.
594  || g1_prime7Val!=0.
595  || g2Val!=0.
596  || g2_primeVal!=0.
597  || g2_prime2Val!=0.
598  || g2_prime3Val!=0.
599  || g2_prime4Val!=0.
600  || g2_prime5Val!=0.
601  || g2_prime6Val!=0.
602  || g2_prime7Val!=0.
603  || g4Val!=0.
604  || g4_primeVal!=0.
605  || g4_prime2Val!=0.
606  || g4_prime3Val!=0.
607  || g4_prime4Val!=0.
608  || g4_prime5Val!=0.
609  || g4_prime6Val!=0.
610  || g4_prime7Val!=0.
611  || g3Val!=0.
612  || g3_primeVal!=0.
613  || g3_prime2Val!=0.
614  || g3_prime3Val!=0.
615  || g3_prime4Val!=0.
616  || g3_prime5Val!=0.
617  || g3_prime6Val!=0.
618  || g3_prime7Val!=0.
619 
620  || g1ValIm!=0.
621  || g1_primeValIm!=0.
622  || g1_prime2ValIm!=0.
623  || g1_prime3ValIm!=0.
624  || g1_prime4ValIm!=0.
625  || g1_prime5ValIm!=0.
626  || g1_prime6ValIm!=0.
627  || g1_prime7ValIm!=0.
628  || g2ValIm!=0.
629  || g2_primeValIm!=0.
630  || g2_prime2ValIm!=0.
631  || g2_prime3ValIm!=0.
632  || g2_prime4ValIm!=0.
633  || g2_prime5ValIm!=0.
634  || g2_prime6ValIm!=0.
635  || g2_prime7ValIm!=0.
636  || g4ValIm!=0.
637  || g4_primeValIm!=0.
638  || g4_prime2ValIm!=0.
639  || g4_prime3ValIm!=0.
640  || g4_prime4ValIm!=0.
641  || g4_prime5ValIm!=0.
642  || g4_prime6ValIm!=0.
643  || g4_prime7ValIm!=0.
644  || g3ValIm!=0.
645  || g3_primeValIm!=0.
646  || g3_prime2ValIm!=0.
647  || g3_prime3ValIm!=0.
648  || g3_prime4ValIm!=0.
649  || g3_prime5ValIm!=0.
650  || g3_prime6ValIm!=0.
651  || g3_prime7ValIm!=0.
652  )
653  )
654  ||
655  (
656  ((VGammaVpmode1==0 && VGammaVpmode2==1) || (VGammaVpmode2==0 && VGammaVpmode1==1)) && (
657  gzgs1_prime2Val!=0.
658  || gzgs2Val!=0.
659  || gzgs4Val!=0.
660  || gzgs3Val!=0.
661 
662  || gzgs1_prime2ValIm!=0.
663  || gzgs2ValIm!=0.
664  || gzgs4ValIm!=0.
665  || gzgs3ValIm!=0.
666  )
667  )
668  ||
669  (
670  (VGammaVpmode1==VGammaVpmode2 && VGammaVpmode1==1) && (
671  ggsgs2Val!=0.
672  || ggsgs4Val!=0.
673  || ggsgs3Val!=0.
674 
675  || ggsgs2ValIm!=0.
676  || ggsgs4ValIm!=0.
677  || ggsgs3ValIm!=0.
678  )
679  )
680  ||
681  (
682  ((VGammaVpmode1==0 && VGammaVpmode2==2) || (VGammaVpmode2==0 && VGammaVpmode1==2)) && (
683  gvvp1Val!=0.
684 
685  || gvvp1ValIm!=0.
686  )
687  )
688  ||
689  (
690  (VGammaVpmode1==VGammaVpmode2 && VGammaVpmode1==2) && (
691  gvpvp1Val!=0.
692 
693  || gvpvp1ValIm!=0.
694  )
695  )
696  );
697 }
698 
699 void RooSpinZero::printParameters() const{
700  MELAout << "g1: (" << g1Val << ", " << g1ValIm << ")" << endl;
701  MELAout << "g2: (" << g2Val << ", " << g2ValIm << ")" << endl;
702  MELAout << "g3: (" << g3Val << ", " << g3ValIm << ")" << endl;
703  MELAout << "g4: (" << g4Val << ", " << g4ValIm << ")" << endl;
704  MELAout << "g1_prime: (" << g1_primeVal << ", " << g1_primeValIm << ")" << endl;
705  MELAout << "g2_prime: (" << g2_primeVal << ", " << g2_primeValIm << ")" << endl;
706  MELAout << "g3_prime: (" << g3_primeVal << ", " << g3_primeValIm << ")" << endl;
707  MELAout << "g4_prime: (" << g4_primeVal << ", " << g4_primeValIm << ")" << endl;
708  MELAout << "g1_prime2: (" << g1_prime2Val << ", " << g1_prime2ValIm << ")" << endl;
709  MELAout << "g2_prime2: (" << g2_prime2Val << ", " << g2_prime2ValIm << ")" << endl;
710  MELAout << "g3_prime2: (" << g3_prime2Val << ", " << g3_prime2ValIm << ")" << endl;
711  MELAout << "g4_prime2: (" << g4_prime2Val << ", " << g4_prime2ValIm << ")" << endl;
712  MELAout << "g1_prime3: (" << g1_prime3Val << ", " << g1_prime3ValIm << ")" << endl;
713  MELAout << "g2_prime3: (" << g2_prime3Val << ", " << g2_prime3ValIm << ")" << endl;
714  MELAout << "g3_prime3: (" << g3_prime3Val << ", " << g3_prime3ValIm << ")" << endl;
715  MELAout << "g4_prime3: (" << g4_prime3Val << ", " << g4_prime3ValIm << ")" << endl;
716  MELAout << "g1_prime4: (" << g1_prime4Val << ", " << g1_prime4ValIm << ")" << endl;
717  MELAout << "g2_prime4: (" << g2_prime4Val << ", " << g2_prime4ValIm << ")" << endl;
718  MELAout << "g3_prime4: (" << g3_prime4Val << ", " << g3_prime4ValIm << ")" << endl;
719  MELAout << "g4_prime4: (" << g4_prime4Val << ", " << g4_prime4ValIm << ")" << endl;
720  MELAout << "g1_prime5: (" << g1_prime5Val << ", " << g1_prime5ValIm << ")" << endl;
721  MELAout << "g2_prime5: (" << g2_prime5Val << ", " << g2_prime5ValIm << ")" << endl;
722  MELAout << "g3_prime5: (" << g3_prime5Val << ", " << g3_prime5ValIm << ")" << endl;
723  MELAout << "g4_prime5: (" << g4_prime5Val << ", " << g4_prime5ValIm << ")" << endl;
724  MELAout << "g1_prime6: (" << g1_prime6Val << ", " << g1_prime6ValIm << ")" << endl;
725  MELAout << "g2_prime6: (" << g2_prime6Val << ", " << g2_prime6ValIm << ")" << endl;
726  MELAout << "g3_prime6: (" << g3_prime6Val << ", " << g3_prime6ValIm << ")" << endl;
727  MELAout << "g4_prime6: (" << g4_prime6Val << ", " << g4_prime6ValIm << ")" << endl;
728  MELAout << "g1_prime7: (" << g1_prime7Val << ", " << g1_prime7ValIm << ")" << endl;
729  MELAout << "g2_prime7: (" << g2_prime7Val << ", " << g2_prime7ValIm << ")" << endl;
730  MELAout << "g3_prime7: (" << g3_prime7Val << ", " << g3_prime7ValIm << ")" << endl;
731  MELAout << "g4_prime7: (" << g4_prime7Val << ", " << g4_prime7ValIm << ")" << endl;
732  MELAout << "gzgs1_prime2: (" << gzgs1_prime2Val << ", " << gzgs1_prime2ValIm << ")" << endl;
733  MELAout << "gzgs2: (" << gzgs2Val << ", " << gzgs2ValIm << ")" << endl;
734  MELAout << "gzgs3: (" << gzgs3Val << ", " << gzgs3ValIm << ")" << endl;
735  MELAout << "gzgs4: (" << gzgs4Val << ", " << gzgs4ValIm << ")" << endl;
736  MELAout << "ggsgs2: (" << ggsgs2Val << ", " << ggsgs2ValIm << ")" << endl;
737  MELAout << "ggsgs3: (" << ggsgs3Val << ", " << ggsgs3ValIm << ")" << endl;
738  MELAout << "ggsgs4: (" << ggsgs4Val << ", " << ggsgs4ValIm << ")" << endl;
739 
740  MELAout << "Lambda: " << Lambda << endl;
741  MELAout << "Lambda_zgs1: " << Lambda_zgs1 << endl;
742  MELAout << "Lambda_z1: " << Lambda_z1 << endl;
743  MELAout << "Lambda_z2: " << Lambda_z2 << endl;
744  MELAout << "Lambda_z3: " << Lambda_z3 << endl;
745  MELAout << "Lambda_z4: " << Lambda_z4 << endl;
746  MELAout << "Lambda_Q: " << Lambda_Q << endl;
747  MELAout << "Lambda_z11: " << Lambda_z11 << endl;
748  MELAout << "Lambda_z21: " << Lambda_z21 << endl;
749  MELAout << "Lambda_z31: " << Lambda_z31 << endl;
750  MELAout << "Lambda_z41: " << Lambda_z41 << endl;
751  MELAout << "Lambda_z12: " << Lambda_z12 << endl;
752  MELAout << "Lambda_z22: " << Lambda_z22 << endl;
753  MELAout << "Lambda_z32: " << Lambda_z32 << endl;
754  MELAout << "Lambda_z42: " << Lambda_z42 << endl;
755  MELAout << "Lambda_z10: " << Lambda_z10 << endl;
756  MELAout << "Lambda_z20: " << Lambda_z20 << endl;
757  MELAout << "Lambda_z30: " << Lambda_z30 << endl;
758  MELAout << "Lambda_z40: " << Lambda_z40 << endl;
759  MELAout << "cz_q1sq: " << cz_q1sq << endl;
760  MELAout << "cz_q2sq: " << cz_q2sq << endl;
761  MELAout << "cz_q12sq: " << cz_q12sq << endl;
762 
763  MELAout << "gvvp1: (" << gvvp1Val << ", " << gvvp1ValIm << ")" << endl;
764  MELAout << "gvpvp1: (" << gvpvp1Val << ", " << gvpvp1ValIm << ")" << endl;
765 
766  RooSpin::printParameters();
767 }
RooSpinZero::g3_prime3ValIm
RooRealProxy g3_prime3ValIm
Definition: RooSpinZero.h:119
RooSpinZero::Lambda_z4
RooRealProxy Lambda_z4
Definition: RooSpinZero.h:152
RooSpin::calculateAmplitudeScale
virtual Double_t calculateAmplitudeScale(int VGammaVpmode1=0, int VGammaVpmode2=0) const
Definition: RooSpin.cc:366
Lambda_z11
double Lambda_z11
Definition: TMCFM.hh:163
RooSpinZero::gzgs2ValIm
RooRealProxy gzgs2ValIm
Definition: RooSpinZero.h:139
RooSpin::printParameters
virtual void printParameters() const
Definition: RooSpin.cc:397
RooSpinZero::g3_prime3Val
RooRealProxy g3_prime3Val
Definition: RooSpinZero.h:77
RooSpinZero::cz_q2sq
RooRealProxy cz_q2sq
Definition: RooSpinZero.h:168
RooSpinZero::g2_prime7ValIm
RooRealProxy g2_prime7ValIm
Definition: RooSpinZero.h:134
RooSpinZero::Lambda_z3
RooRealProxy Lambda_z3
Definition: RooSpinZero.h:151
RooSpinZero::ggsgs2Val
RooRealProxy ggsgs2Val
Definition: RooSpinZero.h:100
RooSpinZero::gzgs2Val
RooRealProxy gzgs2Val
Definition: RooSpinZero.h:97
RooSpinZero::g2Val
RooRealProxy g2Val
Definition: RooSpinZero.h:64
TVar::VerbosityLevel
VerbosityLevel
Definition: TVar.hh:47
Lambda_z41
double Lambda_z41
Definition: TMCFM.hh:163
RooSpinZero::g1_prime2ValIm
RooRealProxy g1_prime2ValIm
Definition: RooSpinZero.h:113
RooSpinZero::g3_prime5Val
RooRealProxy g3_prime5Val
Definition: RooSpinZero.h:85
RooSpinZero::g4_prime3Val
RooRealProxy g4_prime3Val
Definition: RooSpinZero.h:78
RooSpinZero::g4_prime5Val
RooRealProxy g4_prime5Val
Definition: RooSpinZero.h:86
RooSpinZero::g3_prime4Val
RooRealProxy g3_prime4Val
Definition: RooSpinZero.h:81
RooSpinZero::g1_prime4Val
RooRealProxy g1_prime4Val
Definition: RooSpinZero.h:79
RooSpinZero::Lambda_z12
RooRealProxy Lambda_z12
Definition: RooSpinZero.h:159
RooSpinZero::g1_prime7ValIm
RooRealProxy g1_prime7ValIm
Definition: RooSpinZero.h:133
RooSpinZero::g3ValIm
RooRealProxy g3ValIm
Definition: RooSpinZero.h:107
RooSpinZero::g3_prime7Val
RooRealProxy g3_prime7Val
Definition: RooSpinZero.h:93
Lambda_z4
double Lambda_z4
Definition: TMCFM.hh:162
RooSpinZero::g1ValIm
RooRealProxy g1ValIm
Definition: RooSpinZero.h:105
RooSpinZero::g4_prime4ValIm
RooRealProxy g4_prime4ValIm
Definition: RooSpinZero.h:124
RooSpinZero::g4_primeValIm
RooRealProxy g4_primeValIm
Definition: RooSpinZero.h:112
RooSpinZero::g4Val
RooRealProxy g4Val
Definition: RooSpinZero.h:66
RooSpinZero::g1_prime6Val
RooRealProxy g1_prime6Val
Definition: RooSpinZero.h:87
RooSpinZero::g2_prime5Val
RooRealProxy g2_prime5Val
Definition: RooSpinZero.h:84
RooSpinZero::gzgs3ValIm
RooRealProxy gzgs3ValIm
Definition: RooSpinZero.h:140
RooSpin::m12
RooRealProxy m12
Definition: RooSpin.h:113
RooSpinZero::Lambda_z31
RooRealProxy Lambda_z31
Definition: RooSpinZero.h:157
Lambda_z32
double Lambda_z32
Definition: TMCFM.hh:164
RooSpinZero::g2_prime3Val
RooRealProxy g2_prime3Val
Definition: RooSpinZero.h:76
RooSpinZero::g4_prime6ValIm
RooRealProxy g4_prime6ValIm
Definition: RooSpinZero.h:132
RooSpin::calculatePropagator
virtual void calculatePropagator(Double_t &propRe, Double_t &propIm, Double_t mass, Int_t propType=1) const
Definition: RooSpin.cc:133
RooSpinZero::g2_prime7Val
RooRealProxy g2_prime7Val
Definition: RooSpinZero.h:92
RooSpinZero::Lambda_zgs1
RooRealProxy Lambda_zgs1
Definition: RooSpinZero.h:148
RooSpinZero::Lambda_z30
RooRealProxy Lambda_z30
Definition: RooSpinZero.h:165
RooSpinZero::g3_prime6Val
RooRealProxy g3_prime6Val
Definition: RooSpinZero.h:89
RooSpinZero::Lambda_z21
RooRealProxy Lambda_z21
Definition: RooSpinZero.h:156
RooSpinZero::gzgs3Val
RooRealProxy gzgs3Val
Definition: RooSpinZero.h:98
Lambda_z10
double Lambda_z10
Definition: TMCFM.hh:165
RooSpinZero::g3_prime6ValIm
RooRealProxy g3_prime6ValIm
Definition: RooSpinZero.h:131
RooSpinZero::g1_primeValIm
RooRealProxy g1_primeValIm
Definition: RooSpinZero.h:109
RooSpinZero::ggsgs3ValIm
RooRealProxy ggsgs3ValIm
Definition: RooSpinZero.h:143
RooSpinZero::g1_prime3Val
RooRealProxy g1_prime3Val
Definition: RooSpinZero.h:75
RooSpinZero::g3_prime5ValIm
RooRealProxy g3_prime5ValIm
Definition: RooSpinZero.h:127
RooSpinZero::ggsgs3Val
RooRealProxy ggsgs3Val
Definition: RooSpinZero.h:101
RooSpinZero::g4_prime2ValIm
RooRealProxy g4_prime2ValIm
Definition: RooSpinZero.h:116
RooSpinZero::Lambda_Q
RooRealProxy Lambda_Q
Definition: RooSpinZero.h:153
RooSpinZero::RooSpinZero
RooSpinZero()
Definition: RooSpinZero.cc:8
RooSpin::m1
RooRealProxy m1
Definition: RooSpin.h:111
RooSpinZero::computeNeededAmplitude
virtual Bool_t computeNeededAmplitude(int VGammaVpmode1, int VGammaVpmode2) const final
Definition: RooSpinZero.cc:583
RooSpinZero::g4_prime2Val
RooRealProxy g4_prime2Val
Definition: RooSpinZero.h:74
RooSpinZero::g1_primeVal
RooRealProxy g1_primeVal
Definition: RooSpinZero.h:67
RooSpinZero::ggsgs2ValIm
RooRealProxy ggsgs2ValIm
Definition: RooSpinZero.h:142
RooSpinZero::Lambda
RooRealProxy Lambda
Definition: RooSpinZero.h:147
RooSpinZero::g1Val
RooRealProxy g1Val
Definition: RooSpinZero.h:63
RooSpinZero::g3_prime2ValIm
RooRealProxy g3_prime2ValIm
Definition: RooSpinZero.h:115
AnaMelaHelpers::multiplyComplexNumbers
void multiplyComplexNumbers(std::vector< Double_t > const &reals, std::vector< Double_t > const &imags, Double_t &resRe, Double_t &resIm)
Definition: RooSpin.cc:7
RooSpinZero::cz_q1sq
RooRealProxy cz_q1sq
Definition: RooSpinZero.h:167
RooSpinZero::gzgs4ValIm
RooRealProxy gzgs4ValIm
Definition: RooSpinZero.h:141
RooSpinZero::gvpvp1Val
RooRealProxy gvpvp1Val
Definition: RooSpinZero.h:173
RooSpinZero::Lambda_z40
RooRealProxy Lambda_z40
Definition: RooSpinZero.h:166
Lambda_Q
double Lambda_Q
Definition: TMCFM.hh:160
RooSpinZero::g3_primeVal
RooRealProxy g3_primeVal
Definition: RooSpinZero.h:69
MELAStreamHelpers::MELAout
MELAOutputStreamer MELAout
RooSpinZero::g3_prime2Val
RooRealProxy g3_prime2Val
Definition: RooSpinZero.h:73
RooSpinZero::gvpvp1ValIm
RooRealProxy gvpvp1ValIm
Definition: RooSpinZero.h:176
RooSpinZero::g1_prime2Val
RooRealProxy g1_prime2Val
Definition: RooSpinZero.h:71
Lambda_z1
double Lambda_z1
Definition: TMCFM.hh:162
RooSpinZero::g2_prime2Val
RooRealProxy g2_prime2Val
Definition: RooSpinZero.h:72
RooSpinZero::Lambda_z42
RooRealProxy Lambda_z42
Definition: RooSpinZero.h:162
modparameters::kappa
complex(8), public kappa
Definition: mod_Parameters.F90:882
RooSpinZero::g2_prime2ValIm
RooRealProxy g2_prime2ValIm
Definition: RooSpinZero.h:114
RooSpinZero::g4_prime7ValIm
RooRealProxy g4_prime7ValIm
Definition: RooSpinZero.h:136
RooSpinZero::Lambda_z2
RooRealProxy Lambda_z2
Definition: RooSpinZero.h:150
cz_q1sq
int cz_q1sq
Definition: TMCFM.hh:153
RooSpinZero::g4_prime3ValIm
RooRealProxy g4_prime3ValIm
Definition: RooSpinZero.h:120
RooSpin::verbosity
TVar::VerbosityLevel verbosity
Definition: RooSpin.h:139
Lambda_z42
double Lambda_z42
Definition: TMCFM.hh:164
Lambda_z40
double Lambda_z40
Definition: TMCFM.hh:165
RooSpinZero::gvvp1ValIm
RooRealProxy gvvp1ValIm
Definition: RooSpinZero.h:175
RooSpin::Vdecay2
RooSpin::VdecayType Vdecay2
Definition: RooSpin.h:136
RooSpinZero::g4_primeVal
RooRealProxy g4_primeVal
Definition: RooSpinZero.h:70
Lambda_z20
double Lambda_z20
Definition: TMCFM.hh:165
RooSpinZero
Definition: RooSpinZero.h:7
RooSpinZero::ggsgs4ValIm
RooRealProxy ggsgs4ValIm
Definition: RooSpinZero.h:144
RooSpin::kVdecayType_GammaOnshell
@ kVdecayType_GammaOnshell
Definition: RooSpin.h:30
RooSpinZero::g3_primeValIm
RooRealProxy g3_primeValIm
Definition: RooSpinZero.h:111
RooSpinZero::g1_prime3ValIm
RooRealProxy g1_prime3ValIm
Definition: RooSpinZero.h:117
RooSpinZero::g1_prime5ValIm
RooRealProxy g1_prime5ValIm
Definition: RooSpinZero.h:125
RooSpinZero::g4_prime5ValIm
RooRealProxy g4_prime5ValIm
Definition: RooSpinZero.h:128
RooSpinZero::calculateAi
virtual void calculateAi(Double_t &a1Re, Double_t &a1Im, Double_t &a2Re, Double_t &a2Im, Double_t &a3Re, Double_t &a3Im, int VGammaVpmode1=0, int VGammaVpmode2=0) const
Definition: RooSpinZero.cc:288
Lambda_z22
double Lambda_z22
Definition: TMCFM.hh:164
RooSpin
Definition: RooSpin.h:25
Lambda_z3
double Lambda_z3
Definition: TMCFM.hh:162
RooSpinZero::g4ValIm
RooRealProxy g4ValIm
Definition: RooSpinZero.h:108
RooSpinZero::Lambda_z32
RooRealProxy Lambda_z32
Definition: RooSpinZero.h:161
RooSpinZero::g2_primeValIm
RooRealProxy g2_primeValIm
Definition: RooSpinZero.h:110
RooSpinZero.h
RooSpinZero::Lambda_z41
RooRealProxy Lambda_z41
Definition: RooSpinZero.h:158
RooSpinZero::g4_prime6Val
RooRealProxy g4_prime6Val
Definition: RooSpinZero.h:90
RooSpinZero::Lambda_z1
RooRealProxy Lambda_z1
Definition: RooSpinZero.h:149
RooSpinZero::g1_prime6ValIm
RooRealProxy g1_prime6ValIm
Definition: RooSpinZero.h:129
RooSpinZero::modelCouplings
Definition: RooSpinZero.h:10
RooSpinZero::g2_prime4ValIm
RooRealProxy g2_prime4ValIm
Definition: RooSpinZero.h:122
RooSpin::m2
RooRealProxy m2
Definition: RooSpin.h:112
RooSpinZero::gzgs4Val
RooRealProxy gzgs4Val
Definition: RooSpinZero.h:99
RooSpinZero::calculateAmplitudes
virtual void calculateAmplitudes(Double_t &A00Re, Double_t &A00Im, Double_t &AppRe, Double_t &AppIm, Double_t &AmmRe, Double_t &AmmIm, int VGammaVpmode1=0, int VGammaVpmode2=0) const
Definition: RooSpinZero.cc:477
RooSpinZero::g4_prime7Val
RooRealProxy g4_prime7Val
Definition: RooSpinZero.h:94
RooSpinZero::g3Val
RooRealProxy g3Val
Definition: RooSpinZero.h:65
RooSpinZero::ggsgs4Val
RooRealProxy ggsgs4Val
Definition: RooSpinZero.h:102
RooSpinZero::g2ValIm
RooRealProxy g2ValIm
Definition: RooSpinZero.h:106
RooSpinZero::gvvp1Val
RooRealProxy gvvp1Val
Definition: RooSpinZero.h:172
RooSpinZero::g2_prime6Val
RooRealProxy g2_prime6Val
Definition: RooSpinZero.h:88
RooSpinZero::printParameters
virtual void printParameters() const
Definition: RooSpinZero.cc:699
RooSpinZero::g2_prime3ValIm
RooRealProxy g2_prime3ValIm
Definition: RooSpinZero.h:118
Lambda_z30
double Lambda_z30
Definition: TMCFM.hh:165
RooSpinZero::gzgs1_prime2Val
RooRealProxy gzgs1_prime2Val
Definition: RooSpinZero.h:96
Lambda_z31
double Lambda_z31
Definition: TMCFM.hh:163
RooSpinZero::g1_prime5Val
RooRealProxy g1_prime5Val
Definition: RooSpinZero.h:83
Lambda_zgs1
double Lambda_zgs1
Definition: TMCFM.hh:161
cz_q2sq
int cz_q2sq
Definition: TMCFM.hh:153
RooSpinZero::g1_prime7Val
RooRealProxy g1_prime7Val
Definition: RooSpinZero.h:91
MELAStreamHelpers
Definition: MELAStreamHelpers.hh:7
RooSpinZero::g1_prime4ValIm
RooRealProxy g1_prime4ValIm
Definition: RooSpinZero.h:121
RooSpinZero::g2_prime6ValIm
RooRealProxy g2_prime6ValIm
Definition: RooSpinZero.h:130
RooSpinZero::g2_primeVal
RooRealProxy g2_primeVal
Definition: RooSpinZero.h:68
RooSpin::GeVunit
static constexpr Double_t GeVunit
Definition: RooSpin.h:141
RooSpin::getMVGamV
virtual void getMVGamV(Double_t *mV=0, Double_t *gamV=0) const
Definition: RooSpin.cc:305
Lambda_z2
double Lambda_z2
Definition: TMCFM.hh:162
RooSpinZero::cz_q12sq
RooRealProxy cz_q12sq
Definition: RooSpinZero.h:169
RooSpinZero::gzgs1_prime2ValIm
RooRealProxy gzgs1_prime2ValIm
Definition: RooSpinZero.h:138
TVar::DEBUG
@ DEBUG
Definition: TVar.hh:51
RooSpinZero::Lambda_z22
RooRealProxy Lambda_z22
Definition: RooSpinZero.h:160
RooSpin::Vdecay1
RooSpin::VdecayType Vdecay1
Definition: RooSpin.h:135
RooSpinZero::Lambda_z10
RooRealProxy Lambda_z10
Definition: RooSpinZero.h:163
Lambda_z12
double Lambda_z12
Definition: TMCFM.hh:164
RooSpin::modelMeasurables
Definition: RooSpin.h:50
RooSpinZero::g2_prime4Val
RooRealProxy g2_prime4Val
Definition: RooSpinZero.h:80
RooSpinZero::Lambda_z20
RooRealProxy Lambda_z20
Definition: RooSpinZero.h:164
RooSpin::modelParameters
Definition: RooSpin.h:61
RooSpinZero::g2_prime5ValIm
RooRealProxy g2_prime5ValIm
Definition: RooSpinZero.h:126
cz_q12sq
int cz_q12sq
Definition: TMCFM.hh:153
RooSpinZero::Lambda_z11
RooRealProxy Lambda_z11
Definition: RooSpinZero.h:155
RooSpin::VdecayType
VdecayType
Definition: RooSpin.h:28
Lambda_z21
double Lambda_z21
Definition: TMCFM.hh:163
RooSpinZero::g3_prime4ValIm
RooRealProxy g3_prime4ValIm
Definition: RooSpinZero.h:123
RooSpinZero::g3_prime7ValIm
RooRealProxy g3_prime7ValIm
Definition: RooSpinZero.h:135
RooSpinZero::g4_prime4Val
RooRealProxy g4_prime4Val
Definition: RooSpinZero.h:82
Generated by   doxygen 1.8.17