BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDsToPipPipPimPi0.cc
Go to the documentation of this file.
1//--------------------------------------------------------------------------
2// Environment:
3// This software is part of the EvtGen package developed jointly
4// for the BaBar and CLEO collaborations. If you use all or part
5// of it, please give an appropriate acknowledgement.
6//
7// Copyright Information: See EvtGen/COPYRIGHT
8// Copyright (C) 1998 Caltech, UCSB
9//
10// Module: EvtDsToPipPipPimPi0.cc
11//
12// Description: Routine to handle three-body decays of Ds+/Ds-
13//
14// Modification history:
15// Liaoyuan Dong Jan 08 2023 Module created
16//------------------------------------------------------------------------
17#include "../EvtGenBase/EvtPatches.hh"
18
19#include "../EvtGenBase/EvtGenKine.hh"
20#include "../EvtGenBase/EvtPDL.hh"
21#include "../EvtGenBase/EvtParticle.hh"
22#include "../EvtGenBase/EvtReport.hh"
23#include "../EvtGenBase/EvtVector4R.hh"
24#include "EvtDsToPipPipPimPi0.hh"
25#include <math.h>
26#include <stdlib.h>
27#include <string>
28
29EvtDsToPipPipPimPi0::~EvtDsToPipPipPimPi0() {}
30
31void EvtDsToPipPipPimPi0::getName( std::string& model_name ) {
32
33 model_name = "DsToPipPipPimPi0";
34}
35
36EvtDecayBase* EvtDsToPipPipPimPi0::clone() { return new EvtDsToPipPipPimPi0; }
37
38void EvtDsToPipPipPimPi0::init() {
39
40 // std::cout << "EvtDsToPipPipPimPi0 (Oct 17, 2022) ==> Initialization !" << std::endl;
41
42 // check that there are 0 arguments
43 checkNArg( 0 );
44 // check that there are 4 daughters
45 checkNDaug( 4 );
46 // check that parent and daughters are correct
47 if ( getParentId() == EvtPDL::getId( "D_s+" ) )
48 {
49 if ( getDaug( 0 ) != EvtPDL::getId( "pi+" ) )
50 {
51 report( ERROR, "EvtGen" ) << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!" << endl;
52 ::abort();
53 }
54 if ( getDaug( 1 ) != EvtPDL::getId( "pi+" ) )
55 {
56 report( ERROR, "EvtGen" ) << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!" << endl;
57 ::abort();
58 }
59 if ( getDaug( 2 ) != EvtPDL::getId( "pi-" ) )
60 {
61 report( ERROR, "EvtGen" ) << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!" << endl;
62 ::abort();
63 }
64 if ( getDaug( 3 ) != EvtPDL::getId( "pi0" ) )
65 {
66 report( ERROR, "EvtGen" ) << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!" << endl;
67 ::abort();
68 }
69 }
70 else if ( getParentId() == EvtPDL::getId( "D_s-" ) )
71 {
72 if ( getDaug( 0 ) != EvtPDL::getId( "pi-" ) )
73 {
74 report( ERROR, "EvtGen" ) << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!" << endl;
75 ::abort();
76 }
77 if ( getDaug( 1 ) != EvtPDL::getId( "pi-" ) )
78 {
79 report( ERROR, "EvtGen" ) << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!" << endl;
80 ::abort();
81 }
82 if ( getDaug( 2 ) != EvtPDL::getId( "pi+" ) )
83 {
84 report( ERROR, "EvtGen" ) << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!" << endl;
85 ::abort();
86 }
87 if ( getDaug( 3 ) != EvtPDL::getId( "pi0" ) )
88 {
89 report( ERROR, "EvtGen" ) << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!" << endl;
90 ::abort();
91 }
92 }
93 else
94 {
95 report( ERROR, "EvtGen" ) << "In EvtDsToPipPipPimPi0, Mother is Wrong!" << endl;
96 ::abort();
97 }
98}
99
100void EvtDsToPipPipPimPi0::initProbMax() { setProbMax( 10000.0 ); }
101double EvtDsToPipPipPimPi0::mass1[11] = { 0.96500, 0.77526, 1.01946, 0.78200, 1.23000, 1.23000,
102 1.35000, 1.27550, 1.46500, 0.77526, 1.37100 };
103double EvtDsToPipPipPimPi0::mass2[11] = { 0.77526, 0.77526, 0.77526, 0.77526, 0.77526, 0.77526,
104 0.77526, 0.77526, 0.77526, 1.46500, 0.77526 };
105double EvtDsToPipPipPimPi0::width1[11] = { 0.05500, 0.14910, 0.00425, 0.00800,
106 0.42000, 0.42000, 0.26500, 0.18670,
107 0.40000, 0.14910, 0.17000 };
108double EvtDsToPipPipPimPi0::width2[11] = { 0.14910, 0.14910, 0.14910, 0.14910,
109 0.14910, 0.14910, 0.14910, 0.14910,
110 0.14910, 0.40000, 0.14910 };
111double EvtDsToPipPipPimPi0::amp[11] = { 4.0026e-1, 1.2202e-1, 7.4682e-2, 1.9591e-1,
112 5.7255e-1, 3.1005e-1, 1.0, 7.6861e-1,
113 1.1482e+0, 5.7077e-1, 2.7970e-1 };
114double EvtDsToPipPipPimPi0::phase[11] = { -2.2962e+0, 1.0977e+0, 2.8957e+0, 3.2197e+0,
115 3.7810e+0, -1.4586e+0, 0.0, 1.1081e+0,
116 4.3439e-1, -1.6993e+0, 2.2240e+0 };
117double EvtDsToPipPipPimPi0::mD = 1.86486;
118double EvtDsToPipPipPimPi0::metap = 0.95778;
119double EvtDsToPipPipPimPi0::mKst0 = 0.89581;
120double EvtDsToPipPipPimPi0::GKst0 = 0.0474;
121double EvtDsToPipPipPimPi0::mKstp = 0.89166;
122double EvtDsToPipPipPimPi0::Gkstp = 0.0508;
123double EvtDsToPipPipPimPi0::mk0 = 0.497614;
124double EvtDsToPipPipPimPi0::mass_Kaon = 0.49368;
125double EvtDsToPipPipPimPi0::mass_Ks = 0.497611;
126double EvtDsToPipPipPimPi0::mass_Eta = 0.547862;
127double EvtDsToPipPipPimPi0::mass_Pion = 0.13957;
128double EvtDsToPipPipPimPi0::math_pi = 3.1415926;
129double EvtDsToPipPipPimPi0::mass_Pion2 = 0.0194797849;
130double EvtDsToPipPipPimPi0::mass_2Pion = 0.27914;
131double EvtDsToPipPipPimPi0::math_2pi = 6.2831852;
132double EvtDsToPipPipPimPi0::rD2 = 25.0; // 5*5
133double EvtDsToPipPipPimPi0::rRes1 = 9.0; // 3*3
134double EvtDsToPipPipPimPi0::rRes2 = 9.0; // 3*3
135double EvtDsToPipPipPimPi0::g1 = 0.5468;
136double EvtDsToPipPipPimPi0::g2 = 0.23; // K*0(1430)
137double EvtDsToPipPipPimPi0::GS1 = 0.636619783;
138double EvtDsToPipPipPimPi0::GS2 = 0.01860182466;
139double EvtDsToPipPipPimPi0::GS3 = 0.1591549458; // 1/(2*math_2pi)
140double EvtDsToPipPipPimPi0::GS4 = 0.00620060822;
141int EvtDsToPipPipPimPi0::G[4][4] = {
142 { 1, 0, 0, 0 }, { 0, -1, 0, 0 }, { 0, 0, -1, 0 }, { 0, 0, 0, -1 } };
143int EvtDsToPipPipPimPi0::E[4][4][4][4] = {
144 { { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
145 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 1 }, { 0, 0, -1, 0 } },
146 { { 0, 0, 0, 0 }, { 0, 0, 0, -1 }, { 0, 0, 0, 0 }, { 0, 1, 0, 0 } },
147 { { 0, 0, 0, 0 }, { 0, 0, 1, 0 }, { 0, -1, 0, 0 }, { 0, 0, 0, 0 } } },
148 { { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, -1 }, { 0, 0, 1, 0 } },
149 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
150 { { 0, 0, 0, 1 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { -1, 0, 0, 0 } },
151 { { 0, 0, -1, 0 }, { 0, 0, 0, 0 }, { 1, 0, 0, 0 }, { 0, 0, 0, 0 } } },
152 { { { 0, 0, 0, 0 }, { 0, 0, 0, 1 }, { 0, 0, 0, 0 }, { 0, -1, 0, 0 } },
153 { { 0, 0, 0, -1 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 1, 0, 0, 0 } },
154 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
155 { { 0, 1, 0, 0 }, { -1, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } } },
156 { { { 0, 0, 0, 0 }, { 0, 0, -1, 0 }, { 0, 1, 0, 0 }, { 0, 0, 0, 0 } },
157 { { 0, 0, 1, 0 }, { 0, 0, 0, 0 }, { -1, 0, 0, 0 }, { 0, 0, 0, 0 } },
158 { { 0, -1, 0, 0 }, { 1, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
159 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } } } };
160
161void EvtDsToPipPipPimPi0::decay( EvtParticle* p ) {
162
163 p->initializePhaseSpace( getNDaug(), getDaugs() );
164
165 EvtVector4R PPip1 = p->getDaug( 0 )->getP4();
166 EvtVector4R PPip2 = p->getDaug( 1 )->getP4();
167 EvtVector4R PPim = p->getDaug( 2 )->getP4();
168 EvtVector4R PPi0 = p->getDaug( 3 )->getP4();
169
170 double Pip1[4], Pip2[4], Pim[4], Pi0[4];
171 double cof[2], amp_tmp[2], amp_PDF[2], PDF[2];
172 if ( getParentId() == EvtPDL::getId( "D_s+" ) )
173 {
174 Pip1[0] = PPip1.get( 0 );
175 Pip1[1] = PPip1.get( 1 );
176 Pip1[2] = PPip1.get( 2 );
177 Pip1[3] = PPip1.get( 3 );
178 Pip2[0] = PPip2.get( 0 );
179 Pip2[1] = PPip2.get( 1 );
180 Pip2[2] = PPip2.get( 2 );
181 Pip2[3] = PPip2.get( 3 );
182 Pim[0] = PPim.get( 0 );
183 Pim[1] = PPim.get( 1 );
184 Pim[2] = PPim.get( 2 );
185 Pim[3] = PPim.get( 3 );
186 Pi0[0] = PPi0.get( 0 );
187 Pi0[1] = PPi0.get( 1 );
188 Pi0[2] = PPi0.get( 2 );
189 Pi0[3] = PPi0.get( 3 );
190 }
191 else
192 {
193 Pip1[0] = PPip1.get( 0 );
194 Pip1[1] = -( PPip1.get( 1 ) );
195 Pip1[2] = -( PPip1.get( 2 ) );
196 Pip1[3] = -( PPip1.get( 3 ) );
197 Pip2[0] = PPip2.get( 0 );
198 Pip2[1] = -( PPip2.get( 1 ) );
199 Pip2[2] = -( PPip2.get( 2 ) );
200 Pip2[3] = -( PPip2.get( 3 ) );
201 Pim[0] = PPim.get( 0 );
202 Pim[1] = -( PPim.get( 1 ) );
203 Pim[2] = -( PPim.get( 2 ) );
204 Pim[3] = -( PPim.get( 3 ) );
205 Pi0[0] = PPi0.get( 0 );
206 Pi0[1] = -( PPi0.get( 1 ) );
207 Pi0[2] = -( PPi0.get( 2 ) );
208 Pi0[3] = -( PPi0.get( 3 ) );
209 }
210
211 double p12[4], p13[4], p14[4], p23[4], p24[4], p34[4];
212 double p123[4], p124[4], p134[4], p234[4], pD[4];
213 double Sf1, Sf2, Sf3, Sf4, Sf5, Sf6, Sf7, Sf8;
214 Sf1 = spinf1( Pip1[0], Pip1[1], Pip1[2], Pip1[3], Pip2[0], Pip2[1], Pip2[2], Pip2[3], Pim[0],
215 Pim[1], Pim[2], Pim[3], Pi0[0], Pi0[1], Pi0[2], Pi0[3] );
216 Sf2 = spinf2( Pip1[0], Pip1[1], Pip1[2], Pip1[3], Pip2[0], Pip2[1], Pip2[2], Pip2[3], Pim[0],
217 Pim[1], Pim[2], Pim[3], Pi0[0], Pi0[1], Pi0[2], Pi0[3] );
218 Sf3 = spinf3( Pip1[0], Pip1[1], Pip1[2], Pip1[3], Pip2[0], Pip2[1], Pip2[2], Pip2[3], Pim[0],
219 Pim[1], Pim[2], Pim[3], Pi0[0], Pi0[1], Pi0[2], Pi0[3] );
220 Sf4 = spinf4( Pip1[0], Pip1[1], Pip1[2], Pip1[3], Pip2[0], Pip2[1], Pip2[2], Pip2[3], Pim[0],
221 Pim[1], Pim[2], Pim[3], Pi0[0], Pi0[1], Pi0[2], Pi0[3] );
222 Sf5 = spinf5( Pip1[0], Pip1[1], Pip1[2], Pip1[3], Pip2[0], Pip2[1], Pip2[2], Pip2[3], Pim[0],
223 Pim[1], Pim[2], Pim[3], Pi0[0], Pi0[1], Pi0[2], Pi0[3] );
224 Sf6 = spinf6( Pip1[0], Pip1[1], Pip1[2], Pip1[3], Pip2[0], Pip2[1], Pip2[2], Pip2[3], Pim[0],
225 Pim[1], Pim[2], Pim[3], Pi0[0], Pi0[1], Pi0[2], Pi0[3] );
226 Sf7 = spinf7( Pip1[0], Pip1[1], Pip1[2], Pip1[3], Pip2[0], Pip2[1], Pip2[2], Pip2[3], Pim[0],
227 Pim[1], Pim[2], Pim[3], Pi0[0], Pi0[1], Pi0[2], Pi0[3] );
228 Sf8 = spinf8( Pip1[0], Pip1[1], Pip1[2], Pip1[3], Pip2[0], Pip2[1], Pip2[2], Pip2[3], Pim[0],
229 Pim[1], Pim[2], Pim[3], Pi0[0], Pi0[1], Pi0[2], Pi0[3] );
230
231 for ( int i = 0; i != 4; i++ )
232 {
233 p12[i] = Pip1[i] + Pip2[i];
234 p13[i] = Pip1[i] + Pim[i];
235 p14[i] = Pip1[i] + Pi0[i];
236 p23[i] = Pip2[i] + Pim[i];
237 p24[i] = Pip2[i] + Pi0[i];
238 p34[i] = Pim[i] + Pi0[i];
239 p123[i] = Pip1[i] + Pip2[i] + Pim[i];
240 p124[i] = Pip1[i] + Pip2[i] + Pi0[i];
241 p134[i] = Pip1[i] + Pim[i] + Pi0[i];
242 p234[i] = Pip2[i] + Pim[i] + Pi0[i];
243 pD[i] = Pip1[i] + Pip2[i] + Pim[i] + Pi0[i];
244 }
245 double sPip1, sPip2, sPim, sPi0, s12, s13, s14, s23, s24, s34, s123, s124, s134, s234, sD;
246 sPi0 = SCADot( Pi0, Pi0 );
247 sPip1 = SCADot( Pip1, Pip1 );
248 sPip2 = SCADot( Pip2, Pip2 );
249 sPim = SCADot( Pim, Pim );
250
251 s12 = SCADot( p12, p12 );
252 s13 = SCADot( p13, p13 );
253 s14 = SCADot( p14, p14 );
254 s23 = SCADot( p23, p23 );
255 s24 = SCADot( p24, p24 );
256 s34 = SCADot( p34, p34 );
257
258 s123 = SCADot( p123, p123 );
259 s124 = SCADot( p124, p124 );
260 s134 = SCADot( p134, p134 );
261 s234 = SCADot( p234, p234 );
262
263 sD = SCADot( pD, pD );
264
265 double spion12[2] = { mass_Kaon * mass_Kaon, sPip1 };
266 double spion22[2] = { mass_Kaon * mass_Kaon, sPip2 };
267 double spim2[2] = { mass_Kaon * mass_Kaon, sPim };
268
269 double t1_24[4], t1_14[4], t1_13[4], t1_23[4], t1_12[4], t1_34[4];
270 double t1_14_23[4], t1_13_24[4], t1_12_34[4];
271 double t1_134_2[4], t1_234_1[4], t1_123_4[4];
272 double t1_13_4[4], t1_23_4[4], t1_14_3[4], t1_24_3[4], t1_34_1[4], t1_34_2[4];
273
274 double t2_13[4][4], t2_23[4][4];
275 double t2_13_24[4][4], t2_14_23[4][4];
276 double t2_123_4[4][4];
277 double t2_13_2[4][4], t2_23_1[4][4];
278
279 calt1( Pi0, Pip2, t1_24 );
280 calt1( Pi0, Pip1, t1_14 );
281
282 calt1( Pip1, Pim, t1_13 );
283 calt1( Pip2, Pim, t1_23 );
284
285 calt1( Pip1, Pip2, t1_12 );
286 calt1( Pim, Pi0, t1_34 );
287
288 calt1( p14, p23, t1_14_23 );
289 calt1( p24, p13, t1_13_24 );
290 calt1( p12, p34, t1_12_34 );
291
292 calt1( p134, Pip2, t1_134_2 );
293 calt1( p234, Pip1, t1_234_1 );
294 calt1( p123, Pi0, t1_123_4 );
295
296 calt1( p13, Pi0, t1_13_4 );
297 calt1( p23, Pi0, t1_23_4 );
298 calt1( p14, Pim, t1_14_3 );
299 calt1( p24, Pim, t1_24_3 );
300 calt1( p34, Pip1, t1_34_1 );
301 calt1( p34, Pip2, t1_34_2 );
302
303 calt2( Pip1, Pim, t2_13 );
304 calt2( Pip2, Pim, t2_23 );
305
306 calt2( p13, p24, t2_13_24 );
307 calt2( p14, p23, t2_14_23 );
308
309 calt2( p123, Pi0, t2_123_4 );
310
311 calt2( p13, Pip2, t2_13_2 );
312 calt2( p23, Pip1, t2_23_1 );
313
314 amp_PDF[0] = 0;
315 amp_PDF[1] = 0;
316 PDF[0] = 0;
317 PDF[1] = 0;
318
319 double temp_PDF, tmp1, tmp2, tmp3, tmp4, temp_PDF1, amp_tmp1[2], amp_tmp2[2];
320 double pro[2], pro0[2], pro1[2], pro2[2], pro3[2], pro4[2];
321 double t1D[4], B[3], Bc[3];
322 double mass1sq, mass2sq;
323
324 for ( int i = 0; i < 11; i++ )
325 {
326 amp_tmp[0] = 0;
327 amp_tmp[1] = 0;
328 tmp1 = 0;
329 tmp2 = 0;
330 temp_PDF = 0;
331 ;
332 cof[0] = amp[i] * cos( phase[i] );
333 cof[1] = amp[i] * sin( phase[i] );
334 mass1sq = mass1[i] * mass1[i];
335 mass2sq = mass2[i] * mass2[i];
336 temp_PDF = 0;
337 temp_PDF1 = 0;
338 amp_tmp[0] = 0;
339 amp_tmp[1] = 0;
340 amp_tmp1[0] = 0;
341 amp_tmp1[1] = 0;
342 amp_tmp2[0] = 0;
343 amp_tmp2[1] = 0;
344 tmp1 = 0;
345 tmp2 = 0;
346 tmp3 = 0;
347 tmp4 = 0;
348
349 if ( i == 0 )
350 {
351 temp_PDF = 0;
352 propagator980( mass1[i], s23, spion12, spion22, pro0 );
353 propagatorGS( mass2sq, mass2[i], width2[i], s14, sPip1, sPi0, rRes2, pro1 );
354 Com_Multi( pro0, pro1, pro );
355
356 for ( int a = 0; a < 4; a++ ) { temp_PDF += G[a][a] * t1_14_23[a] * t1_14[a]; }
357
358 B[1] = barrier( 1, s14, sPip1, sPi0, rRes2, mass2[i] );
359 B[2] = barrier( 1, sD, s23, s14, rD2, 1.9683 );
360 tmp1 = B[1] * B[2] * temp_PDF;
361 amp_tmp1[0] = tmp1 * pro[0];
362 amp_tmp1[1] = tmp1 * pro[1];
363
364 temp_PDF = 0;
365 propagator980( mass1[i], s13, spion12, spion22, pro0 );
366 propagatorGS( mass2sq, mass2[i], width2[i], s24, sPip2, sPi0, rRes2, pro1 );
367 Com_Multi( pro0, pro1, pro );
368
369 for ( int a = 0; a < 4; a++ ) { temp_PDF += G[a][a] * t1_13_24[a] * t1_24[a]; }
370
371 B[1] = barrier( 1, s24, sPip2, sPi0, rRes2, mass2[i] );
372 B[2] = barrier( 1, sD, s13, s24, rD2, 1.9683 );
373 tmp2 = B[1] * B[2] * temp_PDF;
374 amp_tmp2[0] = tmp2 * pro[0];
375 amp_tmp2[1] = tmp2 * pro[1];
376 }
377
378 if ( i == 1 )
379 {
380 temp_PDF = 0;
381 propagatorGS( mass2sq, mass2[i], width2[i], s14, sPip1, sPi0, rRes2, pro0 );
382 propagatorGS( mass2sq, mass2[i], width2[i], s23, sPip2, sPim, rRes2, pro1 );
383
384 for ( int a = 0; a < 4; a++ ) { temp_PDF += G[a][a] * t1_14[a] * t1_23[a]; }
385
386 B[1] = barrier( 1, s14, sPip1, sPi0, rRes2, mass2[i] );
387 B[2] = barrier( 1, s23, sPip2, sPim, rRes2, mass2[i] );
388 tmp1 = B[1] * B[2] * temp_PDF;
389 amp_tmp1[0] = tmp1 * pro0[0] * pro1[0];
390 amp_tmp1[1] = tmp1 * pro0[1] * pro1[1];
391
392 temp_PDF = 0;
393 propagatorGS( mass2sq, mass2[i], width2[i], s24, sPip2, sPi0, rRes2, pro0 );
394 propagatorGS( mass2sq, mass2[i], width2[i], s13, sPip1, sPim, rRes2, pro1 );
395
396 for ( int a = 0; a < 4; a++ ) { temp_PDF += G[a][a] * t1_24[a] * t1_13[a]; }
397
398 B[1] = barrier( 1, s24, sPip2, sPi0, rRes2, mass2[i] );
399 B[2] = barrier( 1, s13, sPip1, sPim, rRes2, mass2[i] );
400 tmp2 = B[1] * B[2] * temp_PDF;
401 amp_tmp2[0] = tmp2 * pro0[0] * pro1[0];
402 amp_tmp2[1] = tmp2 * pro0[1] * pro1[1];
403 }
404
405 if ( i == 2 )
406 {
407 amp_tmp1[0] = 0;
408 amp_tmp1[1] = 0;
409 amp_tmp2[0] = 0;
410 amp_tmp2[1] = 0;
411
412 temp_PDF = 0;
413 propagatorRBW( mass1sq, mass1[i], width1[i], s134, s14, sPim, rRes1, 1, pro1 );
414 propagatorGS( mass2sq, mass2[i], width2[i], s14, sPip1, sPi0, rRes2, pro0 );
415
416 temp_PDF = Sf3;
417
418 B[0] = barrier( 1, sD, s134, sPip2, rD2, 1.9683 );
419 B[1] = barrier( 1, s134, s14, sPim, rRes1, mass1[i] );
420 B[2] = barrier( 1, s14, sPip1, sPi0, rRes2, mass2[i] );
421 tmp1 = B[0] * B[1] * B[2] * temp_PDF;
422 amp_tmp1[0] += tmp1 * pro0[0] * pro1[0];
423 amp_tmp1[1] += tmp1 * pro0[1] * pro1[1];
424
425 temp_PDF = 0;
426 propagatorRBW( mass1sq, mass1[i], width1[i], s234, s24, sPim, rRes1, 1, pro1 );
427 propagatorGS( mass2sq, mass2[i], width2[i], s24, sPip2, sPi0, rRes2, pro0 );
428
429 temp_PDF = Sf4;
430
431 B[0] = barrier( 1, sD, s234, sPip1, rD2, 1.9683 );
432 B[1] = barrier( 1, s234, s24, sPim, rRes1, mass1[i] );
433 B[2] = barrier( 1, s24, sPip2, sPi0, rRes2, mass2[i] );
434 tmp2 = B[0] * B[1] * B[2] * temp_PDF;
435 amp_tmp2[0] += tmp2 * pro0[0] * pro1[0];
436 amp_tmp2[1] += tmp2 * pro0[1] * pro1[1];
437
438 temp_PDF = 0;
439 propagatorRBW( mass1sq, mass1[i], width1[i], s134, s14, sPim, rRes1, 1, pro1 );
440 propagatorGS( mass2sq, mass2[i], width2[i], s13, sPip1, sPim, rRes2, pro0 );
441
442 temp_PDF = Sf5;
443
444 B[0] = barrier( 1, sD, s134, sPip2, rD2, 1.9683 );
445 B[1] = barrier( 1, s134, s13, sPi0, rRes1, mass1[i] );
446 B[2] = barrier( 1, s13, sPip1, sPim, rRes2, mass2[i] );
447 tmp1 = B[0] * B[1] * B[2] * temp_PDF;
448 amp_tmp1[0] -= tmp1 * pro0[0] * pro1[0];
449 amp_tmp1[1] -= tmp1 * pro0[1] * pro1[1];
450
451 temp_PDF = 0;
452 propagatorRBW( mass1sq, mass1[i], width1[i], s234, s24, sPim, rRes1, 1, pro1 );
453 propagatorGS( mass2sq, mass2[i], width2[i], s23, sPip2, sPim, rRes2, pro0 );
454
455 temp_PDF = Sf6;
456
457 B[0] = barrier( 1, sD, s234, sPip1, rD2, 1.9683 );
458 B[1] = barrier( 1, s234, s23, sPi0, rRes1, mass1[i] );
459 B[2] = barrier( 1, s23, sPip2, sPim, rRes2, mass2[i] );
460 tmp2 = B[0] * B[1] * B[2] * temp_PDF;
461 amp_tmp2[0] -= tmp2 * pro0[0] * pro1[0];
462 amp_tmp2[1] -= tmp2 * pro0[1] * pro1[1];
463
464 temp_PDF = 0;
465 propagatorRBW( mass1sq, mass1[i], width1[i], s134, s34, sPip1, rRes1, 1, pro1 );
466 propagatorGS( mass2sq, mass2[i], width2[i], s34, sPim, sPi0, rRes2, pro0 );
467
468 temp_PDF = Sf7;
469
470 B[0] = barrier( 1, sD, s134, sPip2, rD2, 1.9683 );
471 B[1] = barrier( 1, s134, s34, sPip1, rRes1, mass1[i] );
472 B[2] = barrier( 1, s34, sPim, sPi0, rRes2, mass2[i] );
473 tmp1 = B[0] * B[1] * B[2] * temp_PDF;
474 amp_tmp1[0] -= tmp1 * pro0[0] * pro1[0];
475 amp_tmp1[1] -= tmp1 * pro0[1] * pro1[1];
476
477 temp_PDF = 0;
478 propagatorRBW( mass1sq, mass1[i], width1[i], s234, s34, sPip2, rRes1, 1, pro1 );
479 propagatorGS( mass2sq, mass2[i], width2[i], s34, sPim, sPi0, rRes2, pro0 );
480
481 temp_PDF = Sf8;
482
483 B[0] = barrier( 1, sD, s234, sPip1, rD2, 1.9683 );
484 B[1] = barrier( 1, s234, s34, sPip2, rRes1, mass1[i] );
485 B[2] = barrier( 1, s34, sPim, sPi0, rRes2, mass2[i] );
486 tmp2 = B[0] * B[1] * B[2] * temp_PDF;
487 amp_tmp2[0] -= tmp2 * pro0[0] * pro1[0];
488 amp_tmp2[1] -= tmp2 * pro0[1] * pro1[1];
489 }
490
491 if ( i == 3 )
492 {
493 amp_tmp1[0] = 0;
494 amp_tmp1[1] = 0;
495 amp_tmp2[0] = 0;
496 amp_tmp2[1] = 0;
497
498 temp_PDF = 0;
499 propagatorRBW( mass1sq, mass1[i], width1[i], s134, s14, sPim, rRes1, 1, pro1 );
500 propagatorGS( mass2sq, mass2[i], width2[i], s14, sPip1, sPi0, rRes2, pro0 );
501
502 temp_PDF = Sf3;
503
504 B[0] = barrier( 1, sD, s134, sPip2, rD2, 1.9683 );
505 B[1] = barrier( 1, s134, s14, sPim, rRes1, mass1[i] );
506 B[2] = barrier( 1, s14, sPip1, sPi0, rRes2, mass2[i] );
507 tmp1 = B[0] * B[1] * B[2] * temp_PDF;
508 amp_tmp1[0] += tmp1 * pro0[0] * pro1[0];
509 amp_tmp1[1] += tmp1 * pro0[1] * pro1[1];
510
511 temp_PDF = 0;
512 propagatorRBW( mass1sq, mass1[i], width1[i], s234, s24, sPim, rRes1, 1, pro1 );
513 propagatorGS( mass2sq, mass2[i], width2[i], s24, sPip2, sPi0, rRes2, pro0 );
514
515 temp_PDF = Sf4;
516
517 B[0] = barrier( 1, sD, s234, sPip1, rD2, 1.9683 );
518 B[1] = barrier( 1, s234, s24, sPim, rRes1, mass1[i] );
519 B[2] = barrier( 1, s24, sPip2, sPi0, rRes2, mass2[i] );
520 tmp2 = B[0] * B[1] * B[2] * temp_PDF;
521 amp_tmp2[0] += tmp2 * pro0[0] * pro1[0];
522 amp_tmp2[1] += tmp2 * pro0[1] * pro1[1];
523
524 temp_PDF = 0;
525 propagatorRBW( mass1sq, mass1[i], width1[i], s134, s14, sPim, rRes1, 1, pro1 );
526 propagatorGS( mass2sq, mass2[i], width2[i], s13, sPip1, sPim, rRes2, pro0 );
527
528 temp_PDF = Sf5;
529
530 B[0] = barrier( 1, sD, s134, sPip2, rD2, 1.9683 );
531 B[1] = barrier( 1, s134, s13, sPi0, rRes1, mass1[i] );
532 B[2] = barrier( 1, s13, sPip1, sPim, rRes2, mass2[i] );
533 tmp1 = B[0] * B[1] * B[2] * temp_PDF;
534 amp_tmp1[0] -= tmp1 * pro0[0] * pro1[0];
535 amp_tmp1[1] -= tmp1 * pro0[1] * pro1[1];
536
537 temp_PDF = 0;
538 propagatorRBW( mass1sq, mass1[i], width1[i], s234, s24, sPim, rRes1, 1, pro1 );
539 propagatorGS( mass2sq, mass2[i], width2[i], s23, sPip2, sPim, rRes2, pro0 );
540
541 temp_PDF = Sf6;
542
543 B[0] = barrier( 1, sD, s234, sPip1, rD2, 1.9683 );
544 B[1] = barrier( 1, s234, s23, sPi0, rRes1, mass1[i] );
545 B[2] = barrier( 1, s23, sPip2, sPim, rRes2, mass2[i] );
546 tmp2 = B[0] * B[1] * B[2] * temp_PDF;
547 amp_tmp2[0] -= tmp2 * pro0[0] * pro1[0];
548 amp_tmp2[1] -= tmp2 * pro0[1] * pro1[1];
549
550 temp_PDF = 0;
551 propagatorRBW( mass1sq, mass1[i], width1[i], s134, s34, sPip1, rRes1, 1, pro1 );
552 propagatorGS( mass2sq, mass2[i], width2[i], s34, sPim, sPi0, rRes2, pro0 );
553
554 temp_PDF = Sf7;
555
556 B[0] = barrier( 1, sD, s134, sPip2, rD2, 1.9683 );
557 B[1] = barrier( 1, s134, s34, sPip1, rRes1, mass1[i] );
558 B[2] = barrier( 1, s34, sPim, sPi0, rRes2, mass2[i] );
559 tmp1 = B[0] * B[1] * B[2] * temp_PDF;
560 amp_tmp1[0] -= tmp1 * pro0[0] * pro1[0];
561 amp_tmp1[1] -= tmp1 * pro0[1] * pro1[1];
562
563 temp_PDF = 0;
564 propagatorRBW( mass1sq, mass1[i], width1[i], s234, s34, sPip2, rRes1, 1, pro1 );
565 propagatorGS( mass2sq, mass2[i], width2[i], s34, sPim, sPi0, rRes2, pro0 );
566
567 temp_PDF = Sf8;
568
569 B[0] = barrier( 1, sD, s234, sPip1, rD2, 1.9683 );
570 B[1] = barrier( 1, s234, s34, sPip2, rRes1, mass1[i] );
571 B[2] = barrier( 1, s34, sPim, sPi0, rRes2, mass2[i] );
572 tmp2 = B[0] * B[1] * B[2] * temp_PDF;
573 amp_tmp2[0] -= tmp2 * pro0[0] * pro1[0];
574 amp_tmp2[1] -= tmp2 * pro0[1] * pro1[1];
575 }
576
577 if ( i == 4 )
578 {
579 temp_PDF = 0;
580 // propagatorRBW(mass1sq,mass1[i],width1[i],s123,s13,sPip2,rRes1,0,pro1);
581 propagatorMW( 1.230, s123, 1, pro1 );
582 propagatorGS( mass2sq, mass2[i], width2[i], s13, sPip1, sPim, rRes2, pro0 );
583
584 for ( int a = 0; a < 4; a++ )
585 for ( int b = 0; b < 4; b++ )
586 temp_PDF += G[a][a] * G[b][b] * t1_123_4[a] *
587 ( -G[a][b] + p123[a] * p123[b] / s123 ) * t1_13[b];
588
589 B[0] = barrier( 1, sD, s123, sPi0, rD2, 1.9683 );
590 B[1] = barrier( 1, s13, sPip1, sPim, rRes2, mass2[i] );
591 tmp1 = B[0] * B[1] * temp_PDF;
592 amp_tmp1[0] = tmp1 * pro0[0] * pro1[0];
593 amp_tmp1[1] = tmp1 * pro0[1] * pro1[1];
594
595 temp_PDF = 0;
596 // propagatorRBW(mass1sq,mass1[i],width1[i],s123,s23,sPip1,rRes1,0,pro1);
597 propagatorMW( 1.230, s123, 1, pro1 );
598 propagatorGS( mass2sq, mass2[i], width2[i], s23, sPip2, sPim, rRes2, pro0 );
599
600 for ( int a = 0; a < 4; a++ )
601 for ( int b = 0; b < 4; b++ )
602 temp_PDF += G[a][a] * G[b][b] * t1_123_4[a] *
603 ( -G[a][b] + p123[a] * p123[b] / s123 ) * t1_23[b];
604
605 B[0] = barrier( 1, sD, s123, sPi0, rD2, 1.9683 );
606 B[1] = barrier( 1, s23, sPip2, sPim, rRes2, mass2[i] );
607 tmp2 = B[0] * B[1] * temp_PDF;
608 amp_tmp2[0] = tmp2 * pro0[0] * pro1[0];
609 amp_tmp2[1] = tmp2 * pro0[1] * pro1[1];
610 }
611
612 if ( i == 5 )
613 {
614 amp_tmp1[0] = 0;
615 amp_tmp1[1] = 0;
616 amp_tmp2[0] = 0;
617 amp_tmp2[1] = 0;
618
619 temp_PDF = 0;
620 // propagatorRBW(mass1sq,mass1[i],width1[i],s134,s14,sPim,rRes1,0,pro1);
621 propagatorMW( 1.230, s134, 0, pro1 );
622 propagatorGS( mass2sq, mass2[i], width2[i], s14, sPip1, sPi0, rRes2, pro0 );
623
624 for ( int a = 0; a < 4; a++ )
625 for ( int b = 0; b < 4; b++ )
626 temp_PDF += G[a][a] * G[b][b] * t1_134_2[a] *
627 ( -G[a][b] + p134[a] * p134[b] / s134 ) * t1_14[b];
628
629 B[0] = barrier( 1, sD, s134, sPip2, rD2, 1.9683 );
630 B[1] = barrier( 1, s14, sPip1, sPi0, rRes2, mass2[i] );
631 tmp1 = B[0] * B[1] * temp_PDF;
632 amp_tmp1[0] += tmp1 * pro0[0] * pro1[0] * 1.;
633 amp_tmp1[1] += tmp1 * pro0[1] * pro1[1] * 1.;
634
635 temp_PDF = 0;
636 // propagatorRBW(mass1sq,mass1[i],width1[i],s234,s24,sPim,rRes1,0,pro1);
637 propagatorMW( 1.230, s234, 0, pro1 );
638 propagatorGS( mass2sq, mass2[i], width2[i], s24, sPip2, sPi0, rRes2, pro0 );
639
640 for ( int a = 0; a < 4; a++ )
641 for ( int b = 0; b < 4; b++ )
642 temp_PDF += G[a][a] * G[b][b] * t1_234_1[a] *
643 ( -G[a][b] + p234[a] * p234[b] / s234 ) * t1_24[b];
644
645 B[0] = barrier( 1, sD, s234, sPip1, rD2, 1.9683 );
646 B[1] = barrier( 1, s24, sPip2, sPi0, rRes2, mass2[i] );
647 tmp2 = B[0] * B[1] * temp_PDF;
648 amp_tmp2[0] += tmp2 * pro0[0] * pro1[0] * 1.;
649 amp_tmp2[1] += tmp2 * pro0[1] * pro1[1] * 1.;
650
651 // a10 -> rho0 pi0 violation of C-Parity
652 // temp_PDF=0;
653 // propagatorRBW(mass1sq,mass1[i],width1[i],s134,s14,sPim,rRes1,0,pro1);
654 // propagatorGS(mass2sq,mass2[i],width2[i],s13,sPip1,sPim,rRes2,pro0);
655 //
656 // for(int a=0; a<4; a++)
657 // for(int b=0;b<4;b++)
658 // temp_PDF +=
659 // G[a][a]*G[b][b]*t1_134_2[a]*(-G[a][b]+p134[a]*p134[b]/s134)*t1_13[b];
660 //
661 // B[0] = barrier(1,sD,s134,sPip2,rD2,1.9683);
662 // B[1] = barrier(1,s13,sPip1,sPim,rRes2,mass2[i]);
663 // tmp1 = B[0]*B[1]*temp_PDF;
664 // amp_tmp1[0] += tmp1*pro0[0]*pro1[0];
665 // amp_tmp1[1] += tmp1*pro0[1]*pro1[1];
666 //
667 // temp_PDF=0;
668 // propagatorRBW(mass1sq,mass1[i],width1[i],s234,s24,sPim,rRes1,0,pro1);
669 // propagatorGS(mass2sq,mass2[i],width2[i],s23,sPip2,sPim,rRes2,pro0);
670 //
671 // for(int a=0; a<4; a++)
672 // for(int b=0;b<4;b++)
673 // temp_PDF +=
674 // G[a][a]*G[b][b]*t1_234_1[a]*(-G[a][b]+p234[a]*p234[b]/s234)*t1_23[b];
675 //
676 // B[0] = barrier(1,sD,s234,sPip1,rD2,1.9683);
677 // B[1] = barrier(1,s23,sPip2,sPim,rRes2,mass2[i]);
678 // tmp2 = B[0]*B[1]*temp_PDF;
679 // amp_tmp2[0] += tmp2*pro0[0]*pro1[0];
680 // amp_tmp2[1] += tmp2*pro0[1]*pro1[1];
681
682 temp_PDF = 0;
683 // propagatorRBW(mass1sq,mass1[i],width1[i],s134,s34,sPip1,rRes1,0,pro1);
684 propagatorMW( 1.230, s134, 0, pro1 );
685 propagatorGS( mass2sq, mass2[i], width2[i], s34, sPim, sPi0, rRes2, pro0 );
686
687 for ( int a = 0; a < 4; a++ )
688 for ( int b = 0; b < 4; b++ )
689 temp_PDF += G[a][a] * G[b][b] * t1_134_2[a] *
690 ( -G[a][b] + p134[a] * p134[b] / s134 ) * t1_34[b];
691
692 B[0] = barrier( 1, sD, s134, sPip2, rD2, 1.9683 );
693 B[1] = barrier( 1, s34, sPim, sPi0, rRes2, mass2[i] );
694 tmp1 = B[0] * B[1] * temp_PDF;
695 amp_tmp1[0] += tmp1 * pro0[0] * pro1[0] * -1.;
696 amp_tmp1[1] += tmp1 * pro0[1] * pro1[1] * -1.;
697
698 temp_PDF = 0;
699 // propagatorRBW(mass1sq,mass1[i],width1[i],s234,s34,sPip2,rRes1,0,pro1);
700 propagatorMW( 1.230, s234, 0, pro1 );
701 propagatorGS( mass2sq, mass2[i], width2[i], s34, sPim, sPi0, rRes2, pro0 );
702
703 for ( int a = 0; a < 4; a++ )
704 for ( int b = 0; b < 4; b++ )
705 temp_PDF += G[a][a] * G[b][b] * t1_234_1[a] *
706 ( -G[a][b] + p234[a] * p234[b] / s234 ) * t1_34[b];
707
708 B[0] = barrier( 1, sD, s234, sPip1, rD2, 1.9683 );
709 B[1] = barrier( 1, s34, sPim, sPi0, rRes2, mass2[i] );
710 tmp2 = B[0] * B[1] * temp_PDF;
711 amp_tmp2[0] += tmp2 * pro0[0] * pro1[0] * -1.;
712 amp_tmp2[1] += tmp2 * pro0[1] * pro1[1] * -1.;
713 }
714
715 if ( i == 6 )
716 {
717 temp_PDF = 0;
718 propagatorRBW( mass1sq, mass1[i], width1[i], s23, sPip2, sPim, rRes1, 0, pro0 );
719 propagatorGS( mass2sq, mass2[i], width2[i], s14, sPip1, sPi0, rRes2, pro1 );
720 Com_Multi( pro0, pro1, pro );
721
722 for ( int a = 0; a < 4; a++ ) { temp_PDF += G[a][a] * t1_14_23[a] * t1_14[a]; }
723
724 B[1] = barrier( 1, s14, sPip1, sPi0, rRes2, mass2[i] );
725 B[2] = barrier( 1, sD, s23, s14, rD2, 1.9683 );
726 tmp1 = B[1] * B[2] * temp_PDF;
727 amp_tmp1[0] = tmp1 * pro[0];
728 amp_tmp1[1] = tmp1 * pro[1];
729
730 temp_PDF = 0;
731 propagatorRBW( mass1sq, mass1[i], width1[i], s13, sPip1, sPim, rRes1, 0, pro0 );
732 propagatorGS( mass2sq, mass2[i], width2[i], s24, sPip2, sPi0, rRes2, pro1 );
733 Com_Multi( pro0, pro1, pro );
734
735 for ( int a = 0; a < 4; a++ ) { temp_PDF += G[a][a] * t1_13_24[a] * t1_24[a]; }
736
737 B[1] = barrier( 1, s24, sPip2, sPi0, rRes2, mass2[i] );
738 B[2] = barrier( 1, sD, s13, s24, rD2, 1.9683 );
739 tmp2 = B[1] * B[2] * temp_PDF;
740 amp_tmp2[0] = tmp2 * pro[0];
741 amp_tmp2[1] = tmp2 * pro[1];
742 }
743
744 if ( i == 7 )
745 {
746 temp_PDF = 0;
747 propagatorRBW( mass1sq, mass1[i], width1[i], s23, sPip2, sPim, rRes1, 2, pro0 );
748 propagatorGS( mass2sq, mass2[i], width2[i], s14, sPip1, sPi0, rRes2, pro1 );
749 Com_Multi( pro0, pro1, pro );
750
751 for ( int a = 0; a < 4; a++ )
752 {
753 for ( int b = 0; b < 4; b++ )
754 { temp_PDF += G[a][a] * G[b][b] * t1_14_23[a] * t1_14[b] * t2_23[a][b]; }
755 }
756
757 B[0] = barrier( 2, s23, sPip2, sPim, rRes1, mass1[i] );
758 B[1] = barrier( 1, s14, sPip1, sPi0, rRes2, mass2[i] );
759 B[2] = barrier( 1, sD, s23, s14, rD2, 1.9683 );
760 tmp1 = B[0] * B[1] * B[2] * temp_PDF;
761 amp_tmp1[0] = tmp1 * pro[0];
762 amp_tmp1[1] = tmp1 * pro[1];
763
764 temp_PDF = 0;
765 propagatorRBW( mass1sq, mass1[i], width1[i], s13, sPip1, sPim, rRes1, 2, pro0 );
766 propagatorGS( mass2sq, mass2[i], width2[i], s24, sPip2, sPi0, rRes2, pro1 );
767 Com_Multi( pro0, pro1, pro );
768
769 for ( int a = 0; a < 4; a++ )
770 {
771 for ( int b = 0; b < 4; b++ )
772 { temp_PDF += G[a][a] * G[b][b] * t1_13_24[a] * t1_24[b] * t2_13[a][b]; }
773 }
774
775 B[0] = barrier( 2, s13, sPip1, sPim, rRes1, mass1[i] );
776 B[1] = barrier( 1, s24, sPip2, sPi0, rRes2, mass2[i] );
777 B[2] = barrier( 1, sD, s13, s24, rD2, 1.9683 );
778 tmp2 = B[0] * B[1] * B[2] * temp_PDF;
779 amp_tmp2[0] = tmp2 * pro[0];
780 amp_tmp2[1] = tmp2 * pro[1];
781 }
782
783 if ( i == 8 )
784 {
785 temp_PDF = 0;
786 propagatorRBW( mass1sq, mass1[i], width1[i], s14, sPip1, sPi0, rRes1, 1, pro0 );
787 propagatorGS( mass2sq, mass2[i], width2[i], s23, sPip2, sPim, rRes2, pro1 );
788
789 for ( int a = 0; a < 4; a++ ) { temp_PDF += G[a][a] * t1_14[a] * t1_23[a]; }
790
791 B[1] = barrier( 1, s14, sPip1, sPi0, rRes2, mass2[i] );
792 B[2] = barrier( 1, s23, sPip2, sPim, rRes2, mass2[i] );
793 tmp1 = B[1] * B[2] * temp_PDF;
794 amp_tmp1[0] = tmp1 * pro0[0] * pro1[0];
795 amp_tmp1[1] = tmp1 * pro0[1] * pro1[1];
796
797 temp_PDF = 0;
798 propagatorRBW( mass1sq, mass1[i], width1[i], s24, sPip2, sPi0, rRes1, 1, pro0 );
799 propagatorGS( mass2sq, mass2[i], width2[i], s13, sPip1, sPim, rRes2, pro1 );
800
801 for ( int a = 0; a < 4; a++ ) { temp_PDF += G[a][a] * t1_24[a] * t1_13[a]; }
802
803 B[1] = barrier( 1, s24, sPip2, sPi0, rRes2, mass2[i] );
804 B[2] = barrier( 1, s13, sPip1, sPim, rRes2, mass2[i] );
805 tmp2 = B[1] * B[2] * temp_PDF;
806 amp_tmp2[0] = tmp2 * pro0[0] * pro1[0];
807 amp_tmp2[1] = tmp2 * pro0[1] * pro1[1];
808 }
809
810 if ( i == 9 )
811 {
812 temp_PDF = 0;
813 propagatorGS( mass1sq, mass1[i], width1[i], s14, sPip1, sPi0, rRes1, pro0 );
814 propagatorRBW( mass2sq, mass2[i], width2[i], s23, sPip2, sPim, rRes2, 1, pro1 );
815
816 temp_PDF = Sf1;
817
818 B[0] = barrier( 1, sD, s14, s23, rD2, 1.9683 );
819 B[1] = barrier( 1, s14, sPip1, sPi0, rRes2, mass2[i] );
820 B[2] = barrier( 1, s23, sPip2, sPim, rRes2, mass2[i] );
821 tmp1 = B[1] * B[2] * B[0] * temp_PDF;
822 amp_tmp1[0] = tmp1 * pro0[0] * pro1[0];
823 amp_tmp1[1] = tmp1 * pro0[1] * pro1[1];
824
825 temp_PDF = 0;
826 propagatorGS( mass1sq, mass1[i], width1[i], s24, sPip2, sPi0, rRes1, pro0 );
827 propagatorRBW( mass2sq, mass2[i], width2[i], s13, sPip1, sPim, rRes2, 1, pro1 );
828
829 temp_PDF = Sf2;
830
831 B[0] = barrier( 1, sD, s24, s13, rD2, 1.9683 );
832 B[1] = barrier( 1, s24, sPip2, sPi0, rRes2, mass2[i] );
833 B[2] = barrier( 1, s13, sPip1, sPim, rRes2, mass2[i] );
834 tmp2 = B[1] * B[2] * B[0] * temp_PDF;
835 amp_tmp2[0] = tmp2 * pro0[0] * pro1[0];
836 amp_tmp2[1] = tmp2 * pro0[1] * pro1[1];
837 }
838
839 if ( i == 10 )
840 {
841 amp_tmp1[0] = 0;
842 amp_tmp1[1] = 0;
843 amp_tmp2[0] = 0;
844 amp_tmp2[1] = 0;
845
846 temp_PDF = 0;
847 propagatorRBW( mass1sq, mass1[i], width1[i], s134, s14, sPim, rRes1, 1, pro1 );
848 propagatorGS( mass2sq, mass2[i], width2[i], s14, sPip1, sPi0, rRes2, pro0 );
849
850 for ( int a = 0; a < 4; a++ ) temp_PDF += G[a][a] * t1_134_2[a] * t1_14[a];
851 // temp_PDF += G[a][a]*G[b][b]*t1_134_2[a]*(-G[a][b]+p134[a]*p134[b]/s134)*t1_14[b];
852
853 B[0] = barrier( 1, s134, s14, sPim, rRes1, mass1[i] );
854 B[1] = barrier( 1, s14, sPip1, sPi0, rRes2, mass2[i] );
855 tmp1 = B[0] * B[1] * temp_PDF;
856 amp_tmp1[0] += tmp1 * pro0[0] * pro1[0] * 1.;
857 amp_tmp1[1] += tmp1 * pro0[1] * pro1[1] * 1.;
858
859 temp_PDF = 0;
860 propagatorRBW( mass1sq, mass1[i], width1[i], s234, s24, sPim, rRes1, 1, pro1 );
861 propagatorGS( mass2sq, mass2[i], width2[i], s24, sPip2, sPi0, rRes2, pro0 );
862
863 for ( int a = 0; a < 4; a++ ) temp_PDF += G[a][a] * t1_234_1[a] * t1_24[a];
864 // temp_PDF += G[a][a]*G[b][b]*t1_234_1[a]*(-G[a][b]+p234[a]*p234[b]/s234)*t1_24[b];
865
866 B[0] = barrier( 1, s234, s24, sPim, rRes1, mass1[i] );
867 B[1] = barrier( 1, s24, sPip2, sPi0, rRes2, mass2[i] );
868 tmp2 = B[0] * B[1] * temp_PDF;
869 amp_tmp2[0] += tmp2 * pro0[0] * pro1[0] * 1.;
870 amp_tmp2[1] += tmp2 * pro0[1] * pro1[1] * 1.;
871
872 // pi1300 -> rho0 pi0 violation of C-Parity
873 // temp_PDF=0;
874 // propagatorRBW(mass1sq,mass1[i],width1[i],s134,s14,sPim,rRes1,0,pro1);
875 // propagatorGS(mass2sq,mass2[i],width2[i],s13,sPip1,sPim,rRes2,pro0);
876 //
877 // for(int a=0; a<4; a++)
878 // for(int b=0;b<4;b++)
879 // temp_PDF +=
880 // G[a][a]*G[b][b]*t1_134_2[a]*(-G[a][b]+p134[a]*p134[b]/s134)*t1_13[b];
881 //
882 // B[0] = barrier(1,sD,s134,sPip2,rD2,1.9683);
883 // B[1] = barrier(1,s13,sPip1,sPim,rRes2,mass2[i]);
884 // tmp1 = B[0]*B[1]*temp_PDF;
885 // amp_tmp1[0] += tmp1*pro0[0]*pro1[0];
886 // amp_tmp1[1] += tmp1*pro0[1]*pro1[1];
887 //
888 // temp_PDF=0;
889 // propagatorRBW(mass1sq,mass1[i],width1[i],s234,s24,sPim,rRes1,0,pro1);
890 // propagatorGS(mass2sq,mass2[i],width2[i],s23,sPip2,sPim,rRes2,pro0);
891 //
892 // for(int a=0; a<4; a++)
893 // for(int b=0;b<4;b++)
894 // temp_PDF +=
895 // G[a][a]*G[b][b]*t1_234_1[a]*(-G[a][b]+p234[a]*p234[b]/s234)*t1_23[b];
896 //
897 // B[0] = barrier(1,sD,s234,sPip1,rD2,1.9683);
898 // B[1] = barrier(1,s23,sPip2,sPim,rRes2,mass2[i]);
899 // tmp2 = B[0]*B[1]*temp_PDF;
900 // amp_tmp2[0] += tmp2*pro0[0]*pro1[0];
901 // amp_tmp2[1] += tmp2*pro0[1]*pro1[1];
902
903 temp_PDF = 0;
904 propagatorRBW( mass1sq, mass1[i], width1[i], s134, s34, sPip1, rRes1, 1, pro1 );
905 propagatorGS( mass2sq, mass2[i], width2[i], s34, sPim, sPi0, rRes2, pro0 );
906
907 for ( int a = 0; a < 4; a++ ) temp_PDF += G[a][a] * t1_134_2[a] * t1_34[a];
908 // temp_PDF += G[a][a]*G[b][b]*t1_134_2[a]*(-G[a][b]+p134[a]*p134[b]/s134)*t1_34[b];
909
910 B[0] = barrier( 1, s134, s34, sPip1, rRes1, mass1[i] );
911 B[1] = barrier( 1, s34, sPim, sPi0, rRes2, mass2[i] );
912 tmp1 = B[0] * B[1] * temp_PDF;
913 amp_tmp1[0] += tmp1 * pro0[0] * pro1[0] * -1.;
914 amp_tmp1[1] += tmp1 * pro0[1] * pro1[1] * -1.;
915
916 temp_PDF = 0;
917 propagatorRBW( mass1sq, mass1[i], width1[i], s234, s34, sPip2, rRes1, 1, pro1 );
918 propagatorGS( mass2sq, mass2[i], width2[i], s34, sPim, sPi0, rRes2, pro0 );
919
920 for ( int a = 0; a < 4; a++ ) temp_PDF += G[a][a] * t1_234_1[a] * t1_34[a];
921 // temp_PDF += G[a][a]*G[b][b]*t1_234_1[a]*(-G[a][b]+p234[a]*p234[b]/s234)*t1_34[b];
922
923 B[0] = barrier( 1, s234, s34, sPip2, rRes1, mass1[i] );
924 B[1] = barrier( 1, s34, sPim, sPi0, rRes2, mass2[i] );
925 tmp2 = B[0] * B[1] * temp_PDF;
926 amp_tmp2[0] += tmp2 * pro0[0] * pro1[0] * -1.;
927 amp_tmp2[1] += tmp2 * pro0[1] * pro1[1] * -1.;
928 }
929
930 amp_tmp[0] = amp_tmp1[0] + amp_tmp2[0];
931 amp_tmp[1] = amp_tmp1[1] + amp_tmp2[1];
932 // printf("amp_tmp = %f , %f\n", amp_tmp[0], amp_tmp[1]);
933 // printf("cof = %f , %f\n", cof[0], cof[1]);
934 Com_Multi( amp_tmp, cof, amp_PDF );
935 PDF[0] += amp_PDF[0];
936 PDF[1] += amp_PDF[1];
937 // printf("PDF = %f , %f\n", PDF[0], PDF[1]);
938 }
939
940 double value = PDF[0] * PDF[0] + PDF[1] * PDF[1];
941 // value = value/Gen_eva[id];
942 value = value;
943 if ( value <= 0 ) value = 1e-20;
944 // if(value > 0) printf("max prob = %f\n",value);
945
946 setProb( value );
947
948 return;
949}
950
951double EvtDsToPipPipPimPi0::a10gam_arr[3000] = {
952 0, 0, 0, 0, 0, 0, 0,
953 0, 0, 0, 0, 0, 0, 0,
954 0, 0, 0, 0, 0, 0, 0,
955 0, 0, 0, 0, 0, 0, 0,
956 0, 0, 0, 0, 0, 0, 0,
957 0, 0, 0, 0, 0, 0, 0,
958 0, 0, 0, 0, 0, 0, 0,
959 0, 0, 0, 0, 0, 0, 0,
960 0, 0, 0, 0, 0, 0, 0,
961 0, 0, 0, 0, 0, 0, 0,
962 0, 0, 0, 0, 0, 0, 0,
963 0, 0, 0, 0, 0, 0, 0,
964 0, 0, 0, 0, 0, 0, 0,
965 0, 0, 0, 0, 0, 0, 0,
966 0, 0, 0, 0, 0, 0, 0,
967 0, 0, 0, 0, 0, 0, 0,
968 0, 0, 0, 0, 0, 0, 0,
969 0, 0, 0, 0, 0, 0, 0,
970 0, 0, 0, 0, 0, 0, 0,
971 0, 0, 0, 0, 0, 0, 0,
972 0, 0, 0, 0, 0, 0, 0,
973 0, 0, 0, 0, 0, 0, 0,
974 0, 0, 0, 0, 0, 0, 0,
975 0, 0, 0, 0, 0, 0, 0,
976 0, 0, 0, 0, 0, 0, 0,
977 0, 0, 0, 0, 0, 0, 0,
978 0, 0, 0, 0, 0, 0, 0,
979 0, 0, 0, 0, 0, 0, 0,
980 0, 0, 0, 0, 0, 0, 0,
981 0, 0, 0, 0, 0, 0, 0,
982 0, 0, 0, 0, 0, 0, 0,
983 0, 0, 0, 0, 0, 0, 0,
984 0, 0, 0, 0, 0, 0, 0,
985 0, 0, 0, 0, 0, 0, 0,
986 0, 0, 0, 0, 0, 0, 0,
987 0, 0, 0, 0, 0, 0, 0,
988 0, 0, 0, 0, 0, 0, 0,
989 0, 0, 0, 0, 0, 0, 0,
990 0, 0, 0, 0, 0, 0, 0,
991 0, 0, 0, 0, 0, 0, 0,
992 0, 0, 0, 0, 0, 0, 0,
993 0, 0, 0, 0, 0, 0, 0,
994 0, 0, 0, 0, 0, 0, 0,
995 0, 0, 0, 0, 0, 0, 0,
996 0, 0, 0, 0, 0, 0, 0,
997 0, 0, 0, 0, 0, 0, 0,
998 0, 0, 0, 0, 0, 0, 0,
999 0, 0, 0, 0, 0, 0, 0,
1000 0, 0, 0, 0, 0, 0, 0,
1001 0, 0, 0, 0, 0, 0, 0,
1002 0, 0, 0, 0, 0, 0, 0,
1003 0, 0, 0, 0, 0, 0, 0,
1004 0, 0, 0, 0, 0, 0, 0,
1005 0, 0, 0, 0, 0, 0, 0,
1006 0, 0, 0, 0, 0, 0, 0,
1007 0, 0, 0, 0, 0, 0, 0,
1008 0, 0, 0, 0, 0, 0, 0,
1009 0, 0, 0, 0, 0, 0, 0,
1010 0, 0, 0, 0, 0, 0, 0,
1011 0, 0, 0, 0, 0, 0, 0,
1012 0, 0, 0, 0, 0, 0, 5.10847e-07,
1013 6.61977e-07, 8.41244e-07, 1.04861e-06, 1.29389e-06, 1.5632e-06, 1.881e-06, 2.22767e-06,
1014 2.61443e-06, 3.05012e-06, 3.52012e-06, 4.05306e-06, 4.62497e-06, 5.24606e-06, 5.92555e-06,
1015 6.68573e-06, 7.48203e-06, 8.30439e-06, 9.23754e-06, 1.02432e-05, 1.12879e-05, 1.24557e-05,
1016 1.35884e-05, 1.49144e-05, 1.62366e-05, 1.76005e-05, 1.91314e-05, 2.0737e-05, 2.23708e-05,
1017 2.41359e-05, 2.60124e-05, 2.79615e-05, 3.00748e-05, 3.20951e-05, 3.4398e-05, 3.67148e-05,
1018 3.89229e-05, 4.16021e-05, 4.43174e-05, 4.7257e-05, 4.97924e-05, 5.28634e-05, 5.60937e-05,
1019 5.90265e-05, 6.25579e-05, 6.59154e-05, 6.96911e-05, 7.34966e-05, 7.68773e-05, 8.10852e-05,
1020 8.52922e-05, 8.93417e-05, 9.39913e-05, 9.85651e-05, 0.000103419, 0.000108147, 0.000112791,
1021 0.00011776, 0.000123403, 0.000128711, 0.000134492, 0.000140132, 0.00014662, 0.000152356,
1022 0.000158013, 0.000164076, 0.000170779, 0.000177999, 0.00018501, 0.000192253, 0.000199358,
1023 0.000207475, 0.000214667, 0.000222587, 0.000231224, 0.000238213, 0.000247703, 0.000256315,
1024 0.000265662, 0.000273847, 0.000283963, 0.000293846, 0.000303285, 0.000311502, 0.000322079,
1025 0.000334168, 0.000343724, 0.000354855, 0.000366474, 0.000377177, 0.000388813, 0.00040168,
1026 0.000413493, 0.000424643, 0.000437264, 0.000450301, 0.000463383, 0.000476909, 0.000489312,
1027 0.00050142, 0.000517696, 0.000534774, 0.000547504, 0.000562717, 0.000579844, 0.000591737,
1028 0.000609933, 0.00062713, 0.000642629, 0.000659459, 0.000675847, 0.000692308, 0.00071141,
1029 0.000729495, 0.000744429, 0.000766867, 0.000785325, 0.000802652, 0.00082197, 0.000843976,
1030 0.000863138, 0.000884998, 0.000906708, 0.000929898, 0.000947559, 0.000968399, 0.000996404,
1031 0.00101473, 0.00104141, 0.00106038, 0.00108682, 0.00111357, 0.00113832, 0.00116272,
1032 0.00118918, 0.00121012, 0.00123639, 0.00126683, 0.00129586, 0.00132024, 0.00135298,
1033 0.00137635, 0.00140952, 0.00144097, 0.00146779, 0.00149647, 0.00152659, 0.00155896,
1034 0.00158477, 0.00162065, 0.00165104, 0.00168977, 0.00171834, 0.00175372, 0.00178951,
1035 0.0018246, 0.00185573, 0.00190013, 0.00193404, 0.00197543, 0.00201023, 0.0020423,
1036 0.00208759, 0.00211846, 0.00216766, 0.00220693, 0.00224472, 0.00228621, 0.00232947,
1037 0.00237215, 0.00242014, 0.00246763, 0.00250157, 0.00254408, 0.00259896, 0.00264598,
1038 0.00268984, 0.00274089, 0.0027876, 0.00284181, 0.00289009, 0.00294365, 0.00298791,
1039 0.00303216, 0.00308017, 0.00314244, 0.00319637, 0.00325229, 0.00330926, 0.00335885,
1040 0.00342271, 0.00346827, 0.00353228, 0.00359511, 0.00364882, 0.00371527, 0.00378434,
1041 0.00383893, 0.0039045, 0.00396506, 0.00403849, 0.00409357, 0.00416097, 0.00422666,
1042 0.00429508, 0.00437824, 0.00442343, 0.00451049, 0.00456903, 0.00463377, 0.00471978,
1043 0.00479821, 0.00486691, 0.00495857, 0.00501581, 0.00511116, 0.00517971, 0.00524602,
1044 0.00533189, 0.00541212, 0.00550358, 0.00559244, 0.00566766, 0.00575527, 0.00586013,
1045 0.005914, 0.00601374, 0.00610486, 0.00616553, 0.0062721, 0.00637741, 0.00645871,
1046 0.00657222, 0.00665231, 0.00676846, 0.00685208, 0.00694391, 0.00704976, 0.00716695,
1047 0.00722817, 0.0073624, 0.00745222, 0.00757624, 0.00765055, 0.00778373, 0.00788613,
1048 0.00799569, 0.00810135, 0.00821636, 0.00832424, 0.00842107, 0.00855451, 0.0086506,
1049 0.00879026, 0.00889932, 0.00902619, 0.00915076, 0.00924049, 0.00937964, 0.009491,
1050 0.00962401, 0.00978053, 0.00988756, 0.0100305, 0.0101353, 0.0102867, 0.0104441,
1051 0.0105558, 0.0107242, 0.0108266, 0.0109571, 0.0110897, 0.0112405, 0.0113613,
1052 0.0114894, 0.0116503, 0.0118025, 0.0119203, 0.0121199, 0.0122381, 0.0124294,
1053 0.0125371, 0.0126689, 0.0128119, 0.0130044, 0.0131396, 0.013325, 0.0134785,
1054 0.0136555, 0.0137881, 0.0139677, 0.0141636, 0.0142625, 0.0144934, 0.0146466,
1055 0.0148004, 0.0150156, 0.0151146, 0.0153291, 0.0154632, 0.015703, 0.0158963,
1056 0.0160075, 0.0162392, 0.0164283, 0.0165921, 0.0167503, 0.0170055, 0.0171526,
1057 0.0174251, 0.017558, 0.0177089, 0.017877, 0.0181561, 0.0183833, 0.018488,
1058 0.0186674, 0.0188561, 0.0191999, 0.0193798, 0.0195454, 0.0197736, 0.0199596,
1059 0.0201606, 0.0204273, 0.02067, 0.0207535, 0.0210466, 0.0212628, 0.0213996,
1060 0.0217146, 0.0218709, 0.0221903, 0.0225007, 0.0226146, 0.0229532, 0.0230756,
1061 0.0233669, 0.0235796, 0.0238283, 0.024111, 0.0243377, 0.0245557, 0.0249343,
1062 0.0250247, 0.0252691, 0.0255234, 0.0257978, 0.026107, 0.0263546, 0.0266509,
1063 0.0268604, 0.0271287, 0.0275146, 0.0278669, 0.0280236, 0.0282697, 0.0286572,
1064 0.0289781, 0.0291602, 0.029455, 0.0298135, 0.0299364, 0.0305142, 0.0306114,
1065 0.0309928, 0.0312366, 0.0316771, 0.031846, 0.0321499, 0.0324471, 0.0327513,
1066 0.0330496, 0.0335854, 0.0339567, 0.0340194, 0.0344356, 0.0348833, 0.0351823,
1067 0.035629, 0.0358131, 0.0362275, 0.0364262, 0.036772, 0.0372882, 0.0375451,
1068 0.0379141, 0.0385863, 0.0387748, 0.0394341, 0.0397017, 0.0399935, 0.0402568,
1069 0.0407253, 0.0411851, 0.0415508, 0.0418898, 0.0423181, 0.0429063, 0.0433494,
1070 0.0437551, 0.0440054, 0.0447398, 0.0450333, 0.0454889, 0.0460102, 0.0462613,
1071 0.0467099, 0.0470934, 0.0478124, 0.0483974, 0.0489387, 0.0491559, 0.0494822,
1072 0.0501761, 0.0509158, 0.0514123, 0.0519352, 0.0522839, 0.0529311, 0.0533958,
1073 0.0536812, 0.0545877, 0.0552246, 0.0558031, 0.0563144, 0.0567699, 0.0573668,
1074 0.0582165, 0.0589321, 0.0591083, 0.059916, 0.0604653, 0.0611959, 0.0618002,
1075 0.06273, 0.0631461, 0.0639283, 0.064459, 0.0651378, 0.06585, 0.0664203,
1076 0.0674299, 0.0676451, 0.0689073, 0.0695797, 0.0705521, 0.070957, 0.0719581,
1077 0.0723552, 0.0731694, 0.0738023, 0.0746405, 0.0756375, 0.0768234, 0.0774481,
1078 0.0784771, 0.0786407, 0.0799796, 0.0808581, 0.0818104, 0.0827765, 0.083604,
1079 0.0844563, 0.0857548, 0.086459, 0.0877274, 0.0878173, 0.0897353, 0.0899906,
1080 0.0912604, 0.0921967, 0.093305, 0.094176, 0.0956489, 0.0962456, 0.0973797,
1081 0.0978631, 0.0997736, 0.100602, 0.101641, 0.103275, 0.104253, 0.104856,
1082 0.106542, 0.107601, 0.108657, 0.109335, 0.111208, 0.112411, 0.112519,
1083 0.115234, 0.115867, 0.117292, 0.11831, 0.119598, 0.121313, 0.122108,
1084 0.123445, 0.124304, 0.12644, 0.12744, 0.127965, 0.128983, 0.131425,
1085 0.131891, 0.132958, 0.135072, 0.136372, 0.137859, 0.139469, 0.139895,
1086 0.141239, 0.142312, 0.144194, 0.145061, 0.14797, 0.148421, 0.149148,
1087 0.151263, 0.152417, 0.153852, 0.154373, 0.155729, 0.15729, 0.158637,
1088 0.160362, 0.162451, 0.162713, 0.165762, 0.164869, 0.166916, 0.168694,
1089 0.170652, 0.171882, 0.172525, 0.174622, 0.17715, 0.176453, 0.17879,
1090 0.180351, 0.181534, 0.182003, 0.184618, 0.184791, 0.185818, 0.187754,
1091 0.18921, 0.19104, 0.191221, 0.192176, 0.194501, 0.196039, 0.197454,
1092 0.19889, 0.200492, 0.202103, 0.202012, 0.203932, 0.206371, 0.207377,
1093 0.207238, 0.207741, 0.210689, 0.212003, 0.213854, 0.214701, 0.215246,
1094 0.217599, 0.220232, 0.219966, 0.219449, 0.221044, 0.222879, 0.224735,
1095 0.225572, 0.225291, 0.228986, 0.229253, 0.231717, 0.232168, 0.232043,
1096 0.235528, 0.235789, 0.23709, 0.236679, 0.240255, 0.239063, 0.241123,
1097 0.242824, 0.245739, 0.245582, 0.246838, 0.248066, 0.25014, 0.250404,
1098 0.25013, 0.252638, 0.254109, 0.254707, 0.257212, 0.258817, 0.259441,
1099 0.259976, 0.261818, 0.262706, 0.263087, 0.265652, 0.266072, 0.267088,
1100 0.268155, 0.268635, 0.271209, 0.27077, 0.271821, 0.27319, 0.274764,
1101 0.275037, 0.276091, 0.278019, 0.277085, 0.281357, 0.28164, 0.28116,
1102 0.283233, 0.284078, 0.285207, 0.284973, 0.287791, 0.288662, 0.288786,
1103 0.290711, 0.292558, 0.290947, 0.293373, 0.294503, 0.295762, 0.297076,
1104 0.297324, 0.29907, 0.299182, 0.300861, 0.301929, 0.301915, 0.303147,
1105 0.305086, 0.306043, 0.307639, 0.307859, 0.308634, 0.309576, 0.309839,
1106 0.309856, 0.3122, 0.312669, 0.313976, 0.315044, 0.317234, 0.318186,
1107 0.318364, 0.319883, 0.319694, 0.320653, 0.320546, 0.321426, 0.32232,
1108 0.323426, 0.3238, 0.326689, 0.327256, 0.327845, 0.329071, 0.330415,
1109 0.330224, 0.331656, 0.332284, 0.332526, 0.333567, 0.334426, 0.336384,
1110 0.337273, 0.337489, 0.33854, 0.338955, 0.339142, 0.341662, 0.342302,
1111 0.343481, 0.343899, 0.34338, 0.345925, 0.346405, 0.34765, 0.348484,
1112 0.347483, 0.349126, 0.34929, 0.350908, 0.352362, 0.352584, 0.353779,
1113 0.353871, 0.354532, 0.354155, 0.355416, 0.356756, 0.358541, 0.3579,
1114 0.3592, 0.360009, 0.361232, 0.361467, 0.363298, 0.364186, 0.366374,
1115 0.364917, 0.364327, 0.366702, 0.367645, 0.367933, 0.369071, 0.370738,
1116 0.367952, 0.371982, 0.37071, 0.370969, 0.37351, 0.3725, 0.373201,
1117 0.374315, 0.376762, 0.376296, 0.378849, 0.377269, 0.378157, 0.379463,
1118 0.380044, 0.380239, 0.380691, 0.380712, 0.381581, 0.382821, 0.382769,
1119 0.383597, 0.385241, 0.385573, 0.386094, 0.387948, 0.388176, 0.388829,
1120 0.387262, 0.389284, 0.390265, 0.390807, 0.392308, 0.392241, 0.393039,
1121 0.392672, 0.393517, 0.393756, 0.395187, 0.394738, 0.396592, 0.397667,
1122 0.397501, 0.397686, 0.398778, 0.397756, 0.399491, 0.400856, 0.4008,
1123 0.402307, 0.40258, 0.402486, 0.403773, 0.403612, 0.405249, 0.406007,
1124 0.406182, 0.405619, 0.406478, 0.40722, 0.407959, 0.409399, 0.407149,
1125 0.410364, 0.411445, 0.41226, 0.409334, 0.41238, 0.413182, 0.411058,
1126 0.413268, 0.414802, 0.415712, 0.414306, 0.414652, 0.417959, 0.415588,
1127 0.417994, 0.41751, 0.417698, 0.419487, 0.419246, 0.42, 0.421417,
1128 0.421104, 0.422509, 0.421368, 0.421375, 0.422968, 0.422329, 0.424118,
1129 0.424895, 0.426347, 0.424783, 0.425062, 0.427771, 0.4272, 0.427332,
1130 0.428065, 0.428681, 0.428306, 0.428324, 0.42934, 0.431072, 0.43031,
1131 0.431238, 0.431038, 0.431953, 0.432834, 0.433656, 0.435662, 0.43307,
1132 0.434848, 0.43479, 0.435059, 0.434561, 0.43494, 0.437631, 0.436843,
1133 0.436721, 0.438563, 0.436746, 0.438498, 0.439264, 0.440454, 0.439462,
1134 0.441818, 0.441259, 0.440586, 0.441971, 0.442086, 0.441904, 0.443643,
1135 0.442794, 0.445147, 0.44321, 0.444166, 0.444704, 0.445945, 0.444766,
1136 0.44681, 0.446487, 0.446334, 0.448603, 0.448152, 0.448447, 0.449051,
1137 0.449171, 0.450434, 0.449876, 0.451072, 0.450965, 0.450642, 0.451657,
1138 0.451473, 0.451965, 0.454273, 0.452422, 0.453766, 0.455101, 0.455226,
1139 0.452829, 0.455623, 0.456213, 0.453769, 0.456463, 0.458039, 0.45612,
1140 0.456609, 0.458332, 0.458103, 0.457758, 0.457209, 0.460109, 0.458361,
1141 0.459533, 0.459841, 0.461629, 0.462698, 0.462221, 0.460548, 0.461829,
1142 0.463798, 0.463285, 0.465082, 0.465756, 0.463391, 0.465083, 0.462681,
1143 0.464189, 0.466095, 0.46847, 0.466325, 0.466708, 0.465522, 0.468529,
1144 0.468893, 0.467525, 0.469223, 0.46837, 0.468224, 0.470765, 0.468767,
1145 0.469725, 0.47122, 0.469998, 0.471174, 0.471, 0.470084, 0.472925,
1146 0.473875, 0.47366, 0.475595, 0.472112, 0.474616, 0.475303, 0.474338,
1147 0.475903, 0.475051, 0.475786, 0.475884, 0.475342, 0.477809, 0.476628,
1148 0.478056, 0.477943, 0.479718, 0.480277, 0.479325, 0.478487, 0.478945,
1149 0.481972, 0.482123, 0.482074, 0.481763, 0.48447, 0.484614, 0.481945,
1150 0.483892, 0.487872, 0.484796, 0.485473, 0.489195, 0.485496, 0.486415,
1151 0.486113, 0.486683, 0.489275, 0.488448, 0.488583, 0.488278, 0.492439,
1152 0.491196, 0.490612, 0.492102, 0.493087, 0.495125, 0.492948, 0.493526,
1153 0.49445, 0.494222, 0.495337, 0.495459, 0.495539, 0.494363, 0.496046,
1154 0.495832, 0.496912, 0.498911, 0.499184, 0.498805, 0.499241, 0.500922,
1155 0.501565, 0.499584, 0.501578, 0.502322, 0.501752, 0.502496, 0.502788,
1156 0.505379, 0.503861, 0.506021, 0.506488, 0.506952, 0.508643, 0.506462,
1157 0.507231, 0.507387, 0.509251, 0.50949, 0.508067, 0.508938, 0.508887,
1158 0.509544, 0.51172, 0.509405, 0.511026, 0.512919, 0.511507, 0.512892,
1159 0.51042, 0.513466, 0.514605, 0.51404, 0.515337, 0.514276, 0.514417,
1160 0.514953, 0.514552, 0.518115, 0.518749, 0.516392, 0.517777, 0.520258,
1161 0.51817, 0.517867, 0.519882, 0.51881, 0.519493, 0.518606, 0.520654,
1162 0.522035, 0.518634, 0.520473, 0.521678, 0.522549, 0.520957, 0.521824,
1163 0.523186, 0.521642, 0.52482, 0.523135, 0.525027, 0.525555, 0.525987,
1164 0.523265, 0.526671, 0.524975, 0.524577, 0.527919, 0.525689, 0.528478,
1165 0.528185, 0.52857, 0.529382, 0.527601, 0.529476, 0.527973, 0.531799,
1166 0.529369, 0.529891, 0.529205, 0.528191, 0.530105, 0.529386, 0.531998,
1167 0.529379, 0.534389, 0.531392, 0.53429, 0.534107, 0.532299, 0.535032,
1168 0.530939, 0.533253, 0.53627, 0.534767, 0.536838, 0.532772, 0.535239,
1169 0.536046, 0.534192, 0.537879, 0.536891, 0.53856, 0.539537, 0.538163,
1170 0.539638, 0.536308, 0.541708, 0.539197, 0.538073, 0.53991, 0.537786,
1171 0.542624, 0.538922, 0.540254, 0.539483, 0.53884, 0.539868, 0.540382,
1172 0.542084, 0.544621, 0.542604, 0.543015, 0.542192, 0.544019, 0.544245,
1173 0.543841, 0.542288, 0.545539, 0.544668, 0.542409, 0.543304, 0.54556,
1174 0.544934, 0.547257, 0.541864, 0.543311, 0.545769, 0.545772, 0.545936,
1175 0.546728, 0.547074, 0.5482, 0.548422, 0.549682, 0.549333, 0.549396,
1176 0.549981, 0.546348, 0.549856, 0.547825, 0.548367, 0.547422, 0.550924,
1177 0.551186, 0.547844, 0.550414, 0.551172, 0.54655, 0.549614, 0.55211,
1178 0.553047, 0.554943, 0.552244, 0.551954, 0.551853, 0.549702, 0.553917,
1179 0.55485, 0.551653, 0.553974, 0.553223, 0.553047, 0.552811, 0.5525,
1180 0.554389, 0.55657, 0.554466, 0.555791, 0.556522, 0.554262, 0.556422,
1181 0.558315, 0.557895, 0.558123, 0.554199, 0.55749, 0.557732, 0.556954,
1182 0.557359, 0.556586, 0.557766, 0.557517, 0.554655, 0.559342, 0.562112,
1183 0.559345, 0.557585, 0.55697, 0.559785, 0.562197, 0.56008, 0.558988,
1184 0.559874, 0.561188, 0.564173, 0.560372, 0.56124, 0.561108, 0.563741,
1185 0.558707, 0.561086, 0.563327, 0.562036, 0.564004, 0.561138, 0.563496,
1186 0.565594, 0.562067, 0.562929, 0.564893, 0.567015, 0.562481, 0.562019,
1187 0.563653, 0.561878, 0.566143, 0.5639, 0.564638, 0.563107, 0.566875,
1188 0.564398, 0.566127, 0.566607, 0.567252, 0.566528, 0.566628, 0.56831,
1189 0.564859, 0.569815, 0.566654, 0.564939, 0.568844, 0.568502, 0.565886,
1190 0.568556, 0.566762, 0.568847, 0.568534, 0.566041, 0.566397, 0.569942,
1191 0.568914, 0.567241, 0.570659, 0.57171, 0.56819, 0.57075, 0.568438,
1192 0.567604, 0.569452, 0.570547, 0.571437, 0.573075, 0.571073, 0.571652,
1193 0.568959, 0.569775, 0.570844, 0.570371, 0.572211, 0.570553, 0.571842,
1194 0.570598, 0.569651, 0.572605, 0.574228, 0.572185, 0.572898, 0.572446,
1195 0.574438, 0.57708, 0.57026, 0.574982, 0.574789, 0.570822, 0.573091,
1196 0.571822, 0.57595, 0.574218, 0.576605, 0.573796, 0.576091, 0.576826,
1197 0.575359, 0.575667, 0.575334, 0.574895, 0.572222, 0.572793, 0.576912,
1198 0.577936, 0.57358, 0.57554, 0.578902, 0.578798, 0.574994, 0.57923,
1199 0.57836, 0.57877, 0.578834, 0.574414, 0.575213, 0.571854, 0.578509,
1200 0.57951, 0.578496, 0.580245, 0.57465, 0.579298, 0.577356, 0.577914,
1201 0.579471, 0.58146, 0.578826, 0.578709, 0.578496, 0.577299, 0.57901,
1202 0.577267, 0.580718, 0.580655, 0.58123, 0.578204, 0.577744, 0.580787,
1203 0.579833, 0.578184, 0.585243, 0.580038, 0.579339, 0.578362, 0.5854,
1204 0.580457, 0.580576, 0.580914, 0.584755, 0.582103, 0.581257, 0.58125,
1205 0.58311, 0.580425, 0.58365, 0.579959, 0.585634, 0.58076, 0.585159,
1206 0.58226, 0.584189, 0.582075, 0.580576, 0.58043, 0.580405, 0.582512,
1207 0.585068, 0.58504, 0.583435, 0.584687, 0.584157, 0.582721, 0.585834,
1208 0.58427, 0.584793, 0.585872, 0.58757, 0.582234, 0.584769, 0.586066,
1209 0.584464, 0.586815, 0.58484, 0.58566, 0.588422, 0.588056, 0.584557,
1210 0.584627, 0.585923, 0.587775, 0.586758, 0.587405, 0.587869, 0.585873,
1211 0.585099, 0.582264, 0.584199, 0.590424, 0.588416, 0.589516, 0.58845,
1212 0.587837, 0.586989, 0.589098, 0.588176, 0.585571, 0.588833, 0.588869,
1213 0.590503, 0.585261, 0.590453, 0.587992, 0.586263, 0.589369, 0.589266,
1214 0.586583, 0.588911, 0.585802, 0.585045, 0.588152, 0.587742, 0.590043,
1215 0.591548, 0.589419, 0.591158, 0.588734, 0.591094, 0.592836, 0.590413,
1216 0.589916, 0.591878, 0.587068, 0.590096, 0.589218, 0.592117, 0.590057,
1217 0.590293, 0.59132, 0.593525, 0.588977, 0.589158, 0.588977, 0.589911,
1218 0.588572, 0.592564, 0.591573, 0.591105, 0.590711, 0.590143, 0.59068,
1219 0.592997, 0.591367, 0.594414, 0.592478, 0.590993, 0.592463, 0.590829,
1220 0.592986, 0.592069, 0.594267, 0.591337, 0.591242, 0.590092, 0.59447,
1221 0.588796, 0.594219, 0.592487, 0.591164, 0.589996, 0.593095, 0.592464,
1222 0.593464, 0.591231, 0.591955, 0.59051, 0.588754, 0.595189, 0.594417,
1223 0.598559, 0.597429, 0.591446, 0.597268, 0.592084, 0.593667, 0.594351,
1224 0.595908, 0.593284, 0.594093, 0.594685, 0.590196, 0.595528, 0.593947,
1225 0.596839, 0.593559, 0.595751, 0.595056, 0.593025, 0.595207, 0.599865,
1226 0.595802, 0.597093, 0.596852, 0.600151, 0.594165, 0.59669, 0.593693,
1227 0.597007, 0.59448, 0.59692, 0.598447, 0.58967, 0.594291, 0.596814,
1228 0.601778, 0.595606, 0.593575, 0.599196, 0.596763, 0.592731, 0.593649,
1229 0.596815, 0.596612, 0.598389, 0.600288, 0.601971, 0.593411, 0.595628,
1230 0.597471, 0.596617, 0.597107, 0.596928, 0.601419, 0.600712, 0.60062,
1231 0.597027, 0.598212, 0.598517, 0.598585, 0.598457, 0.595458, 0.601006,
1232 0.597436, 0.598442, 0.59978, 0.598042, 0.600751, 0.604877, 0.599682,
1233 0.602193, 0.600038, 0.602133, 0.601971, 0.603675, 0.598525, 0.597034,
1234 0.600436, 0.599307, 0.600178, 0.597561, 0.601128, 0.598386, 0.604022,
1235 0.600333, 0.596953, 0.596153, 0.601371, 0.59905, 0.598775, 0.598234,
1236 0.603056, 0.598461, 0.599035, 0.60247, 0.598961, 0.600638, 0.602174,
1237 0.603546, 0.602759, 0.595493, 0.599566, 0.603676, 0.601448, 0.601889,
1238 0.601923, 0.59879, 0.602152, 0.602868, 0.603142, 0.599779, 0.602996,
1239 0.603121, 0.604151, 0.601858, 0.599859, 0.605559, 0.600411, 0.605023,
1240 0.59481, 0.60354, 0.601028, 0.605539, 0.599386, 0.605815, 0.601308,
1241 0.60254, 0.605936, 0.598884, 0.602046, 0.604051, 0.604064, 0.600417,
1242 0.600486, 0.603358, 0.597937, 0.601253, 0.601612, 0.603221, 0.60466,
1243 0.604129, 0.604015, 0.605971, 0.604743, 0.603351, 0.602063, 0.603581,
1244 0.602386, 0.606447, 0.602734, 0.603914, 0.602157, 0.608614, 0.607969,
1245 0.604008, 0.607868, 0.607401, 0.605498, 0.604538, 0.605259, 0.605669,
1246 0.605429, 0.604667, 0.604459, 0.603132, 0.607398, 0.607448, 0.606535,
1247 0.606449, 0.603193, 0.603195, 0.605879, 0.604217, 0.602632, 0.603078,
1248 0.607168, 0.608578, 0.607912, 0.606881, 0.605766, 0.604638, 0.608067,
1249 0.606821, 0.606599, 0.603447, 0.608147, 0.608554, 0.608716, 0.613558,
1250 0.614474, 0.610203, 0.605614, 0.60737, 0.608189, 0.606754, 0.610616,
1251 0.603147, 0.603729, 0.611811, 0.606968, 0.60742, 0.606781, 0.604482,
1252 0.6097, 0.606153, 0.609647, 0.607856, 0.607672, 0.604834, 0.606836,
1253 0.611446, 0.608694, 0.611505, 0.609354, 0.608521, 0.60922, 0.60297,
1254 0.607662, 0.608237, 0.6063, 0.608534, 0.609033, 0.606615, 0.61312,
1255 0.613844, 0.608296, 0.614646, 0.61078, 0.607609, 0.606543, 0.607679,
1256 0.60831, 0.612263, 0.610579, 0.612233, 0.612682, 0.609124, 0.61026,
1257 0.606912, 0.609823, 0.608506, 0.610746, 0.610349, 0.611911, 0.608253,
1258 0.610799, 0.60775, 0.609174, 0.609446, 0.613876, 0.606614, 0.610548,
1259 0.607299, 0.617604, 0.611957, 0.608182, 0.61352, 0.610114, 0.607217,
1260 0.615439, 0.610774, 0.61353, 0.607945, 0.611571, 0.614412, 0.607529,
1261 0.612346, 0.615903, 0.613468, 0.615199, 0.612175, 0.612662, 0.60956,
1262 0.608311, 0.612912, 0.611678, 0.613576, 0.61002, 0.610685, 0.60852,
1263 0.609646, 0.612874, 0.611657, 0.608268, 0.609351, 0.613332, 0.610005,
1264 0.616228, 0.609752, 0.610077, 0.610468, 0.612632, 0.611267, 0.614137,
1265 0.617096, 0.61805, 0.615931, 0.610193, 0.606448, 0.610412, 0.614418,
1266 0.613298, 0.616366, 0.610521, 0.609473, 0.617223, 0.614392, 0.609531,
1267 0.611503, 0.619216, 0.612787, 0.616896, 0.613963, 0.613148, 0.611531,
1268 0.612897, 0.613263, 0.611573, 0.616155, 0.615285, 0.61387, 0.613969,
1269 0.61706, 0.615567, 0.614016, 0.611149, 0.611988, 0.613614, 0.613824,
1270 0.609791, 0.613618, 0.609961, 0.620689, 0.615341, 0.610682, 0.614414,
1271 0.61655, 0.615938, 0.615047, 0.61412, 0.614024, 0.611481, 0.613393,
1272 0.616073, 0.61456, 0.614876, 0.617353, 0.617323, 0.614166, 0.618364,
1273 0.615672, 0.615384, 0.616775, 0.614972, 0.615856, 0.612369, 0.613385,
1274 0.615188, 0.621464, 0.621196, 0.617695, 0.617198, 0.615838, 0.620844,
1275 0.614442, 0.616499, 0.611469, 0.61797, 0.613381, 0.612686, 0.617832,
1276 0.616912, 0.615726, 0.613965, 0.617923, 0.616926, 0.613962, 0.616878,
1277 0.616216, 0.615167, 0.615741, 0.621837, 0.622816, 0.618735, 0.614708,
1278 0.616259, 0.615296, 0.617336, 0.613412, 0.617638, 0.619523, 0.617039,
1279 0.614946, 0.620744, 0.618552, 0.616013, 0.617093, 0.616921, 0.618917,
1280 0.61971, 0.616581, 0.614899, 0.619123, 0.613554, 0.6208, 0.618294,
1281 0.619059, 0.620993, 0.619525, 0.621942, 0.614049, 0.616872, 0.617373,
1282 0.621464, 0.615137, 0.621412, 0.617671, 0.619178, 0.617229, 0.620786,
1283 0.618965, 0.619412, 0.624408, 0.620756, 0.622836, 0.619064, 0.622596,
1284 0.618601, 0.612976, 0.614824, 0.621118, 0.621563, 0.621251, 0.619092,
1285 0.62368, 0.620553, 0.625767, 0.621926, 0.618632, 0.616635, 0.616295,
1286 0.625815, 0.620596, 0.616814, 0.619908, 0.620009, 0.618273, 0.618491,
1287 0.621039, 0.619414, 0.618904, 0.617644, 0.622816, 0.620315, 0.622365,
1288 0.62554, 0.62318, 0.621576, 0.622506, 0.622526, 0.621412, 0.617583,
1289 0.627276, 0.619036, 0.618457, 0.620758, 0.618635, 0.618848, 0.616887,
1290 0.625413, 0.622702, 0.626354, 0.619957, 0.621203, 0.622965, 0.620388,
1291 0.624028, 0.619793, 0.619449, 0.619894, 0.622475, 0.625728, 0.620663,
1292 0.625763, 0.620276, 0.6219, 0.622332, 0.62029, 0.623742, 0.620927,
1293 0.616047, 0.62078, 0.619601, 0.625109, 0.623586, 0.618206, 0.624081,
1294 0.622982, 0.617744, 0.622281, 0.621567, 0.62139, 0.62237, 0.625921,
1295 0.625531, 0.625433, 0.621814, 0.624017, 0.622389, 0.623663, 0.623262,
1296 0.627376, 0.623763, 0.621025, 0.625044, 0.622423, 0.619603, 0.622101,
1297 0.61922, 0.626379, 0.622221, 0.625152, 0.6239, 0.626617, 0.617879,
1298 0.626182, 0.622782, 0.624102, 0.625998, 0.629432, 0.626768, 0.624665,
1299 0.623231, 0.62919, 0.620829, 0.620281, 0.624708, 0.626512, 0.622157,
1300 0.626986, 0.622909, 0.619203, 0.623121, 0.631429, 0.623336, 0.624622,
1301 0.623138, 0.62677, 0.623556, 0.629962, 0.626398, 0.622657, 0.62277,
1302 0.626691, 0.622954, 0.623641, 0.621758, 0.625917, 0.626179, 0.624352,
1303 0.629404, 0.624387, 0.623461, 0.627372, 0.624249, 0.62444, 0.625099,
1304 0.624841, 0.62775, 0.625561, 0.628993, 0.624342, 0.621275, 0.629437,
1305 0.632163, 0.623176, 0.624028, 0.624016, 0.624913, 0.62146, 0.625645,
1306 0.632107, 0.62442, 0.626444, 0.627954, 0.628726, 0.629863, 0.626985,
1307 0.628023, 0.62284, 0.634704, 0.622555, 0.628842, 0.633196, 0.624833,
1308 0.630645, 0.62962, 0.627391, 0.625646, 0.624605, 0.623669, 0.631903,
1309 0.627895, 0.631525, 0.623852, 0.635089, 0.632845, 0.626416, 0.630422,
1310 0.626313, 0.628298, 0.632386, 0.623754, 0.627711, 0.627182, 0.631031,
1311 0.632696, 0.6265, 0.626044, 0.631712, 0.631146, 0.630076, 0.6246,
1312 0.626339, 0.626073, 0.62547, 0.624333, 0.627483, 0.629035, 0.623781,
1313 0.626105, 0.63322, 0.628441, 0.625099, 0.631457, 0.631161, 0.628228,
1314 0.629891, 0.62785, 0.627987, 0.627642, 0.63232, 0.629908, 0.633618,
1315 0.631015, 0.631903, 0.627583, 0.627801, 0.631097, 0.63238, 0.634586,
1316 0.627661, 0.625959, 0.629329, 0.625851, 0.632545, 0.631533, 0.628986,
1317 0.631847, 0.632295, 0.641355, 0.627727, 0.634586, 0.625373, 0.627816,
1318 0.633512, 0.63269, 0.627166, 0.632661, 0.628999, 0.631664, 0.633816,
1319 0.631742, 0.633348, 0.630743, 0.63382, 0.635101, 0.632981, 0.6346,
1320 0.632103, 0.632281, 0.63134, 0.634626, 0.630111, 0.631351, 0.630224,
1321 0.630546, 0.632019, 0.632453, 0.631172, 0.632314, 0.631817, 0.62508,
1322 0.632514, 0.634765, 0.629544, 0.635109, 0.630736, 0.63048, 0.635906,
1323 0.635783, 0.624545, 0.630134, 0.628176, 0.630117, 0.63592, 0.634577,
1324 0.632769, 0.634533, 0.63141, 0.628919, 0.636203, 0.638389, 0.63231,
1325 0.632995, 0.629657, 0.630346, 0.632677, 0.634667, 0.634583, 0.629406,
1326 0.63471, 0.635601, 0.631947, 0.635692, 0.638141, 0.630747, 0.631889,
1327 0.635176, 0.633654, 0.62976, 0.636197, 0.633478, 0.634548, 0.63132,
1328 0.634575, 0.633873, 0.632214, 0.63282, 0.631591, 0.629881, 0.635444,
1329 0.636161, 0.634927, 0.637786, 0.636701, 0.630247, 0.6411, 0.640469,
1330 0.633986, 0.632693, 0.633078, 0.631706, 0.628018, 0.640413, 0.63462,
1331 0.631846, 0.639895, 0.634679, 0.629001, 0.641085, 0.636796, 0.636365,
1332 0.634147, 0.639004, 0.633508, 0.636331, 0.642073, 0.634742, 0.635549,
1333 0.639292, 0.633249, 0.629993, 0.639224, 0.635845, 0.635741, 0.631689,
1334 0.633921, 0.636043, 0.636046, 0.634273, 0.633504, 0.635015, 0.63648,
1335 0.634258, 0.631794, 0.637351, 0.637003, 0.635062, 0.641799, 0.630846,
1336 0.636755, 0.637128, 0.636238, 0.635672, 0.642322, 0.63981, 0.638706,
1337 0.640004, 0.640292, 0.640726, 0.640107, 0.637919, 0.63783, 0.642114,
1338 0.632441, 0.637189, 0.635049, 0.633693, 0.638689, 0.632568, 0.646298,
1339 0.631719, 0.633435, 0.633835, 0.636585, 0.637045, 0.637189, 0.634202,
1340 0.639623, 0.640147, 0.634653, 0.640554, 0.639917, 0.64078, 0.638127,
1341 0.63798, 0.637607, 0.639998, 0.63473, 0.643037, 0.636686, 0.63415,
1342 0.645426, 0.639442, 0.639877, 0.636225, 0.636739, 0.628628, 0.640501,
1343 0.646981, 0.637553, 0.63904, 0.639558, 0.632821, 0.641581, 0.634517,
1344 0.639993, 0.640098, 0.640046, 0.634301, 0.63906, 0.644402, 0.643639,
1345 0.638688, 0.643383, 0.63776, 0.639042, 0.639758, 0.640635, 0.641742,
1346 0.641212, 0.639423, 0.636773, 0.635784, 0.639141, 0.637001, 0.638086,
1347 0.639779, 0.633228, 0.641054, 0.636265, 0.642418, 0.644187, 0.640025,
1348 0.637431, 0.64234, 0.639466, 0.639929, 0.643658, 0.645672, 0.63875,
1349 0.638856, 0.643978, 0.647565, 0.647039, 0.642832, 0.640113, 0.642341,
1350 0.644558, 0.637816, 0.638742, 0.641432, 0.644949, 0.641804, 0.641167,
1351 0.646557, 0.64298, 0.642473, 0.641937, 0.644673, 0.638359, 0.639813,
1352 0.639308, 0.636239, 0.640583, 0.641604, 0.647463, 0.643723, 0.638941,
1353 0.640267, 0.641271, 0.648004, 0.647338, 0.640559, 0.640159, 0.639687,
1354 0.647811, 0.644427, 0.639744, 0.645816, 0.640913, 0.647267, 0.644434,
1355 0.645361, 0.640173, 0.647387, 0.643858, 0.638137, 0.642417, 0.647948,
1356 0.638894, 0.646922, 0.650392, 0.643519, 0.645082, 0.63969, 0.643681,
1357 0.647513, 0.641228, 0.642628, 0.644522, 0.645842, 0.647846, 0.64716,
1358 0.643352, 0.651159, 0.644077, 0.650025, 0.644696, 0.643421, 0.641637,
1359 0.647484, 0.640235, 0.643457, 0.642693, 0.645508, 0.642223, 0.639686,
1360 0.646631, 0.641975, 0.643323, 0.647343, 0.64402, 0.648346, 0.649236,
1361 0.647768, 0.645609, 0.641971, 0.641319, 0.641983, 0.644494, 0.644607,
1362 0.649457, 0.643873, 0.647442, 0.646043, 0.64612, 0.648729, 0.646065,
1363 0.652991, 0.645837, 0.650514, 0.651147, 0.644631, 0.645813, 0.643999,
1364 0.647581, 0.646826, 0.64706, 0.650602, 0.654695, 0.648685, 0.651021,
1365 0.651924, 0.649291, 0.652016, 0.654856, 0.644839, 0.644466, 0.648472,
1366 0.645495, 0.648701, 0.644017, 0.648088, 0.647046, 0.647111, 0.654831,
1367 0.649127, 0.643345, 0.643242, 0.651873, 0.649742, 0.648395, 0.643895,
1368 0.642581, 0.643208, 0.651363, 0.653665, 0.65346, 0.654502, 0.64849,
1369 0.643907, 0.653994, 0.648345, 0.646242, 0.647899, 0.651866, 0.648614,
1370 0.654628, 0.649, 0.650335, 0.643417, 0.645178, 0.653398, 0.65154,
1371 0.653001, 0.653917, 0.645493, 0.643492, 0.653624, 0.650601, 0.654645,
1372 0.650887, 0.653012, 0.651437, 0.652385, 0.655021, 0.645354, 0.643753,
1373 0.649658, 0.645741, 0.648658, 0.649114, 0.643689, 0.651045, 0.65228,
1374 0.653334, 0.65179, 0.654058, 0.647741, 0.650629, 0.654112, 0.654751,
1375 0.658177, 0.647574, 0.653494, 0.654006, 0.656714, 0.653971, 0.652711,
1376 0.648534, 0.65007, 0.655662, 0.649423, 0.65662, 0.66087, 0.650238,
1377 0.653608, 0.655941, 0.650294, 0.650623, 0.652413, 0.64931, 0.649623,
1378 0.653893, 0.652993, 0.650144, 0.652416, 0.649424, 0.652945, 0.649755,
1379 0.655409, 0.654311, 0.65301, 0.648886, 0.651531, 0.651158, 0.655668,
1380 0.657696, 0.654918, 0.653121, 0.649591 };
1381
1382double EvtDsToPipPipPimPi0::a1pgam_arr[3000] = {
1383 0, 0, 0, 0, 0, 0, 0,
1384 0, 0, 0, 0, 0, 0, 0,
1385 0, 0, 0, 0, 0, 0, 0,
1386 0, 0, 0, 0, 0, 0, 0,
1387 0, 0, 0, 0, 0, 0, 0,
1388 0, 0, 0, 0, 0, 0, 0,
1389 0, 0, 0, 0, 0, 0, 0,
1390 0, 0, 0, 0, 0, 0, 0,
1391 0, 0, 0, 0, 0, 0, 0,
1392 0, 0, 0, 0, 0, 0, 0,
1393 0, 0, 0, 0, 0, 0, 0,
1394 0, 0, 0, 0, 0, 0, 0,
1395 0, 0, 0, 0, 0, 0, 0,
1396 0, 0, 0, 0, 0, 0, 0,
1397 0, 0, 0, 0, 0, 0, 0,
1398 0, 0, 0, 0, 0, 0, 0,
1399 0, 0, 0, 0, 0, 0, 0,
1400 0, 0, 0, 0, 0, 0, 0,
1401 0, 0, 0, 0, 0, 0, 0,
1402 0, 0, 0, 0, 0, 0, 0,
1403 0, 0, 0, 0, 0, 0, 0,
1404 0, 0, 0, 0, 0, 0, 0,
1405 0, 0, 0, 0, 0, 0, 0,
1406 0, 0, 0, 0, 0, 0, 0,
1407 0, 0, 0, 0, 0, 0, 0,
1408 0, 0, 0, 0, 0, 0, 0,
1409 0, 0, 0, 0, 0, 0, 0,
1410 0, 0, 0, 0, 0, 0, 0,
1411 0, 0, 0, 0, 0, 0, 0,
1412 0, 0, 0, 0, 0, 0, 0,
1413 0, 0, 0, 0, 0, 0, 0,
1414 0, 0, 0, 0, 0, 0, 0,
1415 0, 0, 0, 0, 0, 0, 0,
1416 0, 0, 0, 0, 0, 0, 0,
1417 0, 0, 0, 0, 0, 0, 0,
1418 0, 0, 0, 0, 0, 0, 0,
1419 0, 0, 0, 0, 0, 0, 0,
1420 0, 0, 0, 0, 0, 0, 0,
1421 0, 0, 0, 0, 0, 0, 0,
1422 0, 0, 0, 0, 0, 0, 0,
1423 0, 0, 0, 0, 0, 0, 0,
1424 0, 0, 0, 0, 0, 0, 0,
1425 0, 0, 0, 0, 0, 0, 0,
1426 0, 0, 0, 0, 0, 0, 0,
1427 0, 0, 0, 0, 0, 0, 0,
1428 0, 0, 0, 0, 0, 0, 0,
1429 0, 0, 0, 0, 0, 0, 0,
1430 0, 0, 0, 0, 0, 0, 0,
1431 0, 0, 0, 0, 0, 0, 0,
1432 0, 0, 0, 0, 0, 0, 0,
1433 0, 0, 0, 0, 0, 0, 0,
1434 0, 0, 0, 0, 0, 0, 0,
1435 0, 0, 0, 0, 0, 0, 0,
1436 0, 0, 0, 0, 0, 0, 0,
1437 0, 0, 0, 0, 0, 0, 0,
1438 0, 0, 0, 0, 0, 0, 0,
1439 0, 0, 0, 0, 0, 0, 0,
1440 0, 0, 0, 0, 0, 0, 0,
1441 0, 0, 0, 0, 0, 0, 0,
1442 0, 0, 0, 0, 0, 0, 0,
1443 0, 0, 0, 0, 0, 0, 0,
1444 0, 0, 0, 0, 1.16255e-06, 1.43224e-06, 1.73208e-06,
1445 2.07745e-06, 2.46811e-06, 2.8958e-06, 3.37455e-06, 3.89844e-06, 4.49816e-06, 5.13849e-06,
1446 5.80413e-06, 6.5552e-06, 7.36488e-06, 8.22848e-06, 9.15607e-06, 1.02051e-05, 1.12368e-05,
1447 1.24356e-05, 1.36177e-05, 1.49035e-05, 1.63023e-05, 1.77732e-05, 1.93243e-05, 2.1039e-05,
1448 2.2603e-05, 2.4414e-05, 2.63433e-05, 2.84536e-05, 3.04303e-05, 3.27502e-05, 3.48853e-05,
1449 3.73187e-05, 3.97525e-05, 4.24245e-05, 4.51341e-05, 4.80726e-05, 5.08268e-05, 5.406e-05,
1450 5.73412e-05, 6.06034e-05, 6.37824e-05, 6.74847e-05, 7.10286e-05, 7.49234e-05, 7.90774e-05,
1451 8.2991e-05, 8.74405e-05, 9.17103e-05, 9.59903e-05, 0.000100626, 0.000105334, 0.00011065,
1452 0.000115699, 0.000120608, 0.000126206, 0.000131872, 0.000137576, 0.000143318, 0.000149702,
1453 0.000155464, 0.000161623, 0.000168204, 0.000175364, 0.000182611, 0.000189817, 0.000196914,
1454 0.000203524, 0.000211917, 0.000219516, 0.000227111, 0.000236204, 0.000243656, 0.000253163,
1455 0.000261804, 0.000270058, 0.000279425, 0.000289415, 0.000299098, 0.000308215, 0.000317935,
1456 0.000328688, 0.000339008, 0.000349371, 0.000360037, 0.00037213, 0.000383524, 0.00039489,
1457 0.00040754, 0.000419597, 0.00043143, 0.000444358, 0.000456249, 0.000469994, 0.00048248,
1458 0.000496576, 0.000509287, 0.000523465, 0.000537663, 0.000552615, 0.000565672, 0.000581505,
1459 0.000598695, 0.000614738, 0.000630295, 0.000646497, 0.000661822, 0.000680691, 0.000697636,
1460 0.000710653, 0.000731611, 0.000747991, 0.000768662, 0.000786658, 0.000804289, 0.000825827,
1461 0.000842601, 0.000863005, 0.000884184, 0.000901145, 0.00092432, 0.000947943, 0.000966836,
1462 0.000990153, 0.00101148, 0.00103307, 0.00105758, 0.00107941, 0.00110262, 0.00113216,
1463 0.00114936, 0.00117555, 0.00119969, 0.00122795, 0.00124984, 0.00127791, 0.00130504,
1464 0.00132556, 0.00135835, 0.0013913, 0.00141671, 0.00144362, 0.00147138, 0.00150494,
1465 0.00153143, 0.00156205, 0.00159195, 0.0016225, 0.00165337, 0.00168743, 0.00172087,
1466 0.00175389, 0.00179009, 0.00181904, 0.00185995, 0.00189109, 0.00192425, 0.00196063,
1467 0.00199956, 0.00203658, 0.00207548, 0.00211007, 0.00214713, 0.00218205, 0.00222067,
1468 0.00226647, 0.00230394, 0.00234372, 0.00239147, 0.00243497, 0.00248179, 0.00252284,
1469 0.00258151, 0.00261317, 0.00265228, 0.0027003, 0.00275796, 0.00279158, 0.00284893,
1470 0.0029003, 0.00294198, 0.00299252, 0.00304708, 0.00310524, 0.00314952, 0.0031982,
1471 0.00326313, 0.00329882, 0.00336018, 0.00341997, 0.00346604, 0.00353923, 0.00359243,
1472 0.00364612, 0.00369917, 0.00376501, 0.00382314, 0.00388662, 0.00395138, 0.00401275,
1473 0.00408231, 0.00415988, 0.00421809, 0.00427482, 0.00436948, 0.004428, 0.00449458,
1474 0.00455717, 0.00462989, 0.0047196, 0.00478372, 0.00487112, 0.00494045, 0.00500795,
1475 0.00509995, 0.00517203, 0.00525552, 0.00533371, 0.00542471, 0.00551365, 0.00558177,
1476 0.00568082, 0.00577522, 0.00587262, 0.00594447, 0.00605402, 0.00611928, 0.00620531,
1477 0.00632442, 0.00642466, 0.00651901, 0.00663929, 0.00672864, 0.00682712, 0.00693815,
1478 0.00703404, 0.00715259, 0.00725306, 0.00736991, 0.00743309, 0.00757263, 0.00769897,
1479 0.00782515, 0.00793906, 0.00804642, 0.00813755, 0.00828158, 0.00842992, 0.00853162,
1480 0.0086354, 0.00875561, 0.00891002, 0.0090398, 0.00921414, 0.0093221, 0.00942081,
1481 0.00958837, 0.0097279, 0.0098762, 0.0100076, 0.0101537, 0.0102833, 0.0104462,
1482 0.0106011, 0.0107291, 0.0109382, 0.011029, 0.0112307, 0.0113983, 0.0115593,
1483 0.011673, 0.0119066, 0.0120356, 0.0122478, 0.0124057, 0.0125293, 0.0126917,
1484 0.0129166, 0.0130909, 0.0132769, 0.0134728, 0.0136977, 0.0138256, 0.0140528,
1485 0.0142311, 0.0144146, 0.0146287, 0.0148169, 0.0150073, 0.0151786, 0.0154452,
1486 0.0156354, 0.0158328, 0.016026, 0.0162738, 0.0164946, 0.01674, 0.0169668,
1487 0.0171655, 0.017386, 0.0176114, 0.0178374, 0.0180972, 0.0182908, 0.0185786,
1488 0.0188305, 0.019112, 0.0192812, 0.0195456, 0.0198758, 0.0200246, 0.0203729,
1489 0.0205916, 0.0209285, 0.0211636, 0.0213849, 0.0217342, 0.0219524, 0.0222627,
1490 0.0225478, 0.0228542, 0.0231567, 0.0233594, 0.0237955, 0.0240128, 0.024268,
1491 0.0246971, 0.0250577, 0.0252016, 0.0256289, 0.0258298, 0.0263225, 0.0265181,
1492 0.0269615, 0.0272329, 0.0276034, 0.0278982, 0.02828, 0.0285273, 0.0288389,
1493 0.0292445, 0.0295984, 0.0299843, 0.0303156, 0.030751, 0.0311039, 0.0315635,
1494 0.0318173, 0.032281, 0.0327688, 0.0331304, 0.033465, 0.0338762, 0.0342434,
1495 0.0348077, 0.0351089, 0.0356362, 0.0361068, 0.0363789, 0.0369809, 0.037303,
1496 0.0377374, 0.0383674, 0.0386694, 0.0391393, 0.0396053, 0.040115, 0.040666,
1497 0.0409595, 0.0415084, 0.0419804, 0.0426746, 0.0430245, 0.043641, 0.044243,
1498 0.044666, 0.0454314, 0.0456115, 0.0461765, 0.0468266, 0.0472917, 0.048017,
1499 0.0485655, 0.0490725, 0.0495151, 0.0501254, 0.0508391, 0.0513236, 0.0520693,
1500 0.0526731, 0.0534017, 0.0540798, 0.0544842, 0.0552119, 0.0557412, 0.0563874,
1501 0.0572445, 0.0579557, 0.0587949, 0.0593476, 0.0598343, 0.0607187, 0.061472,
1502 0.0620747, 0.0628726, 0.0636456, 0.0642243, 0.0650331, 0.0657648, 0.0667879,
1503 0.067215, 0.068234, 0.0689635, 0.069516, 0.0703112, 0.0715982, 0.0723831,
1504 0.0730516, 0.0741697, 0.0749361, 0.0757809, 0.0767279, 0.0773934, 0.0787365,
1505 0.0796201, 0.0805431, 0.0810382, 0.0821846, 0.0832118, 0.0839816, 0.0852385,
1506 0.0859977, 0.0872543, 0.0882491, 0.0893764, 0.090282, 0.0914214, 0.0922846,
1507 0.0936055, 0.0945924, 0.0957427, 0.0966392, 0.097442, 0.0990003, 0.100014,
1508 0.101187, 0.10229, 0.103675, 0.105012, 0.105897, 0.107099, 0.108671,
1509 0.109855, 0.110575, 0.111588, 0.113503, 0.114285, 0.116308, 0.11702,
1510 0.11869, 0.119724, 0.120891, 0.123123, 0.123607, 0.125191, 0.127004,
1511 0.128134, 0.129361, 0.130804, 0.13236, 0.133601, 0.13523, 0.136752,
1512 0.138099, 0.138969, 0.141041, 0.142319, 0.14346, 0.144497, 0.146986,
1513 0.147875, 0.149762, 0.151885, 0.153448, 0.154588, 0.155649, 0.157296,
1514 0.158514, 0.1603, 0.161745, 0.16335, 0.164218, 0.167009, 0.168287,
1515 0.16982, 0.170987, 0.172676, 0.174711, 0.176843, 0.177706, 0.178822,
1516 0.180431, 0.182226, 0.183719, 0.184878, 0.186864, 0.188443, 0.191067,
1517 0.190953, 0.193082, 0.195138, 0.197183, 0.198269, 0.198657, 0.200299,
1518 0.202385, 0.203421, 0.205591, 0.207379, 0.208769, 0.209583, 0.211311,
1519 0.212606, 0.213773, 0.215906, 0.216704, 0.219097, 0.219965, 0.222538,
1520 0.22322, 0.225034, 0.225126, 0.227784, 0.229013, 0.230955, 0.231951,
1521 0.232997, 0.234823, 0.237259, 0.237356, 0.239527, 0.239873, 0.241319,
1522 0.24279, 0.24475, 0.245765, 0.247502, 0.247923, 0.249917, 0.251959,
1523 0.252184, 0.253304, 0.254642, 0.256658, 0.257593, 0.259196, 0.26103,
1524 0.261744, 0.262631, 0.263768, 0.264823, 0.266831, 0.268262, 0.269319,
1525 0.270256, 0.271079, 0.27332, 0.274539, 0.27433, 0.276119, 0.278634,
1526 0.279295, 0.279376, 0.281229, 0.280986, 0.283135, 0.284504, 0.284874,
1527 0.287492, 0.287401, 0.289914, 0.290073, 0.291362, 0.29294, 0.29308,
1528 0.293553, 0.294678, 0.296669, 0.296969, 0.297775, 0.298636, 0.300789,
1529 0.300916, 0.303054, 0.303864, 0.305603, 0.305079, 0.305364, 0.306847,
1530 0.309441, 0.309088, 0.31116, 0.313184, 0.312351, 0.313737, 0.315687,
1531 0.315318, 0.316242, 0.317757, 0.31832, 0.318201, 0.319554, 0.320554,
1532 0.321531, 0.321619, 0.321803, 0.324131, 0.325346, 0.326026, 0.325801,
1533 0.32723, 0.329022, 0.330649, 0.331494, 0.332213, 0.332949, 0.332773,
1534 0.332706, 0.333732, 0.334117, 0.336933, 0.33752, 0.337368, 0.33903,
1535 0.336397, 0.339849, 0.339573, 0.342183, 0.341996, 0.342705, 0.342589,
1536 0.344926, 0.345523, 0.344874, 0.345849, 0.347146, 0.348979, 0.347889,
1537 0.350886, 0.351242, 0.35026, 0.350276, 0.352683, 0.353797, 0.353656,
1538 0.354533, 0.354091, 0.355784, 0.35688, 0.35727, 0.358318, 0.36001,
1539 0.358788, 0.362043, 0.360437, 0.360933, 0.360429, 0.362404, 0.363339,
1540 0.363307, 0.362733, 0.364112, 0.365205, 0.366048, 0.367317, 0.367111,
1541 0.367419, 0.367187, 0.369484, 0.368694, 0.37072, 0.369628, 0.371577,
1542 0.371741, 0.371562, 0.372591, 0.373813, 0.374935, 0.375387, 0.374697,
1543 0.375046, 0.376664, 0.378074, 0.377495, 0.376918, 0.378356, 0.376853,
1544 0.379804, 0.380568, 0.380397, 0.38164, 0.380695, 0.383465, 0.379018,
1545 0.382784, 0.386237, 0.383919, 0.382384, 0.385775, 0.386094, 0.385067,
1546 0.386536, 0.388298, 0.387304, 0.387484, 0.390188, 0.388102, 0.388931,
1547 0.391046, 0.390972, 0.392173, 0.392513, 0.393543, 0.390081, 0.390173,
1548 0.389479, 0.393965, 0.394575, 0.394153, 0.395897, 0.395575, 0.395761,
1549 0.396228, 0.39594, 0.398025, 0.397283, 0.399532, 0.398754, 0.39678,
1550 0.396219, 0.400911, 0.397152, 0.402222, 0.398215, 0.399644, 0.400621,
1551 0.403141, 0.402423, 0.402603, 0.4024, 0.402452, 0.400182, 0.405272,
1552 0.40486, 0.405836, 0.40311, 0.406523, 0.4054, 0.408891, 0.407743,
1553 0.405528, 0.408623, 0.407901, 0.409105, 0.407848, 0.408936, 0.40985,
1554 0.410865, 0.410143, 0.406514, 0.412128, 0.410341, 0.413463, 0.413164,
1555 0.410219, 0.412328, 0.411408, 0.411539, 0.409508, 0.4124, 0.41493,
1556 0.416484, 0.413764, 0.414824, 0.416969, 0.41434, 0.418504, 0.417209,
1557 0.416756, 0.414673, 0.419981, 0.419941, 0.417795, 0.418548, 0.420473,
1558 0.418942, 0.421251, 0.420259, 0.420498, 0.420359, 0.42, 0.418842,
1559 0.421591, 0.418179, 0.42091, 0.420125, 0.423146, 0.422455, 0.421435,
1560 0.420877, 0.422758, 0.422608, 0.422679, 0.424559, 0.424062, 0.425715,
1561 0.424857, 0.424602, 0.428664, 0.427413, 0.423959, 0.430009, 0.425758,
1562 0.425193, 0.425283, 0.427155, 0.42827, 0.427879, 0.428337, 0.426313,
1563 0.425814, 0.428359, 0.430518, 0.430796, 0.429357, 0.427775, 0.42792,
1564 0.431126, 0.430628, 0.429625, 0.434041, 0.432825, 0.431701, 0.431599,
1565 0.428025, 0.432971, 0.432064, 0.433326, 0.4315, 0.430761, 0.432701,
1566 0.432203, 0.434129, 0.433326, 0.433766, 0.435136, 0.434353, 0.434856,
1567 0.435385, 0.434678, 0.435788, 0.435205, 0.434885, 0.437207, 0.438295,
1568 0.434515, 0.438352, 0.437057, 0.4365, 0.437515, 0.440111, 0.438409,
1569 0.436411, 0.440059, 0.435848, 0.437692, 0.435946, 0.438257, 0.439022,
1570 0.441926, 0.439565, 0.439481, 0.440565, 0.442036, 0.440238, 0.440979,
1571 0.440926, 0.439638, 0.444596, 0.440835, 0.439819, 0.441838, 0.442397,
1572 0.442598, 0.442112, 0.441093, 0.444941, 0.440799, 0.44294, 0.443923,
1573 0.446154, 0.444274, 0.442494, 0.441571, 0.443644, 0.446548, 0.446601,
1574 0.44691, 0.446974, 0.444332, 0.445558, 0.442494, 0.445232, 0.447516,
1575 0.446635, 0.444198, 0.446976, 0.446334, 0.448763, 0.44505, 0.446623,
1576 0.447772, 0.44823, 0.447627, 0.447718, 0.447844, 0.450088, 0.447984,
1577 0.447244, 0.450647, 0.448979, 0.447463, 0.448813, 0.450568, 0.451572,
1578 0.451958, 0.449579, 0.451486, 0.452818, 0.450341, 0.451603, 0.451527,
1579 0.450674, 0.454367, 0.451096, 0.454005, 0.452618, 0.454843, 0.454316,
1580 0.455046, 0.454825, 0.454406, 0.454624, 0.452211, 0.456073, 0.456859,
1581 0.456543, 0.456608, 0.458317, 0.458728, 0.456571, 0.456989, 0.461559,
1582 0.458669, 0.459518, 0.463071, 0.45999, 0.459771, 0.458386, 0.460108,
1583 0.463027, 0.460509, 0.462204, 0.459846, 0.464034, 0.460146, 0.461821,
1584 0.462694, 0.464422, 0.464904, 0.46283, 0.464172, 0.463687, 0.465088,
1585 0.466588, 0.464877, 0.464343, 0.463757, 0.467875, 0.464723, 0.469002,
1586 0.466809, 0.467942, 0.468663, 0.469546, 0.469821, 0.467963, 0.467366,
1587 0.470452, 0.469731, 0.470751, 0.471091, 0.470409, 0.474515, 0.47055,
1588 0.474283, 0.474409, 0.472911, 0.472518, 0.474577, 0.474959, 0.474488,
1589 0.474743, 0.475597, 0.475335, 0.473161, 0.47643, 0.474537, 0.477407,
1590 0.476401, 0.475578, 0.477679, 0.477219, 0.477499, 0.477062, 0.479224,
1591 0.478395, 0.478, 0.479369, 0.478258, 0.478699, 0.481036, 0.479984,
1592 0.48036, 0.481885, 0.480911, 0.481145, 0.483952, 0.481125, 0.480148,
1593 0.482274, 0.481469, 0.482899, 0.480925, 0.481978, 0.484312, 0.480981,
1594 0.482839, 0.482758, 0.484997, 0.483989, 0.482316, 0.484268, 0.483607,
1595 0.485649, 0.485508, 0.485669, 0.485859, 0.487146, 0.48618, 0.487155,
1596 0.485832, 0.487672, 0.486863, 0.486041, 0.488219, 0.488653, 0.488266,
1597 0.489176, 0.486769, 0.488176, 0.488342, 0.490516, 0.489628, 0.488147,
1598 0.489412, 0.489563, 0.488379, 0.491313, 0.491121, 0.487614, 0.492559,
1599 0.489252, 0.493605, 0.493249, 0.490203, 0.492976, 0.489073, 0.490983,
1600 0.494752, 0.493049, 0.49281, 0.491329, 0.492186, 0.49269, 0.491985,
1601 0.493295, 0.493805, 0.494683, 0.4951, 0.49412, 0.49653, 0.492813,
1602 0.496693, 0.495433, 0.494249, 0.495617, 0.493161, 0.497948, 0.495179,
1603 0.497073, 0.495682, 0.495163, 0.49698, 0.495499, 0.4964, 0.499198,
1604 0.498682, 0.497512, 0.498361, 0.498176, 0.500341, 0.498605, 0.497527,
1605 0.499698, 0.499794, 0.496928, 0.498103, 0.499679, 0.497799, 0.500336,
1606 0.498096, 0.496839, 0.500489, 0.499934, 0.500009, 0.50085, 0.499986,
1607 0.502149, 0.501043, 0.50334, 0.501172, 0.50259, 0.503647, 0.500887,
1608 0.502283, 0.502236, 0.502643, 0.500347, 0.503511, 0.504083, 0.502938,
1609 0.503588, 0.502952, 0.50014, 0.503376, 0.504829, 0.505831, 0.505635,
1610 0.503967, 0.504616, 0.503945, 0.50351, 0.504207, 0.50789, 0.502813,
1611 0.507481, 0.504709, 0.505123, 0.50434, 0.503771, 0.505899, 0.50618,
1612 0.505126, 0.505793, 0.509522, 0.50452, 0.507515, 0.507257, 0.507825,
1613 0.507041, 0.506245, 0.506905, 0.509064, 0.507124, 0.506719, 0.506902,
1614 0.509027, 0.508743, 0.506517, 0.510001, 0.511781, 0.508727, 0.510317,
1615 0.506898, 0.509835, 0.510286, 0.510439, 0.509651, 0.509386, 0.51086,
1616 0.513609, 0.508515, 0.511194, 0.511241, 0.511412, 0.508399, 0.510259,
1617 0.511576, 0.511302, 0.512029, 0.509875, 0.510727, 0.512554, 0.51035,
1618 0.512636, 0.512859, 0.515007, 0.512985, 0.511215, 0.512157, 0.512022,
1619 0.514922, 0.512503, 0.511743, 0.511697, 0.513383, 0.512563, 0.514872,
1620 0.513686, 0.51463, 0.513491, 0.51491, 0.515901, 0.512171, 0.516725,
1621 0.515435, 0.512999, 0.515207, 0.515277, 0.514774, 0.516273, 0.514556,
1622 0.516018, 0.515397, 0.515952, 0.514911, 0.515583, 0.514949, 0.514043,
1623 0.515837, 0.517428, 0.515797, 0.517289, 0.515266, 0.515199, 0.515763,
1624 0.515808, 0.516554, 0.519186, 0.518709, 0.517786, 0.516753, 0.517333,
1625 0.516965, 0.515283, 0.519513, 0.516244, 0.519102, 0.516382, 0.517016,
1626 0.518363, 0.519862, 0.518993, 0.516998, 0.517437, 0.519582, 0.521192,
1627 0.516754, 0.520025, 0.519587, 0.516353, 0.520777, 0.518041, 0.521555,
1628 0.520629, 0.521225, 0.520611, 0.522332, 0.519664, 0.52123, 0.520107,
1629 0.520044, 0.521267, 0.519781, 0.520637, 0.522087, 0.521725, 0.519613,
1630 0.521127, 0.524505, 0.522686, 0.521674, 0.523535, 0.523107, 0.523355,
1631 0.523493, 0.519151, 0.521462, 0.519819, 0.521864, 0.522879, 0.522073,
1632 0.523275, 0.519581, 0.52384, 0.521515, 0.523733, 0.526311, 0.526098,
1633 0.52433, 0.522178, 0.523064, 0.522078, 0.524207, 0.52156, 0.524184,
1634 0.523995, 0.525642, 0.523386, 0.522591, 0.524086, 0.52377, 0.521208,
1635 0.527144, 0.521631, 0.52514, 0.524118, 0.528724, 0.524244, 0.525148,
1636 0.524989, 0.527725, 0.525912, 0.52612, 0.525197, 0.52598, 0.524787,
1637 0.526908, 0.525031, 0.529178, 0.52444, 0.528708, 0.527348, 0.528284,
1638 0.527454, 0.525352, 0.524812, 0.523676, 0.525699, 0.525356, 0.527578,
1639 0.527662, 0.52777, 0.527981, 0.526332, 0.528916, 0.529222, 0.528545,
1640 0.528026, 0.529138, 0.526978, 0.527438, 0.52948, 0.527731, 0.529236,
1641 0.527679, 0.527083, 0.531072, 0.532148, 0.527576, 0.52708, 0.52899,
1642 0.530282, 0.530114, 0.530201, 0.529259, 0.530194, 0.52944, 0.526997,
1643 0.527695, 0.533865, 0.529303, 0.531477, 0.52986, 0.529447, 0.52961,
1644 0.532221, 0.529561, 0.528776, 0.531029, 0.530342, 0.530347, 0.529125,
1645 0.532999, 0.531918, 0.52894, 0.532453, 0.530963, 0.5303, 0.530267,
1646 0.528643, 0.52949, 0.530468, 0.530107, 0.531941, 0.534161, 0.531372,
1647 0.532518, 0.532483, 0.533227, 0.533226, 0.533284, 0.531491, 0.532766,
1648 0.530457, 0.532497, 0.531895, 0.535718, 0.531715, 0.532508, 0.533666,
1649 0.53416, 0.532363, 0.531327, 0.531607, 0.531271, 0.530863, 0.534436,
1650 0.531824, 0.532806, 0.532995, 0.533911, 0.533086, 0.536026, 0.534301,
1651 0.534823, 0.535135, 0.531684, 0.533152, 0.532981, 0.533428, 0.535088,
1652 0.535174, 0.532605, 0.531998, 0.534175, 0.535562, 0.531693, 0.536124,
1653 0.532424, 0.534834, 0.532604, 0.534527, 0.534851, 0.535341, 0.534146,
1654 0.535489, 0.533177, 0.532962, 0.536297, 0.537311, 0.538172, 0.538443,
1655 0.534239, 0.538912, 0.535075, 0.535317, 0.53505, 0.537858, 0.537639,
1656 0.536153, 0.537624, 0.534491, 0.535245, 0.53452, 0.538806, 0.536585,
1657 0.535028, 0.537618, 0.536423, 0.538687, 0.538609, 0.536411, 0.53854,
1658 0.538888, 0.540546, 0.536277, 0.537438, 0.536971, 0.539004, 0.537137,
1659 0.538022, 0.538497, 0.533824, 0.536522, 0.537436, 0.54117, 0.537337,
1660 0.536147, 0.539299, 0.540078, 0.535115, 0.535998, 0.538684, 0.539206,
1661 0.539454, 0.541245, 0.543455, 0.534915, 0.538849, 0.539623, 0.538584,
1662 0.536525, 0.537951, 0.540761, 0.540869, 0.541007, 0.539463, 0.541607,
1663 0.540128, 0.540722, 0.53975, 0.538474, 0.541993, 0.53982, 0.540508,
1664 0.54152, 0.540565, 0.541718, 0.543339, 0.539217, 0.5412, 0.541414,
1665 0.542482, 0.542555, 0.542531, 0.540586, 0.539943, 0.540798, 0.540654,
1666 0.541274, 0.540014, 0.541943, 0.540362, 0.543648, 0.54109, 0.539127,
1667 0.539, 0.542587, 0.5405, 0.541717, 0.541137, 0.542931, 0.541256,
1668 0.540733, 0.543748, 0.541273, 0.541947, 0.542462, 0.544027, 0.541896,
1669 0.538944, 0.541228, 0.544038, 0.542104, 0.545458, 0.544605, 0.540212,
1670 0.543099, 0.544071, 0.543506, 0.542657, 0.54369, 0.544288, 0.547427,
1671 0.54318, 0.539869, 0.545902, 0.54276, 0.544856, 0.540374, 0.543457,
1672 0.542108, 0.546985, 0.540793, 0.545665, 0.542855, 0.544058, 0.546318,
1673 0.541934, 0.542663, 0.545231, 0.545161, 0.541598, 0.542779, 0.543125,
1674 0.539964, 0.544597, 0.545379, 0.543101, 0.546417, 0.547671, 0.545231,
1675 0.547041, 0.54409, 0.543794, 0.545575, 0.544796, 0.545207, 0.546179,
1676 0.544222, 0.545494, 0.544425, 0.549579, 0.547971, 0.546074, 0.546866,
1677 0.547873, 0.546682, 0.546789, 0.546476, 0.546983, 0.548962, 0.549911,
1678 0.546108, 0.545207, 0.549275, 0.545967, 0.546158, 0.546814, 0.54504,
1679 0.546335, 0.546309, 0.544027, 0.546949, 0.547314, 0.547413, 0.549132,
1680 0.548506, 0.546994, 0.54753, 0.546117, 0.54826, 0.549958, 0.548654,
1681 0.544949, 0.547196, 0.550506, 0.549048, 0.551489, 0.553668, 0.550428,
1682 0.548964, 0.547776, 0.546852, 0.547314, 0.5492, 0.546406, 0.546505,
1683 0.551475, 0.549953, 0.548113, 0.549375, 0.546074, 0.549023, 0.546288,
1684 0.551046, 0.549058, 0.548697, 0.547909, 0.548509, 0.552286, 0.548383,
1685 0.552206, 0.551249, 0.549937, 0.552169, 0.547418, 0.549525, 0.549365,
1686 0.548287, 0.548609, 0.548641, 0.545756, 0.554655, 0.551645, 0.551535,
1687 0.552328, 0.551244, 0.550255, 0.552096, 0.550083, 0.549222, 0.553287,
1688 0.549592, 0.552809, 0.552951, 0.549147, 0.549725, 0.550262, 0.552019,
1689 0.549318, 0.553297, 0.550599, 0.5534, 0.551108, 0.553334, 0.549486,
1690 0.551931, 0.552586, 0.554142, 0.550783, 0.552973, 0.549018, 0.555185,
1691 0.552498, 0.550528, 0.554374, 0.551396, 0.55128, 0.556513, 0.553616,
1692 0.552864, 0.552209, 0.552188, 0.554872, 0.549812, 0.552671, 0.556074,
1693 0.554393, 0.556044, 0.552497, 0.552524, 0.551355, 0.550879, 0.552187,
1694 0.553634, 0.555197, 0.55219, 0.553756, 0.551232, 0.553567, 0.555009,
1695 0.552611, 0.552221, 0.551817, 0.554271, 0.551477, 0.555141, 0.553705,
1696 0.553635, 0.553739, 0.555041, 0.553131, 0.557623, 0.556932, 0.558183,
1697 0.559317, 0.553292, 0.550044, 0.552954, 0.555819, 0.554345, 0.559802,
1698 0.553091, 0.5507, 0.558188, 0.554716, 0.550222, 0.55223, 0.557408,
1699 0.553816, 0.557835, 0.556988, 0.554317, 0.554479, 0.554807, 0.555596,
1700 0.555005, 0.556904, 0.558648, 0.554836, 0.556575, 0.556295, 0.558238,
1701 0.556626, 0.554118, 0.554228, 0.55693, 0.556709, 0.552743, 0.556664,
1702 0.553648, 0.561244, 0.556389, 0.553583, 0.557999, 0.558632, 0.560102,
1703 0.55774, 0.553642, 0.555382, 0.556047, 0.557315, 0.557459, 0.559487,
1704 0.55687, 0.558063, 0.557711, 0.557409, 0.560746, 0.558234, 0.556146,
1705 0.558541, 0.558111, 0.558044, 0.556974, 0.556342, 0.55813, 0.560921,
1706 0.561373, 0.55843, 0.558542, 0.558594, 0.562144, 0.557889, 0.558017,
1707 0.554889, 0.559593, 0.557181, 0.556927, 0.559527, 0.559096, 0.557727,
1708 0.55767, 0.558729, 0.559692, 0.556463, 0.561864, 0.558947, 0.55897,
1709 0.559805, 0.563493, 0.563203, 0.560529, 0.559445, 0.560133, 0.559202,
1710 0.559178, 0.556942, 0.558551, 0.561764, 0.560133, 0.557944, 0.561129,
1711 0.557547, 0.558825, 0.557005, 0.561255, 0.563616, 0.564167, 0.559045,
1712 0.558354, 0.561385, 0.55955, 0.561884, 0.559428, 0.560908, 0.561489,
1713 0.561759, 0.562946, 0.556628, 0.559323, 0.561308, 0.563207, 0.559584,
1714 0.56363, 0.561338, 0.56417, 0.561375, 0.561335, 0.562197, 0.561145,
1715 0.56871, 0.56479, 0.56347, 0.562562, 0.56457, 0.563354, 0.558787,
1716 0.559027, 0.562868, 0.563334, 0.561828, 0.561801, 0.565675, 0.564157,
1717 0.568886, 0.564563, 0.564034, 0.560687, 0.561346, 0.567538, 0.565062,
1718 0.559714, 0.562784, 0.561863, 0.561518, 0.566366, 0.561876, 0.562857,
1719 0.56424, 0.56227, 0.56566, 0.566914, 0.566943, 0.567153, 0.56502,
1720 0.564887, 0.564136, 0.565241, 0.564616, 0.562678, 0.569477, 0.561688,
1721 0.563319, 0.562613, 0.563915, 0.559818, 0.563493, 0.568038, 0.565651,
1722 0.567155, 0.563786, 0.563731, 0.56761, 0.566182, 0.56808, 0.564559,
1723 0.564938, 0.564552, 0.567519, 0.565601, 0.565827, 0.56987, 0.563697,
1724 0.56704, 0.565665, 0.563851, 0.569622, 0.563456, 0.563018, 0.564004,
1725 0.56549, 0.568403, 0.568227, 0.562631, 0.566965, 0.564924, 0.565541,
1726 0.565865, 0.565741, 0.56577, 0.566443, 0.571398, 0.568946, 0.568286,
1727 0.566545, 0.56912, 0.565507, 0.568109, 0.565578, 0.571635, 0.566771,
1728 0.566112, 0.568029, 0.566198, 0.566501, 0.56401, 0.563124, 0.571897,
1729 0.565815, 0.56805, 0.564366, 0.570441, 0.563243, 0.570155, 0.566313,
1730 0.568843, 0.569122, 0.570106, 0.56893, 0.567451, 0.569489, 0.573791,
1731 0.56629, 0.563365, 0.567514, 0.571623, 0.566948, 0.572632, 0.568619,
1732 0.566112, 0.567199, 0.573534, 0.571681, 0.566823, 0.565965, 0.569542,
1733 0.567893, 0.574422, 0.571514, 0.567935, 0.570829, 0.571001, 0.571396,
1734 0.568273, 0.568196, 0.569275, 0.57053, 0.568056, 0.571515, 0.571651,
1735 0.570821, 0.573593, 0.570831, 0.571778, 0.567559, 0.570017, 0.57121,
1736 0.571879, 0.574164, 0.569904, 0.567905, 0.572251, 0.576637, 0.568732,
1737 0.572377, 0.568813, 0.569962, 0.566516, 0.570752, 0.574221, 0.572185,
1738 0.571619, 0.57182, 0.573078, 0.574478, 0.571787, 0.57357, 0.568609,
1739 0.577304, 0.571892, 0.572633, 0.575569, 0.56975, 0.574168, 0.574377,
1740 0.571945, 0.571511, 0.572294, 0.568596, 0.57615, 0.574653, 0.577153,
1741 0.571409, 0.578643, 0.57639, 0.571439, 0.575279, 0.572562, 0.572535,
1742 0.577689, 0.57259, 0.574433, 0.571285, 0.575849, 0.576038, 0.57019,
1743 0.573275, 0.576725, 0.575433, 0.575766, 0.571802, 0.57247, 0.571827,
1744 0.573072, 0.572525, 0.577401, 0.574592, 0.570752, 0.575398, 0.577319,
1745 0.575112, 0.567575, 0.576563, 0.576735, 0.572341, 0.579265, 0.57467,
1746 0.574273, 0.576863, 0.57808, 0.574279, 0.576533, 0.57649, 0.576091,
1747 0.57389, 0.57258, 0.575254, 0.576797, 0.578916, 0.575775, 0.572161,
1748 0.575083, 0.574561, 0.580308, 0.575883, 0.574889, 0.576224, 0.57843,
1749 0.587399, 0.579503, 0.578239, 0.571089, 0.575414, 0.579622, 0.580299,
1750 0.575199, 0.577805, 0.577529, 0.578705, 0.580518, 0.577221, 0.580147,
1751 0.577344, 0.579751, 0.578509, 0.58019, 0.579631, 0.575439, 0.57825,
1752 0.582402, 0.578612, 0.578125, 0.575566, 0.578315, 0.575552, 0.575943,
1753 0.576915, 0.580044, 0.579177, 0.577357, 0.573965, 0.582317, 0.577726,
1754 0.577649, 0.583044, 0.576295, 0.576201, 0.580481, 0.581689, 0.575013,
1755 0.579414, 0.574689, 0.577331, 0.583958, 0.578259, 0.580574, 0.582636,
1756 0.577774, 0.574513, 0.583397, 0.584202, 0.57918, 0.577046, 0.578281,
1757 0.578361, 0.57914, 0.580066, 0.582006, 0.578385, 0.580817, 0.581184,
1758 0.578551, 0.581335, 0.582952, 0.579492, 0.581003, 0.582552, 0.579876,
1759 0.579589, 0.580873, 0.581491, 0.580155, 0.579513, 0.582733, 0.579968,
1760 0.580018, 0.579553, 0.581104, 0.57785, 0.583129, 0.583348, 0.580091,
1761 0.58456, 0.584324, 0.577438, 0.590731, 0.583957, 0.581641, 0.578626,
1762 0.580458, 0.579534, 0.578925, 0.583496, 0.583054, 0.581946, 0.585855,
1763 0.582747, 0.578985, 0.583717, 0.586242, 0.585105, 0.585066, 0.584821,
1764 0.582321, 0.584162, 0.586393, 0.582921, 0.583427, 0.587302, 0.584023,
1765 0.57957, 0.583823, 0.585567, 0.587179, 0.580456, 0.583967, 0.585735,
1766 0.579745, 0.581204, 0.580524, 0.583377, 0.58247, 0.582338, 0.578497,
1767 0.586158, 0.58618, 0.587141, 0.589735, 0.578221, 0.585685, 0.584407,
1768 0.583547, 0.583947, 0.590027, 0.588673, 0.586552, 0.588762, 0.590184,
1769 0.587378, 0.587835, 0.586443, 0.583279, 0.59183, 0.582998, 0.587141,
1770 0.586623, 0.584549, 0.588225, 0.581107, 0.592989, 0.583229, 0.582401,
1771 0.585253, 0.585102, 0.585882, 0.586341, 0.584575, 0.587689, 0.586155,
1772 0.583827, 0.590553, 0.590345, 0.590344, 0.586981, 0.585238, 0.587449,
1773 0.587378, 0.582548, 0.590458, 0.59106, 0.584135, 0.59128, 0.586538,
1774 0.588367, 0.585908, 0.584932, 0.582699, 0.589895, 0.594493, 0.585279,
1775 0.587439, 0.58919, 0.581968, 0.591246, 0.584493, 0.590344, 0.587523,
1776 0.590098, 0.583483, 0.586165, 0.593377, 0.592627, 0.587669, 0.592341,
1777 0.584057, 0.588531, 0.58817, 0.587276, 0.591735, 0.590744, 0.590016,
1778 0.58959, 0.585778, 0.589998, 0.589391, 0.590903, 0.589063, 0.585712,
1779 0.59058, 0.586772, 0.590911, 0.593076, 0.590414, 0.588227, 0.590888,
1780 0.591537, 0.592885, 0.593597, 0.597208, 0.589186, 0.590275, 0.592207,
1781 0.595607, 0.59283, 0.594183, 0.589069, 0.592242, 0.594122, 0.589475,
1782 0.590137, 0.589487, 0.591448, 0.591108, 0.589797, 0.594664, 0.593426,
1783 0.591803, 0.591522, 0.592354, 0.588561, 0.591159, 0.591905, 0.586348,
1784 0.589954, 0.591269, 0.597903, 0.593487, 0.59121, 0.593892, 0.594689,
1785 0.600124, 0.596178, 0.593128, 0.591731, 0.589207, 0.599487, 0.592424,
1786 0.593232, 0.597144, 0.594733, 0.598763, 0.596094, 0.593358, 0.593704,
1787 0.596968, 0.592885, 0.590531, 0.593849, 0.595086, 0.589113, 0.59521,
1788 0.600783, 0.595183, 0.594104, 0.592028, 0.598734, 0.595671, 0.590967,
1789 0.592103, 0.594867, 0.593684, 0.598303, 0.598886, 0.594022, 0.599185,
1790 0.594191, 0.598571, 0.59825, 0.596676, 0.592284, 0.597617, 0.590493,
1791 0.592461, 0.594437, 0.598739, 0.596292, 0.592047, 0.595335, 0.593602,
1792 0.596831, 0.597131, 0.595548, 0.601269, 0.599161, 0.598828, 0.592952,
1793 0.595053, 0.597111, 0.592736, 0.592976, 0.59527, 0.598559, 0.593337,
1794 0.600347, 0.599399, 0.598607, 0.59987, 0.597814, 0.601842, 0.598387,
1795 0.603452, 0.600824, 0.595827, 0.598327, 0.595488, 0.60024, 0.598835,
1796 0.597184, 0.598103, 0.606117, 0.602632, 0.602415, 0.601875, 0.60091,
1797 0.602629, 0.60604, 0.598444, 0.59611, 0.597286, 0.596517, 0.601671,
1798 0.597373, 0.596818, 0.598115, 0.597333, 0.604342, 0.600689, 0.59807,
1799 0.596776, 0.602913, 0.600332, 0.597332, 0.596352, 0.593713, 0.595789,
1800 0.60253, 0.607278, 0.606196, 0.605733, 0.600273, 0.597554, 0.607263,
1801 0.602899, 0.600015, 0.598794, 0.601687, 0.601923, 0.606928, 0.601463,
1802 0.602877, 0.596957, 0.600312, 0.601913, 0.60606, 0.603118, 0.606835,
1803 0.597332, 0.599756, 0.601804, 0.606572, 0.607449, 0.602602, 0.605645,
1804 0.605075, 0.602008, 0.606686, 0.601018, 0.597173, 0.601425, 0.599703,
1805 0.599699, 0.60309, 0.594846, 0.601009, 0.604537, 0.606915, 0.607101,
1806 0.605826, 0.604978, 0.604253, 0.607059, 0.606595, 0.606202, 0.600498,
1807 0.60769, 0.608191, 0.60771, 0.60535, 0.603211, 0.601866, 0.602063,
1808 0.60757, 0.601387, 0.603771, 0.609048, 0.607257, 0.605849, 0.607487,
1809 0.600519, 0.601944, 0.605909, 0.601311, 0.603116, 0.607694, 0.607612,
1810 0.601607, 0.608192, 0.604752, 0.606097, 0.604806, 0.609834, 0.608187,
1811 0.606182, 0.604892, 0.604471, 0.608789 };
report
ostream & report(Severity severity, const char *facility)
Definition EvtReport.cc:34
ERROR
@ ERROR
Definition EvtReport.hh:49
sin
double sin(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:185
cos
double cos(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:187
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition EvtDecayBase.cc:261
EvtDecayBase::getParentId
EvtId getParentId()
Definition EvtDecayBase.hh:63
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition EvtDecayBase.cc:465
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:68
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:449
EvtDecayBase::getDaug
EvtId getDaug(int i)
Definition EvtDecayBase.hh:69
EvtDecayProb::setProb
void setProb(double prob)
Definition EvtDecayProb.hh:40
EvtDsToPipPipPimPi0::getName
void getName(std::string &name)
Definition EvtDsToPipPipPimPi0.cc:31
EvtDsToPipPipPimPi0::decay
void decay(EvtParticle *p)
Definition EvtDsToPipPipPimPi0.cc:161
EvtPDL::getId
static EvtId getId(const std::string &name)
Definition EvtPDL.cc:272
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:108
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:80
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1042
EvtVector4R::get
double get(int i) const
Definition EvtVector4R.hh:171