BOSS 7.1.2
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 <stdlib.h>
20#include <math.h>
21#include "EvtGenBase/EvtVector4R.hh"
22#include "EvtGenBase/EvtParticle.hh"
23#include "EvtGenBase/EvtGenKine.hh"
24#include "EvtGenBase/EvtPDL.hh"
25#include "EvtGenModels/EvtDsToPipPipPimPi0.hh"
26#include "EvtGenBase/EvtReport.hh"
27#include <string>
28
29EvtDsToPipPipPimPi0::~EvtDsToPipPipPimPi0() {}
30
31void EvtDsToPipPipPimPi0::getName(std::string& model_name){
32
33 model_name="DsToPipPipPimPi0";
34
35}
36
37
43
44void EvtDsToPipPipPimPi0::init(){
45
46 //std::cout << "EvtDsToPipPipPimPi0 (Oct 17, 2022) ==> Initialization !" << std::endl;
47
48 // check that there are 0 arguments
49 checkNArg(0);
50 // check that there are 4 daughters
51 checkNDaug(4);
52 // check that parent and daughters are correct
53 if(getParentId()==EvtPDL::getId("D_s+")){
54 if(getDaug(0)!=EvtPDL::getId("pi+")){
55 report(ERROR,"EvtGen") << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!"<<endl;
56 ::abort();
57 }
58 if(getDaug(1)!=EvtPDL::getId("pi+")){
59 report(ERROR,"EvtGen") << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!"<<endl;
60 ::abort();
61 }
62 if(getDaug(2)!=EvtPDL::getId("pi-")){
63 report(ERROR,"EvtGen") << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!"<<endl;
64 ::abort();
65 }
66 if(getDaug(3)!=EvtPDL::getId("pi0")){
67 report(ERROR,"EvtGen") << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!"<<endl;
68 ::abort();
69 }
70 }else if(getParentId()==EvtPDL::getId("D_s-")){
71 if(getDaug(0)!=EvtPDL::getId("pi-")){
72 report(ERROR,"EvtGen") << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!"<<endl;
73 ::abort();
74 }
75 if(getDaug(1)!=EvtPDL::getId("pi-")){
76 report(ERROR,"EvtGen") << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!"<<endl;
77 ::abort();
78 }
79 if(getDaug(2)!=EvtPDL::getId("pi+")){
80 report(ERROR,"EvtGen") << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!"<<endl;
81 ::abort();
82 }
83 if(getDaug(3)!=EvtPDL::getId("pi0")){
84 report(ERROR,"EvtGen") << "In EvtDsToPipPipPimPi0, Decay Chain is Wrong!"<<endl;
85 ::abort();
86 }
87 }else{
88 report(ERROR,"EvtGen") << "In EvtDsToPipPipPimPi0, Mother is Wrong!"<<endl;
89 ::abort();
90 }
91
92}
93
94
95void EvtDsToPipPipPimPi0::initProbMax() {
96
97 setProbMax(10000.0);
98
99}
100double EvtDsToPipPipPimPi0::mass1[11] ={0.96500,0.77526,1.01946,0.78200,1.23000,1.23000,1.35000,1.27550,1.46500,0.77526,1.37100};
101double EvtDsToPipPipPimPi0::mass2[11] ={0.77526,0.77526,0.77526,0.77526,0.77526,0.77526,0.77526,0.77526,0.77526,1.46500,0.77526};
102double EvtDsToPipPipPimPi0::width1[11]={0.05500,0.14910,0.00425,0.00800,0.42000,0.42000,0.26500,0.18670,0.40000,0.14910,0.17000};
103double EvtDsToPipPipPimPi0::width2[11]={0.14910,0.14910,0.14910,0.14910,0.14910,0.14910,0.14910,0.14910,0.14910,0.40000,0.14910};
104double EvtDsToPipPipPimPi0::amp[11]={ 4.0026e-1,1.2202e-1, 7.4682e-2,1.9591e-1,5.7255e-1,3.1005e-1,1.0,7.6861e-1,1.1482e+0,5.7077e-1,2.7970e-1};
105double EvtDsToPipPipPimPi0::phase[11]={-2.2962e+0,1.0977e+0,2.8957e+0,3.2197e+0,3.7810e+0,-1.4586e+0,0.0,1.1081e+0,4.3439e-1,-1.6993e+0,2.2240e+0};
106double EvtDsToPipPipPimPi0::mD = 1.86486;
107double EvtDsToPipPipPimPi0::metap = 0.95778;
108double EvtDsToPipPipPimPi0::mKst0 = 0.89581;
109double EvtDsToPipPipPimPi0::GKst0 = 0.0474;
110double EvtDsToPipPipPimPi0::mKstp = 0.89166;
111double EvtDsToPipPipPimPi0::Gkstp = 0.0508;
112double EvtDsToPipPipPimPi0::mk0 = 0.497614;
113double EvtDsToPipPipPimPi0::mass_Kaon = 0.49368;
114double EvtDsToPipPipPimPi0::mass_Ks = 0.497611;
115double EvtDsToPipPipPimPi0::mass_Eta = 0.547862;
116double EvtDsToPipPipPimPi0::mass_Pion = 0.13957;
117double EvtDsToPipPipPimPi0::math_pi = 3.1415926;
118double EvtDsToPipPipPimPi0::mass_Pion2 = 0.0194797849;
119double EvtDsToPipPipPimPi0::mass_2Pion = 0.27914;
120double EvtDsToPipPipPimPi0::math_2pi = 6.2831852;
121double EvtDsToPipPipPimPi0::rD2 = 25.0; // 5*5
122double EvtDsToPipPipPimPi0::rRes1 = 9.0; // 3*3
123double EvtDsToPipPipPimPi0::rRes2 = 9.0; // 3*3
124double EvtDsToPipPipPimPi0::g1 = 0.5468;
125double EvtDsToPipPipPimPi0::g2 = 0.23; //K*0(1430)
126double EvtDsToPipPipPimPi0::GS1 = 0.636619783;
127double EvtDsToPipPipPimPi0::GS2 = 0.01860182466;
128double EvtDsToPipPipPimPi0::GS3 = 0.1591549458; // 1/(2*math_2pi)
129double EvtDsToPipPipPimPi0::GS4 = 0.00620060822;
130int EvtDsToPipPipPimPi0::G[4][4] = { {1,0,0,0}, {0,-1,0,0}, {0,0,-1,0}, {0,0,0,-1} };
131int EvtDsToPipPipPimPi0::E[4][4][4][4] =
132{ { {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
133 {{0,0,0,0}, {0,0,0,0}, {0,0,0,1}, {0,0,-1,0}},
134 {{0,0,0,0}, {0,0,0,-1}, {0,0,0,0}, {0,1,0,0} },
135 {{0,0,0,0}, {0,0,1,0}, {0,-1,0,0}, {0,0,0,0} } },
136{ {{0,0,0,0}, {0,0,0,0}, {0,0,0,-1}, {0,0,1,0} },
137 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
138 {{0,0,0,1}, {0,0,0,0}, {0,0,0,0}, {-1,0,0,0}},
139 {{0,0,-1,0}, {0,0,0,0}, {1,0,0,0}, {0,0,0,0} } },
140{ {{0,0,0,0}, {0,0,0,1}, {0,0,0,0}, {0,-1,0,0}},
141 {{0,0,0,-1}, {0,0,0,0}, {0,0,0,0}, {1,0,0,0} },
142 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} },
143 {{0,1,0,0}, {-1,0,0,0}, {0,0,0,0}, {0,0,0,0} } },
144{ {{0,0,0,0}, {0,0,-1,0}, {0,1,0,0}, {0,0,0,0} },
145 {{0,0,1,0}, {0,0,0,0}, {-1,0,0,0}, {0,0,0,0} },
146 {{0,-1,0,0}, {1,0,0,0}, {0,0,0,0}, {0,0,0,0} },
147 {{0,0,0,0}, {0,0,0,0}, {0,0,0,0}, {0,0,0,0} } } };
148
149
150void EvtDsToPipPipPimPi0::decay( EvtParticle *p){
151
152 p->initializePhaseSpace(getNDaug(),getDaugs());
153
154 EvtVector4R PPip1 = p->getDaug(0)->getP4();
155 EvtVector4R PPip2 = p->getDaug(1)->getP4();
156 EvtVector4R PPim = p->getDaug(2)->getP4();
157 EvtVector4R PPi0 = p->getDaug(3)->getP4();
158
159 double Pip1[4], Pip2[4], Pim[4], Pi0[4];
160 double cof[2],amp_tmp[2],amp_PDF[2], PDF[2];
161 if(getParentId()==EvtPDL::getId("D_s+")){
162 Pip1[0] = PPip1.get(0); Pip1[1] = PPip1.get(1); Pip1[2] = PPip1.get(2); Pip1[3] = PPip1.get(3);
163 Pip2[0] = PPip2.get(0); Pip2[1] = PPip2.get(1); Pip2[2] = PPip2.get(2); Pip2[3] = PPip2.get(3);
164 Pim[0] = PPim.get(0); Pim[1] = PPim.get(1); Pim[2] = PPim.get(2); Pim[3] = PPim.get(3);
165 Pi0[0] = PPi0.get(0); Pi0[1] = PPi0.get(1); Pi0[2] = PPi0.get(2); Pi0[3] = PPi0.get(3);
166 }else{
167 Pip1[0] = PPip1.get(0); Pip1[1] = -(PPip1.get(1)); Pip1[2] = -(PPip1.get(2)); Pip1[3] = -(PPip1.get(3));
168 Pip2[0] = PPip2.get(0); Pip2[1] = -(PPip2.get(1)); Pip2[2] = -(PPip2.get(2)); Pip2[3] = -(PPip2.get(3));
169 Pim[0] = PPim.get(0); Pim[1] = -(PPim.get(1)); Pim[2] = -(PPim.get(2)); Pim[3] = -(PPim.get(3));
170 Pi0[0] = PPi0.get(0); Pi0[1] = -(PPi0.get(1)); Pi0[2] = -(PPi0.get(2)); Pi0[3] = -(PPi0.get(3));
171 }
172
173 double p12[4],p13[4],p14[4],p23[4],p24[4],p34[4];
174 double p123[4],p124[4],p134[4],p234[4],pD[4];
175 double Sf1,Sf2,Sf3,Sf4,Sf5,Sf6,Sf7,Sf8;
176 Sf1 = spinf1( Pip1[0],Pip1[1],Pip1[2],Pip1[3],
177 Pip2[0],Pip2[1],Pip2[2],Pip2[3],
178 Pim[0] ,Pim[1] ,Pim[2] ,Pim[3] ,
179 Pi0[0] ,Pi0[1] ,Pi0[2] ,Pi0[3]);
180 Sf2 = spinf2( Pip1[0],Pip1[1],Pip1[2],Pip1[3],
181 Pip2[0],Pip2[1],Pip2[2],Pip2[3],
182 Pim[0] ,Pim[1] ,Pim[2] ,Pim[3] ,
183 Pi0[0] ,Pi0[1] ,Pi0[2] ,Pi0[3]);
184 Sf3 = spinf3( Pip1[0],Pip1[1],Pip1[2],Pip1[3],
185 Pip2[0],Pip2[1],Pip2[2],Pip2[3],
186 Pim[0] ,Pim[1] ,Pim[2] ,Pim[3] ,
187 Pi0[0] ,Pi0[1] ,Pi0[2] ,Pi0[3]);
188 Sf4 = spinf4( Pip1[0],Pip1[1],Pip1[2],Pip1[3],
189 Pip2[0],Pip2[1],Pip2[2],Pip2[3],
190 Pim[0] ,Pim[1] ,Pim[2] ,Pim[3] ,
191 Pi0[0] ,Pi0[1] ,Pi0[2] ,Pi0[3]);
192 Sf5 = spinf5( Pip1[0],Pip1[1],Pip1[2],Pip1[3],
193 Pip2[0],Pip2[1],Pip2[2],Pip2[3],
194 Pim[0] ,Pim[1] ,Pim[2] ,Pim[3] ,
195 Pi0[0] ,Pi0[1] ,Pi0[2] ,Pi0[3]);
196 Sf6 = spinf6( Pip1[0],Pip1[1],Pip1[2],Pip1[3],
197 Pip2[0],Pip2[1],Pip2[2],Pip2[3],
198 Pim[0] ,Pim[1] ,Pim[2] ,Pim[3] ,
199 Pi0[0] ,Pi0[1] ,Pi0[2] ,Pi0[3]);
200 Sf7 = spinf7( Pip1[0],Pip1[1],Pip1[2],Pip1[3],
201 Pip2[0],Pip2[1],Pip2[2],Pip2[3],
202 Pim[0] ,Pim[1] ,Pim[2] ,Pim[3] ,
203 Pi0[0] ,Pi0[1] ,Pi0[2] ,Pi0[3]);
204 Sf8 = spinf8( Pip1[0],Pip1[1],Pip1[2],Pip1[3],
205 Pip2[0],Pip2[1],Pip2[2],Pip2[3],
206 Pim[0] ,Pim[1] ,Pim[2] ,Pim[3] ,
207 Pi0[0] ,Pi0[1] ,Pi0[2] ,Pi0[3]);
208
209
210 for(int i=0;i!=4;i++){
211 p12[i]=Pip1[i]+Pip2[i];
212 p13[i]=Pip1[i]+Pim[i];
213 p14[i]=Pip1[i]+Pi0[i];
214 p23[i]=Pip2[i]+Pim[i];
215 p24[i]=Pip2[i]+Pi0[i];
216 p34[i]=Pim[i]+Pi0[i];
217 p123[i]=Pip1[i]+Pip2[i]+Pim[i];
218 p124[i]=Pip1[i]+Pip2[i]+Pi0[i];
219 p134[i]=Pip1[i]+Pim[i]+Pi0[i];
220 p234[i]=Pip2[i]+Pim[i]+Pi0[i];
221 pD[i]=Pip1[i]+Pip2[i]+Pim[i]+Pi0[i];
222 }
223 double sPip1,sPip2,sPim,sPi0,s12,s13,s14,s23,s24,s34,s123,s124,s134,s234,sD;
224 sPi0 = SCADot(Pi0,Pi0);
225 sPip1 = SCADot(Pip1,Pip1);
226 sPip2 = SCADot(Pip2,Pip2);
227 sPim = SCADot(Pim,Pim);
228
229 s12 = SCADot(p12,p12);
230 s13 = SCADot(p13,p13);
231 s14 = SCADot(p14,p14);
232 s23 = SCADot(p23,p23);
233 s24 = SCADot(p24,p24);
234 s34 = SCADot(p34,p34);
235
236 s123 = SCADot(p123,p123);
237 s124 = SCADot(p124,p124);
238 s134 = SCADot(p134,p134);
239 s234 = SCADot(p234,p234);
240
241 sD = SCADot(pD,pD);
242
243 double spion12[2]={mass_Kaon*mass_Kaon,sPip1};
244 double spion22[2]={mass_Kaon*mass_Kaon,sPip2};
245 double spim2[2]={mass_Kaon*mass_Kaon,sPim};
246
247 double t1_24[4],t1_14[4],t1_13[4],t1_23[4],t1_12[4],t1_34[4];
248 double t1_14_23[4],t1_13_24[4],t1_12_34[4];
249 double t1_134_2[4],t1_234_1[4],t1_123_4[4];
250 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];
251
252 double t2_13[4][4],t2_23[4][4];
253 double t2_13_24[4][4],t2_14_23[4][4];
254 double t2_123_4[4][4];
255 double t2_13_2[4][4],t2_23_1[4][4];
256
257 calt1(Pi0,Pip2,t1_24);
258 calt1(Pi0,Pip1,t1_14);
259
260 calt1(Pip1,Pim,t1_13);
261 calt1(Pip2,Pim,t1_23);
262
263 calt1(Pip1,Pip2,t1_12);
264 calt1(Pim,Pi0,t1_34);
265
266 calt1(p14,p23,t1_14_23);
267 calt1(p24,p13,t1_13_24);
268 calt1(p12,p34,t1_12_34);
269
270 calt1(p134,Pip2,t1_134_2);
271 calt1(p234,Pip1,t1_234_1);
272 calt1(p123,Pi0,t1_123_4);
273
274 calt1(p13,Pi0,t1_13_4);
275 calt1(p23,Pi0,t1_23_4);
276 calt1(p14,Pim,t1_14_3);
277 calt1(p24,Pim,t1_24_3);
278 calt1(p34,Pip1,t1_34_1);
279 calt1(p34,Pip2,t1_34_2);
280
281 calt2(Pip1,Pim,t2_13);
282 calt2(Pip2,Pim,t2_23);
283
284 calt2(p13,p24,t2_13_24);
285 calt2(p14,p23,t2_14_23);
286
287 calt2(p123,Pi0,t2_123_4);
288
289 calt2(p13,Pip2,t2_13_2);
290 calt2(p23,Pip1,t2_23_1);
291
292
293 amp_PDF[0] = 0;
294 amp_PDF[1] = 0;
295 PDF[0] = 0;
296 PDF[1] = 0;
297
298 double temp_PDF, tmp1, tmp2, tmp3,tmp4,temp_PDF1,amp_tmp1[2],amp_tmp2[2];
299 double pro[2], pro0[2], pro1[2],pro2[2],pro3[2],pro4[2];
300 double t1D[4],B[3], Bc[3];
301 double mass1sq, mass2sq;
302
303 for(int i=0; i<11; i++){
304 amp_tmp[0] = 0;
305 amp_tmp[1] = 0;
306 tmp1 = 0;tmp2 = 0;temp_PDF = 0;;
307 cof[0] = amp[i]*cos(phase[i]);
308 cof[1] = amp[i]*sin(phase[i]);
309 mass1sq = mass1[i]*mass1[i];
310 mass2sq = mass2[i]*mass2[i];
311 temp_PDF = 0;
312 temp_PDF1 = 0;
313 amp_tmp[0] = 0;
314 amp_tmp[1] = 0;
315 amp_tmp1[0] = 0;
316 amp_tmp1[1] = 0;
317 amp_tmp2[0] = 0;
318 amp_tmp2[1] = 0;
319 tmp1 = 0;
320 tmp2 = 0;
321 tmp3 = 0;
322 tmp4 = 0;
323
324 if(i==0){
325 temp_PDF=0;
326 propagator980(mass1[i],s23,spion12,spion22,pro0);
327 propagatorGS(mass2sq,mass2[i],width2[i],s14,sPip1,sPi0,rRes2,pro1);
328 Com_Multi(pro0,pro1,pro);
329
330 for(int a=0; a<4; a++){
331 temp_PDF += G[a][a]*t1_14_23[a]*t1_14[a];
332 }
333
334 B[1] = barrier(1,s14,sPip1,sPi0,rRes2,mass2[i]);
335 B[2] = barrier(1,sD,s23,s14,rD2,1.9683);
336 tmp1 = B[1]*B[2]*temp_PDF;
337 amp_tmp1[0] = tmp1*pro[0];
338 amp_tmp1[1] = tmp1*pro[1];
339
340 temp_PDF=0;
341 propagator980(mass1[i],s13,spion12,spion22,pro0);
342 propagatorGS(mass2sq,mass2[i],width2[i],s24,sPip2,sPi0,rRes2,pro1);
343 Com_Multi(pro0,pro1,pro);
344
345 for(int a=0; a<4; a++){
346 temp_PDF += G[a][a]*t1_13_24[a]*t1_24[a];
347 }
348
349 B[1] = barrier(1,s24,sPip2,sPi0,rRes2,mass2[i]);
350 B[2] = barrier(1,sD,s13,s24,rD2,1.9683);
351 tmp2 = B[1]*B[2]*temp_PDF;
352 amp_tmp2[0] = tmp2*pro[0];
353 amp_tmp2[1] = tmp2*pro[1];
354 }
355
356 if(i==1){
357 temp_PDF=0;
358 propagatorGS(mass2sq,mass2[i],width2[i],s14,sPip1,sPi0,rRes2,pro0);
359 propagatorGS(mass2sq,mass2[i],width2[i],s23,sPip2,sPim,rRes2,pro1);
360
361 for(int a=0; a<4; a++){
362 temp_PDF += G[a][a]*t1_14[a]*t1_23[a];
363 }
364
365 B[1] = barrier(1,s14,sPip1,sPi0,rRes2,mass2[i]);
366 B[2] = barrier(1,s23,sPip2,sPim,rRes2,mass2[i]);
367 tmp1 = B[1]*B[2]*temp_PDF;
368 amp_tmp1[0] = tmp1*pro0[0]*pro1[0];
369 amp_tmp1[1] = tmp1*pro0[1]*pro1[1];
370
371 temp_PDF=0;
372 propagatorGS(mass2sq,mass2[i],width2[i],s24,sPip2,sPi0,rRes2,pro0);
373 propagatorGS(mass2sq,mass2[i],width2[i],s13,sPip1,sPim,rRes2,pro1);
374
375 for(int a=0; a<4; a++){
376 temp_PDF += G[a][a]*t1_24[a]*t1_13[a];
377 }
378
379 B[1] = barrier(1,s24,sPip2,sPi0,rRes2,mass2[i]);
380 B[2] = barrier(1,s13,sPip1,sPim,rRes2,mass2[i]);
381 tmp2 = B[1]*B[2]*temp_PDF;
382 amp_tmp2[0] = tmp2*pro0[0]*pro1[0];
383 amp_tmp2[1] = tmp2*pro0[1]*pro1[1];
384 }
385
386 if(i==2){
387 amp_tmp1[0]=0;
388 amp_tmp1[1]=0;
389 amp_tmp2[0]=0;
390 amp_tmp2[1]=0;
391
392 temp_PDF=0;
393 propagatorRBW(mass1sq,mass1[i],width1[i],s134,s14,sPim,rRes1,1,pro1);
394 propagatorGS(mass2sq,mass2[i],width2[i],s14,sPip1,sPi0,rRes2,pro0);
395
396 temp_PDF=Sf3;
397
398 B[0] = barrier(1,sD,s134,sPip2,rD2,1.9683);
399 B[1] = barrier(1,s134,s14,sPim,rRes1,mass1[i]);
400 B[2] = barrier(1,s14,sPip1,sPi0,rRes2,mass2[i]);
401 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
402 amp_tmp1[0] += tmp1*pro0[0]*pro1[0];
403 amp_tmp1[1] += tmp1*pro0[1]*pro1[1];
404
405 temp_PDF=0;
406 propagatorRBW(mass1sq,mass1[i],width1[i],s234,s24,sPim,rRes1,1,pro1);
407 propagatorGS(mass2sq,mass2[i],width2[i],s24,sPip2,sPi0,rRes2,pro0);
408
409 temp_PDF=Sf4;
410
411 B[0] = barrier(1,sD,s234,sPip1,rD2,1.9683);
412 B[1] = barrier(1,s234,s24,sPim,rRes1,mass1[i]);
413 B[2] = barrier(1,s24,sPip2,sPi0,rRes2,mass2[i]);
414 tmp2 = B[0]*B[1]*B[2]*temp_PDF;
415 amp_tmp2[0] += tmp2*pro0[0]*pro1[0];
416 amp_tmp2[1] += tmp2*pro0[1]*pro1[1];
417
418 temp_PDF=0;
419 propagatorRBW(mass1sq,mass1[i],width1[i],s134,s14,sPim,rRes1,1,pro1);
420 propagatorGS(mass2sq,mass2[i],width2[i],s13,sPip1,sPim,rRes2,pro0);
421
422 temp_PDF=Sf5;
423
424 B[0] = barrier(1,sD,s134,sPip2,rD2,1.9683);
425 B[1] = barrier(1,s134,s13,sPi0,rRes1,mass1[i]);
426 B[2] = barrier(1,s13,sPip1,sPim,rRes2,mass2[i]);
427 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
428 amp_tmp1[0] -= tmp1*pro0[0]*pro1[0];
429 amp_tmp1[1] -= tmp1*pro0[1]*pro1[1];
430
431 temp_PDF=0;
432 propagatorRBW(mass1sq,mass1[i],width1[i],s234,s24,sPim,rRes1,1,pro1);
433 propagatorGS(mass2sq,mass2[i],width2[i],s23,sPip2,sPim,rRes2,pro0);
434
435 temp_PDF=Sf6;
436
437 B[0] = barrier(1,sD,s234,sPip1,rD2,1.9683);
438 B[1] = barrier(1,s234,s23,sPi0,rRes1,mass1[i]);
439 B[2] = barrier(1,s23,sPip2,sPim,rRes2,mass2[i]);
440 tmp2 = B[0]*B[1]*B[2]*temp_PDF;
441 amp_tmp2[0] -= tmp2*pro0[0]*pro1[0];
442 amp_tmp2[1] -= tmp2*pro0[1]*pro1[1];
443
444 temp_PDF=0;
445 propagatorRBW(mass1sq,mass1[i],width1[i],s134,s34,sPip1,rRes1,1,pro1);
446 propagatorGS(mass2sq,mass2[i],width2[i],s34,sPim,sPi0,rRes2,pro0);
447
448 temp_PDF=Sf7;
449
450 B[0] = barrier(1,sD,s134,sPip2,rD2,1.9683);
451 B[1] = barrier(1,s134,s34,sPip1,rRes1,mass1[i]);
452 B[2] = barrier(1,s34,sPim,sPi0,rRes2,mass2[i]);
453 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
454 amp_tmp1[0] -= tmp1*pro0[0]*pro1[0];
455 amp_tmp1[1] -= tmp1*pro0[1]*pro1[1];
456
457 temp_PDF=0;
458 propagatorRBW(mass1sq,mass1[i],width1[i],s234,s34,sPip2,rRes1,1,pro1);
459 propagatorGS(mass2sq,mass2[i],width2[i],s34,sPim,sPi0,rRes2,pro0);
460
461 temp_PDF=Sf8;
462
463 B[0] = barrier(1,sD,s234,sPip1,rD2,1.9683);
464 B[1] = barrier(1,s234,s34,sPip2,rRes1,mass1[i]);
465 B[2] = barrier(1,s34,sPim,sPi0,rRes2,mass2[i]);
466 tmp2 = B[0]*B[1]*B[2]*temp_PDF;
467 amp_tmp2[0] -= tmp2*pro0[0]*pro1[0];
468 amp_tmp2[1] -= tmp2*pro0[1]*pro1[1];
469 }
470
471 if(i==3){
472 amp_tmp1[0]=0;
473 amp_tmp1[1]=0;
474 amp_tmp2[0]=0;
475 amp_tmp2[1]=0;
476
477 temp_PDF=0;
478 propagatorRBW(mass1sq,mass1[i],width1[i],s134,s14,sPim,rRes1,1,pro1);
479 propagatorGS(mass2sq,mass2[i],width2[i],s14,sPip1,sPi0,rRes2,pro0);
480
481 temp_PDF=Sf3;
482
483 B[0] = barrier(1,sD,s134,sPip2,rD2,1.9683);
484 B[1] = barrier(1,s134,s14,sPim,rRes1,mass1[i]);
485 B[2] = barrier(1,s14,sPip1,sPi0,rRes2,mass2[i]);
486 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
487 amp_tmp1[0] += tmp1*pro0[0]*pro1[0];
488 amp_tmp1[1] += tmp1*pro0[1]*pro1[1];
489
490 temp_PDF=0;
491 propagatorRBW(mass1sq,mass1[i],width1[i],s234,s24,sPim,rRes1,1,pro1);
492 propagatorGS(mass2sq,mass2[i],width2[i],s24,sPip2,sPi0,rRes2,pro0);
493
494 temp_PDF=Sf4;
495
496 B[0] = barrier(1,sD,s234,sPip1,rD2,1.9683);
497 B[1] = barrier(1,s234,s24,sPim,rRes1,mass1[i]);
498 B[2] = barrier(1,s24,sPip2,sPi0,rRes2,mass2[i]);
499 tmp2 = B[0]*B[1]*B[2]*temp_PDF;
500 amp_tmp2[0] += tmp2*pro0[0]*pro1[0];
501 amp_tmp2[1] += tmp2*pro0[1]*pro1[1];
502
503 temp_PDF=0;
504 propagatorRBW(mass1sq,mass1[i],width1[i],s134,s14,sPim,rRes1,1,pro1);
505 propagatorGS(mass2sq,mass2[i],width2[i],s13,sPip1,sPim,rRes2,pro0);
506
507 temp_PDF=Sf5;
508
509 B[0] = barrier(1,sD,s134,sPip2,rD2,1.9683);
510 B[1] = barrier(1,s134,s13,sPi0,rRes1,mass1[i]);
511 B[2] = barrier(1,s13,sPip1,sPim,rRes2,mass2[i]);
512 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
513 amp_tmp1[0] -= tmp1*pro0[0]*pro1[0];
514 amp_tmp1[1] -= tmp1*pro0[1]*pro1[1];
515
516 temp_PDF=0;
517 propagatorRBW(mass1sq,mass1[i],width1[i],s234,s24,sPim,rRes1,1,pro1);
518 propagatorGS(mass2sq,mass2[i],width2[i],s23,sPip2,sPim,rRes2,pro0);
519
520 temp_PDF=Sf6;
521
522 B[0] = barrier(1,sD,s234,sPip1,rD2,1.9683);
523 B[1] = barrier(1,s234,s23,sPi0,rRes1,mass1[i]);
524 B[2] = barrier(1,s23,sPip2,sPim,rRes2,mass2[i]);
525 tmp2 = B[0]*B[1]*B[2]*temp_PDF;
526 amp_tmp2[0] -= tmp2*pro0[0]*pro1[0];
527 amp_tmp2[1] -= tmp2*pro0[1]*pro1[1];
528
529 temp_PDF=0;
530 propagatorRBW(mass1sq,mass1[i],width1[i],s134,s34,sPip1,rRes1,1,pro1);
531 propagatorGS(mass2sq,mass2[i],width2[i],s34,sPim,sPi0,rRes2,pro0);
532
533 temp_PDF=Sf7;
534
535 B[0] = barrier(1,sD,s134,sPip2,rD2,1.9683);
536 B[1] = barrier(1,s134,s34,sPip1,rRes1,mass1[i]);
537 B[2] = barrier(1,s34,sPim,sPi0,rRes2,mass2[i]);
538 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
539 amp_tmp1[0] -= tmp1*pro0[0]*pro1[0];
540 amp_tmp1[1] -= tmp1*pro0[1]*pro1[1];
541
542 temp_PDF=0;
543 propagatorRBW(mass1sq,mass1[i],width1[i],s234,s34,sPip2,rRes1,1,pro1);
544 propagatorGS(mass2sq,mass2[i],width2[i],s34,sPim,sPi0,rRes2,pro0);
545
546 temp_PDF=Sf8;
547
548 B[0] = barrier(1,sD,s234,sPip1,rD2,1.9683);
549 B[1] = barrier(1,s234,s34,sPip2,rRes1,mass1[i]);
550 B[2] = barrier(1,s34,sPim,sPi0,rRes2,mass2[i]);
551 tmp2 = B[0]*B[1]*B[2]*temp_PDF;
552 amp_tmp2[0] -= tmp2*pro0[0]*pro1[0];
553 amp_tmp2[1] -= tmp2*pro0[1]*pro1[1];
554 }
555
556 if(i==4){
557 temp_PDF=0;
558 //propagatorRBW(mass1sq,mass1[i],width1[i],s123,s13,sPip2,rRes1,0,pro1);
559 propagatorMW(1.230,s123,1,pro1);
560 propagatorGS(mass2sq,mass2[i],width2[i],s13,sPip1,sPim,rRes2,pro0);
561
562 for(int a=0; a<4; a++)
563 for(int b=0;b<4;b++)
564 temp_PDF += G[a][a]*G[b][b]*t1_123_4[a]*(-G[a][b]+p123[a]*p123[b]/s123)*t1_13[b];
565
566 B[0] = barrier(1,sD,s123,sPi0,rD2,1.9683);
567 B[1] = barrier(1,s13,sPip1,sPim,rRes2,mass2[i]);
568 tmp1 = B[0]*B[1]*temp_PDF;
569 amp_tmp1[0] = tmp1*pro0[0]*pro1[0];
570 amp_tmp1[1] = tmp1*pro0[1]*pro1[1];
571
572 temp_PDF=0;
573 //propagatorRBW(mass1sq,mass1[i],width1[i],s123,s23,sPip1,rRes1,0,pro1);
574 propagatorMW(1.230,s123,1,pro1);
575 propagatorGS(mass2sq,mass2[i],width2[i],s23,sPip2,sPim,rRes2,pro0);
576
577 for(int a=0; a<4; a++)
578 for(int b=0;b<4;b++)
579 temp_PDF += G[a][a]*G[b][b]*t1_123_4[a]*(-G[a][b]+p123[a]*p123[b]/s123)*t1_23[b];
580
581 B[0] = barrier(1,sD,s123,sPi0,rD2,1.9683);
582 B[1] = barrier(1,s23,sPip2,sPim,rRes2,mass2[i]);
583 tmp2 = B[0]*B[1]*temp_PDF;
584 amp_tmp2[0] = tmp2*pro0[0]*pro1[0];
585 amp_tmp2[1] = tmp2*pro0[1]*pro1[1];
586 }
587
588 if(i==5){
589 amp_tmp1[0]=0;
590 amp_tmp1[1]=0;
591 amp_tmp2[0]=0;
592 amp_tmp2[1]=0;
593
594 temp_PDF=0;
595 //propagatorRBW(mass1sq,mass1[i],width1[i],s134,s14,sPim,rRes1,0,pro1);
596 propagatorMW(1.230,s134,0,pro1);
597 propagatorGS(mass2sq,mass2[i],width2[i],s14,sPip1,sPi0,rRes2,pro0);
598
599 for(int a=0; a<4; a++)
600 for(int b=0;b<4;b++)
601 temp_PDF += G[a][a]*G[b][b]*t1_134_2[a]*(-G[a][b]+p134[a]*p134[b]/s134)*t1_14[b];
602
603 B[0] = barrier(1,sD,s134,sPip2,rD2,1.9683);
604 B[1] = barrier(1,s14,sPip1,sPi0,rRes2,mass2[i]);
605 tmp1 = B[0]*B[1]*temp_PDF;
606 amp_tmp1[0] += tmp1*pro0[0]*pro1[0]*1.;
607 amp_tmp1[1] += tmp1*pro0[1]*pro1[1]*1.;
608
609 temp_PDF=0;
610 //propagatorRBW(mass1sq,mass1[i],width1[i],s234,s24,sPim,rRes1,0,pro1);
611 propagatorMW(1.230,s234,0,pro1);
612 propagatorGS(mass2sq,mass2[i],width2[i],s24,sPip2,sPi0,rRes2,pro0);
613
614 for(int a=0; a<4; a++)
615 for(int b=0;b<4;b++)
616 temp_PDF += G[a][a]*G[b][b]*t1_234_1[a]*(-G[a][b]+p234[a]*p234[b]/s234)*t1_24[b];
617
618 B[0] = barrier(1,sD,s234,sPip1,rD2,1.9683);
619 B[1] = barrier(1,s24,sPip2,sPi0,rRes2,mass2[i]);
620 tmp2 = B[0]*B[1]*temp_PDF;
621 amp_tmp2[0] += tmp2*pro0[0]*pro1[0]*1.;
622 amp_tmp2[1] += tmp2*pro0[1]*pro1[1]*1.;
623
624 // a10 -> rho0 pi0 violation of C-Parity
625 //temp_PDF=0;
626 //propagatorRBW(mass1sq,mass1[i],width1[i],s134,s14,sPim,rRes1,0,pro1);
627 //propagatorGS(mass2sq,mass2[i],width2[i],s13,sPip1,sPim,rRes2,pro0);
628 //
629 //for(int a=0; a<4; a++)
630 // for(int b=0;b<4;b++)
631 // temp_PDF += G[a][a]*G[b][b]*t1_134_2[a]*(-G[a][b]+p134[a]*p134[b]/s134)*t1_13[b];
632 //
633 //B[0] = barrier(1,sD,s134,sPip2,rD2,1.9683);
634 //B[1] = barrier(1,s13,sPip1,sPim,rRes2,mass2[i]);
635 //tmp1 = B[0]*B[1]*temp_PDF;
636 //amp_tmp1[0] += tmp1*pro0[0]*pro1[0];
637 //amp_tmp1[1] += tmp1*pro0[1]*pro1[1];
638 //
639 //temp_PDF=0;
640 //propagatorRBW(mass1sq,mass1[i],width1[i],s234,s24,sPim,rRes1,0,pro1);
641 //propagatorGS(mass2sq,mass2[i],width2[i],s23,sPip2,sPim,rRes2,pro0);
642 //
643 //for(int a=0; a<4; a++)
644 // for(int b=0;b<4;b++)
645 // temp_PDF += G[a][a]*G[b][b]*t1_234_1[a]*(-G[a][b]+p234[a]*p234[b]/s234)*t1_23[b];
646 //
647 //B[0] = barrier(1,sD,s234,sPip1,rD2,1.9683);
648 //B[1] = barrier(1,s23,sPip2,sPim,rRes2,mass2[i]);
649 //tmp2 = B[0]*B[1]*temp_PDF;
650 //amp_tmp2[0] += tmp2*pro0[0]*pro1[0];
651 //amp_tmp2[1] += tmp2*pro0[1]*pro1[1];
652
653 temp_PDF=0;
654 //propagatorRBW(mass1sq,mass1[i],width1[i],s134,s34,sPip1,rRes1,0,pro1);
655 propagatorMW(1.230,s134,0,pro1);
656 propagatorGS(mass2sq,mass2[i],width2[i],s34,sPim,sPi0,rRes2,pro0);
657
658 for(int a=0; a<4; a++)
659 for(int b=0;b<4;b++)
660 temp_PDF += G[a][a]*G[b][b]*t1_134_2[a]*(-G[a][b]+p134[a]*p134[b]/s134)*t1_34[b];
661
662 B[0] = barrier(1,sD,s134,sPip2,rD2,1.9683);
663 B[1] = barrier(1,s34,sPim,sPi0,rRes2,mass2[i]);
664 tmp1 = B[0]*B[1]*temp_PDF;
665 amp_tmp1[0] += tmp1*pro0[0]*pro1[0]*-1.;
666 amp_tmp1[1] += tmp1*pro0[1]*pro1[1]*-1.;
667
668 temp_PDF=0;
669 //propagatorRBW(mass1sq,mass1[i],width1[i],s234,s34,sPip2,rRes1,0,pro1);
670 propagatorMW(1.230,s234,0,pro1);
671 propagatorGS(mass2sq,mass2[i],width2[i],s34,sPim,sPi0,rRes2,pro0);
672
673 for(int a=0; a<4; a++)
674 for(int b=0;b<4;b++)
675 temp_PDF += G[a][a]*G[b][b]*t1_234_1[a]*(-G[a][b]+p234[a]*p234[b]/s234)*t1_34[b];
676
677 B[0] = barrier(1,sD,s234,sPip1,rD2,1.9683);
678 B[1] = barrier(1,s34,sPim,sPi0,rRes2,mass2[i]);
679 tmp2 = B[0]*B[1]*temp_PDF;
680 amp_tmp2[0] += tmp2*pro0[0]*pro1[0]*-1.;
681 amp_tmp2[1] += tmp2*pro0[1]*pro1[1]*-1.;
682 }
683
684 if(i==6){
685 temp_PDF=0;
686 propagatorRBW(mass1sq,mass1[i],width1[i],s23,sPip2,sPim,rRes1,0,pro0);
687 propagatorGS(mass2sq,mass2[i],width2[i],s14,sPip1,sPi0,rRes2,pro1);
688 Com_Multi(pro0,pro1,pro);
689
690 for(int a=0; a<4; a++){
691 temp_PDF += G[a][a]*t1_14_23[a]*t1_14[a];
692 }
693
694 B[1] = barrier(1,s14,sPip1,sPi0,rRes2,mass2[i]);
695 B[2] = barrier(1,sD,s23,s14,rD2,1.9683);
696 tmp1 = B[1]*B[2]*temp_PDF;
697 amp_tmp1[0] = tmp1*pro[0];
698 amp_tmp1[1] = tmp1*pro[1];
699
700 temp_PDF=0;
701 propagatorRBW(mass1sq,mass1[i],width1[i],s13,sPip1,sPim,rRes1,0,pro0);
702 propagatorGS(mass2sq,mass2[i],width2[i],s24,sPip2,sPi0,rRes2,pro1);
703 Com_Multi(pro0,pro1,pro);
704
705 for(int a=0; a<4; a++){
706 temp_PDF += G[a][a]*t1_13_24[a]*t1_24[a];
707 }
708
709 B[1] = barrier(1,s24,sPip2,sPi0,rRes2,mass2[i]);
710 B[2] = barrier(1,sD,s13,s24,rD2,1.9683);
711 tmp2 = B[1]*B[2]*temp_PDF;
712 amp_tmp2[0] = tmp2*pro[0];
713 amp_tmp2[1] = tmp2*pro[1];
714 }
715
716 if(i==7){
717 temp_PDF=0;
718 propagatorRBW(mass1sq,mass1[i],width1[i],s23,sPip2,sPim,rRes1,2,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++){
723 for(int b=0; b<4; b++){
724 temp_PDF += G[a][a]*G[b][b]*t1_14_23[a]*t1_14[b]*t2_23[a][b];
725 }
726 }
727
728 B[0] = barrier(2,s23,sPip2,sPim,rRes1,mass1[i]);
729 B[1] = barrier(1,s14,sPip1,sPi0,rRes2,mass2[i]);
730 B[2] = barrier(1,sD,s23,s14,rD2,1.9683);
731 tmp1 = B[0]*B[1]*B[2]*temp_PDF;
732 amp_tmp1[0] = tmp1*pro[0];
733 amp_tmp1[1] = tmp1*pro[1];
734
735 temp_PDF=0;
736 propagatorRBW(mass1sq,mass1[i],width1[i],s13,sPip1,sPim,rRes1,2,pro0);
737 propagatorGS(mass2sq,mass2[i],width2[i],s24,sPip2,sPi0,rRes2,pro1);
738 Com_Multi(pro0,pro1,pro);
739
740 for(int a=0; a<4; a++){
741 for(int b=0; b<4; b++){
742 temp_PDF += G[a][a]*G[b][b]*t1_13_24[a]*t1_24[b]*t2_13[a][b];
743 }
744 }
745
746 B[0] = barrier(2,s13,sPip1,sPim,rRes1,mass1[i]);
747 B[1] = barrier(1,s24,sPip2,sPi0,rRes2,mass2[i]);
748 B[2] = barrier(1,sD,s13,s24,rD2,1.9683);
749 tmp2 = B[0]*B[1]*B[2]*temp_PDF;
750 amp_tmp2[0] = tmp2*pro[0];
751 amp_tmp2[1] = tmp2*pro[1];
752 }
753
754 if(i==8){
755 temp_PDF=0;
756 propagatorRBW(mass1sq,mass1[i],width1[i],s14,sPip1,sPi0,rRes1,1,pro0);
757 propagatorGS(mass2sq,mass2[i],width2[i],s23,sPip2,sPim,rRes2,pro1);
758
759 for(int a=0; a<4; a++){
760 temp_PDF += G[a][a]*t1_14[a]*t1_23[a];
761 }
762
763 B[1] = barrier(1,s14,sPip1,sPi0,rRes2,mass2[i]);
764 B[2] = barrier(1,s23,sPip2,sPim,rRes2,mass2[i]);
765 tmp1 = B[1]*B[2]*temp_PDF;
766 amp_tmp1[0] = tmp1*pro0[0]*pro1[0];
767 amp_tmp1[1] = tmp1*pro0[1]*pro1[1];
768
769 temp_PDF=0;
770 propagatorRBW(mass1sq,mass1[i],width1[i],s24,sPip2,sPi0,rRes1,1,pro0);
771 propagatorGS(mass2sq,mass2[i],width2[i],s13,sPip1,sPim,rRes2,pro1);
772
773 for(int a=0; a<4; a++){
774 temp_PDF += G[a][a]*t1_24[a]*t1_13[a];
775 }
776
777 B[1] = barrier(1,s24,sPip2,sPi0,rRes2,mass2[i]);
778 B[2] = barrier(1,s13,sPip1,sPim,rRes2,mass2[i]);
779 tmp2 = B[1]*B[2]*temp_PDF;
780 amp_tmp2[0] = tmp2*pro0[0]*pro1[0];
781 amp_tmp2[1] = tmp2*pro0[1]*pro1[1];
782 }
783
784 if(i==9){
785 temp_PDF=0;
786 propagatorGS(mass1sq,mass1[i],width1[i],s14,sPip1,sPi0,rRes1,pro0);
787 propagatorRBW(mass2sq,mass2[i],width2[i],s23,sPip2,sPim,rRes2,1,pro1);
788
789 temp_PDF=Sf1;
790
791 B[0] = barrier(1,sD,s14,s23,rD2,1.9683);
792 B[1] = barrier(1,s14,sPip1,sPi0,rRes2,mass2[i]);
793 B[2] = barrier(1,s23,sPip2,sPim,rRes2,mass2[i]);
794 tmp1 = B[1]*B[2]*B[0]*temp_PDF;
795 amp_tmp1[0] = tmp1*pro0[0]*pro1[0];
796 amp_tmp1[1] = tmp1*pro0[1]*pro1[1];
797
798 temp_PDF=0;
799 propagatorGS(mass1sq,mass1[i],width1[i],s24,sPip2,sPi0,rRes1,pro0);
800 propagatorRBW(mass2sq,mass2[i],width2[i],s13,sPip1,sPim,rRes2,1,pro1);
801
802 temp_PDF=Sf2;
803
804 B[0] = barrier(1,sD,s24,s13,rD2,1.9683);
805 B[1] = barrier(1,s24,sPip2,sPi0,rRes2,mass2[i]);
806 B[2] = barrier(1,s13,sPip1,sPim,rRes2,mass2[i]);
807 tmp2 = B[1]*B[2]*B[0]*temp_PDF;
808 amp_tmp2[0] = tmp2*pro0[0]*pro1[0];
809 amp_tmp2[1] = tmp2*pro0[1]*pro1[1];
810 }
811
812 if(i==10){
813 amp_tmp1[0]=0;
814 amp_tmp1[1]=0;
815 amp_tmp2[0]=0;
816 amp_tmp2[1]=0;
817
818 temp_PDF=0;
819 propagatorRBW(mass1sq,mass1[i],width1[i],s134,s14,sPim,rRes1,1,pro1);
820 propagatorGS(mass2sq,mass2[i],width2[i],s14,sPip1,sPi0,rRes2,pro0);
821
822 for(int a=0; a<4; a++)
823 temp_PDF += G[a][a]*t1_134_2[a]*t1_14[a];
824 //temp_PDF += G[a][a]*G[b][b]*t1_134_2[a]*(-G[a][b]+p134[a]*p134[b]/s134)*t1_14[b];
825
826 B[0] = barrier(1,s134,s14,sPim,rRes1,mass1[i]);
827 B[1] = barrier(1,s14,sPip1,sPi0,rRes2,mass2[i]);
828 tmp1 = B[0]*B[1]*temp_PDF;
829 amp_tmp1[0] += tmp1*pro0[0]*pro1[0]*1.;
830 amp_tmp1[1] += tmp1*pro0[1]*pro1[1]*1.;
831
832 temp_PDF=0;
833 propagatorRBW(mass1sq,mass1[i],width1[i],s234,s24,sPim,rRes1,1,pro1);
834 propagatorGS(mass2sq,mass2[i],width2[i],s24,sPip2,sPi0,rRes2,pro0);
835
836 for(int a=0; a<4; a++)
837 temp_PDF += G[a][a]*t1_234_1[a]*t1_24[a];
838 //temp_PDF += G[a][a]*G[b][b]*t1_234_1[a]*(-G[a][b]+p234[a]*p234[b]/s234)*t1_24[b];
839
840 B[0] = barrier(1,s234,s24,sPim,rRes1,mass1[i]);
841 B[1] = barrier(1,s24,sPip2,sPi0,rRes2,mass2[i]);
842 tmp2 = B[0]*B[1]*temp_PDF;
843 amp_tmp2[0] += tmp2*pro0[0]*pro1[0]*1.;
844 amp_tmp2[1] += tmp2*pro0[1]*pro1[1]*1.;
845
846 // pi1300 -> rho0 pi0 violation of C-Parity
847 //temp_PDF=0;
848 //propagatorRBW(mass1sq,mass1[i],width1[i],s134,s14,sPim,rRes1,0,pro1);
849 //propagatorGS(mass2sq,mass2[i],width2[i],s13,sPip1,sPim,rRes2,pro0);
850 //
851 //for(int a=0; a<4; a++)
852 // for(int b=0;b<4;b++)
853 // temp_PDF += G[a][a]*G[b][b]*t1_134_2[a]*(-G[a][b]+p134[a]*p134[b]/s134)*t1_13[b];
854 //
855 //B[0] = barrier(1,sD,s134,sPip2,rD2,1.9683);
856 //B[1] = barrier(1,s13,sPip1,sPim,rRes2,mass2[i]);
857 //tmp1 = B[0]*B[1]*temp_PDF;
858 //amp_tmp1[0] += tmp1*pro0[0]*pro1[0];
859 //amp_tmp1[1] += tmp1*pro0[1]*pro1[1];
860 //
861 //temp_PDF=0;
862 //propagatorRBW(mass1sq,mass1[i],width1[i],s234,s24,sPim,rRes1,0,pro1);
863 //propagatorGS(mass2sq,mass2[i],width2[i],s23,sPip2,sPim,rRes2,pro0);
864 //
865 //for(int a=0; a<4; a++)
866 // for(int b=0;b<4;b++)
867 // temp_PDF += G[a][a]*G[b][b]*t1_234_1[a]*(-G[a][b]+p234[a]*p234[b]/s234)*t1_23[b];
868 //
869 //B[0] = barrier(1,sD,s234,sPip1,rD2,1.9683);
870 //B[1] = barrier(1,s23,sPip2,sPim,rRes2,mass2[i]);
871 //tmp2 = B[0]*B[1]*temp_PDF;
872 //amp_tmp2[0] += tmp2*pro0[0]*pro1[0];
873 //amp_tmp2[1] += tmp2*pro0[1]*pro1[1];
874
875 temp_PDF=0;
876 propagatorRBW(mass1sq,mass1[i],width1[i],s134,s34,sPip1,rRes1,1,pro1);
877 propagatorGS(mass2sq,mass2[i],width2[i],s34,sPim,sPi0,rRes2,pro0);
878
879 for(int a=0; a<4; a++)
880 temp_PDF += G[a][a]*t1_134_2[a]*t1_34[a];
881 //temp_PDF += G[a][a]*G[b][b]*t1_134_2[a]*(-G[a][b]+p134[a]*p134[b]/s134)*t1_34[b];
882
883 B[0] = barrier(1,s134,s34,sPip1,rRes1,mass1[i]);
884 B[1] = barrier(1,s34,sPim,sPi0,rRes2,mass2[i]);
885 tmp1 = B[0]*B[1]*temp_PDF;
886 amp_tmp1[0] += tmp1*pro0[0]*pro1[0]*-1.;
887 amp_tmp1[1] += tmp1*pro0[1]*pro1[1]*-1.;
888
889 temp_PDF=0;
890 propagatorRBW(mass1sq,mass1[i],width1[i],s234,s34,sPip2,rRes1,1,pro1);
891 propagatorGS(mass2sq,mass2[i],width2[i],s34,sPim,sPi0,rRes2,pro0);
892
893 for(int a=0; a<4; a++)
894 temp_PDF += G[a][a]*t1_234_1[a]*t1_34[a];
895 //temp_PDF += G[a][a]*G[b][b]*t1_234_1[a]*(-G[a][b]+p234[a]*p234[b]/s234)*t1_34[b];
896
897 B[0] = barrier(1,s234,s34,sPip2,rRes1,mass1[i]);
898 B[1] = barrier(1,s34,sPim,sPi0,rRes2,mass2[i]);
899 tmp2 = B[0]*B[1]*temp_PDF;
900 amp_tmp2[0] += tmp2*pro0[0]*pro1[0]*-1.;
901 amp_tmp2[1] += tmp2*pro0[1]*pro1[1]*-1.;
902 }
903
904 amp_tmp[0] = amp_tmp1[0]+amp_tmp2[0];
905 amp_tmp[1] = amp_tmp1[1]+amp_tmp2[1];
906 //printf("amp_tmp = %f , %f\n", amp_tmp[0], amp_tmp[1]);
907 //printf("cof = %f , %f\n", cof[0], cof[1]);
908 Com_Multi(amp_tmp,cof,amp_PDF);
909 PDF[0] += amp_PDF[0];
910 PDF[1] += amp_PDF[1];
911 //printf("PDF = %f , %f\n", PDF[0], PDF[1]);
912 }
913
914 double value = PDF[0]*PDF[0] + PDF[1]*PDF[1];
915 //value = value/Gen_eva[id];
916 value = value;
917 if(value <=0) value = 1e-20;
918 //if(value > 0) printf("max prob = %f\n",value);
919
920 setProb(value);
921
922 return ;
923
924}
925
926
927
928
929
930
931
932double EvtDsToPipPipPimPi0::a10gam_arr[3000]={
9330
934,0
935,0
936,0
937,0
938,0
939,0
940,0
941,0
942,0
943,0
944,0
945,0
946,0
947,0
948,0
949,0
950,0
951,0
952,0
953,0
954,0
955,0
956,0
957,0
958,0
959,0
960,0
961,0
962,0
963,0
964,0
965,0
966,0
967,0
968,0
969,0
970,0
971,0
972,0
973,0
974,0
975,0
976,0
977,0
978,0
979,0
980,0
981,0
982,0
983,0
984,0
985,0
986,0
987,0
988,0
989,0
990,0
991,0
992,0
993,0
994,0
995,0
996,0
997,0
998,0
999,0
1000,0
1001,0
1002,0
1003,0
1004,0
1005,0
1006,0
1007,0
1008,0
1009,0
1010,0
1011,0
1012,0
1013,0
1014,0
1015,0
1016,0
1017,0
1018,0
1019,0
1020,0
1021,0
1022,0
1023,0
1024,0
1025,0
1026,0
1027,0
1028,0
1029,0
1030,0
1031,0
1032,0
1033,0
1034,0
1035,0
1036,0
1037,0
1038,0
1039,0
1040,0
1041,0
1042,0
1043,0
1044,0
1045,0
1046,0
1047,0
1048,0
1049,0
1050,0
1051,0
1052,0
1053,0
1054,0
1055,0
1056,0
1057,0
1058,0
1059,0
1060,0
1061,0
1062,0
1063,0
1064,0
1065,0
1066,0
1067,0
1068,0
1069,0
1070,0
1071,0
1072,0
1073,0
1074,0
1075,0
1076,0
1077,0
1078,0
1079,0
1080,0
1081,0
1082,0
1083,0
1084,0
1085,0
1086,0
1087,0
1088,0
1089,0
1090,0
1091,0
1092,0
1093,0
1094,0
1095,0
1096,0
1097,0
1098,0
1099,0
1100,0
1101,0
1102,0
1103,0
1104,0
1105,0
1106,0
1107,0
1108,0
1109,0
1110,0
1111,0
1112,0
1113,0
1114,0
1115,0
1116,0
1117,0
1118,0
1119,0
1120,0
1121,0
1122,0
1123,0
1124,0
1125,0
1126,0
1127,0
1128,0
1129,0
1130,0
1131,0
1132,0
1133,0
1134,0
1135,0
1136,0
1137,0
1138,0
1139,0
1140,0
1141,0
1142,0
1143,0
1144,0
1145,0
1146,0
1147,0
1148,0
1149,0
1150,0
1151,0
1152,0
1153,0
1154,0
1155,0
1156,0
1157,0
1158,0
1159,0
1160,0
1161,0
1162,0
1163,0
1164,0
1165,0
1166,0
1167,0
1168,0
1169,0
1170,0
1171,0
1172,0
1173,0
1174,0
1175,0
1176,0
1177,0
1178,0
1179,0
1180,0
1181,0
1182,0
1183,0
1184,0
1185,0
1186,0
1187,0
1188,0
1189,0
1190,0
1191,0
1192,0
1193,0
1194,0
1195,0
1196,0
1197,0
1198,0
1199,0
1200,0
1201,0
1202,0
1203,0
1204,0
1205,0
1206,0
1207,0
1208,0
1209,0
1210,0
1211,0
1212,0
1213,0
1214,0
1215,0
1216,0
1217,0
1218,0
1219,0
1220,0
1221,0
1222,0
1223,0
1224,0
1225,0
1226,0
1227,0
1228,0
1229,0
1230,0
1231,0
1232,0
1233,0
1234,0
1235,0
1236,0
1237,0
1238,0
1239,0
1240,0
1241,0
1242,0
1243,0
1244,0
1245,0
1246,0
1247,0
1248,0
1249,0
1250,0
1251,0
1252,0
1253,0
1254,0
1255,0
1256,0
1257,0
1258,0
1259,0
1260,0
1261,0
1262,0
1263,0
1264,0
1265,0
1266,0
1267,0
1268,0
1269,0
1270,0
1271,0
1272,0
1273,0
1274,0
1275,0
1276,0
1277,0
1278,0
1279,0
1280,0
1281,0
1282,0
1283,0
1284,0
1285,0
1286,0
1287,0
1288,0
1289,0
1290,0
1291,0
1292,0
1293,0
1294,0
1295,0
1296,0
1297,0
1298,0
1299,0
1300,0
1301,0
1302,0
1303,0
1304,0
1305,0
1306,0
1307,0
1308,0
1309,0
1310,0
1311,0
1312,0
1313,0
1314,0
1315,0
1316,0
1317,0
1318,0
1319,0
1320,0
1321,0
1322,0
1323,0
1324,0
1325,0
1326,0
1327,0
1328,0
1329,0
1330,0
1331,0
1332,0
1333,0
1334,0
1335,0
1336,0
1337,0
1338,0
1339,0
1340,0
1341,0
1342,0
1343,0
1344,0
1345,0
1346,0
1347,0
1348,0
1349,0
1350,0
1351,0
1352,0
1353,0
1354,0
1355,0
1356,0
1357,0
1358,0
1359,5.10847e-07
1360,6.61977e-07
1361,8.41244e-07
1362,1.04861e-06
1363,1.29389e-06
1364,1.5632e-06
1365,1.881e-06
1366,2.22767e-06
1367,2.61443e-06
1368,3.05012e-06
1369,3.52012e-06
1370,4.05306e-06
1371,4.62497e-06
1372,5.24606e-06
1373,5.92555e-06
1374,6.68573e-06
1375,7.48203e-06
1376,8.30439e-06
1377,9.23754e-06
1378,1.02432e-05
1379,1.12879e-05
1380,1.24557e-05
1381,1.35884e-05
1382,1.49144e-05
1383,1.62366e-05
1384,1.76005e-05
1385,1.91314e-05
1386,2.0737e-05
1387,2.23708e-05
1388,2.41359e-05
1389,2.60124e-05
1390,2.79615e-05
1391,3.00748e-05
1392,3.20951e-05
1393,3.4398e-05
1394,3.67148e-05
1395,3.89229e-05
1396,4.16021e-05
1397,4.43174e-05
1398,4.7257e-05
1399,4.97924e-05
1400,5.28634e-05
1401,5.60937e-05
1402,5.90265e-05
1403,6.25579e-05
1404,6.59154e-05
1405,6.96911e-05
1406,7.34966e-05
1407,7.68773e-05
1408,8.10852e-05
1409,8.52922e-05
1410,8.93417e-05
1411,9.39913e-05
1412,9.85651e-05
1413,0.000103419
1414,0.000108147
1415,0.000112791
1416,0.00011776
1417,0.000123403
1418,0.000128711
1419,0.000134492
1420,0.000140132
1421,0.00014662
1422,0.000152356
1423,0.000158013
1424,0.000164076
1425,0.000170779
1426,0.000177999
1427,0.00018501
1428,0.000192253
1429,0.000199358
1430,0.000207475
1431,0.000214667
1432,0.000222587
1433,0.000231224
1434,0.000238213
1435,0.000247703
1436,0.000256315
1437,0.000265662
1438,0.000273847
1439,0.000283963
1440,0.000293846
1441,0.000303285
1442,0.000311502
1443,0.000322079
1444,0.000334168
1445,0.000343724
1446,0.000354855
1447,0.000366474
1448,0.000377177
1449,0.000388813
1450,0.00040168
1451,0.000413493
1452,0.000424643
1453,0.000437264
1454,0.000450301
1455,0.000463383
1456,0.000476909
1457,0.000489312
1458,0.00050142
1459,0.000517696
1460,0.000534774
1461,0.000547504
1462,0.000562717
1463,0.000579844
1464,0.000591737
1465,0.000609933
1466,0.00062713
1467,0.000642629
1468,0.000659459
1469,0.000675847
1470,0.000692308
1471,0.00071141
1472,0.000729495
1473,0.000744429
1474,0.000766867
1475,0.000785325
1476,0.000802652
1477,0.00082197
1478,0.000843976
1479,0.000863138
1480,0.000884998
1481,0.000906708
1482,0.000929898
1483,0.000947559
1484,0.000968399
1485,0.000996404
1486,0.00101473
1487,0.00104141
1488,0.00106038
1489,0.00108682
1490,0.00111357
1491,0.00113832
1492,0.00116272
1493,0.00118918
1494,0.00121012
1495,0.00123639
1496,0.00126683
1497,0.00129586
1498,0.00132024
1499,0.00135298
1500,0.00137635
1501,0.00140952
1502,0.00144097
1503,0.00146779
1504,0.00149647
1505,0.00152659
1506,0.00155896
1507,0.00158477
1508,0.00162065
1509,0.00165104
1510,0.00168977
1511,0.00171834
1512,0.00175372
1513,0.00178951
1514,0.0018246
1515,0.00185573
1516,0.00190013
1517,0.00193404
1518,0.00197543
1519,0.00201023
1520,0.0020423
1521,0.00208759
1522,0.00211846
1523,0.00216766
1524,0.00220693
1525,0.00224472
1526,0.00228621
1527,0.00232947
1528,0.00237215
1529,0.00242014
1530,0.00246763
1531,0.00250157
1532,0.00254408
1533,0.00259896
1534,0.00264598
1535,0.00268984
1536,0.00274089
1537,0.0027876
1538,0.00284181
1539,0.00289009
1540,0.00294365
1541,0.00298791
1542,0.00303216
1543,0.00308017
1544,0.00314244
1545,0.00319637
1546,0.00325229
1547,0.00330926
1548,0.00335885
1549,0.00342271
1550,0.00346827
1551,0.00353228
1552,0.00359511
1553,0.00364882
1554,0.00371527
1555,0.00378434
1556,0.00383893
1557,0.0039045
1558,0.00396506
1559,0.00403849
1560,0.00409357
1561,0.00416097
1562,0.00422666
1563,0.00429508
1564,0.00437824
1565,0.00442343
1566,0.00451049
1567,0.00456903
1568,0.00463377
1569,0.00471978
1570,0.00479821
1571,0.00486691
1572,0.00495857
1573,0.00501581
1574,0.00511116
1575,0.00517971
1576,0.00524602
1577,0.00533189
1578,0.00541212
1579,0.00550358
1580,0.00559244
1581,0.00566766
1582,0.00575527
1583,0.00586013
1584,0.005914
1585,0.00601374
1586,0.00610486
1587,0.00616553
1588,0.0062721
1589,0.00637741
1590,0.00645871
1591,0.00657222
1592,0.00665231
1593,0.00676846
1594,0.00685208
1595,0.00694391
1596,0.00704976
1597,0.00716695
1598,0.00722817
1599,0.0073624
1600,0.00745222
1601,0.00757624
1602,0.00765055
1603,0.00778373
1604,0.00788613
1605,0.00799569
1606,0.00810135
1607,0.00821636
1608,0.00832424
1609,0.00842107
1610,0.00855451
1611,0.0086506
1612,0.00879026
1613,0.00889932
1614,0.00902619
1615,0.00915076
1616,0.00924049
1617,0.00937964
1618,0.009491
1619,0.00962401
1620,0.00978053
1621,0.00988756
1622,0.0100305
1623,0.0101353
1624,0.0102867
1625,0.0104441
1626,0.0105558
1627,0.0107242
1628,0.0108266
1629,0.0109571
1630,0.0110897
1631,0.0112405
1632,0.0113613
1633,0.0114894
1634,0.0116503
1635,0.0118025
1636,0.0119203
1637,0.0121199
1638,0.0122381
1639,0.0124294
1640,0.0125371
1641,0.0126689
1642,0.0128119
1643,0.0130044
1644,0.0131396
1645,0.013325
1646,0.0134785
1647,0.0136555
1648,0.0137881
1649,0.0139677
1650,0.0141636
1651,0.0142625
1652,0.0144934
1653,0.0146466
1654,0.0148004
1655,0.0150156
1656,0.0151146
1657,0.0153291
1658,0.0154632
1659,0.015703
1660,0.0158963
1661,0.0160075
1662,0.0162392
1663,0.0164283
1664,0.0165921
1665,0.0167503
1666,0.0170055
1667,0.0171526
1668,0.0174251
1669,0.017558
1670,0.0177089
1671,0.017877
1672,0.0181561
1673,0.0183833
1674,0.018488
1675,0.0186674
1676,0.0188561
1677,0.0191999
1678,0.0193798
1679,0.0195454
1680,0.0197736
1681,0.0199596
1682,0.0201606
1683,0.0204273
1684,0.02067
1685,0.0207535
1686,0.0210466
1687,0.0212628
1688,0.0213996
1689,0.0217146
1690,0.0218709
1691,0.0221903
1692,0.0225007
1693,0.0226146
1694,0.0229532
1695,0.0230756
1696,0.0233669
1697,0.0235796
1698,0.0238283
1699,0.024111
1700,0.0243377
1701,0.0245557
1702,0.0249343
1703,0.0250247
1704,0.0252691
1705,0.0255234
1706,0.0257978
1707,0.026107
1708,0.0263546
1709,0.0266509
1710,0.0268604
1711,0.0271287
1712,0.0275146
1713,0.0278669
1714,0.0280236
1715,0.0282697
1716,0.0286572
1717,0.0289781
1718,0.0291602
1719,0.029455
1720,0.0298135
1721,0.0299364
1722,0.0305142
1723,0.0306114
1724,0.0309928
1725,0.0312366
1726,0.0316771
1727,0.031846
1728,0.0321499
1729,0.0324471
1730,0.0327513
1731,0.0330496
1732,0.0335854
1733,0.0339567
1734,0.0340194
1735,0.0344356
1736,0.0348833
1737,0.0351823
1738,0.035629
1739,0.0358131
1740,0.0362275
1741,0.0364262
1742,0.036772
1743,0.0372882
1744,0.0375451
1745,0.0379141
1746,0.0385863
1747,0.0387748
1748,0.0394341
1749,0.0397017
1750,0.0399935
1751,0.0402568
1752,0.0407253
1753,0.0411851
1754,0.0415508
1755,0.0418898
1756,0.0423181
1757,0.0429063
1758,0.0433494
1759,0.0437551
1760,0.0440054
1761,0.0447398
1762,0.0450333
1763,0.0454889
1764,0.0460102
1765,0.0462613
1766,0.0467099
1767,0.0470934
1768,0.0478124
1769,0.0483974
1770,0.0489387
1771,0.0491559
1772,0.0494822
1773,0.0501761
1774,0.0509158
1775,0.0514123
1776,0.0519352
1777,0.0522839
1778,0.0529311
1779,0.0533958
1780,0.0536812
1781,0.0545877
1782,0.0552246
1783,0.0558031
1784,0.0563144
1785,0.0567699
1786,0.0573668
1787,0.0582165
1788,0.0589321
1789,0.0591083
1790,0.059916
1791,0.0604653
1792,0.0611959
1793,0.0618002
1794,0.06273
1795,0.0631461
1796,0.0639283
1797,0.064459
1798,0.0651378
1799,0.06585
1800,0.0664203
1801,0.0674299
1802,0.0676451
1803,0.0689073
1804,0.0695797
1805,0.0705521
1806,0.070957
1807,0.0719581
1808,0.0723552
1809,0.0731694
1810,0.0738023
1811,0.0746405
1812,0.0756375
1813,0.0768234
1814,0.0774481
1815,0.0784771
1816,0.0786407
1817,0.0799796
1818,0.0808581
1819,0.0818104
1820,0.0827765
1821,0.083604
1822,0.0844563
1823,0.0857548
1824,0.086459
1825,0.0877274
1826,0.0878173
1827,0.0897353
1828,0.0899906
1829,0.0912604
1830,0.0921967
1831,0.093305
1832,0.094176
1833,0.0956489
1834,0.0962456
1835,0.0973797
1836,0.0978631
1837,0.0997736
1838,0.100602
1839,0.101641
1840,0.103275
1841,0.104253
1842,0.104856
1843,0.106542
1844,0.107601
1845,0.108657
1846,0.109335
1847,0.111208
1848,0.112411
1849,0.112519
1850,0.115234
1851,0.115867
1852,0.117292
1853,0.11831
1854,0.119598
1855,0.121313
1856,0.122108
1857,0.123445
1858,0.124304
1859,0.12644
1860,0.12744
1861,0.127965
1862,0.128983
1863,0.131425
1864,0.131891
1865,0.132958
1866,0.135072
1867,0.136372
1868,0.137859
1869,0.139469
1870,0.139895
1871,0.141239
1872,0.142312
1873,0.144194
1874,0.145061
1875,0.14797
1876,0.148421
1877,0.149148
1878,0.151263
1879,0.152417
1880,0.153852
1881,0.154373
1882,0.155729
1883,0.15729
1884,0.158637
1885,0.160362
1886,0.162451
1887,0.162713
1888,0.165762
1889,0.164869
1890,0.166916
1891,0.168694
1892,0.170652
1893,0.171882
1894,0.172525
1895,0.174622
1896,0.17715
1897,0.176453
1898,0.17879
1899,0.180351
1900,0.181534
1901,0.182003
1902,0.184618
1903,0.184791
1904,0.185818
1905,0.187754
1906,0.18921
1907,0.19104
1908,0.191221
1909,0.192176
1910,0.194501
1911,0.196039
1912,0.197454
1913,0.19889
1914,0.200492
1915,0.202103
1916,0.202012
1917,0.203932
1918,0.206371
1919,0.207377
1920,0.207238
1921,0.207741
1922,0.210689
1923,0.212003
1924,0.213854
1925,0.214701
1926,0.215246
1927,0.217599
1928,0.220232
1929,0.219966
1930,0.219449
1931,0.221044
1932,0.222879
1933,0.224735
1934,0.225572
1935,0.225291
1936,0.228986
1937,0.229253
1938,0.231717
1939,0.232168
1940,0.232043
1941,0.235528
1942,0.235789
1943,0.23709
1944,0.236679
1945,0.240255
1946,0.239063
1947,0.241123
1948,0.242824
1949,0.245739
1950,0.245582
1951,0.246838
1952,0.248066
1953,0.25014
1954,0.250404
1955,0.25013
1956,0.252638
1957,0.254109
1958,0.254707
1959,0.257212
1960,0.258817
1961,0.259441
1962,0.259976
1963,0.261818
1964,0.262706
1965,0.263087
1966,0.265652
1967,0.266072
1968,0.267088
1969,0.268155
1970,0.268635
1971,0.271209
1972,0.27077
1973,0.271821
1974,0.27319
1975,0.274764
1976,0.275037
1977,0.276091
1978,0.278019
1979,0.277085
1980,0.281357
1981,0.28164
1982,0.28116
1983,0.283233
1984,0.284078
1985,0.285207
1986,0.284973
1987,0.287791
1988,0.288662
1989,0.288786
1990,0.290711
1991,0.292558
1992,0.290947
1993,0.293373
1994,0.294503
1995,0.295762
1996,0.297076
1997,0.297324
1998,0.29907
1999,0.299182
2000,0.300861
2001,0.301929
2002,0.301915
2003,0.303147
2004,0.305086
2005,0.306043
2006,0.307639
2007,0.307859
2008,0.308634
2009,0.309576
2010,0.309839
2011,0.309856
2012,0.3122
2013,0.312669
2014,0.313976
2015,0.315044
2016,0.317234
2017,0.318186
2018,0.318364
2019,0.319883
2020,0.319694
2021,0.320653
2022,0.320546
2023,0.321426
2024,0.32232
2025,0.323426
2026,0.3238
2027,0.326689
2028,0.327256
2029,0.327845
2030,0.329071
2031,0.330415
2032,0.330224
2033,0.331656
2034,0.332284
2035,0.332526
2036,0.333567
2037,0.334426
2038,0.336384
2039,0.337273
2040,0.337489
2041,0.33854
2042,0.338955
2043,0.339142
2044,0.341662
2045,0.342302
2046,0.343481
2047,0.343899
2048,0.34338
2049,0.345925
2050,0.346405
2051,0.34765
2052,0.348484
2053,0.347483
2054,0.349126
2055,0.34929
2056,0.350908
2057,0.352362
2058,0.352584
2059,0.353779
2060,0.353871
2061,0.354532
2062,0.354155
2063,0.355416
2064,0.356756
2065,0.358541
2066,0.3579
2067,0.3592
2068,0.360009
2069,0.361232
2070,0.361467
2071,0.363298
2072,0.364186
2073,0.366374
2074,0.364917
2075,0.364327
2076,0.366702
2077,0.367645
2078,0.367933
2079,0.369071
2080,0.370738
2081,0.367952
2082,0.371982
2083,0.37071
2084,0.370969
2085,0.37351
2086,0.3725
2087,0.373201
2088,0.374315
2089,0.376762
2090,0.376296
2091,0.378849
2092,0.377269
2093,0.378157
2094,0.379463
2095,0.380044
2096,0.380239
2097,0.380691
2098,0.380712
2099,0.381581
2100,0.382821
2101,0.382769
2102,0.383597
2103,0.385241
2104,0.385573
2105,0.386094
2106,0.387948
2107,0.388176
2108,0.388829
2109,0.387262
2110,0.389284
2111,0.390265
2112,0.390807
2113,0.392308
2114,0.392241
2115,0.393039
2116,0.392672
2117,0.393517
2118,0.393756
2119,0.395187
2120,0.394738
2121,0.396592
2122,0.397667
2123,0.397501
2124,0.397686
2125,0.398778
2126,0.397756
2127,0.399491
2128,0.400856
2129,0.4008
2130,0.402307
2131,0.40258
2132,0.402486
2133,0.403773
2134,0.403612
2135,0.405249
2136,0.406007
2137,0.406182
2138,0.405619
2139,0.406478
2140,0.40722
2141,0.407959
2142,0.409399
2143,0.407149
2144,0.410364
2145,0.411445
2146,0.41226
2147,0.409334
2148,0.41238
2149,0.413182
2150,0.411058
2151,0.413268
2152,0.414802
2153,0.415712
2154,0.414306
2155,0.414652
2156,0.417959
2157,0.415588
2158,0.417994
2159,0.41751
2160,0.417698
2161,0.419487
2162,0.419246
2163,0.42
2164,0.421417
2165,0.421104
2166,0.422509
2167,0.421368
2168,0.421375
2169,0.422968
2170,0.422329
2171,0.424118
2172,0.424895
2173,0.426347
2174,0.424783
2175,0.425062
2176,0.427771
2177,0.4272
2178,0.427332
2179,0.428065
2180,0.428681
2181,0.428306
2182,0.428324
2183,0.42934
2184,0.431072
2185,0.43031
2186,0.431238
2187,0.431038
2188,0.431953
2189,0.432834
2190,0.433656
2191,0.435662
2192,0.43307
2193,0.434848
2194,0.43479
2195,0.435059
2196,0.434561
2197,0.43494
2198,0.437631
2199,0.436843
2200,0.436721
2201,0.438563
2202,0.436746
2203,0.438498
2204,0.439264
2205,0.440454
2206,0.439462
2207,0.441818
2208,0.441259
2209,0.440586
2210,0.441971
2211,0.442086
2212,0.441904
2213,0.443643
2214,0.442794
2215,0.445147
2216,0.44321
2217,0.444166
2218,0.444704
2219,0.445945
2220,0.444766
2221,0.44681
2222,0.446487
2223,0.446334
2224,0.448603
2225,0.448152
2226,0.448447
2227,0.449051
2228,0.449171
2229,0.450434
2230,0.449876
2231,0.451072
2232,0.450965
2233,0.450642
2234,0.451657
2235,0.451473
2236,0.451965
2237,0.454273
2238,0.452422
2239,0.453766
2240,0.455101
2241,0.455226
2242,0.452829
2243,0.455623
2244,0.456213
2245,0.453769
2246,0.456463
2247,0.458039
2248,0.45612
2249,0.456609
2250,0.458332
2251,0.458103
2252,0.457758
2253,0.457209
2254,0.460109
2255,0.458361
2256,0.459533
2257,0.459841
2258,0.461629
2259,0.462698
2260,0.462221
2261,0.460548
2262,0.461829
2263,0.463798
2264,0.463285
2265,0.465082
2266,0.465756
2267,0.463391
2268,0.465083
2269,0.462681
2270,0.464189
2271,0.466095
2272,0.46847
2273,0.466325
2274,0.466708
2275,0.465522
2276,0.468529
2277,0.468893
2278,0.467525
2279,0.469223
2280,0.46837
2281,0.468224
2282,0.470765
2283,0.468767
2284,0.469725
2285,0.47122
2286,0.469998
2287,0.471174
2288,0.471
2289,0.470084
2290,0.472925
2291,0.473875
2292,0.47366
2293,0.475595
2294,0.472112
2295,0.474616
2296,0.475303
2297,0.474338
2298,0.475903
2299,0.475051
2300,0.475786
2301,0.475884
2302,0.475342
2303,0.477809
2304,0.476628
2305,0.478056
2306,0.477943
2307,0.479718
2308,0.480277
2309,0.479325
2310,0.478487
2311,0.478945
2312,0.481972
2313,0.482123
2314,0.482074
2315,0.481763
2316,0.48447
2317,0.484614
2318,0.481945
2319,0.483892
2320,0.487872
2321,0.484796
2322,0.485473
2323,0.489195
2324,0.485496
2325,0.486415
2326,0.486113
2327,0.486683
2328,0.489275
2329,0.488448
2330,0.488583
2331,0.488278
2332,0.492439
2333,0.491196
2334,0.490612
2335,0.492102
2336,0.493087
2337,0.495125
2338,0.492948
2339,0.493526
2340,0.49445
2341,0.494222
2342,0.495337
2343,0.495459
2344,0.495539
2345,0.494363
2346,0.496046
2347,0.495832
2348,0.496912
2349,0.498911
2350,0.499184
2351,0.498805
2352,0.499241
2353,0.500922
2354,0.501565
2355,0.499584
2356,0.501578
2357,0.502322
2358,0.501752
2359,0.502496
2360,0.502788
2361,0.505379
2362,0.503861
2363,0.506021
2364,0.506488
2365,0.506952
2366,0.508643
2367,0.506462
2368,0.507231
2369,0.507387
2370,0.509251
2371,0.50949
2372,0.508067
2373,0.508938
2374,0.508887
2375,0.509544
2376,0.51172
2377,0.509405
2378,0.511026
2379,0.512919
2380,0.511507
2381,0.512892
2382,0.51042
2383,0.513466
2384,0.514605
2385,0.51404
2386,0.515337
2387,0.514276
2388,0.514417
2389,0.514953
2390,0.514552
2391,0.518115
2392,0.518749
2393,0.516392
2394,0.517777
2395,0.520258
2396,0.51817
2397,0.517867
2398,0.519882
2399,0.51881
2400,0.519493
2401,0.518606
2402,0.520654
2403,0.522035
2404,0.518634
2405,0.520473
2406,0.521678
2407,0.522549
2408,0.520957
2409,0.521824
2410,0.523186
2411,0.521642
2412,0.52482
2413,0.523135
2414,0.525027
2415,0.525555
2416,0.525987
2417,0.523265
2418,0.526671
2419,0.524975
2420,0.524577
2421,0.527919
2422,0.525689
2423,0.528478
2424,0.528185
2425,0.52857
2426,0.529382
2427,0.527601
2428,0.529476
2429,0.527973
2430,0.531799
2431,0.529369
2432,0.529891
2433,0.529205
2434,0.528191
2435,0.530105
2436,0.529386
2437,0.531998
2438,0.529379
2439,0.534389
2440,0.531392
2441,0.53429
2442,0.534107
2443,0.532299
2444,0.535032
2445,0.530939
2446,0.533253
2447,0.53627
2448,0.534767
2449,0.536838
2450,0.532772
2451,0.535239
2452,0.536046
2453,0.534192
2454,0.537879
2455,0.536891
2456,0.53856
2457,0.539537
2458,0.538163
2459,0.539638
2460,0.536308
2461,0.541708
2462,0.539197
2463,0.538073
2464,0.53991
2465,0.537786
2466,0.542624
2467,0.538922
2468,0.540254
2469,0.539483
2470,0.53884
2471,0.539868
2472,0.540382
2473,0.542084
2474,0.544621
2475,0.542604
2476,0.543015
2477,0.542192
2478,0.544019
2479,0.544245
2480,0.543841
2481,0.542288
2482,0.545539
2483,0.544668
2484,0.542409
2485,0.543304
2486,0.54556
2487,0.544934
2488,0.547257
2489,0.541864
2490,0.543311
2491,0.545769
2492,0.545772
2493,0.545936
2494,0.546728
2495,0.547074
2496,0.5482
2497,0.548422
2498,0.549682
2499,0.549333
2500,0.549396
2501,0.549981
2502,0.546348
2503,0.549856
2504,0.547825
2505,0.548367
2506,0.547422
2507,0.550924
2508,0.551186
2509,0.547844
2510,0.550414
2511,0.551172
2512,0.54655
2513,0.549614
2514,0.55211
2515,0.553047
2516,0.554943
2517,0.552244
2518,0.551954
2519,0.551853
2520,0.549702
2521,0.553917
2522,0.55485
2523,0.551653
2524,0.553974
2525,0.553223
2526,0.553047
2527,0.552811
2528,0.5525
2529,0.554389
2530,0.55657
2531,0.554466
2532,0.555791
2533,0.556522
2534,0.554262
2535,0.556422
2536,0.558315
2537,0.557895
2538,0.558123
2539,0.554199
2540,0.55749
2541,0.557732
2542,0.556954
2543,0.557359
2544,0.556586
2545,0.557766
2546,0.557517
2547,0.554655
2548,0.559342
2549,0.562112
2550,0.559345
2551,0.557585
2552,0.55697
2553,0.559785
2554,0.562197
2555,0.56008
2556,0.558988
2557,0.559874
2558,0.561188
2559,0.564173
2560,0.560372
2561,0.56124
2562,0.561108
2563,0.563741
2564,0.558707
2565,0.561086
2566,0.563327
2567,0.562036
2568,0.564004
2569,0.561138
2570,0.563496
2571,0.565594
2572,0.562067
2573,0.562929
2574,0.564893
2575,0.567015
2576,0.562481
2577,0.562019
2578,0.563653
2579,0.561878
2580,0.566143
2581,0.5639
2582,0.564638
2583,0.563107
2584,0.566875
2585,0.564398
2586,0.566127
2587,0.566607
2588,0.567252
2589,0.566528
2590,0.566628
2591,0.56831
2592,0.564859
2593,0.569815
2594,0.566654
2595,0.564939
2596,0.568844
2597,0.568502
2598,0.565886
2599,0.568556
2600,0.566762
2601,0.568847
2602,0.568534
2603,0.566041
2604,0.566397
2605,0.569942
2606,0.568914
2607,0.567241
2608,0.570659
2609,0.57171
2610,0.56819
2611,0.57075
2612,0.568438
2613,0.567604
2614,0.569452
2615,0.570547
2616,0.571437
2617,0.573075
2618,0.571073
2619,0.571652
2620,0.568959
2621,0.569775
2622,0.570844
2623,0.570371
2624,0.572211
2625,0.570553
2626,0.571842
2627,0.570598
2628,0.569651
2629,0.572605
2630,0.574228
2631,0.572185
2632,0.572898
2633,0.572446
2634,0.574438
2635,0.57708
2636,0.57026
2637,0.574982
2638,0.574789
2639,0.570822
2640,0.573091
2641,0.571822
2642,0.57595
2643,0.574218
2644,0.576605
2645,0.573796
2646,0.576091
2647,0.576826
2648,0.575359
2649,0.575667
2650,0.575334
2651,0.574895
2652,0.572222
2653,0.572793
2654,0.576912
2655,0.577936
2656,0.57358
2657,0.57554
2658,0.578902
2659,0.578798
2660,0.574994
2661,0.57923
2662,0.57836
2663,0.57877
2664,0.578834
2665,0.574414
2666,0.575213
2667,0.571854
2668,0.578509
2669,0.57951
2670,0.578496
2671,0.580245
2672,0.57465
2673,0.579298
2674,0.577356
2675,0.577914
2676,0.579471
2677,0.58146
2678,0.578826
2679,0.578709
2680,0.578496
2681,0.577299
2682,0.57901
2683,0.577267
2684,0.580718
2685,0.580655
2686,0.58123
2687,0.578204
2688,0.577744
2689,0.580787
2690,0.579833
2691,0.578184
2692,0.585243
2693,0.580038
2694,0.579339
2695,0.578362
2696,0.5854
2697,0.580457
2698,0.580576
2699,0.580914
2700,0.584755
2701,0.582103
2702,0.581257
2703,0.58125
2704,0.58311
2705,0.580425
2706,0.58365
2707,0.579959
2708,0.585634
2709,0.58076
2710,0.585159
2711,0.58226
2712,0.584189
2713,0.582075
2714,0.580576
2715,0.58043
2716,0.580405
2717,0.582512
2718,0.585068
2719,0.58504
2720,0.583435
2721,0.584687
2722,0.584157
2723,0.582721
2724,0.585834
2725,0.58427
2726,0.584793
2727,0.585872
2728,0.58757
2729,0.582234
2730,0.584769
2731,0.586066
2732,0.584464
2733,0.586815
2734,0.58484
2735,0.58566
2736,0.588422
2737,0.588056
2738,0.584557
2739,0.584627
2740,0.585923
2741,0.587775
2742,0.586758
2743,0.587405
2744,0.587869
2745,0.585873
2746,0.585099
2747,0.582264
2748,0.584199
2749,0.590424
2750,0.588416
2751,0.589516
2752,0.58845
2753,0.587837
2754,0.586989
2755,0.589098
2756,0.588176
2757,0.585571
2758,0.588833
2759,0.588869
2760,0.590503
2761,0.585261
2762,0.590453
2763,0.587992
2764,0.586263
2765,0.589369
2766,0.589266
2767,0.586583
2768,0.588911
2769,0.585802
2770,0.585045
2771,0.588152
2772,0.587742
2773,0.590043
2774,0.591548
2775,0.589419
2776,0.591158
2777,0.588734
2778,0.591094
2779,0.592836
2780,0.590413
2781,0.589916
2782,0.591878
2783,0.587068
2784,0.590096
2785,0.589218
2786,0.592117
2787,0.590057
2788,0.590293
2789,0.59132
2790,0.593525
2791,0.588977
2792,0.589158
2793,0.588977
2794,0.589911
2795,0.588572
2796,0.592564
2797,0.591573
2798,0.591105
2799,0.590711
2800,0.590143
2801,0.59068
2802,0.592997
2803,0.591367
2804,0.594414
2805,0.592478
2806,0.590993
2807,0.592463
2808,0.590829
2809,0.592986
2810,0.592069
2811,0.594267
2812,0.591337
2813,0.591242
2814,0.590092
2815,0.59447
2816,0.588796
2817,0.594219
2818,0.592487
2819,0.591164
2820,0.589996
2821,0.593095
2822,0.592464
2823,0.593464
2824,0.591231
2825,0.591955
2826,0.59051
2827,0.588754
2828,0.595189
2829,0.594417
2830,0.598559
2831,0.597429
2832,0.591446
2833,0.597268
2834,0.592084
2835,0.593667
2836,0.594351
2837,0.595908
2838,0.593284
2839,0.594093
2840,0.594685
2841,0.590196
2842,0.595528
2843,0.593947
2844,0.596839
2845,0.593559
2846,0.595751
2847,0.595056
2848,0.593025
2849,0.595207
2850,0.599865
2851,0.595802
2852,0.597093
2853,0.596852
2854,0.600151
2855,0.594165
2856,0.59669
2857,0.593693
2858,0.597007
2859,0.59448
2860,0.59692
2861,0.598447
2862,0.58967
2863,0.594291
2864,0.596814
2865,0.601778
2866,0.595606
2867,0.593575
2868,0.599196
2869,0.596763
2870,0.592731
2871,0.593649
2872,0.596815
2873,0.596612
2874,0.598389
2875,0.600288
2876,0.601971
2877,0.593411
2878,0.595628
2879,0.597471
2880,0.596617
2881,0.597107
2882,0.596928
2883,0.601419
2884,0.600712
2885,0.60062
2886,0.597027
2887,0.598212
2888,0.598517
2889,0.598585
2890,0.598457
2891,0.595458
2892,0.601006
2893,0.597436
2894,0.598442
2895,0.59978
2896,0.598042
2897,0.600751
2898,0.604877
2899,0.599682
2900,0.602193
2901,0.600038
2902,0.602133
2903,0.601971
2904,0.603675
2905,0.598525
2906,0.597034
2907,0.600436
2908,0.599307
2909,0.600178
2910,0.597561
2911,0.601128
2912,0.598386
2913,0.604022
2914,0.600333
2915,0.596953
2916,0.596153
2917,0.601371
2918,0.59905
2919,0.598775
2920,0.598234
2921,0.603056
2922,0.598461
2923,0.599035
2924,0.60247
2925,0.598961
2926,0.600638
2927,0.602174
2928,0.603546
2929,0.602759
2930,0.595493
2931,0.599566
2932,0.603676
2933,0.601448
2934,0.601889
2935,0.601923
2936,0.59879
2937,0.602152
2938,0.602868
2939,0.603142
2940,0.599779
2941,0.602996
2942,0.603121
2943,0.604151
2944,0.601858
2945,0.599859
2946,0.605559
2947,0.600411
2948,0.605023
2949,0.59481
2950,0.60354
2951,0.601028
2952,0.605539
2953,0.599386
2954,0.605815
2955,0.601308
2956,0.60254
2957,0.605936
2958,0.598884
2959,0.602046
2960,0.604051
2961,0.604064
2962,0.600417
2963,0.600486
2964,0.603358
2965,0.597937
2966,0.601253
2967,0.601612
2968,0.603221
2969,0.60466
2970,0.604129
2971,0.604015
2972,0.605971
2973,0.604743
2974,0.603351
2975,0.602063
2976,0.603581
2977,0.602386
2978,0.606447
2979,0.602734
2980,0.603914
2981,0.602157
2982,0.608614
2983,0.607969
2984,0.604008
2985,0.607868
2986,0.607401
2987,0.605498
2988,0.604538
2989,0.605259
2990,0.605669
2991,0.605429
2992,0.604667
2993,0.604459
2994,0.603132
2995,0.607398
2996,0.607448
2997,0.606535
2998,0.606449
2999,0.603193
3000,0.603195
3001,0.605879
3002,0.604217
3003,0.602632
3004,0.603078
3005,0.607168
3006,0.608578
3007,0.607912
3008,0.606881
3009,0.605766
3010,0.604638
3011,0.608067
3012,0.606821
3013,0.606599
3014,0.603447
3015,0.608147
3016,0.608554
3017,0.608716
3018,0.613558
3019,0.614474
3020,0.610203
3021,0.605614
3022,0.60737
3023,0.608189
3024,0.606754
3025,0.610616
3026,0.603147
3027,0.603729
3028,0.611811
3029,0.606968
3030,0.60742
3031,0.606781
3032,0.604482
3033,0.6097
3034,0.606153
3035,0.609647
3036,0.607856
3037,0.607672
3038,0.604834
3039,0.606836
3040,0.611446
3041,0.608694
3042,0.611505
3043,0.609354
3044,0.608521
3045,0.60922
3046,0.60297
3047,0.607662
3048,0.608237
3049,0.6063
3050,0.608534
3051,0.609033
3052,0.606615
3053,0.61312
3054,0.613844
3055,0.608296
3056,0.614646
3057,0.61078
3058,0.607609
3059,0.606543
3060,0.607679
3061,0.60831
3062,0.612263
3063,0.610579
3064,0.612233
3065,0.612682
3066,0.609124
3067,0.61026
3068,0.606912
3069,0.609823
3070,0.608506
3071,0.610746
3072,0.610349
3073,0.611911
3074,0.608253
3075,0.610799
3076,0.60775
3077,0.609174
3078,0.609446
3079,0.613876
3080,0.606614
3081,0.610548
3082,0.607299
3083,0.617604
3084,0.611957
3085,0.608182
3086,0.61352
3087,0.610114
3088,0.607217
3089,0.615439
3090,0.610774
3091,0.61353
3092,0.607945
3093,0.611571
3094,0.614412
3095,0.607529
3096,0.612346
3097,0.615903
3098,0.613468
3099,0.615199
3100,0.612175
3101,0.612662
3102,0.60956
3103,0.608311
3104,0.612912
3105,0.611678
3106,0.613576
3107,0.61002
3108,0.610685
3109,0.60852
3110,0.609646
3111,0.612874
3112,0.611657
3113,0.608268
3114,0.609351
3115,0.613332
3116,0.610005
3117,0.616228
3118,0.609752
3119,0.610077
3120,0.610468
3121,0.612632
3122,0.611267
3123,0.614137
3124,0.617096
3125,0.61805
3126,0.615931
3127,0.610193
3128,0.606448
3129,0.610412
3130,0.614418
3131,0.613298
3132,0.616366
3133,0.610521
3134,0.609473
3135,0.617223
3136,0.614392
3137,0.609531
3138,0.611503
3139,0.619216
3140,0.612787
3141,0.616896
3142,0.613963
3143,0.613148
3144,0.611531
3145,0.612897
3146,0.613263
3147,0.611573
3148,0.616155
3149,0.615285
3150,0.61387
3151,0.613969
3152,0.61706
3153,0.615567
3154,0.614016
3155,0.611149
3156,0.611988
3157,0.613614
3158,0.613824
3159,0.609791
3160,0.613618
3161,0.609961
3162,0.620689
3163,0.615341
3164,0.610682
3165,0.614414
3166,0.61655
3167,0.615938
3168,0.615047
3169,0.61412
3170,0.614024
3171,0.611481
3172,0.613393
3173,0.616073
3174,0.61456
3175,0.614876
3176,0.617353
3177,0.617323
3178,0.614166
3179,0.618364
3180,0.615672
3181,0.615384
3182,0.616775
3183,0.614972
3184,0.615856
3185,0.612369
3186,0.613385
3187,0.615188
3188,0.621464
3189,0.621196
3190,0.617695
3191,0.617198
3192,0.615838
3193,0.620844
3194,0.614442
3195,0.616499
3196,0.611469
3197,0.61797
3198,0.613381
3199,0.612686
3200,0.617832
3201,0.616912
3202,0.615726
3203,0.613965
3204,0.617923
3205,0.616926
3206,0.613962
3207,0.616878
3208,0.616216
3209,0.615167
3210,0.615741
3211,0.621837
3212,0.622816
3213,0.618735
3214,0.614708
3215,0.616259
3216,0.615296
3217,0.617336
3218,0.613412
3219,0.617638
3220,0.619523
3221,0.617039
3222,0.614946
3223,0.620744
3224,0.618552
3225,0.616013
3226,0.617093
3227,0.616921
3228,0.618917
3229,0.61971
3230,0.616581
3231,0.614899
3232,0.619123
3233,0.613554
3234,0.6208
3235,0.618294
3236,0.619059
3237,0.620993
3238,0.619525
3239,0.621942
3240,0.614049
3241,0.616872
3242,0.617373
3243,0.621464
3244,0.615137
3245,0.621412
3246,0.617671
3247,0.619178
3248,0.617229
3249,0.620786
3250,0.618965
3251,0.619412
3252,0.624408
3253,0.620756
3254,0.622836
3255,0.619064
3256,0.622596
3257,0.618601
3258,0.612976
3259,0.614824
3260,0.621118
3261,0.621563
3262,0.621251
3263,0.619092
3264,0.62368
3265,0.620553
3266,0.625767
3267,0.621926
3268,0.618632
3269,0.616635
3270,0.616295
3271,0.625815
3272,0.620596
3273,0.616814
3274,0.619908
3275,0.620009
3276,0.618273
3277,0.618491
3278,0.621039
3279,0.619414
3280,0.618904
3281,0.617644
3282,0.622816
3283,0.620315
3284,0.622365
3285,0.62554
3286,0.62318
3287,0.621576
3288,0.622506
3289,0.622526
3290,0.621412
3291,0.617583
3292,0.627276
3293,0.619036
3294,0.618457
3295,0.620758
3296,0.618635
3297,0.618848
3298,0.616887
3299,0.625413
3300,0.622702
3301,0.626354
3302,0.619957
3303,0.621203
3304,0.622965
3305,0.620388
3306,0.624028
3307,0.619793
3308,0.619449
3309,0.619894
3310,0.622475
3311,0.625728
3312,0.620663
3313,0.625763
3314,0.620276
3315,0.6219
3316,0.622332
3317,0.62029
3318,0.623742
3319,0.620927
3320,0.616047
3321,0.62078
3322,0.619601
3323,0.625109
3324,0.623586
3325,0.618206
3326,0.624081
3327,0.622982
3328,0.617744
3329,0.622281
3330,0.621567
3331,0.62139
3332,0.62237
3333,0.625921
3334,0.625531
3335,0.625433
3336,0.621814
3337,0.624017
3338,0.622389
3339,0.623663
3340,0.623262
3341,0.627376
3342,0.623763
3343,0.621025
3344,0.625044
3345,0.622423
3346,0.619603
3347,0.622101
3348,0.61922
3349,0.626379
3350,0.622221
3351,0.625152
3352,0.6239
3353,0.626617
3354,0.617879
3355,0.626182
3356,0.622782
3357,0.624102
3358,0.625998
3359,0.629432
3360,0.626768
3361,0.624665
3362,0.623231
3363,0.62919
3364,0.620829
3365,0.620281
3366,0.624708
3367,0.626512
3368,0.622157
3369,0.626986
3370,0.622909
3371,0.619203
3372,0.623121
3373,0.631429
3374,0.623336
3375,0.624622
3376,0.623138
3377,0.62677
3378,0.623556
3379,0.629962
3380,0.626398
3381,0.622657
3382,0.62277
3383,0.626691
3384,0.622954
3385,0.623641
3386,0.621758
3387,0.625917
3388,0.626179
3389,0.624352
3390,0.629404
3391,0.624387
3392,0.623461
3393,0.627372
3394,0.624249
3395,0.62444
3396,0.625099
3397,0.624841
3398,0.62775
3399,0.625561
3400,0.628993
3401,0.624342
3402,0.621275
3403,0.629437
3404,0.632163
3405,0.623176
3406,0.624028
3407,0.624016
3408,0.624913
3409,0.62146
3410,0.625645
3411,0.632107
3412,0.62442
3413,0.626444
3414,0.627954
3415,0.628726
3416,0.629863
3417,0.626985
3418,0.628023
3419,0.62284
3420,0.634704
3421,0.622555
3422,0.628842
3423,0.633196
3424,0.624833
3425,0.630645
3426,0.62962
3427,0.627391
3428,0.625646
3429,0.624605
3430,0.623669
3431,0.631903
3432,0.627895
3433,0.631525
3434,0.623852
3435,0.635089
3436,0.632845
3437,0.626416
3438,0.630422
3439,0.626313
3440,0.628298
3441,0.632386
3442,0.623754
3443,0.627711
3444,0.627182
3445,0.631031
3446,0.632696
3447,0.6265
3448,0.626044
3449,0.631712
3450,0.631146
3451,0.630076
3452,0.6246
3453,0.626339
3454,0.626073
3455,0.62547
3456,0.624333
3457,0.627483
3458,0.629035
3459,0.623781
3460,0.626105
3461,0.63322
3462,0.628441
3463,0.625099
3464,0.631457
3465,0.631161
3466,0.628228
3467,0.629891
3468,0.62785
3469,0.627987
3470,0.627642
3471,0.63232
3472,0.629908
3473,0.633618
3474,0.631015
3475,0.631903
3476,0.627583
3477,0.627801
3478,0.631097
3479,0.63238
3480,0.634586
3481,0.627661
3482,0.625959
3483,0.629329
3484,0.625851
3485,0.632545
3486,0.631533
3487,0.628986
3488,0.631847
3489,0.632295
3490,0.641355
3491,0.627727
3492,0.634586
3493,0.625373
3494,0.627816
3495,0.633512
3496,0.63269
3497,0.627166
3498,0.632661
3499,0.628999
3500,0.631664
3501,0.633816
3502,0.631742
3503,0.633348
3504,0.630743
3505,0.63382
3506,0.635101
3507,0.632981
3508,0.6346
3509,0.632103
3510,0.632281
3511,0.63134
3512,0.634626
3513,0.630111
3514,0.631351
3515,0.630224
3516,0.630546
3517,0.632019
3518,0.632453
3519,0.631172
3520,0.632314
3521,0.631817
3522,0.62508
3523,0.632514
3524,0.634765
3525,0.629544
3526,0.635109
3527,0.630736
3528,0.63048
3529,0.635906
3530,0.635783
3531,0.624545
3532,0.630134
3533,0.628176
3534,0.630117
3535,0.63592
3536,0.634577
3537,0.632769
3538,0.634533
3539,0.63141
3540,0.628919
3541,0.636203
3542,0.638389
3543,0.63231
3544,0.632995
3545,0.629657
3546,0.630346
3547,0.632677
3548,0.634667
3549,0.634583
3550,0.629406
3551,0.63471
3552,0.635601
3553,0.631947
3554,0.635692
3555,0.638141
3556,0.630747
3557,0.631889
3558,0.635176
3559,0.633654
3560,0.62976
3561,0.636197
3562,0.633478
3563,0.634548
3564,0.63132
3565,0.634575
3566,0.633873
3567,0.632214
3568,0.63282
3569,0.631591
3570,0.629881
3571,0.635444
3572,0.636161
3573,0.634927
3574,0.637786
3575,0.636701
3576,0.630247
3577,0.6411
3578,0.640469
3579,0.633986
3580,0.632693
3581,0.633078
3582,0.631706
3583,0.628018
3584,0.640413
3585,0.63462
3586,0.631846
3587,0.639895
3588,0.634679
3589,0.629001
3590,0.641085
3591,0.636796
3592,0.636365
3593,0.634147
3594,0.639004
3595,0.633508
3596,0.636331
3597,0.642073
3598,0.634742
3599,0.635549
3600,0.639292
3601,0.633249
3602,0.629993
3603,0.639224
3604,0.635845
3605,0.635741
3606,0.631689
3607,0.633921
3608,0.636043
3609,0.636046
3610,0.634273
3611,0.633504
3612,0.635015
3613,0.63648
3614,0.634258
3615,0.631794
3616,0.637351
3617,0.637003
3618,0.635062
3619,0.641799
3620,0.630846
3621,0.636755
3622,0.637128
3623,0.636238
3624,0.635672
3625,0.642322
3626,0.63981
3627,0.638706
3628,0.640004
3629,0.640292
3630,0.640726
3631,0.640107
3632,0.637919
3633,0.63783
3634,0.642114
3635,0.632441
3636,0.637189
3637,0.635049
3638,0.633693
3639,0.638689
3640,0.632568
3641,0.646298
3642,0.631719
3643,0.633435
3644,0.633835
3645,0.636585
3646,0.637045
3647,0.637189
3648,0.634202
3649,0.639623
3650,0.640147
3651,0.634653
3652,0.640554
3653,0.639917
3654,0.64078
3655,0.638127
3656,0.63798
3657,0.637607
3658,0.639998
3659,0.63473
3660,0.643037
3661,0.636686
3662,0.63415
3663,0.645426
3664,0.639442
3665,0.639877
3666,0.636225
3667,0.636739
3668,0.628628
3669,0.640501
3670,0.646981
3671,0.637553
3672,0.63904
3673,0.639558
3674,0.632821
3675,0.641581
3676,0.634517
3677,0.639993
3678,0.640098
3679,0.640046
3680,0.634301
3681,0.63906
3682,0.644402
3683,0.643639
3684,0.638688
3685,0.643383
3686,0.63776
3687,0.639042
3688,0.639758
3689,0.640635
3690,0.641742
3691,0.641212
3692,0.639423
3693,0.636773
3694,0.635784
3695,0.639141
3696,0.637001
3697,0.638086
3698,0.639779
3699,0.633228
3700,0.641054
3701,0.636265
3702,0.642418
3703,0.644187
3704,0.640025
3705,0.637431
3706,0.64234
3707,0.639466
3708,0.639929
3709,0.643658
3710,0.645672
3711,0.63875
3712,0.638856
3713,0.643978
3714,0.647565
3715,0.647039
3716,0.642832
3717,0.640113
3718,0.642341
3719,0.644558
3720,0.637816
3721,0.638742
3722,0.641432
3723,0.644949
3724,0.641804
3725,0.641167
3726,0.646557
3727,0.64298
3728,0.642473
3729,0.641937
3730,0.644673
3731,0.638359
3732,0.639813
3733,0.639308
3734,0.636239
3735,0.640583
3736,0.641604
3737,0.647463
3738,0.643723
3739,0.638941
3740,0.640267
3741,0.641271
3742,0.648004
3743,0.647338
3744,0.640559
3745,0.640159
3746,0.639687
3747,0.647811
3748,0.644427
3749,0.639744
3750,0.645816
3751,0.640913
3752,0.647267
3753,0.644434
3754,0.645361
3755,0.640173
3756,0.647387
3757,0.643858
3758,0.638137
3759,0.642417
3760,0.647948
3761,0.638894
3762,0.646922
3763,0.650392
3764,0.643519
3765,0.645082
3766,0.63969
3767,0.643681
3768,0.647513
3769,0.641228
3770,0.642628
3771,0.644522
3772,0.645842
3773,0.647846
3774,0.64716
3775,0.643352
3776,0.651159
3777,0.644077
3778,0.650025
3779,0.644696
3780,0.643421
3781,0.641637
3782,0.647484
3783,0.640235
3784,0.643457
3785,0.642693
3786,0.645508
3787,0.642223
3788,0.639686
3789,0.646631
3790,0.641975
3791,0.643323
3792,0.647343
3793,0.64402
3794,0.648346
3795,0.649236
3796,0.647768
3797,0.645609
3798,0.641971
3799,0.641319
3800,0.641983
3801,0.644494
3802,0.644607
3803,0.649457
3804,0.643873
3805,0.647442
3806,0.646043
3807,0.64612
3808,0.648729
3809,0.646065
3810,0.652991
3811,0.645837
3812,0.650514
3813,0.651147
3814,0.644631
3815,0.645813
3816,0.643999
3817,0.647581
3818,0.646826
3819,0.64706
3820,0.650602
3821,0.654695
3822,0.648685
3823,0.651021
3824,0.651924
3825,0.649291
3826,0.652016
3827,0.654856
3828,0.644839
3829,0.644466
3830,0.648472
3831,0.645495
3832,0.648701
3833,0.644017
3834,0.648088
3835,0.647046
3836,0.647111
3837,0.654831
3838,0.649127
3839,0.643345
3840,0.643242
3841,0.651873
3842,0.649742
3843,0.648395
3844,0.643895
3845,0.642581
3846,0.643208
3847,0.651363
3848,0.653665
3849,0.65346
3850,0.654502
3851,0.64849
3852,0.643907
3853,0.653994
3854,0.648345
3855,0.646242
3856,0.647899
3857,0.651866
3858,0.648614
3859,0.654628
3860,0.649
3861,0.650335
3862,0.643417
3863,0.645178
3864,0.653398
3865,0.65154
3866,0.653001
3867,0.653917
3868,0.645493
3869,0.643492
3870,0.653624
3871,0.650601
3872,0.654645
3873,0.650887
3874,0.653012
3875,0.651437
3876,0.652385
3877,0.655021
3878,0.645354
3879,0.643753
3880,0.649658
3881,0.645741
3882,0.648658
3883,0.649114
3884,0.643689
3885,0.651045
3886,0.65228
3887,0.653334
3888,0.65179
3889,0.654058
3890,0.647741
3891,0.650629
3892,0.654112
3893,0.654751
3894,0.658177
3895,0.647574
3896,0.653494
3897,0.654006
3898,0.656714
3899,0.653971
3900,0.652711
3901,0.648534
3902,0.65007
3903,0.655662
3904,0.649423
3905,0.65662
3906,0.66087
3907,0.650238
3908,0.653608
3909,0.655941
3910,0.650294
3911,0.650623
3912,0.652413
3913,0.64931
3914,0.649623
3915,0.653893
3916,0.652993
3917,0.650144
3918,0.652416
3919,0.649424
3920,0.652945
3921,0.649755
3922,0.655409
3923,0.654311
3924,0.65301
3925,0.648886
3926,0.651531
3927,0.651158
3928,0.655668
3929,0.657696
3930,0.654918
3931,0.653121
3932,0.649591
3933};
3934
3935double EvtDsToPipPipPimPi0::a1pgam_arr[3000]={
39360
3937,0
3938,0
3939,0
3940,0
3941,0
3942,0
3943,0
3944,0
3945,0
3946,0
3947,0
3948,0
3949,0
3950,0
3951,0
3952,0
3953,0
3954,0
3955,0
3956,0
3957,0
3958,0
3959,0
3960,0
3961,0
3962,0
3963,0
3964,0
3965,0
3966,0
3967,0
3968,0
3969,0
3970,0
3971,0
3972,0
3973,0
3974,0
3975,0
3976,0
3977,0
3978,0
3979,0
3980,0
3981,0
3982,0
3983,0
3984,0
3985,0
3986,0
3987,0
3988,0
3989,0
3990,0
3991,0
3992,0
3993,0
3994,0
3995,0
3996,0
3997,0
3998,0
3999,0
4000,0
4001,0
4002,0
4003,0
4004,0
4005,0
4006,0
4007,0
4008,0
4009,0
4010,0
4011,0
4012,0
4013,0
4014,0
4015,0
4016,0
4017,0
4018,0
4019,0
4020,0
4021,0
4022,0
4023,0
4024,0
4025,0
4026,0
4027,0
4028,0
4029,0
4030,0
4031,0
4032,0
4033,0
4034,0
4035,0
4036,0
4037,0
4038,0
4039,0
4040,0
4041,0
4042,0
4043,0
4044,0
4045,0
4046,0
4047,0
4048,0
4049,0
4050,0
4051,0
4052,0
4053,0
4054,0
4055,0
4056,0
4057,0
4058,0
4059,0
4060,0
4061,0
4062,0
4063,0
4064,0
4065,0
4066,0
4067,0
4068,0
4069,0
4070,0
4071,0
4072,0
4073,0
4074,0
4075,0
4076,0
4077,0
4078,0
4079,0
4080,0
4081,0
4082,0
4083,0
4084,0
4085,0
4086,0
4087,0
4088,0
4089,0
4090,0
4091,0
4092,0
4093,0
4094,0
4095,0
4096,0
4097,0
4098,0
4099,0
4100,0
4101,0
4102,0
4103,0
4104,0
4105,0
4106,0
4107,0
4108,0
4109,0
4110,0
4111,0
4112,0
4113,0
4114,0
4115,0
4116,0
4117,0
4118,0
4119,0
4120,0
4121,0
4122,0
4123,0
4124,0
4125,0
4126,0
4127,0
4128,0
4129,0
4130,0
4131,0
4132,0
4133,0
4134,0
4135,0
4136,0
4137,0
4138,0
4139,0
4140,0
4141,0
4142,0
4143,0
4144,0
4145,0
4146,0
4147,0
4148,0
4149,0
4150,0
4151,0
4152,0
4153,0
4154,0
4155,0
4156,0
4157,0
4158,0
4159,0
4160,0
4161,0
4162,0
4163,0
4164,0
4165,0
4166,0
4167,0
4168,0
4169,0
4170,0
4171,0
4172,0
4173,0
4174,0
4175,0
4176,0
4177,0
4178,0
4179,0
4180,0
4181,0
4182,0
4183,0
4184,0
4185,0
4186,0
4187,0
4188,0
4189,0
4190,0
4191,0
4192,0
4193,0
4194,0
4195,0
4196,0
4197,0
4198,0
4199,0
4200,0
4201,0
4202,0
4203,0
4204,0
4205,0
4206,0
4207,0
4208,0
4209,0
4210,0
4211,0
4212,0
4213,0
4214,0
4215,0
4216,0
4217,0
4218,0
4219,0
4220,0
4221,0
4222,0
4223,0
4224,0
4225,0
4226,0
4227,0
4228,0
4229,0
4230,0
4231,0
4232,0
4233,0
4234,0
4235,0
4236,0
4237,0
4238,0
4239,0
4240,0
4241,0
4242,0
4243,0
4244,0
4245,0
4246,0
4247,0
4248,0
4249,0
4250,0
4251,0
4252,0
4253,0
4254,0
4255,0
4256,0
4257,0
4258,0
4259,0
4260,0
4261,0
4262,0
4263,0
4264,0
4265,0
4266,0
4267,0
4268,0
4269,0
4270,0
4271,0
4272,0
4273,0
4274,0
4275,0
4276,0
4277,0
4278,0
4279,0
4280,0
4281,0
4282,0
4283,0
4284,0
4285,0
4286,0
4287,0
4288,0
4289,0
4290,0
4291,0
4292,0
4293,0
4294,0
4295,0
4296,0
4297,0
4298,0
4299,0
4300,0
4301,0
4302,0
4303,0
4304,0
4305,0
4306,0
4307,0
4308,0
4309,0
4310,0
4311,0
4312,0
4313,0
4314,0
4315,0
4316,0
4317,0
4318,0
4319,0
4320,0
4321,0
4322,0
4323,0
4324,0
4325,0
4326,0
4327,0
4328,0
4329,0
4330,0
4331,0
4332,0
4333,0
4334,0
4335,0
4336,0
4337,0
4338,0
4339,0
4340,0
4341,0
4342,0
4343,0
4344,0
4345,0
4346,0
4347,0
4348,0
4349,0
4350,0
4351,0
4352,0
4353,0
4354,0
4355,0
4356,0
4357,0
4358,0
4359,0
4360,0
4361,0
4362,0
4363,0
4364,0
4365,0
4366,0
4367,1.16255e-06
4368,1.43224e-06
4369,1.73208e-06
4370,2.07745e-06
4371,2.46811e-06
4372,2.8958e-06
4373,3.37455e-06
4374,3.89844e-06
4375,4.49816e-06
4376,5.13849e-06
4377,5.80413e-06
4378,6.5552e-06
4379,7.36488e-06
4380,8.22848e-06
4381,9.15607e-06
4382,1.02051e-05
4383,1.12368e-05
4384,1.24356e-05
4385,1.36177e-05
4386,1.49035e-05
4387,1.63023e-05
4388,1.77732e-05
4389,1.93243e-05
4390,2.1039e-05
4391,2.2603e-05
4392,2.4414e-05
4393,2.63433e-05
4394,2.84536e-05
4395,3.04303e-05
4396,3.27502e-05
4397,3.48853e-05
4398,3.73187e-05
4399,3.97525e-05
4400,4.24245e-05
4401,4.51341e-05
4402,4.80726e-05
4403,5.08268e-05
4404,5.406e-05
4405,5.73412e-05
4406,6.06034e-05
4407,6.37824e-05
4408,6.74847e-05
4409,7.10286e-05
4410,7.49234e-05
4411,7.90774e-05
4412,8.2991e-05
4413,8.74405e-05
4414,9.17103e-05
4415,9.59903e-05
4416,0.000100626
4417,0.000105334
4418,0.00011065
4419,0.000115699
4420,0.000120608
4421,0.000126206
4422,0.000131872
4423,0.000137576
4424,0.000143318
4425,0.000149702
4426,0.000155464
4427,0.000161623
4428,0.000168204
4429,0.000175364
4430,0.000182611
4431,0.000189817
4432,0.000196914
4433,0.000203524
4434,0.000211917
4435,0.000219516
4436,0.000227111
4437,0.000236204
4438,0.000243656
4439,0.000253163
4440,0.000261804
4441,0.000270058
4442,0.000279425
4443,0.000289415
4444,0.000299098
4445,0.000308215
4446,0.000317935
4447,0.000328688
4448,0.000339008
4449,0.000349371
4450,0.000360037
4451,0.00037213
4452,0.000383524
4453,0.00039489
4454,0.00040754
4455,0.000419597
4456,0.00043143
4457,0.000444358
4458,0.000456249
4459,0.000469994
4460,0.00048248
4461,0.000496576
4462,0.000509287
4463,0.000523465
4464,0.000537663
4465,0.000552615
4466,0.000565672
4467,0.000581505
4468,0.000598695
4469,0.000614738
4470,0.000630295
4471,0.000646497
4472,0.000661822
4473,0.000680691
4474,0.000697636
4475,0.000710653
4476,0.000731611
4477,0.000747991
4478,0.000768662
4479,0.000786658
4480,0.000804289
4481,0.000825827
4482,0.000842601
4483,0.000863005
4484,0.000884184
4485,0.000901145
4486,0.00092432
4487,0.000947943
4488,0.000966836
4489,0.000990153
4490,0.00101148
4491,0.00103307
4492,0.00105758
4493,0.00107941
4494,0.00110262
4495,0.00113216
4496,0.00114936
4497,0.00117555
4498,0.00119969
4499,0.00122795
4500,0.00124984
4501,0.00127791
4502,0.00130504
4503,0.00132556
4504,0.00135835
4505,0.0013913
4506,0.00141671
4507,0.00144362
4508,0.00147138
4509,0.00150494
4510,0.00153143
4511,0.00156205
4512,0.00159195
4513,0.0016225
4514,0.00165337
4515,0.00168743
4516,0.00172087
4517,0.00175389
4518,0.00179009
4519,0.00181904
4520,0.00185995
4521,0.00189109
4522,0.00192425
4523,0.00196063
4524,0.00199956
4525,0.00203658
4526,0.00207548
4527,0.00211007
4528,0.00214713
4529,0.00218205
4530,0.00222067
4531,0.00226647
4532,0.00230394
4533,0.00234372
4534,0.00239147
4535,0.00243497
4536,0.00248179
4537,0.00252284
4538,0.00258151
4539,0.00261317
4540,0.00265228
4541,0.0027003
4542,0.00275796
4543,0.00279158
4544,0.00284893
4545,0.0029003
4546,0.00294198
4547,0.00299252
4548,0.00304708
4549,0.00310524
4550,0.00314952
4551,0.0031982
4552,0.00326313
4553,0.00329882
4554,0.00336018
4555,0.00341997
4556,0.00346604
4557,0.00353923
4558,0.00359243
4559,0.00364612
4560,0.00369917
4561,0.00376501
4562,0.00382314
4563,0.00388662
4564,0.00395138
4565,0.00401275
4566,0.00408231
4567,0.00415988
4568,0.00421809
4569,0.00427482
4570,0.00436948
4571,0.004428
4572,0.00449458
4573,0.00455717
4574,0.00462989
4575,0.0047196
4576,0.00478372
4577,0.00487112
4578,0.00494045
4579,0.00500795
4580,0.00509995
4581,0.00517203
4582,0.00525552
4583,0.00533371
4584,0.00542471
4585,0.00551365
4586,0.00558177
4587,0.00568082
4588,0.00577522
4589,0.00587262
4590,0.00594447
4591,0.00605402
4592,0.00611928
4593,0.00620531
4594,0.00632442
4595,0.00642466
4596,0.00651901
4597,0.00663929
4598,0.00672864
4599,0.00682712
4600,0.00693815
4601,0.00703404
4602,0.00715259
4603,0.00725306
4604,0.00736991
4605,0.00743309
4606,0.00757263
4607,0.00769897
4608,0.00782515
4609,0.00793906
4610,0.00804642
4611,0.00813755
4612,0.00828158
4613,0.00842992
4614,0.00853162
4615,0.0086354
4616,0.00875561
4617,0.00891002
4618,0.0090398
4619,0.00921414
4620,0.0093221
4621,0.00942081
4622,0.00958837
4623,0.0097279
4624,0.0098762
4625,0.0100076
4626,0.0101537
4627,0.0102833
4628,0.0104462
4629,0.0106011
4630,0.0107291
4631,0.0109382
4632,0.011029
4633,0.0112307
4634,0.0113983
4635,0.0115593
4636,0.011673
4637,0.0119066
4638,0.0120356
4639,0.0122478
4640,0.0124057
4641,0.0125293
4642,0.0126917
4643,0.0129166
4644,0.0130909
4645,0.0132769
4646,0.0134728
4647,0.0136977
4648,0.0138256
4649,0.0140528
4650,0.0142311
4651,0.0144146
4652,0.0146287
4653,0.0148169
4654,0.0150073
4655,0.0151786
4656,0.0154452
4657,0.0156354
4658,0.0158328
4659,0.016026
4660,0.0162738
4661,0.0164946
4662,0.01674
4663,0.0169668
4664,0.0171655
4665,0.017386
4666,0.0176114
4667,0.0178374
4668,0.0180972
4669,0.0182908
4670,0.0185786
4671,0.0188305
4672,0.019112
4673,0.0192812
4674,0.0195456
4675,0.0198758
4676,0.0200246
4677,0.0203729
4678,0.0205916
4679,0.0209285
4680,0.0211636
4681,0.0213849
4682,0.0217342
4683,0.0219524
4684,0.0222627
4685,0.0225478
4686,0.0228542
4687,0.0231567
4688,0.0233594
4689,0.0237955
4690,0.0240128
4691,0.024268
4692,0.0246971
4693,0.0250577
4694,0.0252016
4695,0.0256289
4696,0.0258298
4697,0.0263225
4698,0.0265181
4699,0.0269615
4700,0.0272329
4701,0.0276034
4702,0.0278982
4703,0.02828
4704,0.0285273
4705,0.0288389
4706,0.0292445
4707,0.0295984
4708,0.0299843
4709,0.0303156
4710,0.030751
4711,0.0311039
4712,0.0315635
4713,0.0318173
4714,0.032281
4715,0.0327688
4716,0.0331304
4717,0.033465
4718,0.0338762
4719,0.0342434
4720,0.0348077
4721,0.0351089
4722,0.0356362
4723,0.0361068
4724,0.0363789
4725,0.0369809
4726,0.037303
4727,0.0377374
4728,0.0383674
4729,0.0386694
4730,0.0391393
4731,0.0396053
4732,0.040115
4733,0.040666
4734,0.0409595
4735,0.0415084
4736,0.0419804
4737,0.0426746
4738,0.0430245
4739,0.043641
4740,0.044243
4741,0.044666
4742,0.0454314
4743,0.0456115
4744,0.0461765
4745,0.0468266
4746,0.0472917
4747,0.048017
4748,0.0485655
4749,0.0490725
4750,0.0495151
4751,0.0501254
4752,0.0508391
4753,0.0513236
4754,0.0520693
4755,0.0526731
4756,0.0534017
4757,0.0540798
4758,0.0544842
4759,0.0552119
4760,0.0557412
4761,0.0563874
4762,0.0572445
4763,0.0579557
4764,0.0587949
4765,0.0593476
4766,0.0598343
4767,0.0607187
4768,0.061472
4769,0.0620747
4770,0.0628726
4771,0.0636456
4772,0.0642243
4773,0.0650331
4774,0.0657648
4775,0.0667879
4776,0.067215
4777,0.068234
4778,0.0689635
4779,0.069516
4780,0.0703112
4781,0.0715982
4782,0.0723831
4783,0.0730516
4784,0.0741697
4785,0.0749361
4786,0.0757809
4787,0.0767279
4788,0.0773934
4789,0.0787365
4790,0.0796201
4791,0.0805431
4792,0.0810382
4793,0.0821846
4794,0.0832118
4795,0.0839816
4796,0.0852385
4797,0.0859977
4798,0.0872543
4799,0.0882491
4800,0.0893764
4801,0.090282
4802,0.0914214
4803,0.0922846
4804,0.0936055
4805,0.0945924
4806,0.0957427
4807,0.0966392
4808,0.097442
4809,0.0990003
4810,0.100014
4811,0.101187
4812,0.10229
4813,0.103675
4814,0.105012
4815,0.105897
4816,0.107099
4817,0.108671
4818,0.109855
4819,0.110575
4820,0.111588
4821,0.113503
4822,0.114285
4823,0.116308
4824,0.11702
4825,0.11869
4826,0.119724
4827,0.120891
4828,0.123123
4829,0.123607
4830,0.125191
4831,0.127004
4832,0.128134
4833,0.129361
4834,0.130804
4835,0.13236
4836,0.133601
4837,0.13523
4838,0.136752
4839,0.138099
4840,0.138969
4841,0.141041
4842,0.142319
4843,0.14346
4844,0.144497
4845,0.146986
4846,0.147875
4847,0.149762
4848,0.151885
4849,0.153448
4850,0.154588
4851,0.155649
4852,0.157296
4853,0.158514
4854,0.1603
4855,0.161745
4856,0.16335
4857,0.164218
4858,0.167009
4859,0.168287
4860,0.16982
4861,0.170987
4862,0.172676
4863,0.174711
4864,0.176843
4865,0.177706
4866,0.178822
4867,0.180431
4868,0.182226
4869,0.183719
4870,0.184878
4871,0.186864
4872,0.188443
4873,0.191067
4874,0.190953
4875,0.193082
4876,0.195138
4877,0.197183
4878,0.198269
4879,0.198657
4880,0.200299
4881,0.202385
4882,0.203421
4883,0.205591
4884,0.207379
4885,0.208769
4886,0.209583
4887,0.211311
4888,0.212606
4889,0.213773
4890,0.215906
4891,0.216704
4892,0.219097
4893,0.219965
4894,0.222538
4895,0.22322
4896,0.225034
4897,0.225126
4898,0.227784
4899,0.229013
4900,0.230955
4901,0.231951
4902,0.232997
4903,0.234823
4904,0.237259
4905,0.237356
4906,0.239527
4907,0.239873
4908,0.241319
4909,0.24279
4910,0.24475
4911,0.245765
4912,0.247502
4913,0.247923
4914,0.249917
4915,0.251959
4916,0.252184
4917,0.253304
4918,0.254642
4919,0.256658
4920,0.257593
4921,0.259196
4922,0.26103
4923,0.261744
4924,0.262631
4925,0.263768
4926,0.264823
4927,0.266831
4928,0.268262
4929,0.269319
4930,0.270256
4931,0.271079
4932,0.27332
4933,0.274539
4934,0.27433
4935,0.276119
4936,0.278634
4937,0.279295
4938,0.279376
4939,0.281229
4940,0.280986
4941,0.283135
4942,0.284504
4943,0.284874
4944,0.287492
4945,0.287401
4946,0.289914
4947,0.290073
4948,0.291362
4949,0.29294
4950,0.29308
4951,0.293553
4952,0.294678
4953,0.296669
4954,0.296969
4955,0.297775
4956,0.298636
4957,0.300789
4958,0.300916
4959,0.303054
4960,0.303864
4961,0.305603
4962,0.305079
4963,0.305364
4964,0.306847
4965,0.309441
4966,0.309088
4967,0.31116
4968,0.313184
4969,0.312351
4970,0.313737
4971,0.315687
4972,0.315318
4973,0.316242
4974,0.317757
4975,0.31832
4976,0.318201
4977,0.319554
4978,0.320554
4979,0.321531
4980,0.321619
4981,0.321803
4982,0.324131
4983,0.325346
4984,0.326026
4985,0.325801
4986,0.32723
4987,0.329022
4988,0.330649
4989,0.331494
4990,0.332213
4991,0.332949
4992,0.332773
4993,0.332706
4994,0.333732
4995,0.334117
4996,0.336933
4997,0.33752
4998,0.337368
4999,0.33903
5000,0.336397
5001,0.339849
5002,0.339573
5003,0.342183
5004,0.341996
5005,0.342705
5006,0.342589
5007,0.344926
5008,0.345523
5009,0.344874
5010,0.345849
5011,0.347146
5012,0.348979
5013,0.347889
5014,0.350886
5015,0.351242
5016,0.35026
5017,0.350276
5018,0.352683
5019,0.353797
5020,0.353656
5021,0.354533
5022,0.354091
5023,0.355784
5024,0.35688
5025,0.35727
5026,0.358318
5027,0.36001
5028,0.358788
5029,0.362043
5030,0.360437
5031,0.360933
5032,0.360429
5033,0.362404
5034,0.363339
5035,0.363307
5036,0.362733
5037,0.364112
5038,0.365205
5039,0.366048
5040,0.367317
5041,0.367111
5042,0.367419
5043,0.367187
5044,0.369484
5045,0.368694
5046,0.37072
5047,0.369628
5048,0.371577
5049,0.371741
5050,0.371562
5051,0.372591
5052,0.373813
5053,0.374935
5054,0.375387
5055,0.374697
5056,0.375046
5057,0.376664
5058,0.378074
5059,0.377495
5060,0.376918
5061,0.378356
5062,0.376853
5063,0.379804
5064,0.380568
5065,0.380397
5066,0.38164
5067,0.380695
5068,0.383465
5069,0.379018
5070,0.382784
5071,0.386237
5072,0.383919
5073,0.382384
5074,0.385775
5075,0.386094
5076,0.385067
5077,0.386536
5078,0.388298
5079,0.387304
5080,0.387484
5081,0.390188
5082,0.388102
5083,0.388931
5084,0.391046
5085,0.390972
5086,0.392173
5087,0.392513
5088,0.393543
5089,0.390081
5090,0.390173
5091,0.389479
5092,0.393965
5093,0.394575
5094,0.394153
5095,0.395897
5096,0.395575
5097,0.395761
5098,0.396228
5099,0.39594
5100,0.398025
5101,0.397283
5102,0.399532
5103,0.398754
5104,0.39678
5105,0.396219
5106,0.400911
5107,0.397152
5108,0.402222
5109,0.398215
5110,0.399644
5111,0.400621
5112,0.403141
5113,0.402423
5114,0.402603
5115,0.4024
5116,0.402452
5117,0.400182
5118,0.405272
5119,0.40486
5120,0.405836
5121,0.40311
5122,0.406523
5123,0.4054
5124,0.408891
5125,0.407743
5126,0.405528
5127,0.408623
5128,0.407901
5129,0.409105
5130,0.407848
5131,0.408936
5132,0.40985
5133,0.410865
5134,0.410143
5135,0.406514
5136,0.412128
5137,0.410341
5138,0.413463
5139,0.413164
5140,0.410219
5141,0.412328
5142,0.411408
5143,0.411539
5144,0.409508
5145,0.4124
5146,0.41493
5147,0.416484
5148,0.413764
5149,0.414824
5150,0.416969
5151,0.41434
5152,0.418504
5153,0.417209
5154,0.416756
5155,0.414673
5156,0.419981
5157,0.419941
5158,0.417795
5159,0.418548
5160,0.420473
5161,0.418942
5162,0.421251
5163,0.420259
5164,0.420498
5165,0.420359
5166,0.42
5167,0.418842
5168,0.421591
5169,0.418179
5170,0.42091
5171,0.420125
5172,0.423146
5173,0.422455
5174,0.421435
5175,0.420877
5176,0.422758
5177,0.422608
5178,0.422679
5179,0.424559
5180,0.424062
5181,0.425715
5182,0.424857
5183,0.424602
5184,0.428664
5185,0.427413
5186,0.423959
5187,0.430009
5188,0.425758
5189,0.425193
5190,0.425283
5191,0.427155
5192,0.42827
5193,0.427879
5194,0.428337
5195,0.426313
5196,0.425814
5197,0.428359
5198,0.430518
5199,0.430796
5200,0.429357
5201,0.427775
5202,0.42792
5203,0.431126
5204,0.430628
5205,0.429625
5206,0.434041
5207,0.432825
5208,0.431701
5209,0.431599
5210,0.428025
5211,0.432971
5212,0.432064
5213,0.433326
5214,0.4315
5215,0.430761
5216,0.432701
5217,0.432203
5218,0.434129
5219,0.433326
5220,0.433766
5221,0.435136
5222,0.434353
5223,0.434856
5224,0.435385
5225,0.434678
5226,0.435788
5227,0.435205
5228,0.434885
5229,0.437207
5230,0.438295
5231,0.434515
5232,0.438352
5233,0.437057
5234,0.4365
5235,0.437515
5236,0.440111
5237,0.438409
5238,0.436411
5239,0.440059
5240,0.435848
5241,0.437692
5242,0.435946
5243,0.438257
5244,0.439022
5245,0.441926
5246,0.439565
5247,0.439481
5248,0.440565
5249,0.442036
5250,0.440238
5251,0.440979
5252,0.440926
5253,0.439638
5254,0.444596
5255,0.440835
5256,0.439819
5257,0.441838
5258,0.442397
5259,0.442598
5260,0.442112
5261,0.441093
5262,0.444941
5263,0.440799
5264,0.44294
5265,0.443923
5266,0.446154
5267,0.444274
5268,0.442494
5269,0.441571
5270,0.443644
5271,0.446548
5272,0.446601
5273,0.44691
5274,0.446974
5275,0.444332
5276,0.445558
5277,0.442494
5278,0.445232
5279,0.447516
5280,0.446635
5281,0.444198
5282,0.446976
5283,0.446334
5284,0.448763
5285,0.44505
5286,0.446623
5287,0.447772
5288,0.44823
5289,0.447627
5290,0.447718
5291,0.447844
5292,0.450088
5293,0.447984
5294,0.447244
5295,0.450647
5296,0.448979
5297,0.447463
5298,0.448813
5299,0.450568
5300,0.451572
5301,0.451958
5302,0.449579
5303,0.451486
5304,0.452818
5305,0.450341
5306,0.451603
5307,0.451527
5308,0.450674
5309,0.454367
5310,0.451096
5311,0.454005
5312,0.452618
5313,0.454843
5314,0.454316
5315,0.455046
5316,0.454825
5317,0.454406
5318,0.454624
5319,0.452211
5320,0.456073
5321,0.456859
5322,0.456543
5323,0.456608
5324,0.458317
5325,0.458728
5326,0.456571
5327,0.456989
5328,0.461559
5329,0.458669
5330,0.459518
5331,0.463071
5332,0.45999
5333,0.459771
5334,0.458386
5335,0.460108
5336,0.463027
5337,0.460509
5338,0.462204
5339,0.459846
5340,0.464034
5341,0.460146
5342,0.461821
5343,0.462694
5344,0.464422
5345,0.464904
5346,0.46283
5347,0.464172
5348,0.463687
5349,0.465088
5350,0.466588
5351,0.464877
5352,0.464343
5353,0.463757
5354,0.467875
5355,0.464723
5356,0.469002
5357,0.466809
5358,0.467942
5359,0.468663
5360,0.469546
5361,0.469821
5362,0.467963
5363,0.467366
5364,0.470452
5365,0.469731
5366,0.470751
5367,0.471091
5368,0.470409
5369,0.474515
5370,0.47055
5371,0.474283
5372,0.474409
5373,0.472911
5374,0.472518
5375,0.474577
5376,0.474959
5377,0.474488
5378,0.474743
5379,0.475597
5380,0.475335
5381,0.473161
5382,0.47643
5383,0.474537
5384,0.477407
5385,0.476401
5386,0.475578
5387,0.477679
5388,0.477219
5389,0.477499
5390,0.477062
5391,0.479224
5392,0.478395
5393,0.478
5394,0.479369
5395,0.478258
5396,0.478699
5397,0.481036
5398,0.479984
5399,0.48036
5400,0.481885
5401,0.480911
5402,0.481145
5403,0.483952
5404,0.481125
5405,0.480148
5406,0.482274
5407,0.481469
5408,0.482899
5409,0.480925
5410,0.481978
5411,0.484312
5412,0.480981
5413,0.482839
5414,0.482758
5415,0.484997
5416,0.483989
5417,0.482316
5418,0.484268
5419,0.483607
5420,0.485649
5421,0.485508
5422,0.485669
5423,0.485859
5424,0.487146
5425,0.48618
5426,0.487155
5427,0.485832
5428,0.487672
5429,0.486863
5430,0.486041
5431,0.488219
5432,0.488653
5433,0.488266
5434,0.489176
5435,0.486769
5436,0.488176
5437,0.488342
5438,0.490516
5439,0.489628
5440,0.488147
5441,0.489412
5442,0.489563
5443,0.488379
5444,0.491313
5445,0.491121
5446,0.487614
5447,0.492559
5448,0.489252
5449,0.493605
5450,0.493249
5451,0.490203
5452,0.492976
5453,0.489073
5454,0.490983
5455,0.494752
5456,0.493049
5457,0.49281
5458,0.491329
5459,0.492186
5460,0.49269
5461,0.491985
5462,0.493295
5463,0.493805
5464,0.494683
5465,0.4951
5466,0.49412
5467,0.49653
5468,0.492813
5469,0.496693
5470,0.495433
5471,0.494249
5472,0.495617
5473,0.493161
5474,0.497948
5475,0.495179
5476,0.497073
5477,0.495682
5478,0.495163
5479,0.49698
5480,0.495499
5481,0.4964
5482,0.499198
5483,0.498682
5484,0.497512
5485,0.498361
5486,0.498176
5487,0.500341
5488,0.498605
5489,0.497527
5490,0.499698
5491,0.499794
5492,0.496928
5493,0.498103
5494,0.499679
5495,0.497799
5496,0.500336
5497,0.498096
5498,0.496839
5499,0.500489
5500,0.499934
5501,0.500009
5502,0.50085
5503,0.499986
5504,0.502149
5505,0.501043
5506,0.50334
5507,0.501172
5508,0.50259
5509,0.503647
5510,0.500887
5511,0.502283
5512,0.502236
5513,0.502643
5514,0.500347
5515,0.503511
5516,0.504083
5517,0.502938
5518,0.503588
5519,0.502952
5520,0.50014
5521,0.503376
5522,0.504829
5523,0.505831
5524,0.505635
5525,0.503967
5526,0.504616
5527,0.503945
5528,0.50351
5529,0.504207
5530,0.50789
5531,0.502813
5532,0.507481
5533,0.504709
5534,0.505123
5535,0.50434
5536,0.503771
5537,0.505899
5538,0.50618
5539,0.505126
5540,0.505793
5541,0.509522
5542,0.50452
5543,0.507515
5544,0.507257
5545,0.507825
5546,0.507041
5547,0.506245
5548,0.506905
5549,0.509064
5550,0.507124
5551,0.506719
5552,0.506902
5553,0.509027
5554,0.508743
5555,0.506517
5556,0.510001
5557,0.511781
5558,0.508727
5559,0.510317
5560,0.506898
5561,0.509835
5562,0.510286
5563,0.510439
5564,0.509651
5565,0.509386
5566,0.51086
5567,0.513609
5568,0.508515
5569,0.511194
5570,0.511241
5571,0.511412
5572,0.508399
5573,0.510259
5574,0.511576
5575,0.511302
5576,0.512029
5577,0.509875
5578,0.510727
5579,0.512554
5580,0.51035
5581,0.512636
5582,0.512859
5583,0.515007
5584,0.512985
5585,0.511215
5586,0.512157
5587,0.512022
5588,0.514922
5589,0.512503
5590,0.511743
5591,0.511697
5592,0.513383
5593,0.512563
5594,0.514872
5595,0.513686
5596,0.51463
5597,0.513491
5598,0.51491
5599,0.515901
5600,0.512171
5601,0.516725
5602,0.515435
5603,0.512999
5604,0.515207
5605,0.515277
5606,0.514774
5607,0.516273
5608,0.514556
5609,0.516018
5610,0.515397
5611,0.515952
5612,0.514911
5613,0.515583
5614,0.514949
5615,0.514043
5616,0.515837
5617,0.517428
5618,0.515797
5619,0.517289
5620,0.515266
5621,0.515199
5622,0.515763
5623,0.515808
5624,0.516554
5625,0.519186
5626,0.518709
5627,0.517786
5628,0.516753
5629,0.517333
5630,0.516965
5631,0.515283
5632,0.519513
5633,0.516244
5634,0.519102
5635,0.516382
5636,0.517016
5637,0.518363
5638,0.519862
5639,0.518993
5640,0.516998
5641,0.517437
5642,0.519582
5643,0.521192
5644,0.516754
5645,0.520025
5646,0.519587
5647,0.516353
5648,0.520777
5649,0.518041
5650,0.521555
5651,0.520629
5652,0.521225
5653,0.520611
5654,0.522332
5655,0.519664
5656,0.52123
5657,0.520107
5658,0.520044
5659,0.521267
5660,0.519781
5661,0.520637
5662,0.522087
5663,0.521725
5664,0.519613
5665,0.521127
5666,0.524505
5667,0.522686
5668,0.521674
5669,0.523535
5670,0.523107
5671,0.523355
5672,0.523493
5673,0.519151
5674,0.521462
5675,0.519819
5676,0.521864
5677,0.522879
5678,0.522073
5679,0.523275
5680,0.519581
5681,0.52384
5682,0.521515
5683,0.523733
5684,0.526311
5685,0.526098
5686,0.52433
5687,0.522178
5688,0.523064
5689,0.522078
5690,0.524207
5691,0.52156
5692,0.524184
5693,0.523995
5694,0.525642
5695,0.523386
5696,0.522591
5697,0.524086
5698,0.52377
5699,0.521208
5700,0.527144
5701,0.521631
5702,0.52514
5703,0.524118
5704,0.528724
5705,0.524244
5706,0.525148
5707,0.524989
5708,0.527725
5709,0.525912
5710,0.52612
5711,0.525197
5712,0.52598
5713,0.524787
5714,0.526908
5715,0.525031
5716,0.529178
5717,0.52444
5718,0.528708
5719,0.527348
5720,0.528284
5721,0.527454
5722,0.525352
5723,0.524812
5724,0.523676
5725,0.525699
5726,0.525356
5727,0.527578
5728,0.527662
5729,0.52777
5730,0.527981
5731,0.526332
5732,0.528916
5733,0.529222
5734,0.528545
5735,0.528026
5736,0.529138
5737,0.526978
5738,0.527438
5739,0.52948
5740,0.527731
5741,0.529236
5742,0.527679
5743,0.527083
5744,0.531072
5745,0.532148
5746,0.527576
5747,0.52708
5748,0.52899
5749,0.530282
5750,0.530114
5751,0.530201
5752,0.529259
5753,0.530194
5754,0.52944
5755,0.526997
5756,0.527695
5757,0.533865
5758,0.529303
5759,0.531477
5760,0.52986
5761,0.529447
5762,0.52961
5763,0.532221
5764,0.529561
5765,0.528776
5766,0.531029
5767,0.530342
5768,0.530347
5769,0.529125
5770,0.532999
5771,0.531918
5772,0.52894
5773,0.532453
5774,0.530963
5775,0.5303
5776,0.530267
5777,0.528643
5778,0.52949
5779,0.530468
5780,0.530107
5781,0.531941
5782,0.534161
5783,0.531372
5784,0.532518
5785,0.532483
5786,0.533227
5787,0.533226
5788,0.533284
5789,0.531491
5790,0.532766
5791,0.530457
5792,0.532497
5793,0.531895
5794,0.535718
5795,0.531715
5796,0.532508
5797,0.533666
5798,0.53416
5799,0.532363
5800,0.531327
5801,0.531607
5802,0.531271
5803,0.530863
5804,0.534436
5805,0.531824
5806,0.532806
5807,0.532995
5808,0.533911
5809,0.533086
5810,0.536026
5811,0.534301
5812,0.534823
5813,0.535135
5814,0.531684
5815,0.533152
5816,0.532981
5817,0.533428
5818,0.535088
5819,0.535174
5820,0.532605
5821,0.531998
5822,0.534175
5823,0.535562
5824,0.531693
5825,0.536124
5826,0.532424
5827,0.534834
5828,0.532604
5829,0.534527
5830,0.534851
5831,0.535341
5832,0.534146
5833,0.535489
5834,0.533177
5835,0.532962
5836,0.536297
5837,0.537311
5838,0.538172
5839,0.538443
5840,0.534239
5841,0.538912
5842,0.535075
5843,0.535317
5844,0.53505
5845,0.537858
5846,0.537639
5847,0.536153
5848,0.537624
5849,0.534491
5850,0.535245
5851,0.53452
5852,0.538806
5853,0.536585
5854,0.535028
5855,0.537618
5856,0.536423
5857,0.538687
5858,0.538609
5859,0.536411
5860,0.53854
5861,0.538888
5862,0.540546
5863,0.536277
5864,0.537438
5865,0.536971
5866,0.539004
5867,0.537137
5868,0.538022
5869,0.538497
5870,0.533824
5871,0.536522
5872,0.537436
5873,0.54117
5874,0.537337
5875,0.536147
5876,0.539299
5877,0.540078
5878,0.535115
5879,0.535998
5880,0.538684
5881,0.539206
5882,0.539454
5883,0.541245
5884,0.543455
5885,0.534915
5886,0.538849
5887,0.539623
5888,0.538584
5889,0.536525
5890,0.537951
5891,0.540761
5892,0.540869
5893,0.541007
5894,0.539463
5895,0.541607
5896,0.540128
5897,0.540722
5898,0.53975
5899,0.538474
5900,0.541993
5901,0.53982
5902,0.540508
5903,0.54152
5904,0.540565
5905,0.541718
5906,0.543339
5907,0.539217
5908,0.5412
5909,0.541414
5910,0.542482
5911,0.542555
5912,0.542531
5913,0.540586
5914,0.539943
5915,0.540798
5916,0.540654
5917,0.541274
5918,0.540014
5919,0.541943
5920,0.540362
5921,0.543648
5922,0.54109
5923,0.539127
5924,0.539
5925,0.542587
5926,0.5405
5927,0.541717
5928,0.541137
5929,0.542931
5930,0.541256
5931,0.540733
5932,0.543748
5933,0.541273
5934,0.541947
5935,0.542462
5936,0.544027
5937,0.541896
5938,0.538944
5939,0.541228
5940,0.544038
5941,0.542104
5942,0.545458
5943,0.544605
5944,0.540212
5945,0.543099
5946,0.544071
5947,0.543506
5948,0.542657
5949,0.54369
5950,0.544288
5951,0.547427
5952,0.54318
5953,0.539869
5954,0.545902
5955,0.54276
5956,0.544856
5957,0.540374
5958,0.543457
5959,0.542108
5960,0.546985
5961,0.540793
5962,0.545665
5963,0.542855
5964,0.544058
5965,0.546318
5966,0.541934
5967,0.542663
5968,0.545231
5969,0.545161
5970,0.541598
5971,0.542779
5972,0.543125
5973,0.539964
5974,0.544597
5975,0.545379
5976,0.543101
5977,0.546417
5978,0.547671
5979,0.545231
5980,0.547041
5981,0.54409
5982,0.543794
5983,0.545575
5984,0.544796
5985,0.545207
5986,0.546179
5987,0.544222
5988,0.545494
5989,0.544425
5990,0.549579
5991,0.547971
5992,0.546074
5993,0.546866
5994,0.547873
5995,0.546682
5996,0.546789
5997,0.546476
5998,0.546983
5999,0.548962
6000,0.549911
6001,0.546108
6002,0.545207
6003,0.549275
6004,0.545967
6005,0.546158
6006,0.546814
6007,0.54504
6008,0.546335
6009,0.546309
6010,0.544027
6011,0.546949
6012,0.547314
6013,0.547413
6014,0.549132
6015,0.548506
6016,0.546994
6017,0.54753
6018,0.546117
6019,0.54826
6020,0.549958
6021,0.548654
6022,0.544949
6023,0.547196
6024,0.550506
6025,0.549048
6026,0.551489
6027,0.553668
6028,0.550428
6029,0.548964
6030,0.547776
6031,0.546852
6032,0.547314
6033,0.5492
6034,0.546406
6035,0.546505
6036,0.551475
6037,0.549953
6038,0.548113
6039,0.549375
6040,0.546074
6041,0.549023
6042,0.546288
6043,0.551046
6044,0.549058
6045,0.548697
6046,0.547909
6047,0.548509
6048,0.552286
6049,0.548383
6050,0.552206
6051,0.551249
6052,0.549937
6053,0.552169
6054,0.547418
6055,0.549525
6056,0.549365
6057,0.548287
6058,0.548609
6059,0.548641
6060,0.545756
6061,0.554655
6062,0.551645
6063,0.551535
6064,0.552328
6065,0.551244
6066,0.550255
6067,0.552096
6068,0.550083
6069,0.549222
6070,0.553287
6071,0.549592
6072,0.552809
6073,0.552951
6074,0.549147
6075,0.549725
6076,0.550262
6077,0.552019
6078,0.549318
6079,0.553297
6080,0.550599
6081,0.5534
6082,0.551108
6083,0.553334
6084,0.549486
6085,0.551931
6086,0.552586
6087,0.554142
6088,0.550783
6089,0.552973
6090,0.549018
6091,0.555185
6092,0.552498
6093,0.550528
6094,0.554374
6095,0.551396
6096,0.55128
6097,0.556513
6098,0.553616
6099,0.552864
6100,0.552209
6101,0.552188
6102,0.554872
6103,0.549812
6104,0.552671
6105,0.556074
6106,0.554393
6107,0.556044
6108,0.552497
6109,0.552524
6110,0.551355
6111,0.550879
6112,0.552187
6113,0.553634
6114,0.555197
6115,0.55219
6116,0.553756
6117,0.551232
6118,0.553567
6119,0.555009
6120,0.552611
6121,0.552221
6122,0.551817
6123,0.554271
6124,0.551477
6125,0.555141
6126,0.553705
6127,0.553635
6128,0.553739
6129,0.555041
6130,0.553131
6131,0.557623
6132,0.556932
6133,0.558183
6134,0.559317
6135,0.553292
6136,0.550044
6137,0.552954
6138,0.555819
6139,0.554345
6140,0.559802
6141,0.553091
6142,0.5507
6143,0.558188
6144,0.554716
6145,0.550222
6146,0.55223
6147,0.557408
6148,0.553816
6149,0.557835
6150,0.556988
6151,0.554317
6152,0.554479
6153,0.554807
6154,0.555596
6155,0.555005
6156,0.556904
6157,0.558648
6158,0.554836
6159,0.556575
6160,0.556295
6161,0.558238
6162,0.556626
6163,0.554118
6164,0.554228
6165,0.55693
6166,0.556709
6167,0.552743
6168,0.556664
6169,0.553648
6170,0.561244
6171,0.556389
6172,0.553583
6173,0.557999
6174,0.558632
6175,0.560102
6176,0.55774
6177,0.553642
6178,0.555382
6179,0.556047
6180,0.557315
6181,0.557459
6182,0.559487
6183,0.55687
6184,0.558063
6185,0.557711
6186,0.557409
6187,0.560746
6188,0.558234
6189,0.556146
6190,0.558541
6191,0.558111
6192,0.558044
6193,0.556974
6194,0.556342
6195,0.55813
6196,0.560921
6197,0.561373
6198,0.55843
6199,0.558542
6200,0.558594
6201,0.562144
6202,0.557889
6203,0.558017
6204,0.554889
6205,0.559593
6206,0.557181
6207,0.556927
6208,0.559527
6209,0.559096
6210,0.557727
6211,0.55767
6212,0.558729
6213,0.559692
6214,0.556463
6215,0.561864
6216,0.558947
6217,0.55897
6218,0.559805
6219,0.563493
6220,0.563203
6221,0.560529
6222,0.559445
6223,0.560133
6224,0.559202
6225,0.559178
6226,0.556942
6227,0.558551
6228,0.561764
6229,0.560133
6230,0.557944
6231,0.561129
6232,0.557547
6233,0.558825
6234,0.557005
6235,0.561255
6236,0.563616
6237,0.564167
6238,0.559045
6239,0.558354
6240,0.561385
6241,0.55955
6242,0.561884
6243,0.559428
6244,0.560908
6245,0.561489
6246,0.561759
6247,0.562946
6248,0.556628
6249,0.559323
6250,0.561308
6251,0.563207
6252,0.559584
6253,0.56363
6254,0.561338
6255,0.56417
6256,0.561375
6257,0.561335
6258,0.562197
6259,0.561145
6260,0.56871
6261,0.56479
6262,0.56347
6263,0.562562
6264,0.56457
6265,0.563354
6266,0.558787
6267,0.559027
6268,0.562868
6269,0.563334
6270,0.561828
6271,0.561801
6272,0.565675
6273,0.564157
6274,0.568886
6275,0.564563
6276,0.564034
6277,0.560687
6278,0.561346
6279,0.567538
6280,0.565062
6281,0.559714
6282,0.562784
6283,0.561863
6284,0.561518
6285,0.566366
6286,0.561876
6287,0.562857
6288,0.56424
6289,0.56227
6290,0.56566
6291,0.566914
6292,0.566943
6293,0.567153
6294,0.56502
6295,0.564887
6296,0.564136
6297,0.565241
6298,0.564616
6299,0.562678
6300,0.569477
6301,0.561688
6302,0.563319
6303,0.562613
6304,0.563915
6305,0.559818
6306,0.563493
6307,0.568038
6308,0.565651
6309,0.567155
6310,0.563786
6311,0.563731
6312,0.56761
6313,0.566182
6314,0.56808
6315,0.564559
6316,0.564938
6317,0.564552
6318,0.567519
6319,0.565601
6320,0.565827
6321,0.56987
6322,0.563697
6323,0.56704
6324,0.565665
6325,0.563851
6326,0.569622
6327,0.563456
6328,0.563018
6329,0.564004
6330,0.56549
6331,0.568403
6332,0.568227
6333,0.562631
6334,0.566965
6335,0.564924
6336,0.565541
6337,0.565865
6338,0.565741
6339,0.56577
6340,0.566443
6341,0.571398
6342,0.568946
6343,0.568286
6344,0.566545
6345,0.56912
6346,0.565507
6347,0.568109
6348,0.565578
6349,0.571635
6350,0.566771
6351,0.566112
6352,0.568029
6353,0.566198
6354,0.566501
6355,0.56401
6356,0.563124
6357,0.571897
6358,0.565815
6359,0.56805
6360,0.564366
6361,0.570441
6362,0.563243
6363,0.570155
6364,0.566313
6365,0.568843
6366,0.569122
6367,0.570106
6368,0.56893
6369,0.567451
6370,0.569489
6371,0.573791
6372,0.56629
6373,0.563365
6374,0.567514
6375,0.571623
6376,0.566948
6377,0.572632
6378,0.568619
6379,0.566112
6380,0.567199
6381,0.573534
6382,0.571681
6383,0.566823
6384,0.565965
6385,0.569542
6386,0.567893
6387,0.574422
6388,0.571514
6389,0.567935
6390,0.570829
6391,0.571001
6392,0.571396
6393,0.568273
6394,0.568196
6395,0.569275
6396,0.57053
6397,0.568056
6398,0.571515
6399,0.571651
6400,0.570821
6401,0.573593
6402,0.570831
6403,0.571778
6404,0.567559
6405,0.570017
6406,0.57121
6407,0.571879
6408,0.574164
6409,0.569904
6410,0.567905
6411,0.572251
6412,0.576637
6413,0.568732
6414,0.572377
6415,0.568813
6416,0.569962
6417,0.566516
6418,0.570752
6419,0.574221
6420,0.572185
6421,0.571619
6422,0.57182
6423,0.573078
6424,0.574478
6425,0.571787
6426,0.57357
6427,0.568609
6428,0.577304
6429,0.571892
6430,0.572633
6431,0.575569
6432,0.56975
6433,0.574168
6434,0.574377
6435,0.571945
6436,0.571511
6437,0.572294
6438,0.568596
6439,0.57615
6440,0.574653
6441,0.577153
6442,0.571409
6443,0.578643
6444,0.57639
6445,0.571439
6446,0.575279
6447,0.572562
6448,0.572535
6449,0.577689
6450,0.57259
6451,0.574433
6452,0.571285
6453,0.575849
6454,0.576038
6455,0.57019
6456,0.573275
6457,0.576725
6458,0.575433
6459,0.575766
6460,0.571802
6461,0.57247
6462,0.571827
6463,0.573072
6464,0.572525
6465,0.577401
6466,0.574592
6467,0.570752
6468,0.575398
6469,0.577319
6470,0.575112
6471,0.567575
6472,0.576563
6473,0.576735
6474,0.572341
6475,0.579265
6476,0.57467
6477,0.574273
6478,0.576863
6479,0.57808
6480,0.574279
6481,0.576533
6482,0.57649
6483,0.576091
6484,0.57389
6485,0.57258
6486,0.575254
6487,0.576797
6488,0.578916
6489,0.575775
6490,0.572161
6491,0.575083
6492,0.574561
6493,0.580308
6494,0.575883
6495,0.574889
6496,0.576224
6497,0.57843
6498,0.587399
6499,0.579503
6500,0.578239
6501,0.571089
6502,0.575414
6503,0.579622
6504,0.580299
6505,0.575199
6506,0.577805
6507,0.577529
6508,0.578705
6509,0.580518
6510,0.577221
6511,0.580147
6512,0.577344
6513,0.579751
6514,0.578509
6515,0.58019
6516,0.579631
6517,0.575439
6518,0.57825
6519,0.582402
6520,0.578612
6521,0.578125
6522,0.575566
6523,0.578315
6524,0.575552
6525,0.575943
6526,0.576915
6527,0.580044
6528,0.579177
6529,0.577357
6530,0.573965
6531,0.582317
6532,0.577726
6533,0.577649
6534,0.583044
6535,0.576295
6536,0.576201
6537,0.580481
6538,0.581689
6539,0.575013
6540,0.579414
6541,0.574689
6542,0.577331
6543,0.583958
6544,0.578259
6545,0.580574
6546,0.582636
6547,0.577774
6548,0.574513
6549,0.583397
6550,0.584202
6551,0.57918
6552,0.577046
6553,0.578281
6554,0.578361
6555,0.57914
6556,0.580066
6557,0.582006
6558,0.578385
6559,0.580817
6560,0.581184
6561,0.578551
6562,0.581335
6563,0.582952
6564,0.579492
6565,0.581003
6566,0.582552
6567,0.579876
6568,0.579589
6569,0.580873
6570,0.581491
6571,0.580155
6572,0.579513
6573,0.582733
6574,0.579968
6575,0.580018
6576,0.579553
6577,0.581104
6578,0.57785
6579,0.583129
6580,0.583348
6581,0.580091
6582,0.58456
6583,0.584324
6584,0.577438
6585,0.590731
6586,0.583957
6587,0.581641
6588,0.578626
6589,0.580458
6590,0.579534
6591,0.578925
6592,0.583496
6593,0.583054
6594,0.581946
6595,0.585855
6596,0.582747
6597,0.578985
6598,0.583717
6599,0.586242
6600,0.585105
6601,0.585066
6602,0.584821
6603,0.582321
6604,0.584162
6605,0.586393
6606,0.582921
6607,0.583427
6608,0.587302
6609,0.584023
6610,0.57957
6611,0.583823
6612,0.585567
6613,0.587179
6614,0.580456
6615,0.583967
6616,0.585735
6617,0.579745
6618,0.581204
6619,0.580524
6620,0.583377
6621,0.58247
6622,0.582338
6623,0.578497
6624,0.586158
6625,0.58618
6626,0.587141
6627,0.589735
6628,0.578221
6629,0.585685
6630,0.584407
6631,0.583547
6632,0.583947
6633,0.590027
6634,0.588673
6635,0.586552
6636,0.588762
6637,0.590184
6638,0.587378
6639,0.587835
6640,0.586443
6641,0.583279
6642,0.59183
6643,0.582998
6644,0.587141
6645,0.586623
6646,0.584549
6647,0.588225
6648,0.581107
6649,0.592989
6650,0.583229
6651,0.582401
6652,0.585253
6653,0.585102
6654,0.585882
6655,0.586341
6656,0.584575
6657,0.587689
6658,0.586155
6659,0.583827
6660,0.590553
6661,0.590345
6662,0.590344
6663,0.586981
6664,0.585238
6665,0.587449
6666,0.587378
6667,0.582548
6668,0.590458
6669,0.59106
6670,0.584135
6671,0.59128
6672,0.586538
6673,0.588367
6674,0.585908
6675,0.584932
6676,0.582699
6677,0.589895
6678,0.594493
6679,0.585279
6680,0.587439
6681,0.58919
6682,0.581968
6683,0.591246
6684,0.584493
6685,0.590344
6686,0.587523
6687,0.590098
6688,0.583483
6689,0.586165
6690,0.593377
6691,0.592627
6692,0.587669
6693,0.592341
6694,0.584057
6695,0.588531
6696,0.58817
6697,0.587276
6698,0.591735
6699,0.590744
6700,0.590016
6701,0.58959
6702,0.585778
6703,0.589998
6704,0.589391
6705,0.590903
6706,0.589063
6707,0.585712
6708,0.59058
6709,0.586772
6710,0.590911
6711,0.593076
6712,0.590414
6713,0.588227
6714,0.590888
6715,0.591537
6716,0.592885
6717,0.593597
6718,0.597208
6719,0.589186
6720,0.590275
6721,0.592207
6722,0.595607
6723,0.59283
6724,0.594183
6725,0.589069
6726,0.592242
6727,0.594122
6728,0.589475
6729,0.590137
6730,0.589487
6731,0.591448
6732,0.591108
6733,0.589797
6734,0.594664
6735,0.593426
6736,0.591803
6737,0.591522
6738,0.592354
6739,0.588561
6740,0.591159
6741,0.591905
6742,0.586348
6743,0.589954
6744,0.591269
6745,0.597903
6746,0.593487
6747,0.59121
6748,0.593892
6749,0.594689
6750,0.600124
6751,0.596178
6752,0.593128
6753,0.591731
6754,0.589207
6755,0.599487
6756,0.592424
6757,0.593232
6758,0.597144
6759,0.594733
6760,0.598763
6761,0.596094
6762,0.593358
6763,0.593704
6764,0.596968
6765,0.592885
6766,0.590531
6767,0.593849
6768,0.595086
6769,0.589113
6770,0.59521
6771,0.600783
6772,0.595183
6773,0.594104
6774,0.592028
6775,0.598734
6776,0.595671
6777,0.590967
6778,0.592103
6779,0.594867
6780,0.593684
6781,0.598303
6782,0.598886
6783,0.594022
6784,0.599185
6785,0.594191
6786,0.598571
6787,0.59825
6788,0.596676
6789,0.592284
6790,0.597617
6791,0.590493
6792,0.592461
6793,0.594437
6794,0.598739
6795,0.596292
6796,0.592047
6797,0.595335
6798,0.593602
6799,0.596831
6800,0.597131
6801,0.595548
6802,0.601269
6803,0.599161
6804,0.598828
6805,0.592952
6806,0.595053
6807,0.597111
6808,0.592736
6809,0.592976
6810,0.59527
6811,0.598559
6812,0.593337
6813,0.600347
6814,0.599399
6815,0.598607
6816,0.59987
6817,0.597814
6818,0.601842
6819,0.598387
6820,0.603452
6821,0.600824
6822,0.595827
6823,0.598327
6824,0.595488
6825,0.60024
6826,0.598835
6827,0.597184
6828,0.598103
6829,0.606117
6830,0.602632
6831,0.602415
6832,0.601875
6833,0.60091
6834,0.602629
6835,0.60604
6836,0.598444
6837,0.59611
6838,0.597286
6839,0.596517
6840,0.601671
6841,0.597373
6842,0.596818
6843,0.598115
6844,0.597333
6845,0.604342
6846,0.600689
6847,0.59807
6848,0.596776
6849,0.602913
6850,0.600332
6851,0.597332
6852,0.596352
6853,0.593713
6854,0.595789
6855,0.60253
6856,0.607278
6857,0.606196
6858,0.605733
6859,0.600273
6860,0.597554
6861,0.607263
6862,0.602899
6863,0.600015
6864,0.598794
6865,0.601687
6866,0.601923
6867,0.606928
6868,0.601463
6869,0.602877
6870,0.596957
6871,0.600312
6872,0.601913
6873,0.60606
6874,0.603118
6875,0.606835
6876,0.597332
6877,0.599756
6878,0.601804
6879,0.606572
6880,0.607449
6881,0.602602
6882,0.605645
6883,0.605075
6884,0.602008
6885,0.606686
6886,0.601018
6887,0.597173
6888,0.601425
6889,0.599703
6890,0.599699
6891,0.60309
6892,0.594846
6893,0.601009
6894,0.604537
6895,0.606915
6896,0.607101
6897,0.605826
6898,0.604978
6899,0.604253
6900,0.607059
6901,0.606595
6902,0.606202
6903,0.600498
6904,0.60769
6905,0.608191
6906,0.60771
6907,0.60535
6908,0.603211
6909,0.601866
6910,0.602063
6911,0.60757
6912,0.601387
6913,0.603771
6914,0.609048
6915,0.607257
6916,0.605849
6917,0.607487
6918,0.600519
6919,0.601944
6920,0.605909
6921,0.601311
6922,0.603116
6923,0.607694
6924,0.607612
6925,0.601607
6926,0.608192
6927,0.604752
6928,0.606097
6929,0.604806
6930,0.609834
6931,0.608187
6932,0.606182
6933,0.604892
6934,0.604471
6935,0.608789
6936};
sin
double sin(const BesAngle a)
Definition BesAngle.h:210
cos
double cos(const BesAngle a)
Definition BesAngle.h:213
report
ostream & report(Severity severity, const char *facility)
Definition EvtReport.cc:36
ERROR
@ ERROR
Definition EvtReport.hh:49
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition EvtDecayBase.cc:297
EvtDecayBase::getParentId
EvtId getParentId()
Definition EvtDecayBase.hh:60
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition EvtDecayBase.cc:504
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:65
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:482
EvtDecayBase::getDaug
EvtId getDaug(int i)
Definition EvtDecayBase.hh:66
EvtDecayProb::setProb
void setProb(double prob)
Definition EvtDecayProb.hh:34
EvtDsToPipPipPimPi0::getName
void getName(std::string &name)
Definition EvtDsToPipPipPimPi0.cc:31
EvtDsToPipPipPimPi0::decay
void decay(EvtParticle *p)
Definition EvtDsToPipPipPimPi0.cc:150
EvtPDL::getId
static EvtId getId(const std::string &name)
Definition EvtPDL.cc:287
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:120
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:84
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1070
EvtVector4R::get
double get(int i) const
Definition EvtVector4R.hh:179
s12
double double double double * s12
Definition qcdloop1.h:77
s23
double double double double double * s23
Definition qcdloop1.h:77
b
const double b
Definition slope.cxx:9