BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtD0ToKSKppimpi0.cc
Go to the documentation of this file.
1//--------------------------------------------------------------------------
2// Environment:
3// This software is part of models developed at BES collaboration
4// based on the EvtGen framework. If you use all or part
5// of it, please give an appropriate acknowledgement.
6//
7// Copyright Information: See EvtGen/BesCopyright
8// Copyright (A) 2006 Ping Rong-Gang @IHEP
9//
10// Module: EvtD0ToKSKppimpi0.cc
11// the necessary file: EvtD0ToKSKppimpi0.hh
12//
13// Description: D0 -> KS K+ pi- pi0
14//
15// Modification history:
16//
17// Zhenxuan Li Jun. 15, 2023 Module created
18//------------------------------------------------------------------------
19#include "EvtD0ToKSKppimpi0.hh"
20#include "../EvtGenBase/EvtComplex.hh"
21#include "../EvtGenBase/EvtDecayTable.hh"
22#include "../EvtGenBase/EvtGenKine.hh"
23#include "../EvtGenBase/EvtPDL.hh"
24#include "../EvtGenBase/EvtParticle.hh"
25#include "../EvtGenBase/EvtPatches.hh"
26#include "../EvtGenBase/EvtReport.hh"
27#include <cmath>
28#include <iostream>
29#include <stdlib.h>
30using namespace std;
31
32EvtD0ToKSKppimpi0::~EvtD0ToKSKppimpi0() {}
33
34void EvtD0ToKSKppimpi0::getName( std::string& model_name ) { model_name = "D0ToKSKppimpi0"; }
35
36EvtDecayBase* EvtD0ToKSKppimpi0::clone() { return new EvtD0ToKSKppimpi0; }
37
38void EvtD0ToKSKppimpi0::init() {
39 checkNArg( 0 );
40 checkNDaug( 4 );
41 checkSpinParent( EvtSpinType::SCALAR );
42 // cout << "Initializing EvtD0ToKSKppimpi0" << endl;
43
44 mD = 1.86484;
45 mPi = 0.13957;
46 mKa = 0.493677;
47 math_pi = 3.1415926;
48 mKstr0 = 0.89555;
49 mrhop = 0.77511;
50 GKstr0 = 0.0473;
51 Grhop = 0.1491;
52 mKstrp = 0.89167;
53 mrho0 = 0.77526;
54 GKstrp = 0.0514;
55 Grho0 = 0.1474;
56 ma980 = 0.980;
57 Ga980 = 0.0972;
58 ma1_1260 = 1.230;
59 mK1_1270 = 1.272;
60 Ga1_1260 = 0.420;
61 GK1_1270 = 0.087;
62 meta1475 = 1.475;
63 mK1_1400 = 1.403;
64 Geta1475 = 0.09;
65 GK1_1400 = 0.174;
66
67 rho[1] = 1.0; // 2S
68 phi[1] = 0.0;
69 rho[0] = 1.1688955687; // 1S
70 phi[0] = 5.4715577597;
71 rho[2] = 1.2403434033; // 2P
72 phi[2] = -4.2080261661;
73 rho[3] = 4.9828746114; // 16S
74 phi[3] = 5.6470624082;
75 rho[4] = 2.8452580959; // 201
76 phi[4] = 0.1889705326;
77 rho[5] = 2.8452580959; // 201
78 phi[5] = 0.1889705326;
79 rho[6] = 8.7553053301; // 101
80 phi[6] = -5.9775594415;
81
82 modetype[0] = 1001;
83 modetype[1] = 1002;
84 modetype[2] = 1002;
85 modetype[3] = 2008;
86 modetype[4] = 4001;
87 modetype[5] = 4002;
88 modetype[6] = 3001;
89
90 int GG[4][4] = { { 1, 0, 0, 0 }, { 0, -1, 0, 0 }, { 0, 0, -1, 0 }, { 0, 0, 0, -1 } };
91 int EE[4][4][4][4] = {
92 { { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
93 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 1 }, { 0, 0, -1, 0 } },
94 { { 0, 0, 0, 0 }, { 0, 0, 0, -1 }, { 0, 0, 0, 0 }, { 0, 1, 0, 0 } },
95 { { 0, 0, 0, 0 }, { 0, 0, 1, 0 }, { 0, -1, 0, 0 }, { 0, 0, 0, 0 } } },
96 { { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, -1 }, { 0, 0, 1, 0 } },
97 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
98 { { 0, 0, 0, 1 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { -1, 0, 0, 0 } },
99 { { 0, 0, -1, 0 }, { 0, 0, 0, 0 }, { 1, 0, 0, 0 }, { 0, 0, 0, 0 } } },
100 { { { 0, 0, 0, 0 }, { 0, 0, 0, 1 }, { 0, 0, 0, 0 }, { 0, -1, 0, 0 } },
101 { { 0, 0, 0, -1 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 1, 0, 0, 0 } },
102 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
103 { { 0, 1, 0, 0 }, { -1, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } } },
104 { { { 0, 0, 0, 0 }, { 0, 0, -1, 0 }, { 0, 1, 0, 0 }, { 0, 0, 0, 0 } },
105 { { 0, 0, 1, 0 }, { 0, 0, 0, 0 }, { -1, 0, 0, 0 }, { 0, 0, 0, 0 } },
106 { { 0, -1, 0, 0 }, { 1, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } },
107 { { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 }, { 0, 0, 0, 0 } } } };
108 for ( int i = 0; i < 4; i++ )
109 {
110 for ( int j = 0; j < 4; j++ )
111 {
112 G[i][j] = GG[i][j];
113 for ( int k = 0; k < 4; k++ )
114 {
115 for ( int l = 0; l < 4; l++ ) { E[i][j][k][l] = EE[i][j][k][l]; }
116 }
117 }
118 }
119}
120
121void EvtD0ToKSKppimpi0::initProbMax() { setProbMax( 172500 ); }
122
123void EvtD0ToKSKppimpi0::decay( EvtParticle* p ) {
124 //-----------for max value------------------
125 /* double maxprob = 0.0;
126 for(int ir=0;ir<=60000000;ir++){
127 p->initializePhaseSpace(getNDaug(),getDaugs());
128 EvtVector4R ks = p->getDaug(0)->getP4();
129 EvtVector4R kp = p->getDaug(1)->getP4();
130 EvtVector4R pim = p->getDaug(2)->getP4();
131 EvtVector4R pi0 = p->getDaug(3)->getP4();
132 mother_c=EvtPDL::getStdHep(p->getId());
133 int cc;
134 if(mother_c==421){
135 cc=1;
136 } else if(mother_c==-421){
137 cc=-1;
138 }
139 double value;
140 double KS[4],Kp[4],Pim[4],Pi0[4];
141 KS[0] = ks.get(0); Kp[0] = kp.get(0); Pim[0] = pim.get(0); Pi0[0] =
142 pi0.get(0); KS[1] = cc*ks.get(1); Kp[1] = cc*kp.get(1); Pim[1] = cc*pim.get(1); Pi0[1] =
143 cc*pi0.get(1); KS[2] = cc*ks.get(2); Kp[2] = cc*kp.get(2); Pim[2] = cc*pim.get(2); Pi0[2]
144 = cc*pi0.get(2); KS[3] = cc*ks.get(3); Kp[3] = cc*kp.get(3); Pim[3] = cc*pim.get(3);
145 Pi0[3] = cc*pi0.get(3); calPDF(KS, Kp, Pim, Pi0, value); if(value>maxprob) {
146 maxprob=value;
147 std::cout << "Max PDF = " << ir << " prob= " << value << std::endl;
148 }
149 }
150 std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
151 return;*/
152 //-----------------------------------------------
153 p->initializePhaseSpace( getNDaug(), getDaugs() );
154 EvtVector4R ks = p->getDaug( 0 )->getP4();
155 EvtVector4R kp = p->getDaug( 1 )->getP4();
156 EvtVector4R pim = p->getDaug( 2 )->getP4();
157 EvtVector4R pi0 = p->getDaug( 3 )->getP4();
158
159 mother_c = EvtPDL::getStdHep( p->getId() );
160 int cc;
161 if ( mother_c == 421 ) { cc = 1; }
162 else if ( mother_c == -421 ) { cc = -1; }
163
164 double KS[4], Kp[4], Pim[4], Pi0[4];
165 KS[0] = ks.get( 0 );
166 Kp[0] = kp.get( 0 );
167 Pim[0] = pim.get( 0 );
168 Pi0[0] = pi0.get( 0 );
169 KS[1] = cc * ks.get( 1 );
170 Kp[1] = cc * kp.get( 1 );
171 Pim[1] = cc * pim.get( 1 );
172 Pi0[1] = cc * pi0.get( 1 );
173 KS[2] = cc * ks.get( 2 );
174 Kp[2] = cc * kp.get( 2 );
175 Pim[2] = cc * pim.get( 2 );
176 Pi0[2] = cc * pi0.get( 2 );
177 KS[3] = cc * ks.get( 3 );
178 Kp[3] = cc * kp.get( 3 );
179 Pim[3] = cc * pim.get( 3 );
180 Pi0[3] = cc * pi0.get( 3 );
181 // KS[0] = 0.581505; Kp[0] = 0.654953; Pim[0] = 0.344256; Pi0[0] = 0.297861;
182 // KS[1] = 0.302649; Kp[1] =-0.369050; Pim[1] =-0.164686; Pi0[1] = 0.027881;
183 // KS[2] = 0.038261; Kp[2] =-0.221134; Pim[2] =-0.247290; Pi0[2] = 0.203629;
184 // KS[3] =-0.032410; Kp[3] =-0.012204; Pim[3] =-0.103723; Pi0[3] = 0.168109;
185
186 double value;
187 calPDF( KS, Kp, Pim, Pi0, value );
188 setProb( value );
189 return;
190}
191
192void EvtD0ToKSKppimpi0::calPDF( double KS[], double Kp[], double Pim[], double Pi0[],
193 double& Result ) {
194 double cof[2], amp_tmp[2], amp_PDF[2], PDF[2];
195 double flag[3], mass_R[2], width_R[2];
196 double sa[3], sb[3], sc[3], B[3];
197 double t1D[4], t1V[4]; //, t1A[4];
198 double pS[4], pV[4], pD[4], pA[4], pP[4];
199 double pro1[2], pro2[2], proKPi_S[2], pro[2];
200 double temp_PDF = 0;
201
202 double rD2 = 25.0;
203 double rRes2 = 9.0;
204 double mass1[7] = { mKstr0, mKstrp, mKstrp, mrhop, mKstrp, mKstrp, ma980 };
205 double mass2[7] = { mKstr0, mKstrp, mKstrp, mK1_1270, meta1475, meta1475, mrhop };
206 double width1[7] = { GKstr0, GKstrp, GKstrp, Grhop, GKstrp, GKstrp, Ga980 };
207 double width2[7] = { GKstr0, GKstrp, GKstrp, GK1_1270, Geta1475, Geta1475, Grhop };
208 double g0[7] = { 1, 1, 1, 1, 1, 1, 1 };
209 double g1[7] = { 1, 1, 1, 1, 1, 1, 1 };
210 double g2[7] = { 0, 0, 1, 0, 0, 0, 1 };
211
212 PDF[0] = 0;
213 PDF[1] = 0;
214
215 for ( int i = 0; i < 7; i++ )
216 {
217 flag[0] = g0[i];
218 flag[1] = g1[i];
219 flag[2] = g2[i];
220 mass_R[0] = mass1[i];
221 mass_R[1] = mass2[i];
222 width_R[0] = width1[i];
223 width_R[1] = width2[i];
224
225 amp_tmp[0] = 0;
226 amp_tmp[1] = 0;
227
228 cof[0] = rho[i] * cos( phi[i] );
229 cof[1] = rho[i] * sin( phi[i] );
230
231 if ( modetype[i] == 1001 )
232 { // D0 -> K*0 K*0
233 double t1V1[4], t1V2[4], t2D[4][4];
234 temp_PDF = 0;
235 double pV1[4], pV2[4];
236 for ( int ii = 0; ii != 4; ii++ )
237 {
238 pV1[ii] = KS[ii] + Pi0[ii];
239 pV2[ii] = Kp[ii] + Pim[ii];
240 pD[ii] = pV1[ii] + pV2[ii];
241 }
242 sa[0] = SCADot( pV1, pV1 );
243 sb[0] = SCADot( KS, KS );
244 sc[0] = SCADot( Pi0, Pi0 );
245 sa[1] = SCADot( pV2, pV2 );
246 sb[1] = SCADot( Kp, Kp );
247 sc[1] = SCADot( Pim, Pim );
248 sa[2] = SCADot( pD, pD );
249 sb[2] = sa[0];
250 sc[2] = sa[1];
251 if ( flag[0] == 1 )
252 {
253 propagatorRBW( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0],
254 rRes2, 1, pro1 );
255 }
256 else if ( flag[0] == 0 )
257 {
258 pro1[0] = 1;
259 pro1[1] = 0;
260 }
261 if ( flag[1] == 1 )
262 {
263 propagatorRBW( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
264 rRes2, 1, pro2 );
265 }
266 else if ( flag[1] == 0 )
267 {
268 pro2[0] = 1;
269 pro2[1] = 0;
270 }
271 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
272 B[1] = Barrier( mass_R[1] * mass_R[1], 1, sa[1], sb[1], sc[1], rRes2 );
273 calt1( KS, Pi0, t1V1 );
274 calt1( Kp, Pim, t1V2 );
275 if ( flag[2] == 0 )
276 {
277 for ( int ii = 0; ii != 4; ii++ ) { temp_PDF += ( G[ii][ii] ) * t1V1[ii] * t1V2[ii]; }
278 B[2] = 1;
279 }
280 if ( flag[2] == 1 )
281 {
282 calt1( pV1, pV2, t1D );
283 for ( int ii = 0; ii != 4; ii++ )
284 {
285 for ( int j = 0; j != 4; j++ )
286 {
287 for ( int k = 0; k != 4; k++ )
288 {
289 for ( int l = 0; l != 4; l++ )
290 {
291 temp_PDF += E[ii][j][k][l] * pD[ii] * t1D[j] * t1V1[k] * t1V2[l] *
292 ( G[ii][ii] ) * ( G[j][j] ) * ( G[l][l] ) * ( G[k][k] );
293 }
294 }
295 }
296 }
297 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
298 }
299 if ( flag[2] == 2 )
300 {
301 calt2( pV1, pV2, t2D );
302 for ( int ii = 0; ii != 4; ii++ )
303 {
304 for ( int j = 0; j != 4; j++ )
305 { temp_PDF += t2D[ii][j] * t1V1[ii] * t1V2[j] * ( G[ii][ii] ) * ( G[j][j] ); }
306 }
307 B[2] = Barrier( mD * mD, 2, sa[2], sb[2], sc[2], rD2 );
308 }
309 Com_Multi( pro1, pro2, pro );
310 amp_tmp[0] = temp_PDF * B[0] * B[1] * B[2] * pro[0];
311 amp_tmp[1] = temp_PDF * B[0] * B[1] * B[2] * pro[1];
312 }
313 else if ( modetype[i] == 1002 )
314 { // D0 -> K*- K*+
315 double t1V1[4], t1V2[4], t2D[4][4];
316 temp_PDF = 0;
317 double pV1[4], pV2[4];
318 for ( int ii = 0; ii != 4; ii++ )
319 {
320 pV1[ii] = KS[ii] + Pim[ii];
321 pV2[ii] = Kp[ii] + Pi0[ii];
322 pD[ii] = pV1[ii] + pV2[ii];
323 }
324 sa[0] = SCADot( pV1, pV1 );
325 sb[0] = SCADot( KS, KS );
326 sc[0] = SCADot( Pim, Pim );
327 sa[1] = SCADot( pV2, pV2 );
328 sb[1] = SCADot( Kp, Kp );
329 sc[1] = SCADot( Pi0, Pi0 );
330 sa[2] = SCADot( pD, pD );
331 sb[2] = sa[0];
332 sc[2] = sa[1];
333 if ( flag[0] == 1 )
334 {
335 propagatorRBW( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0],
336 rRes2, 1, pro1 );
337 }
338 else if ( flag[0] == 0 )
339 {
340 pro1[0] = 1;
341 pro1[1] = 0;
342 }
343 if ( flag[1] == 1 )
344 {
345 propagatorRBW( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
346 rRes2, 1, pro2 );
347 }
348 else if ( flag[1] == 0 )
349 {
350 pro2[0] = 1;
351 pro2[1] = 0;
352 }
353 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
354 B[1] = Barrier( mass_R[1] * mass_R[1], 1, sa[1], sb[1], sc[1], rRes2 );
355 calt1( KS, Pim, t1V1 );
356 calt1( Kp, Pi0, t1V2 );
357 if ( flag[2] == 0 )
358 {
359 for ( int ii = 0; ii != 4; ii++ ) { temp_PDF += ( G[ii][ii] ) * t1V1[ii] * t1V2[ii]; }
360 B[2] = 1;
361 }
362 if ( flag[2] == 1 )
363 {
364 calt1( pV1, pV2, t1D );
365 for ( int ii = 0; ii != 4; ii++ )
366 {
367 for ( int j = 0; j != 4; j++ )
368 {
369 for ( int k = 0; k != 4; k++ )
370 {
371 for ( int l = 0; l != 4; l++ )
372 {
373 temp_PDF += E[ii][j][k][l] * pD[ii] * t1D[j] * t1V1[k] * t1V2[l] *
374 ( G[ii][ii] ) * ( G[j][j] ) * ( G[l][l] ) * ( G[k][k] );
375 }
376 }
377 }
378 }
379 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
380 }
381 if ( flag[2] == 2 )
382 {
383 calt2( pV1, pV2, t2D );
384 for ( int ii = 0; ii != 4; ii++ )
385 {
386 for ( int j = 0; j != 4; j++ )
387 { temp_PDF += t2D[ii][j] * t1V1[ii] * t1V2[j] * ( G[ii][ii] ) * ( G[j][j] ); }
388 }
389 B[2] = Barrier( mD * mD, 2, sa[2], sb[2], sc[2], rD2 );
390 }
391 Com_Multi( pro1, pro2, pro );
392 amp_tmp[0] = temp_PDF * B[0] * B[1] * B[2] * pro[0];
393 amp_tmp[1] = temp_PDF * B[0] * B[1] * B[2] * pro[1];
394 }
395 else if ( modetype[i] == 2001 )
396 { // D0 -> a1()+ pi-, a1()+ -> K*+ KS
397 temp_PDF = 0;
398 double t2A[4][4];
399 for ( int ii = 0; ii != 4; ii++ )
400 {
401 pV[ii] = Kp[ii] + Pi0[ii];
402 pA[ii] = pV[ii] + KS[ii];
403 pD[ii] = pA[ii] + Pim[ii];
404 }
405 sa[0] = SCADot( pV, pV );
406 sb[0] = SCADot( Kp, Kp );
407 sc[0] = SCADot( Pi0, Pi0 );
408 sa[1] = SCADot( pA, pA );
409 sb[1] = sa[0];
410 sc[1] = SCADot( KS, KS );
411 sa[2] = SCADot( pD, pD );
412 sb[2] = sa[1];
413 sc[2] = SCADot( Pim, Pim );
414
415 if ( flag[0] == 1 )
416 {
417 propagatorRBW( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0],
418 rRes2, 1, pro1 ); // K*
419 }
420 else if ( flag[0] == 0 )
421 {
422 pro1[0] = 1;
423 pro1[1] = 0;
424 }
425
426 if ( flag[1] == 1 )
427 {
428 propagatorRBW( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
429 rRes2, flag[2], pro2 );
430 }
431 else if ( flag[1] == 0 )
432 {
433 pro2[0] = 1;
434 pro2[1] = 0;
435 }
436
437 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
438 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
439 calt1( Kp, Pi0, t1V );
440 calt1( pA, Pim, t1D );
441 if ( flag[2] == 0 )
442 {
443 for ( int ii = 0; ii != 4; ii++ )
444 {
445 for ( int j = 0; j != 4; j++ )
446 {
447 temp_PDF +=
448 t1D[ii] * ( G[ii][j] - pA[ii] * pA[j] / sa[1] ) * t1V[j] * G[ii][ii] * G[j][j];
449 }
450 }
451 B[1] = 1;
452 }
453 else if ( flag[2] == 2 )
454 {
455 calt2( pV, KS, t2A );
456 for ( int ii = 0; ii != 4; ii++ )
457 {
458 for ( int j = 0; j != 4; j++ )
459 { temp_PDF += t1D[ii] * t2A[ii][j] * t1V[j] * G[ii][ii] * G[j][j]; }
460 }
461 B[1] = Barrier( mass_R[1] * mass_R[1], 2, sa[1], sb[1], sc[1], rRes2 );
462 }
463 Com_Multi( pro1, pro2, pro );
464 amp_tmp[0] = temp_PDF * B[0] * B[1] * B[2] * pro[0];
465 amp_tmp[1] = temp_PDF * B[0] * B[1] * B[2] * pro[1];
466 }
467 else if ( modetype[i] == 2002 )
468 { // D0 -> a1()+ pi-, a1()+ -> K*0 K+
469 temp_PDF = 0;
470 double t2A[4][4];
471 for ( int ii = 0; ii != 4; ii++ )
472 {
473 pV[ii] = KS[ii] + Pi0[ii];
474 pA[ii] = pV[ii] + Kp[ii];
475 pD[ii] = pA[ii] + Pim[ii];
476 }
477 sa[0] = SCADot( pV, pV );
478 sb[0] = SCADot( KS, KS );
479 sc[0] = SCADot( Pi0, Pi0 );
480 sa[1] = SCADot( pA, pA );
481 sb[1] = sa[0];
482 sc[1] = SCADot( Kp, Kp );
483 sa[2] = SCADot( pD, pD );
484 sb[2] = sa[1];
485 sc[2] = SCADot( Pim, Pim );
486
487 if ( flag[0] == 1 )
488 {
489 propagatorRBW( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0],
490 rRes2, 1, pro1 ); // K*
491 }
492 else if ( flag[0] == 0 )
493 {
494 pro1[0] = 1;
495 pro1[1] = 0;
496 }
497
498 if ( flag[1] == 1 )
499 {
500 propagatorRBW( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
501 rRes2, flag[2], pro2 );
502 }
503 else if ( flag[1] == 0 )
504 {
505 pro2[0] = 1;
506 pro2[1] = 0;
507 }
508
509 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
510 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
511 calt1( KS, Pi0, t1V );
512 calt1( pA, Pim, t1D );
513 if ( flag[2] == 0 )
514 {
515 for ( int ii = 0; ii != 4; ii++ )
516 {
517 for ( int j = 0; j != 4; j++ )
518 {
519 temp_PDF +=
520 t1D[ii] * ( G[ii][j] - pA[ii] * pA[j] / sa[1] ) * t1V[j] * G[ii][ii] * G[j][j];
521 }
522 }
523 B[1] = 1;
524 }
525 else if ( flag[2] == 2 )
526 {
527 calt2( pV, Kp, t2A );
528 for ( int ii = 0; ii != 4; ii++ )
529 {
530 for ( int j = 0; j != 4; j++ )
531 { temp_PDF += t1D[ii] * t2A[ii][j] * t1V[j] * G[ii][ii] * G[j][j]; }
532 }
533 B[1] = Barrier( mass_R[1] * mass_R[1], 2, sa[1], sb[1], sc[1], rRes2 );
534 }
535 Com_Multi( pro1, pro2, pro );
536 amp_tmp[0] = temp_PDF * B[0] * B[1] * B[2] * pro[0];
537 amp_tmp[1] = temp_PDF * B[0] * B[1] * B[2] * pro[1];
538 }
539 else if ( modetype[i] == 2003 )
540 { // D0 -> a1()0 pi0, a1()0 -> K*0 KS
541 temp_PDF = 0;
542 double t2A[4][4];
543 for ( int ii = 0; ii != 4; ii++ )
544 {
545 pV[ii] = Kp[ii] + Pim[ii];
546 pA[ii] = pV[ii] + KS[ii];
547 pD[ii] = pA[ii] + Pi0[ii];
548 }
549 sa[0] = SCADot( pV, pV );
550 sb[0] = SCADot( Kp, Kp );
551 sc[0] = SCADot( Pim, Pim );
552 sa[1] = SCADot( pA, pA );
553 sb[1] = sa[0];
554 sc[1] = SCADot( KS, KS );
555 sa[2] = SCADot( pD, pD );
556 sb[2] = sa[1];
557 sc[2] = SCADot( Pi0, Pi0 );
558
559 if ( flag[0] == 1 )
560 {
561 propagatorRBW( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0],
562 rRes2, 1, pro1 ); // K*
563 }
564 else if ( flag[0] == 0 )
565 {
566 pro1[0] = 1;
567 pro1[1] = 0;
568 }
569
570 if ( flag[1] == 1 )
571 {
572 propagatorRBW( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
573 rRes2, flag[2], pro2 );
574 }
575 else if ( flag[1] == 0 )
576 {
577 pro2[0] = 1;
578 pro2[1] = 0;
579 }
580
581 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
582 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
583 calt1( Kp, Pim, t1V );
584 calt1( pA, Pi0, t1D );
585 if ( flag[2] == 0 )
586 {
587 for ( int ii = 0; ii != 4; ii++ )
588 {
589 for ( int j = 0; j != 4; j++ )
590 {
591 temp_PDF +=
592 t1D[ii] * ( G[ii][j] - pA[ii] * pA[j] / sa[1] ) * t1V[j] * G[ii][ii] * G[j][j];
593 }
594 }
595 B[1] = 1;
596 }
597 else if ( flag[2] == 2 )
598 {
599 calt2( pV, KS, t2A );
600 for ( int ii = 0; ii != 4; ii++ )
601 {
602 for ( int j = 0; j != 4; j++ )
603 { temp_PDF += t1D[ii] * t2A[ii][j] * t1V[j] * G[ii][ii] * G[j][j]; }
604 }
605 B[1] = Barrier( mass_R[1] * mass_R[1], 2, sa[1], sb[1], sc[1], rRes2 );
606 }
607 Com_Multi( pro1, pro2, pro );
608 amp_tmp[0] = temp_PDF * B[0] * B[1] * B[2] * pro[0];
609 amp_tmp[1] = temp_PDF * B[0] * B[1] * B[2] * pro[1];
610 }
611 else if ( modetype[i] == 2004 )
612 { // D0 -> a1()0 pi0, a1()0 -> K*- K+
613 temp_PDF = 0;
614 double t2A[4][4];
615 for ( int ii = 0; ii != 4; ii++ )
616 {
617 pV[ii] = KS[ii] + Pim[ii];
618 pA[ii] = pV[ii] + Kp[ii];
619 pD[ii] = pA[ii] + Pi0[ii];
620 }
621 sa[0] = SCADot( pV, pV );
622 sb[0] = SCADot( KS, KS );
623 sc[0] = SCADot( Pim, Pim );
624 sa[1] = SCADot( pA, pA );
625 sb[1] = sa[0];
626 sc[1] = SCADot( Kp, Kp );
627 sa[2] = SCADot( pD, pD );
628 sb[2] = sa[1];
629 sc[2] = SCADot( Pi0, Pi0 );
630
631 if ( flag[0] == 1 )
632 {
633 propagatorRBW( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0],
634 rRes2, 1, pro1 ); // K*
635 }
636 else if ( flag[0] == 0 )
637 {
638 pro1[0] = 1;
639 pro1[1] = 0;
640 }
641
642 if ( flag[1] == 1 )
643 {
644 propagatorRBW( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
645 rRes2, flag[2], pro2 );
646 }
647 else if ( flag[1] == 0 )
648 {
649 pro2[0] = 1;
650 pro2[1] = 0;
651 }
652
653 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
654 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
655 calt1( KS, Pim, t1V );
656 calt1( pA, Pi0, t1D );
657 if ( flag[2] == 0 )
658 {
659 for ( int ii = 0; ii != 4; ii++ )
660 {
661 for ( int j = 0; j != 4; j++ )
662 {
663 temp_PDF +=
664 t1D[ii] * ( G[ii][j] - pA[ii] * pA[j] / sa[1] ) * t1V[j] * G[ii][ii] * G[j][j];
665 }
666 }
667 B[1] = 1;
668 }
669 else if ( flag[2] == 2 )
670 {
671 calt2( pV, Kp, t2A );
672 for ( int ii = 0; ii != 4; ii++ )
673 {
674 for ( int j = 0; j != 4; j++ )
675 { temp_PDF += t1D[ii] * t2A[ii][j] * t1V[j] * G[ii][ii] * G[j][j]; }
676 }
677 B[1] = Barrier( mass_R[1] * mass_R[1], 2, sa[1], sb[1], sc[1], rRes2 );
678 }
679 Com_Multi( pro1, pro2, pro );
680 amp_tmp[0] = temp_PDF * B[0] * B[1] * B[2] * pro[0];
681 amp_tmp[1] = temp_PDF * B[0] * B[1] * B[2] * pro[1];
682 }
683 else if ( modetype[i] == 2005 )
684 { // D0 -> K1(1270)- K+, K1()- -> KS rho-
685 temp_PDF = 0;
686 double t2A[4][4];
687 for ( int ii = 0; ii != 4; ii++ )
688 {
689 pV[ii] = Pim[ii] + Pi0[ii];
690 pA[ii] = pV[ii] + KS[ii];
691 pD[ii] = pA[ii] + Kp[ii];
692 }
693 sa[0] = SCADot( pV, pV );
694 sb[0] = SCADot( Pim, Pim );
695 sc[0] = SCADot( Pi0, Pi0 );
696 sa[1] = SCADot( pA, pA );
697 sb[1] = sa[0];
698 sc[1] = SCADot( KS, KS );
699 sa[2] = SCADot( pD, pD );
700 sb[2] = sa[1];
701 sc[2] = SCADot( Kp, Kp );
702
703 if ( flag[0] == 1 )
704 {
705 propagatorGS( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2,
706 pro1 ); // rho
707 }
708 else if ( flag[0] == 0 )
709 {
710 pro1[0] = 1;
711 pro1[1] = 0;
712 }
713
714 if ( flag[1] == 1 )
715 {
716 propagatorRBW_K1( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
717 rRes2, flag[2], pro2 );
718 }
719 else if ( flag[1] == 0 )
720 {
721 pro2[0] = 1;
722 pro2[1] = 0;
723 }
724
725 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
726 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
727 calt1( Pim, Pi0, t1V );
728 calt1( pA, Kp, t1D );
729 if ( flag[2] == 0 )
730 {
731 for ( int ii = 0; ii != 4; ii++ )
732 {
733 for ( int j = 0; j != 4; j++ )
734 {
735 temp_PDF +=
736 t1D[ii] * ( G[ii][j] - pA[ii] * pA[j] / sa[1] ) * t1V[j] * G[ii][ii] * G[j][j];
737 }
738 }
739 B[1] = 1;
740 }
741 else if ( flag[2] == 2 )
742 {
743 calt2( pV, KS, t2A );
744 for ( int ii = 0; ii != 4; ii++ )
745 {
746 for ( int j = 0; j != 4; j++ )
747 { temp_PDF += t1D[ii] * t2A[ii][j] * t1V[j] * G[ii][ii] * G[j][j]; }
748 }
749 B[1] = Barrier( mass_R[1] * mass_R[1], 2, sa[1], sb[1], sc[1], rRes2 );
750 }
751 Com_Multi( pro1, pro2, pro );
752 amp_tmp[0] = temp_PDF * B[0] * B[1] * B[2] * pro[0];
753 amp_tmp[1] = temp_PDF * B[0] * B[1] * B[2] * pro[1];
754 }
755 else if ( modetype[i] == 2006 )
756 { // D0 -> K1(1270/1400)- K+, K1()- -> K*0 pi-
757 temp_PDF = 0;
758 double t2A[4][4];
759 for ( int ii = 0; ii != 4; ii++ )
760 {
761 pV[ii] = KS[ii] + Pi0[ii];
762 pA[ii] = pV[ii] + Pim[ii];
763 pD[ii] = pA[ii] + Kp[ii];
764 }
765 sa[0] = SCADot( pV, pV );
766 sb[0] = SCADot( KS, KS );
767 sc[0] = SCADot( Pi0, Pi0 );
768 sa[1] = SCADot( pA, pA );
769 sb[1] = sa[0];
770 sc[1] = SCADot( Pim, Pim );
771 sa[2] = SCADot( pD, pD );
772 sb[2] = sa[1];
773 sc[2] = SCADot( Kp, Kp );
774
775 if ( flag[0] == 1 )
776 {
777 // propagatorGS(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2,
778 // pro1);//rho
779 propagatorRBW( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0],
780 rRes2, 1, pro1 ); // K*
781 }
782 else if ( flag[0] == 0 )
783 {
784 pro1[0] = 1;
785 pro1[1] = 0;
786 }
787
788 if ( flag[1] == 1 )
789 {
790 propagatorRBW_K1( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
791 rRes2, flag[2], pro2 );
792 }
793 else if ( flag[1] == 0 )
794 {
795 pro2[0] = 1;
796 pro2[1] = 0;
797 }
798
799 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
800 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
801 calt1( KS, Pi0, t1V );
802 calt1( pA, Kp, t1D );
803 if ( flag[2] == 0 )
804 {
805 for ( int ii = 0; ii != 4; ii++ )
806 {
807 for ( int j = 0; j != 4; j++ )
808 {
809 temp_PDF +=
810 t1D[ii] * ( G[ii][j] - pA[ii] * pA[j] / sa[1] ) * t1V[j] * G[ii][ii] * G[j][j];
811 }
812 }
813 B[1] = 1;
814 }
815 else if ( flag[2] == 2 )
816 {
817 calt2( pV, Pim, t2A );
818 for ( int ii = 0; ii != 4; ii++ )
819 {
820 for ( int j = 0; j != 4; j++ )
821 { temp_PDF += t1D[ii] * t2A[ii][j] * t1V[j] * G[ii][ii] * G[j][j]; }
822 }
823 B[1] = Barrier( mass_R[1] * mass_R[1], 2, sa[1], sb[1], sc[1], rRes2 );
824 }
825 Com_Multi( pro1, pro2, pro );
826 amp_tmp[0] = temp_PDF * B[0] * B[1] * B[2] * pro[0];
827 amp_tmp[1] = temp_PDF * B[0] * B[1] * B[2] * pro[1];
828 }
829 else if ( modetype[i] == 2007 )
830 { // D0 -> K1(1270/1400)- K+, K1()- -> K*- pi0
831 temp_PDF = 0;
832 double t2A[4][4];
833 for ( int ii = 0; ii != 4; ii++ )
834 {
835 pV[ii] = KS[ii] + Pim[ii];
836 pA[ii] = pV[ii] + Pi0[ii];
837 pD[ii] = pA[ii] + Kp[ii];
838 }
839 sa[0] = SCADot( pV, pV );
840 sb[0] = SCADot( KS, KS );
841 sc[0] = SCADot( Pim, Pim );
842 sa[1] = SCADot( pA, pA );
843 sb[1] = sa[0];
844 sc[1] = SCADot( Pi0, Pi0 );
845 sa[2] = SCADot( pD, pD );
846 sb[2] = sa[1];
847 sc[2] = SCADot( Kp, Kp );
848
849 if ( flag[0] == 1 )
850 {
851 // propagatorGS(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2,
852 // pro1);//rho
853 propagatorRBW( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0],
854 rRes2, 1, pro1 ); // K*
855 }
856 else if ( flag[0] == 0 )
857 {
858 pro1[0] = 1;
859 pro1[1] = 0;
860 }
861
862 if ( flag[1] == 1 )
863 {
864 propagatorRBW_K1( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
865 rRes2, flag[2], pro2 );
866 }
867 else if ( flag[1] == 0 )
868 {
869 pro2[0] = 1;
870 pro2[1] = 0;
871 }
872
873 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
874 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
875 calt1( KS, Pim, t1V );
876 calt1( pA, Kp, t1D );
877 if ( flag[2] == 0 )
878 {
879 for ( int ii = 0; ii != 4; ii++ )
880 {
881 for ( int j = 0; j != 4; j++ )
882 {
883 temp_PDF +=
884 t1D[ii] * ( G[ii][j] - pA[ii] * pA[j] / sa[1] ) * t1V[j] * G[ii][ii] * G[j][j];
885 }
886 }
887 B[1] = 1;
888 }
889 else if ( flag[2] == 2 )
890 {
891 calt2( pV, Pi0, t2A );
892 for ( int ii = 0; ii != 4; ii++ )
893 {
894 for ( int j = 0; j != 4; j++ )
895 { temp_PDF += t1D[ii] * t2A[ii][j] * t1V[j] * G[ii][ii] * G[j][j]; }
896 }
897 B[1] = Barrier( mass_R[1] * mass_R[1], 2, sa[1], sb[1], sc[1], rRes2 );
898 }
899 Com_Multi( pro1, pro2, pro );
900 amp_tmp[0] = -1.0 * temp_PDF * B[0] * B[1] * B[2] * pro[0];
901 amp_tmp[1] = -1.0 * temp_PDF * B[0] * B[1] * B[2] * pro[1];
902 }
903 else if ( modetype[i] == 2008 )
904 { // D0 -> K1(1270)0 KS, K1()0 -> K+ rho-
905 temp_PDF = 0;
906 double t2A[4][4];
907 for ( int ii = 0; ii != 4; ii++ )
908 {
909 pV[ii] = Pim[ii] + Pi0[ii];
910 pA[ii] = pV[ii] + Kp[ii];
911 pD[ii] = pA[ii] + KS[ii];
912 }
913 sa[0] = SCADot( pV, pV );
914 sb[0] = SCADot( Pim, Pim );
915 sc[0] = SCADot( Pi0, Pi0 );
916 sa[1] = SCADot( pA, pA );
917 sb[1] = sa[0];
918 sc[1] = SCADot( Kp, Kp );
919 sa[2] = SCADot( pD, pD );
920 sb[2] = sa[1];
921 sc[2] = SCADot( KS, KS );
922
923 if ( flag[0] == 1 )
924 {
925 propagatorGS( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2,
926 pro1 ); // rho
927 // propagatorRBW(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2,
928 // 1, pro1);//K*
929 }
930 else if ( flag[0] == 0 )
931 {
932 pro1[0] = 1;
933 pro1[1] = 0;
934 }
935
936 if ( flag[1] == 1 )
937 {
938 propagatorRBW_K1( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
939 rRes2, flag[2], pro2 );
940 }
941 else if ( flag[1] == 0 )
942 {
943 pro2[0] = 1;
944 pro2[1] = 0;
945 }
946
947 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
948 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
949 calt1( Pim, Pi0, t1V );
950 calt1( pA, KS, t1D );
951 if ( flag[2] == 0 )
952 {
953 for ( int ii = 0; ii != 4; ii++ )
954 {
955 for ( int j = 0; j != 4; j++ )
956 {
957 temp_PDF +=
958 t1D[ii] * ( G[ii][j] - pA[ii] * pA[j] / sa[1] ) * t1V[j] * G[ii][ii] * G[j][j];
959 }
960 }
961 B[1] = 1;
962 }
963 else if ( flag[2] == 2 )
964 {
965 calt2( pV, Kp, t2A );
966 for ( int ii = 0; ii != 4; ii++ )
967 {
968 for ( int j = 0; j != 4; j++ )
969 { temp_PDF += t1D[ii] * t2A[ii][j] * t1V[j] * G[ii][ii] * G[j][j]; }
970 }
971 B[1] = Barrier( mass_R[1] * mass_R[1], 2, sa[1], sb[1], sc[1], rRes2 );
972 }
973 Com_Multi( pro1, pro2, pro );
974 amp_tmp[0] = temp_PDF * B[0] * B[1] * B[2] * pro[0];
975 amp_tmp[1] = temp_PDF * B[0] * B[1] * B[2] * pro[1];
976 }
977 else if ( modetype[i] == 2009 )
978 { // D0 -> K1(1270/1400)0 KS, K1()0 -> K*+ pi-
979 temp_PDF = 0;
980 double t2A[4][4];
981 for ( int ii = 0; ii != 4; ii++ )
982 {
983 pV[ii] = Kp[ii] + Pi0[ii];
984 pA[ii] = pV[ii] + Pim[ii];
985 pD[ii] = pA[ii] + KS[ii];
986 }
987 sa[0] = SCADot( pV, pV );
988 sb[0] = SCADot( Kp, Kp );
989 sc[0] = SCADot( Pi0, Pi0 );
990 sa[1] = SCADot( pA, pA );
991 sb[1] = sa[0];
992 sc[1] = SCADot( Pim, Pim );
993 sa[2] = SCADot( pD, pD );
994 sb[2] = sa[1];
995 sc[2] = SCADot( KS, KS );
996
997 if ( flag[0] == 1 )
998 {
999 // propagatorGS(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2,
1000 // pro1);//rho
1001 propagatorRBW( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0],
1002 rRes2, 1, pro1 ); // K*
1003 }
1004 else if ( flag[0] == 0 )
1005 {
1006 pro1[0] = 1;
1007 pro1[1] = 0;
1008 }
1009
1010 if ( flag[1] == 1 )
1011 {
1012 propagatorRBW_K1( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
1013 rRes2, flag[2], pro2 );
1014 }
1015 else if ( flag[1] == 0 )
1016 {
1017 pro2[0] = 1;
1018 pro2[1] = 0;
1019 }
1020
1021 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
1022 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
1023 calt1( Kp, Pi0, t1V );
1024 calt1( pA, KS, t1D );
1025 if ( flag[2] == 0 )
1026 {
1027 for ( int ii = 0; ii != 4; ii++ )
1028 {
1029 for ( int j = 0; j != 4; j++ )
1030 {
1031 temp_PDF +=
1032 t1D[ii] * ( G[ii][j] - pA[ii] * pA[j] / sa[1] ) * t1V[j] * G[ii][ii] * G[j][j];
1033 }
1034 }
1035 B[1] = 1;
1036 }
1037 else if ( flag[2] == 2 )
1038 {
1039 calt2( pV, Pim, t2A );
1040 for ( int ii = 0; ii != 4; ii++ )
1041 {
1042 for ( int j = 0; j != 4; j++ )
1043 { temp_PDF += t1D[ii] * t2A[ii][j] * t1V[j] * G[ii][ii] * G[j][j]; }
1044 }
1045 B[1] = Barrier( mass_R[1] * mass_R[1], 2, sa[1], sb[1], sc[1], rRes2 );
1046 }
1047 Com_Multi( pro1, pro2, pro );
1048 amp_tmp[0] = temp_PDF * B[0] * B[1] * B[2] * pro[0];
1049 amp_tmp[1] = temp_PDF * B[0] * B[1] * B[2] * pro[1];
1050 }
1051 else if ( modetype[i] == 2010 )
1052 { // D0 -> K1(1270/1400)0 KS, K1()0 -> K*0 pi0
1053 temp_PDF = 0;
1054 double t2A[4][4];
1055 for ( int ii = 0; ii != 4; ii++ )
1056 {
1057 pV[ii] = Kp[ii] + Pim[ii];
1058 pA[ii] = pV[ii] + Pi0[ii];
1059 pD[ii] = pA[ii] + KS[ii];
1060 }
1061 sa[0] = SCADot( pV, pV );
1062 sb[0] = SCADot( Kp, Kp );
1063 sc[0] = SCADot( Pim, Pim );
1064 sa[1] = SCADot( pA, pA );
1065 sb[1] = sa[0];
1066 sc[1] = SCADot( Pi0, Pi0 );
1067 sa[2] = SCADot( pD, pD );
1068 sb[2] = sa[1];
1069 sc[2] = SCADot( KS, KS );
1070
1071 if ( flag[0] == 1 )
1072 {
1073 // propagatorGS(mass_R[0]*mass_R[0],mass_R[0], width_R[0], sa[0], sb[0], sc[0], rRes2,
1074 // pro1);//rho
1075 propagatorRBW( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0],
1076 rRes2, 1, pro1 ); // K*
1077 }
1078 else if ( flag[0] == 0 )
1079 {
1080 pro1[0] = 1;
1081 pro1[1] = 0;
1082 }
1083
1084 if ( flag[1] == 1 )
1085 {
1086 propagatorRBW_K1( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
1087 rRes2, flag[2], pro2 );
1088 }
1089 else if ( flag[1] == 0 )
1090 {
1091 pro2[0] = 1;
1092 pro2[1] = 0;
1093 }
1094
1095 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
1096 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
1097 calt1( Kp, Pim, t1V );
1098 calt1( pA, KS, t1D );
1099 if ( flag[2] == 0 )
1100 {
1101 for ( int ii = 0; ii != 4; ii++ )
1102 {
1103 for ( int j = 0; j != 4; j++ )
1104 {
1105 temp_PDF +=
1106 t1D[ii] * ( G[ii][j] - pA[ii] * pA[j] / sa[1] ) * t1V[j] * G[ii][ii] * G[j][j];
1107 }
1108 }
1109 B[1] = 1;
1110 }
1111 else if ( flag[2] == 2 )
1112 {
1113 calt2( pV, Pi0, t2A );
1114 for ( int ii = 0; ii != 4; ii++ )
1115 {
1116 for ( int j = 0; j != 4; j++ )
1117 { temp_PDF += t1D[ii] * t2A[ii][j] * t1V[j] * G[ii][ii] * G[j][j]; }
1118 }
1119 B[1] = Barrier( mass_R[1] * mass_R[1], 2, sa[1], sb[1], sc[1], rRes2 );
1120 }
1121 Com_Multi( pro1, pro2, pro );
1122 amp_tmp[0] = -1.0 * temp_PDF * B[0] * B[1] * B[2] * pro[0];
1123 amp_tmp[1] = -1.0 * temp_PDF * B[0] * B[1] * B[2] * pro[1];
1124 }
1125 else if ( modetype[i] == 3001 )
1126 { // D0 to a0(980)+ rho-
1127 temp_PDF = 0;
1128 for ( int ii = 0; ii != 4; ii++ )
1129 {
1130 pS[ii] = KS[ii] + Kp[ii];
1131 pV[ii] = Pim[ii] + Pi0[ii];
1132 pD[ii] = pS[ii] + pV[ii];
1133 }
1134 sa[0] = SCADot( pS, pS );
1135 sb[0] = SCADot( KS, KS );
1136 sc[0] = SCADot( Kp, Kp );
1137 sa[1] = SCADot( pV, pV );
1138 sb[1] = SCADot( Pim, Pim );
1139 sc[1] = SCADot( Pi0, Pi0 );
1140 sa[2] = SCADot( pD, pD );
1141 sb[2] = sa[0];
1142 sc[2] = sa[1];
1143 if ( flag[0] == 1 ) { propagatorFlatte( mass_R[0], width_R[0], sa[0], pro1 ); }
1144 else if ( flag[0] == 0 )
1145 {
1146 pro1[0] = 1;
1147 pro1[1] = 0;
1148 }
1149
1150 if ( flag[1] == 1 )
1151 {
1152 propagatorGS( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1], rRes2,
1153 pro2 ); // rho
1154 }
1155 else if ( flag[1] == 0 )
1156 {
1157 pro2[0] = 1;
1158 pro2[1] = 0;
1159 }
1160
1161 B[1] = Barrier( mass_R[1] * mass_R[1], 1, sa[1], sb[1], sc[1], rRes2 );
1162 B[2] = Barrier( mD * mD, 1, sa[2], sb[2], sc[2], rD2 );
1163
1164 calt1( Pim, Pi0, t1V );
1165 calt1( pS, pV, t1D );
1166 for ( int ii = 0; ii != 4; ii++ ) { temp_PDF += G[ii][ii] * t1D[ii] * t1V[ii]; }
1167 Com_Multi( pro1, pro2, pro );
1168 amp_tmp[0] = temp_PDF * B[1] * B[2] * pro[0];
1169 amp_tmp[1] = temp_PDF * B[1] * B[2] * pro[1];
1170 }
1171 else if ( modetype[i] == 4001 )
1172 { // D0 to eta() pi0, eta() -> K*- K+
1173 temp_PDF = 0;
1174 for ( int ii = 0; ii != 4; ii++ )
1175 {
1176 pV[ii] = KS[ii] + Pim[ii];
1177 pP[ii] = pV[ii] + Kp[ii];
1178 pD[ii] = pP[ii] + Pi0[ii];
1179 }
1180 sa[0] = SCADot( pV, pV );
1181 sb[0] = SCADot( KS, KS );
1182 sc[0] = SCADot( Pim, Pim );
1183 sa[1] = SCADot( pP, pP );
1184 sb[1] = sa[0];
1185 sc[1] = SCADot( Kp, Kp );
1186 sa[2] = SCADot( pD, pD );
1187 sb[2] = sa[1];
1188 sc[2] = SCADot( Pi0, Pi0 );
1189 ;
1190 if ( flag[0] == 1 )
1191 {
1192 propagatorRBW( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0],
1193 rRes2, 1, pro1 ); // K*
1194 }
1195 else if ( flag[0] == 0 )
1196 {
1197 pro1[0] = 1;
1198 pro1[1] = 0;
1199 }
1200
1201 if ( flag[1] == 1 )
1202 {
1203 propagatorRBW( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
1204 rRes2, 1, pro2 ); // eta
1205 }
1206 else if ( flag[1] == 0 )
1207 {
1208 pro2[0] = 1;
1209 pro2[1] = 0;
1210 }
1211
1212 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
1213 B[1] = Barrier( mass_R[1] * mass_R[1], 1, sa[1], sb[1], sc[1], rRes2 );
1214
1215 calt1( KS, Pim, t1V );
1216 for ( int ii = 0; ii != 4; ii++ ) { temp_PDF += Kp[ii] * t1V[ii] * ( G[ii][ii] ); }
1217 Com_Multi( pro1, pro2, pro );
1218 amp_tmp[0] = temp_PDF * B[0] * B[1] * pro[0];
1219 amp_tmp[1] = temp_PDF * B[0] * B[1] * pro[1];
1220 }
1221 else if ( modetype[i] == 4002 )
1222 { // D0 to eta() pi0, eta() -> K*0 KS
1223 temp_PDF = 0;
1224 for ( int ii = 0; ii != 4; ii++ )
1225 {
1226 pV[ii] = Kp[ii] + Pim[ii];
1227 pP[ii] = pV[ii] + KS[ii];
1228 pD[ii] = pP[ii] + Pi0[ii];
1229 }
1230 sa[0] = SCADot( pV, pV );
1231 sb[0] = SCADot( Kp, Kp );
1232 sc[0] = SCADot( Pim, Pim );
1233 sa[1] = SCADot( pP, pP );
1234 sb[1] = sa[0];
1235 sc[1] = SCADot( KS, KS );
1236 sa[2] = SCADot( pD, pD );
1237 sb[2] = sa[1];
1238 sc[2] = SCADot( Pi0, Pi0 );
1239 ;
1240 if ( flag[0] == 1 )
1241 {
1242 propagatorRBW( mass_R[0] * mass_R[0], mass_R[0], width_R[0], sa[0], sb[0], sc[0],
1243 rRes2, 1, pro1 ); // K*
1244 }
1245 else if ( flag[0] == 0 )
1246 {
1247 pro1[0] = 1;
1248 pro1[1] = 0;
1249 }
1250
1251 if ( flag[1] == 1 )
1252 {
1253 propagatorRBW( mass_R[1] * mass_R[1], mass_R[1], width_R[1], sa[1], sb[1], sc[1],
1254 rRes2, 1, pro2 ); // eta
1255 }
1256 else if ( flag[1] == 0 )
1257 {
1258 pro2[0] = 1;
1259 pro2[1] = 0;
1260 }
1261
1262 B[0] = Barrier( mass_R[0] * mass_R[0], 1, sa[0], sb[0], sc[0], rRes2 );
1263 B[1] = Barrier( mass_R[1] * mass_R[1], 1, sa[1], sb[1], sc[1], rRes2 );
1264
1265 calt1( Kp, Pim, t1V );
1266 for ( int ii = 0; ii != 4; ii++ ) { temp_PDF += KS[ii] * t1V[ii] * ( G[ii][ii] ); }
1267 Com_Multi( pro1, pro2, pro );
1268 amp_tmp[0] = temp_PDF * B[0] * B[1] * pro[0];
1269 amp_tmp[1] = temp_PDF * B[0] * B[1] * pro[1];
1270 }
1271
1272 Com_Multi( amp_tmp, cof, amp_PDF );
1273 // printf("pdf: %d %.10f %.10f\n",i,amp_PDF[0],amp_PDF[1]);
1274 PDF[0] += amp_PDF[0];
1275 PDF[1] += amp_PDF[1];
1276 }
1277 double value = PDF[0] * PDF[0] + PDF[1] * PDF[1];
1278 Result = value;
1279 // printf("PDF[0] : %.10f \n",PDF[0]);
1280 // printf("PDF[1] : %.10f \n",PDF[1]);
1281 // printf("Value : %.10f \n",value);
1282}
1283
1284void EvtD0ToKSKppimpi0::Com_Multi( double a1[2], double a2[2], double res[2] ) {
1285 res[0] = a1[0] * a2[0] - a1[1] * a2[1];
1286 res[1] = a1[1] * a2[0] + a1[0] * a2[1];
1287}
1288void EvtD0ToKSKppimpi0::Com_Divide( double a1[2], double a2[2], double res[2] ) {
1289 double tmp = a2[0] * a2[0] + a2[1] * a2[1];
1290 res[0] = ( a1[0] * a2[0] + a1[1] * a2[1] ) / tmp;
1291 res[1] = ( a1[1] * a2[0] - a1[0] * a2[1] ) / tmp;
1292}
1293double EvtD0ToKSKppimpi0::SCADot( double a1[4], double a2[4] ) {
1294 double _cal = a1[0] * a2[0] - a1[1] * a2[1] - a1[2] * a2[2] - a1[3] * a2[3];
1295 return _cal;
1296}
1297double EvtD0ToKSKppimpi0::Barrier( double mass2, int l, double sa, double sb, double sc,
1298 double r2 ) {
1299 double F;
1300 double tmp = sa + sb - sc;
1301 double q = fabs( 0.25 * tmp * tmp / sa - sb );
1302 // if (q < 0) q = 1e-16;
1303 double tmp2 = mass2 + sb - sc;
1304 double q0 = fabs( 0.25 * tmp2 * tmp2 / mass2 - sb );
1305 // if (q0 < 0) q0 = 1e-16;
1306 double z = q * r2;
1307 double z0 = q0 * r2;
1308 if ( l == 1 ) { F = sqrt( ( 1.0 + z0 ) / ( 1.0 + z ) ); }
1309 else if ( l == 2 )
1310 {
1311 double z2 = z * z;
1312 double z02 = z0 * z0;
1313 F = sqrt( ( 9.0 + 3.0 * z0 + z02 ) / ( 9.0 + 3.0 * z + z2 ) );
1314 }
1315 else { F = 1.0; }
1316 return F;
1317}
1318void EvtD0ToKSKppimpi0::calt1( double daug1[4], double daug2[4], double t1[4] ) {
1319 double p, pq, tmp;
1320 double pa[4], qa[4];
1321 for ( int i = 0; i < 4; i++ )
1322 {
1323 pa[i] = daug1[i] + daug2[i];
1324 qa[i] = daug1[i] - daug2[i];
1325 }
1326 p = SCADot( pa, pa );
1327 pq = SCADot( pa, qa );
1328 tmp = pq / p;
1329 for ( int i = 0; i < 4; i++ ) { t1[i] = qa[i] - tmp * pa[i]; }
1330}
1331void EvtD0ToKSKppimpi0::calt2( double daug1[4], double daug2[4], double t2[4][4] ) {
1332 double p, r;
1333 double pa[4], t1[4];
1334 calt1( daug1, daug2, t1 );
1335 r = SCADot( t1, t1 ) / 3.0;
1336 for ( int i = 0; i < 4; i++ ) { pa[i] = daug1[i] + daug2[i]; }
1337 p = SCADot( pa, pa );
1338 for ( int i = 0; i < 4; i++ )
1339 {
1340 for ( int j = 0; j < 4; j++ )
1341 { t2[i][j] = t1[i] * t1[j] - r * ( G[i][j] - pa[i] * pa[j] / p ); }
1342 }
1343}
1344void EvtD0ToKSKppimpi0::propagator( double mass2, double mass, double width, double sx,
1345 double prop[2] ) {
1346 double a[2], b[2];
1347 a[0] = 1;
1348 a[1] = 0;
1349 b[0] = mass2 - sx;
1350 b[1] = -mass * width;
1351 Com_Divide( a, b, prop );
1352}
1353double EvtD0ToKSKppimpi0::wid( double mass2, double mass, double sa, double sb, double sc,
1354 double r2, int l ) {
1355 double widm = 0.;
1356 double m = sqrt( sa );
1357 double tmp = sb - sc;
1358 double tmp1 = sa + tmp;
1359 double q = fabs( 0.25 * tmp1 * tmp1 / sa - sb );
1360 // if(q<0) q = 1e-16;
1361 double tmp2 = mass2 + tmp;
1362 double q0 = fabs( 0.25 * tmp2 * tmp2 / mass2 - sb );
1363 // if(q0<0) q0 = 1e-16;
1364 double z = q * r2;
1365 double z0 = q0 * r2;
1366 double t = q / q0;
1367 if ( l == 0 ) { widm = sqrt( t ) * mass / m; }
1368 else if ( l == 1 ) { widm = t * sqrt( t ) * mass / m * ( 1 + z0 ) / ( 1 + z ); }
1369 else if ( l == 2 )
1370 { widm = t * t * sqrt( t ) * mass / m * ( 9 + 3 * z0 + z0 * z0 ) / ( 9 + 3 * z + z * z ); }
1371 return widm;
1372}
1373double EvtD0ToKSKppimpi0::widl1( double mass2, double mass, double sa, double sb, double sc,
1374 double r2 ) {
1375 double widm = 0.;
1376 double m = sqrt( sa );
1377 double tmp = sb - sc;
1378 double tmp1 = sa + tmp;
1379 double q = fabs( 0.25 * tmp1 * tmp1 / sa - sb );
1380 // if(q<0) q = 1e-16;
1381 double tmp2 = mass2 + tmp;
1382 double q0 = fabs( 0.25 * tmp2 * tmp2 / mass2 - sb );
1383 // if(q0<0) q0 = 1e-16;
1384 double z = q * r2;
1385 double z0 = q0 * r2;
1386 double F = ( 1 + z0 ) / ( 1 + z );
1387 double t = q / q0;
1388 widm = t * sqrt( t ) * mass / m * F;
1389 return widm;
1390}
1391void EvtD0ToKSKppimpi0::propagatorRBW( double mass2, double mass, double width, double sa,
1392 double sb, double sc, double r2, int l,
1393 double prop[2] ) {
1394 double a[2], b[2];
1395 a[0] = 1;
1396 a[1] = 0;
1397 b[0] = mass2 - sa;
1398 b[1] = -mass * width * wid( mass2, mass, sa, sb, sc, r2, l );
1399 Com_Divide( a, b, prop );
1400}
1401void EvtD0ToKSKppimpi0::propagatorRBWl1( double mass2, double mass, double width, double sa,
1402 double sb, double sc, double r2, double prop[2] ) {
1403 double a[2], b[2];
1404 a[0] = 1;
1405 a[1] = 0;
1406 b[0] = mass2 - sa;
1407 b[1] = -mass * width * widl1( mass2, mass, sa, sb, sc, r2 );
1408 Com_Divide( a, b, prop );
1409}
1410void EvtD0ToKSKppimpi0::propagatorRBW_a1( double mass2, double mass, double width, double sa,
1411 double sb, double sc, double r2, int l,
1412 double prop[2] ) {
1413 double a[2], b[2];
1414 int iii = int( sa * 1000 / 1.55 ) - 115;
1415 if ( iii < 0 ) { iii = 0; }
1416 double a1_width[2000] = {
1417 // 1885 bin, s = (0.17825, 3.1)
1418 0.0000000172, 0.0000000487, 0.0000001049, 0.0000001927, 0.0000003185, 0.0000004888,
1419 0.0000007099, 0.0000009878, 0.0000013285, 0.0000017378, 0.0000022212, 0.0000027845,
1420 0.0000034328, 0.0000041717, 0.0000050061, 0.0000059414, 0.0000069824, 0.0000081342,
1421 0.0000094014, 0.0000107890, 0.0000123016, 0.0000139437, 0.0000157200, 0.0000176350,
1422 0.0000196931, 0.0000218986, 0.0000242559, 0.0000267692, 0.0000294428, 0.0000322808,
1423 0.0000352875, 0.0000384668, 0.0000418228, 0.0000453596, 0.0000490811, 0.0000529913,
1424 0.0000570941, 0.0000613934, 0.0000658931, 0.0000705969, 0.0000755087, 0.0000806324,
1425 0.0000859716, 0.0000915300, 0.0000973115, 0.0001033197, 0.0001095583, 0.0001160310,
1426 0.0001227415, 0.0001296933, 0.0001368901, 0.0001443355, 0.0001520332, 0.0001599867,
1427 0.0001681995, 0.0001766753, 0.0001854177, 0.0001944301, 0.0002037162, 0.0002132795,
1428 0.0002231235, 0.0002332517, 0.0002436677, 0.0002543750, 0.0002653771, 0.0002766775,
1429 0.0002882799, 0.0003001876, 0.0003124041, 0.0003249331, 0.0003377780, 0.0003509423,
1430 0.0003644296, 0.0003782433, 0.0003923870, 0.0004068643, 0.0004216785, 0.0004368333,
1431 0.0004523323, 0.0004681788, 0.0004843765, 0.0005009290, 0.0005178397, 0.0005351122,
1432 0.0005527502, 0.0005707571, 0.0005891365, 0.0006078922, 0.0006270275, 0.0006465462,
1433 0.0006664518, 0.0006867480, 0.0007074384, 0.0007285267, 0.0007500165, 0.0007719114,
1434 0.0007942151, 0.0008169314, 0.0008400639, 0.0008636163, 0.0008875923, 0.0009119957,
1435 0.0009368303, 0.0009620997, 0.0009878078, 0.0010139580, 0.0010405550, 0.0010676020,
1436 0.0010951030, 0.0011230610, 0.0011514820, 0.0011803670, 0.0012097230, 0.0012395510,
1437 0.0012698570, 0.0013006440, 0.0013319160, 0.0013636770, 0.0013959320, 0.0014286840,
1438 0.0014619370, 0.0014956950, 0.0015299630, 0.0015647450, 0.0016000440, 0.0016358650,
1439 0.0016722130, 0.0017090900, 0.0017465020, 0.0017844530, 0.0018229470, 0.0018619890,
1440 0.0019015820, 0.0019417310, 0.0019824400, 0.0020237150, 0.0020655590, 0.0021079760,
1441 0.0021509730, 0.0021945510, 0.0022387180, 0.0022834760, 0.0023288310, 0.0023747880,
1442 0.0024213500, 0.0024685230, 0.0025163120, 0.0025647210, 0.0026137560, 0.0026634200,
1443 0.0027137200, 0.0027646590, 0.0028162430, 0.0028684780, 0.0029213670, 0.0029749170,
1444 0.0030291320, 0.0030840170, 0.0031395780, 0.0031958200, 0.0032527490, 0.0033103700,
1445 0.0033686880, 0.0034277080, 0.0034874370, 0.0035478800, 0.0036090420, 0.0036709290,
1446 0.0037335480, 0.0037969030, 0.0038610000, 0.0039258460, 0.0039914470, 0.0040578080,
1447 0.0041249360, 0.0041928370, 0.0042615160, 0.0043309810, 0.0044012380, 0.0044722930,
1448 0.0045441520, 0.0046168220, 0.0046903100, 0.0047646230, 0.0048397670, 0.0049157480,
1449 0.0049925750, 0.0050702540, 0.0051487920, 0.0052281970, 0.0053084750, 0.0053896340,
1450 0.0054716810, 0.0055546240, 0.0056384710, 0.0057232280, 0.0058089050, 0.0058955090,
1451 0.0059830470, 0.0060715290, 0.0061609620, 0.0062513540, 0.0063427140, 0.0064350500,
1452 0.0065283720, 0.0066226870, 0.0067180050, 0.0068143350, 0.0069116850, 0.0070100650,
1453 0.0071094850, 0.0072099530, 0.0073114800, 0.0074140740, 0.0075177470, 0.0076225070,
1454 0.0077283660, 0.0078353330, 0.0079434180, 0.0080526330, 0.0081629880, 0.0082744940,
1455 0.0083871620, 0.0085010030, 0.0086160290, 0.0087322510, 0.0088496800, 0.0089683290,
1456 0.0090882100, 0.0092093350, 0.0093317160, 0.0094553660, 0.0095802980, 0.0097065250,
1457 0.0098340590, 0.0099629140, 0.0100931000, 0.0102246400, 0.0103575500, 0.0104918200,
1458 0.0106274900, 0.0107645700, 0.0109030700, 0.0110430000, 0.0111843800, 0.0113272300,
1459 0.0114715700, 0.0116174000, 0.0117647500, 0.0119136300, 0.0120640600, 0.0122160600,
1460 0.0123696300, 0.0125248100, 0.0126816100, 0.0128400500, 0.0130001400, 0.0131619100,
1461 0.0133253800, 0.0134905500, 0.0136574600, 0.0138261300, 0.0139965600, 0.0141688000,
1462 0.0143428500, 0.0145187400, 0.0146964900, 0.0148761200, 0.0150576500, 0.0152411200,
1463 0.0154265300, 0.0156139200, 0.0158033100, 0.0159947200, 0.0161881800, 0.0163837200,
1464 0.0165813500, 0.0167811100, 0.0169830300, 0.0171871200, 0.0173934200, 0.0176019600,
1465 0.0178127600, 0.0180258500, 0.0182412700, 0.0184590400, 0.0186791900, 0.0189017500,
1466 0.0191267600, 0.0193542500, 0.0195842500, 0.0198167900, 0.0200519000, 0.0202896300,
1467 0.0205300000, 0.0207730500, 0.0210188200, 0.0212673300, 0.0215186400, 0.0217727700,
1468 0.0220297700, 0.0222896700, 0.0225525200, 0.0228183400, 0.0230871900, 0.0233591000,
1469 0.0236341200, 0.0239122800, 0.0241936400, 0.0244782300, 0.0247661000, 0.0250572900,
1470 0.0253518600, 0.0256498400, 0.0259512800, 0.0262562300, 0.0265647400, 0.0268768600,
1471 0.0271926300, 0.0275121100, 0.0278353500, 0.0281624000, 0.0284933100, 0.0288281300,
1472 0.0291669200, 0.0295097300, 0.0298566200, 0.0302076400, 0.0305628500, 0.0309223000,
1473 0.0312860500, 0.0316541600, 0.0320266900, 0.0324036900, 0.0327852300, 0.0331713600,
1474 0.0335621500, 0.0339576500, 0.0343579300, 0.0347630500, 0.0351730600, 0.0355880400,
1475 0.0360080500, 0.0364331400, 0.0368633900, 0.0372988500, 0.0377395900, 0.0381856600,
1476 0.0386371500, 0.0390941000, 0.0395565900, 0.0400246700, 0.0404984100, 0.0409778800,
1477 0.0414631300, 0.0419542300, 0.0424512400, 0.0429542200, 0.0434632400, 0.0439783500,
1478 0.0444996200, 0.0450271000, 0.0455608500, 0.0461009400, 0.0466474000, 0.0472003200,
1479 0.0477597200, 0.0483256700, 0.0488982300, 0.0494774300, 0.0500633300, 0.0506559700,
1480 0.0512554000, 0.0518616600, 0.0524748000, 0.0530948400, 0.0537218300, 0.0543558000,
1481 0.0549967800, 0.0556448000, 0.0562998900, 0.0569620700, 0.0576313600, 0.0583077800,
1482 0.0589913500, 0.0596820600, 0.0603799400, 0.0610849900, 0.0617972100, 0.0625166000,
1483 0.0632431500, 0.0639768500, 0.0647176900, 0.0654656600, 0.0662207300, 0.0669828800,
1484 0.0677520900, 0.0685283100, 0.0693115200, 0.0701016800, 0.0708987400, 0.0717026500,
1485 0.0725133700, 0.0733308400, 0.0741550000, 0.0749858000, 0.0758231500, 0.0766670000,
1486 0.0775172700, 0.0783738800, 0.0792367600, 0.0801058200, 0.0809809800, 0.0818621400,
1487 0.0827492100, 0.0836421000, 0.0845407100, 0.0854449400, 0.0863546800, 0.0872698400,
1488 0.0881902900, 0.0891159400, 0.0900466700, 0.0909823700, 0.0919229100, 0.0928682000,
1489 0.0938181000, 0.0947725000, 0.0957312800, 0.0966943100, 0.0976614900, 0.0986326700,
1490 0.0996077500, 0.1005866000, 0.1015691000, 0.1025551000, 0.1035445000, 0.1045372000,
1491 0.1055331000, 0.1065320000, 0.1075338000, 0.1085384000, 0.1095456000, 0.1105555000,
1492 0.1115677000, 0.1125823000, 0.1135991000, 0.1146180000, 0.1156389000, 0.1166616000,
1493 0.1176862000, 0.1187123000, 0.1197401000, 0.1207692000, 0.1217997000, 0.1228315000,
1494 0.1238644000, 0.1248983000, 0.1259332000, 0.1269690000, 0.1280055000, 0.1290427000,
1495 0.1300805000, 0.1311188000, 0.1321575000, 0.1331966000, 0.1342359000, 0.1352754000,
1496 0.1363150000, 0.1373546000, 0.1383942000, 0.1394337000, 0.1404730000, 0.1415120000,
1497 0.1425508000, 0.1435891000, 0.1446271000, 0.1456645000, 0.1467014000, 0.1477377000,
1498 0.1487733000, 0.1498082000, 0.1508423000, 0.1518756000, 0.1529081000, 0.1539396000,
1499 0.1549702000, 0.1559998000, 0.1570284000, 0.1580558000, 0.1590822000, 0.1601074000,
1500 0.1611314000, 0.1621541000, 0.1631757000, 0.1641959000, 0.1652148000, 0.1662323000,
1501 0.1672484000, 0.1682632000, 0.1692765000, 0.1702883000, 0.1712986000, 0.1723075000,
1502 0.1733147000, 0.1743205000, 0.1753246000, 0.1763271000, 0.1773281000, 0.1783273000,
1503 0.1793250000, 0.1803209000, 0.1813152000, 0.1823078000, 0.1832986000, 0.1842877000,
1504 0.1852751000, 0.1862607000, 0.1872445000, 0.1882266000, 0.1892068000, 0.1901853000,
1505 0.1911619000, 0.1921368000, 0.1931098000, 0.1940809000, 0.1950503000, 0.1960177000,
1506 0.1969834000, 0.1979471000, 0.1989090000, 0.1998690000, 0.2008271000, 0.2017834000,
1507 0.2027378000, 0.2036903000, 0.2046408000, 0.2055895000, 0.2065363000, 0.2074812000,
1508 0.2084242000, 0.2093653000, 0.2103045000, 0.2112418000, 0.2121772000, 0.2131106000,
1509 0.2140422000, 0.2149719000, 0.2158996000, 0.2168255000, 0.2177494000, 0.2186715000,
1510 0.2195916000, 0.2205098000, 0.2214262000, 0.2223406000, 0.2232532000, 0.2241638000,
1511 0.2250726000, 0.2259795000, 0.2268844000, 0.2277875000, 0.2286888000, 0.2295881000,
1512 0.2304856000, 0.2313812000, 0.2322749000, 0.2331668000, 0.2340568000, 0.2349449000,
1513 0.2358312000, 0.2367157000, 0.2375983000, 0.2384790000, 0.2393580000, 0.2402351000,
1514 0.2411103000, 0.2419837000, 0.2428554000, 0.2437252000, 0.2445931000, 0.2454593000,
1515 0.2463237000, 0.2471863000, 0.2480471000, 0.2489061000, 0.2497633000, 0.2506187000,
1516 0.2514724000, 0.2523243000, 0.2531744000, 0.2540228000, 0.2548694000, 0.2557143000,
1517 0.2565574000, 0.2573988000, 0.2582385000, 0.2590764000, 0.2599126000, 0.2607471000,
1518 0.2615799000, 0.2624110000, 0.2632404000, 0.2640681000, 0.2648941000, 0.2657184000,
1519 0.2665410000, 0.2673620000, 0.2681813000, 0.2689989000, 0.2698149000, 0.2706292000,
1520 0.2714419000, 0.2722529000, 0.2730623000, 0.2738701000, 0.2746762000, 0.2754807000,
1521 0.2762836000, 0.2770849000, 0.2778846000, 0.2786827000, 0.2794792000, 0.2802741000,
1522 0.2810675000, 0.2818592000, 0.2826494000, 0.2834380000, 0.2842251000, 0.2850106000,
1523 0.2857946000, 0.2865770000, 0.2873578000, 0.2881372000, 0.2889150000, 0.2896913000,
1524 0.2904661000, 0.2912393000, 0.2920111000, 0.2927813000, 0.2935501000, 0.2943173000,
1525 0.2950831000, 0.2958474000, 0.2966102000, 0.2973715000, 0.2981314000, 0.2988898000,
1526 0.2996468000, 0.3004023000, 0.3011564000, 0.3019090000, 0.3026602000, 0.3034099000,
1527 0.3041583000, 0.3049052000, 0.3056507000, 0.3063948000, 0.3071374000, 0.3078787000,
1528 0.3086186000, 0.3093571000, 0.3100942000, 0.3108299000, 0.3115643000, 0.3122973000,
1529 0.3130289000, 0.3137591000, 0.3144880000, 0.3152155000, 0.3159417000, 0.3166666000,
1530 0.3173901000, 0.3181122000, 0.3188331000, 0.3195526000, 0.3202708000, 0.3209877000,
1531 0.3217033000, 0.3224175000, 0.3231305000, 0.3238421000, 0.3245525000, 0.3252616000,
1532 0.3259694000, 0.3266759000, 0.3273811000, 0.3280851000, 0.3287878000, 0.3294893000,
1533 0.3301894000, 0.3308884000, 0.3315861000, 0.3322825000, 0.3329777000, 0.3336716000,
1534 0.3343644000, 0.3350559000, 0.3357461000, 0.3364352000, 0.3371230000, 0.3378096000,
1535 0.3384951000, 0.3391793000, 0.3398623000, 0.3405441000, 0.3412247000, 0.3419042000,
1536 0.3425824000, 0.3432595000, 0.3439354000, 0.3446101000, 0.3452837000, 0.3459561000,
1537 0.3466273000, 0.3472974000, 0.3479663000, 0.3486341000, 0.3493008000, 0.3499663000,
1538 0.3506306000, 0.3512939000, 0.3519560000, 0.3526169000, 0.3532768000, 0.3539355000,
1539 0.3545931000, 0.3552497000, 0.3559051000, 0.3565594000, 0.3572126000, 0.3578647000,
1540 0.3585157000, 0.3591656000, 0.3598144000, 0.3604622000, 0.3611089000, 0.3617545000,
1541 0.3623990000, 0.3630425000, 0.3636849000, 0.3643263000, 0.3649665000, 0.3656058000,
1542 0.3662440000, 0.3668811000, 0.3675172000, 0.3681523000, 0.3687863000, 0.3694193000,
1543 0.3700512000, 0.3706822000, 0.3713121000, 0.3719409000, 0.3725688000, 0.3731957000,
1544 0.3738215000, 0.3744464000, 0.3750702000, 0.3756930000, 0.3763149000, 0.3769357000,
1545 0.3775556000, 0.3781745000, 0.3787924000, 0.3794093000, 0.3800252000, 0.3806401000,
1546 0.3812541000, 0.3818671000, 0.3824792000, 0.3830902000, 0.3837004000, 0.3843095000,
1547 0.3849178000, 0.3855250000, 0.3861313000, 0.3867367000, 0.3873411000, 0.3879446000,
1548 0.3885472000, 0.3891488000, 0.3897495000, 0.3903492000, 0.3909481000, 0.3915460000,
1549 0.3921430000, 0.3927391000, 0.3933343000, 0.3939285000, 0.3945219000, 0.3951143000,
1550 0.3957059000, 0.3962966000, 0.3968863000, 0.3974752000, 0.3980631000, 0.3986502000,
1551 0.3992364000, 0.3998218000, 0.4004062000, 0.4009898000, 0.4015725000, 0.4021543000,
1552 0.4027352000, 0.4033153000, 0.4038946000, 0.4044729000, 0.4050504000, 0.4056271000,
1553 0.4062029000, 0.4067778000, 0.4073519000, 0.4079252000, 0.4084976000, 0.4090692000,
1554 0.4096399000, 0.4102098000, 0.4107789000, 0.4113472000, 0.4119146000, 0.4124812000,
1555 0.4130470000, 0.4136119000, 0.4141760000, 0.4147394000, 0.4153019000, 0.4158636000,
1556 0.4164245000, 0.4169846000, 0.4175439000, 0.4181024000, 0.4186601000, 0.4192170000,
1557 0.4197731000, 0.4203284000, 0.4208830000, 0.4214367000, 0.4219897000, 0.4225419000,
1558 0.4230933000, 0.4236439000, 0.4241938000, 0.4247429000, 0.4252912000, 0.4258388000,
1559 0.4263856000, 0.4269317000, 0.4274769000, 0.4280215000, 0.4285652000, 0.4291083000,
1560 0.4296505000, 0.4301921000, 0.4307329000, 0.4312729000, 0.4318122000, 0.4323508000,
1561 0.4328886000, 0.4334257000, 0.4339621000, 0.4344977000, 0.4350326000, 0.4355668000,
1562 0.4361003000, 0.4366331000, 0.4371651000, 0.4376964000, 0.4382270000, 0.4387569000,
1563 0.4392861000, 0.4398146000, 0.4403424000, 0.4408695000, 0.4413959000, 0.4419216000,
1564 0.4424466000, 0.4429709000, 0.4434945000, 0.4440174000, 0.4445396000, 0.4450612000,
1565 0.4455821000, 0.4461023000, 0.4466218000, 0.4471407000, 0.4476589000, 0.4481764000,
1566 0.4486932000, 0.4492094000, 0.4497249000, 0.4502398000, 0.4507540000, 0.4512675000,
1567 0.4517804000, 0.4522926000, 0.4528042000, 0.4533152000, 0.4538255000, 0.4543351000,
1568 0.4548442000, 0.4553525000, 0.4558603000, 0.4563674000, 0.4568739000, 0.4573797000,
1569 0.4578849000, 0.4583895000, 0.4588935000, 0.4593969000, 0.4598996000, 0.4604018000,
1570 0.4609033000, 0.4614042000, 0.4619045000, 0.4624042000, 0.4629032000, 0.4634017000,
1571 0.4638996000, 0.4643969000, 0.4648936000, 0.4653897000, 0.4658852000, 0.4663802000,
1572 0.4668745000, 0.4673683000, 0.4678614000, 0.4683541000, 0.4688461000, 0.4693375000,
1573 0.4698284000, 0.4703188000, 0.4708085000, 0.4712977000, 0.4717863000, 0.4722744000,
1574 0.4727619000, 0.4732489000, 0.4737353000, 0.4742212000, 0.4747065000, 0.4751913000,
1575 0.4756756000, 0.4761593000, 0.4766424000, 0.4771251000, 0.4776072000, 0.4780888000,
1576 0.4785699000, 0.4790504000, 0.4795304000, 0.4800099000, 0.4804889000, 0.4809674000,
1577 0.4814454000, 0.4819229000, 0.4823999000, 0.4828764000, 0.4833524000, 0.4838279000,
1578 0.4843029000, 0.4847774000, 0.4852515000, 0.4857250000, 0.4861981000, 0.4866708000,
1579 0.4871429000, 0.4876146000, 0.4880858000, 0.4885566000, 0.4890269000, 0.4894968000,
1580 0.4899662000, 0.4904351000, 0.4909037000, 0.4913718000, 0.4918394000, 0.4923067000,
1581 0.4927735000, 0.4932398000, 0.4937058000, 0.4941713000, 0.4946365000, 0.4951012000,
1582 0.4955655000, 0.4960295000, 0.4964930000, 0.4969562000, 0.4974189000, 0.4978813000,
1583 0.4983434000, 0.4988050000, 0.4992663000, 0.4997272000, 0.5001878000, 0.5006480000,
1584 0.5011079000, 0.5015674000, 0.5020266000, 0.5024855000, 0.5029440000, 0.5034023000,
1585 0.5038602000, 0.5043178000, 0.5047751000, 0.5052321000, 0.5056889000, 0.5061453000,
1586 0.5066015000, 0.5070574000, 0.5075131000, 0.5079685000, 0.5084237000, 0.5088786000,
1587 0.5093333000, 0.5097878000, 0.5102420000, 0.5106961000, 0.5111499000, 0.5116036000,
1588 0.5120571000, 0.5125104000, 0.5129636000, 0.5134166000, 0.5138694000, 0.5143221000,
1589 0.5147747000, 0.5152272000, 0.5156796000, 0.5161319000, 0.5165842000, 0.5170363000,
1590 0.5174884000, 0.5179405000, 0.5183925000, 0.5188445000, 0.5192965000, 0.5197485000,
1591 0.5202005000, 0.5206526000, 0.5211048000, 0.5215570000, 0.5220092000, 0.5224616000,
1592 0.5229141000, 0.5233668000, 0.5238195000, 0.5242725000, 0.5247257000, 0.5251790000,
1593 0.5256326000, 0.5260864000, 0.5265405000, 0.5269949000, 0.5274496000, 0.5279047000,
1594 0.5283601000, 0.5288159000, 0.5292721000, 0.5297287000, 0.5301858000, 0.5306434000,
1595 0.5311015000, 0.5315601000, 0.5320194000, 0.5324792000, 0.5329396000, 0.5334008000,
1596 0.5338626000, 0.5343252000, 0.5347885000, 0.5352526000, 0.5357176000, 0.5361834000,
1597 0.5366502000, 0.5371179000, 0.5375866000, 0.5380563000, 0.5385271000, 0.5389989000,
1598 0.5394720000, 0.5399462000, 0.5404217000, 0.5408984000, 0.5413764000, 0.5418558000,
1599 0.5423367000, 0.5428189000, 0.5433026000, 0.5437879000, 0.5442747000, 0.5447631000,
1600 0.5452532000, 0.5457450000, 0.5462385000, 0.5467337000, 0.5472307000, 0.5477296000,
1601 0.5482303000, 0.5487328000, 0.5492373000, 0.5497437000, 0.5502520000, 0.5507623000,
1602 0.5512745000, 0.5517887000, 0.5523048000, 0.5528230000, 0.5533430000, 0.5538651000,
1603 0.5543890000, 0.5549149000, 0.5554426000, 0.5559722000, 0.5565036000, 0.5570367000,
1604 0.5575716000, 0.5581082000, 0.5586464000, 0.5591862000, 0.5597274000, 0.5602702000,
1605 0.5608143000, 0.5613597000, 0.5619063000, 0.5624541000, 0.5630030000, 0.5635529000,
1606 0.5641037000, 0.5646553000, 0.5652078000, 0.5657609000, 0.5663146000, 0.5668689000,
1607 0.5674236000, 0.5679787000, 0.5685341000, 0.5690898000, 0.5696456000, 0.5702015000,
1608 0.5707575000, 0.5713134000, 0.5718692000, 0.5724249000, 0.5729803000, 0.5735355000,
1609 0.5740903000, 0.5746448000, 0.5751989000, 0.5757525000, 0.5763056000, 0.5768581000,
1610 0.5774101000, 0.5779615000, 0.5785121000, 0.5790622000, 0.5796114000, 0.5801600000,
1611 0.5807078000, 0.5812548000, 0.5818009000, 0.5823462000, 0.5828907000, 0.5834343000,
1612 0.5839770000, 0.5845188000, 0.5850597000, 0.5855996000, 0.5861386000, 0.5866767000,
1613 0.5872137000, 0.5877498000, 0.5882850000, 0.5888191000, 0.5893522000, 0.5898844000,
1614 0.5904155000, 0.5909457000, 0.5914748000, 0.5920029000, 0.5925300000, 0.5930560000,
1615 0.5935811000, 0.5941051000, 0.5946282000, 0.5951502000, 0.5956711000, 0.5961911000,
1616 0.5967100000, 0.5972280000, 0.5977449000, 0.5982608000, 0.5987757000, 0.5992896000,
1617 0.5998025000, 0.6003144000, 0.6008253000, 0.6013352000, 0.6018441000, 0.6023520000,
1618 0.6028589000, 0.6033649000, 0.6038699000, 0.6043739000, 0.6048769000, 0.6053790000,
1619 0.6058801000, 0.6063803000, 0.6068795000, 0.6073777000, 0.6078750000, 0.6083714000,
1620 0.6088668000, 0.6093613000, 0.6098549000, 0.6103476000, 0.6108393000, 0.6113301000,
1621 0.6118201000, 0.6123091000, 0.6127972000, 0.6132844000, 0.6137707000, 0.6142562000,
1622 0.6147407000, 0.6152244000, 0.6157072000, 0.6161892000, 0.6166703000, 0.6171505000,
1623 0.6176298000, 0.6181084000, 0.6185860000, 0.6190629000, 0.6195388000, 0.6200140000,
1624 0.6204883000, 0.6209618000, 0.6214345000, 0.6219064000, 0.6223774000, 0.6228476000,
1625 0.6233171000, 0.6237857000, 0.6242536000, 0.6247206000, 0.6251869000, 0.6256523000,
1626 0.6261170000, 0.6265809000, 0.6270441000, 0.6275064000, 0.6279681000, 0.6284289000,
1627 0.6288890000, 0.6293483000, 0.6298069000, 0.6302647000, 0.6307218000, 0.6311782000,
1628 0.6316338000, 0.6320887000, 0.6325428000, 0.6329963000, 0.6334490000, 0.6339009000,
1629 0.6343522000, 0.6348028000, 0.6352526000, 0.6357018000, 0.6361502000, 0.6365979000,
1630 0.6370450000, 0.6374913000, 0.6379370000, 0.6383820000, 0.6388262000, 0.6392699000,
1631 0.6397128000, 0.6401550000, 0.6405966000, 0.6410375000, 0.6414778000, 0.6419174000,
1632 0.6423563000, 0.6427946000, 0.6432322000, 0.6436692000, 0.6441055000, 0.6445412000,
1633 0.6449762000, 0.6454106000, 0.6458443000, 0.6462774000, 0.6467099000, 0.6471418000,
1634 0.6475730000, 0.6480036000, 0.6484336000, 0.6488629000, 0.6492917000, 0.6497198000,
1635 0.6501473000, 0.6505742000, 0.6510005000, 0.6514262000, 0.6518513000, 0.6522758000,
1636 0.6526997000, 0.6531230000, 0.6535457000, 0.6539678000, 0.6543893000, 0.6548103000,
1637 0.6552306000, 0.6556504000, 0.6560696000, 0.6564882000, 0.6569062000, 0.6573237000,
1638 0.6577406000, 0.6581569000, 0.6585727000, 0.6589879000, 0.6594025000, 0.6598166000,
1639 0.6602301000, 0.6606431000, 0.6610555000, 0.6614673000, 0.6618787000, 0.6622894000,
1640 0.6626996000, 0.6631093000, 0.6635184000, 0.6639270000, 0.6643350000, 0.6647425000,
1641 0.6651495000, 0.6655560000, 0.6659619000, 0.6663672000, 0.6667721000, 0.6671764000,
1642 0.6675802000, 0.6679835000, 0.6683863000, 0.6687885000, 0.6691902000, 0.6695914000,
1643 0.6699921000, 0.6703923000, 0.6707920000, 0.6711912000, 0.6715898000, 0.6719880000,
1644 0.6723856000, 0.6727828000, 0.6731794000, 0.6735756000, 0.6739712000, 0.6743664000,
1645 0.6747610000, 0.6751552000, 0.6755489000, 0.6759421000, 0.6763348000, 0.6767270000,
1646 0.6771187000, 0.6775100000, 0.6779008000, 0.6782910000, 0.6786809000, 0.6790702000,
1647 0.6794590000, 0.6798474000, 0.6802353000, 0.6806228000, 0.6810097000, 0.6813962000,
1648 0.6817823000, 0.6821678000, 0.6825530000, 0.6829376000, 0.6833218000, 0.6837055000,
1649 0.6840888000, 0.6844716000, 0.6848539000, 0.6852358000, 0.6856172000, 0.6859982000,
1650 0.6863788000, 0.6867588000, 0.6871385000, 0.6875177000, 0.6878964000, 0.6882747000,
1651 0.6886526000, 0.6890300000, 0.6894070000, 0.6897835000, 0.6901596000, 0.6905353000,
1652 0.6909105000, 0.6912853000, 0.6916596000, 0.6920336000, 0.6924071000, 0.6927801000,
1653 0.6931528000, 0.6935250000, 0.6938967000, 0.6942681000, 0.6946390000, 0.6950095000,
1654 0.6953796000, 0.6957493000, 0.6961185000, 0.6964874000, 0.6968558000, 0.6972238000,
1655 0.6975914000, 0.6979585000, 0.6983253000, 0.6986916000, 0.6990576000, 0.6994231000,
1656 0.6997882000, 0.7001529000, 0.7005172000, 0.7008811000, 0.7012446000, 0.7016077000,
1657 0.7019704000, 0.7023326000, 0.7026945000, 0.7030560000, 0.7034171000, 0.7037778000,
1658 0.7041381000, 0.7044980000, 0.7048575000, 0.7052166000, 0.7055753000, 0.7059336000,
1659 0.7062916000, 0.7066491000, 0.7070063000, 0.7073630000, 0.7077194000, 0.7080754000,
1660 0.7084310000, 0.7087863000, 0.7091411000, 0.7094956000, 0.7098497000, 0.7102034000,
1661 0.7105567000, 0.7109096000, 0.7112622000, 0.7116144000, 0.7119662000, 0.7123177000,
1662 0.7126687000, 0.7130194000, 0.7133697000, 0.7137197000, 0.7140693000, 0.7144185000,
1663 0.7147673000, 0.7151158000, 0.7154639000, 0.7158117000, 0.7161591000, 0.7165061000,
1664 0.7168527000, 0.7171990000, 0.7175450000, 0.7178905000, 0.7182357000, 0.7185806000,
1665 0.7189251000, 0.7192692000, 0.7196130000, 0.7199564000, 0.7202995000, 0.7206422000,
1666 0.7209845000, 0.7213266000, 0.7216682000, 0.7220095000, 0.7223505000, 0.7226911000,
1667 0.7230313000, 0.7233712000, 0.7237108000, 0.7240500000, 0.7243889000, 0.7247274000,
1668 0.7250656000, 0.7254034000, 0.7257409000, 0.7260780000, 0.7264148000, 0.7267513000,
1669 0.7270874000, 0.7274232000, 0.7277587000, 0.7280938000, 0.7284286000, 0.7287630000,
1670 0.7290971000, 0.7294309000, 0.7297643000, 0.7300974000, 0.7304302000, 0.7307626000,
1671 0.7310947000, 0.7314265000, 0.7317579000, 0.7320891000, 0.7324198000, 0.7327503000,
1672 0.7330804000, 0.7334102000, 0.7337397000, 0.7340689000, 0.7343977000, 0.7347262000,
1673 0.7350544000, 0.7353823000, 0.7357098000, 0.7360370000, 0.7363639000, 0.7366905000,
1674 0.7370168000, 0.7373427000, 0.7376683000, 0.7379936000, 0.7383186000, 0.7386433000,
1675 0.7389677000, 0.7392917000, 0.7396155000, 0.7399389000, 0.7402620000, 0.7405848000,
1676 0.7409073000, 0.7412294000, 0.7415513000, 0.7418728000, 0.7421941000, 0.7425150000,
1677 0.7428357000, 0.7431560000, 0.7434760000, 0.7437957000, 0.7441151000, 0.7444342000,
1678 0.7447530000, 0.7450715000, 0.7453897000, 0.7457076000, 0.7460251000, 0.7463424000,
1679 0.7466594000, 0.7469761000, 0.7472925000, 0.7476085000, 0.7479243000, 0.7482398000,
1680 0.7485550000, 0.7488699000, 0.7491845000, 0.7494988000, 0.7498128000, 0.7501265000,
1681 0.7504399000, 0.7507531000, 0.7510659000, 0.7513784000, 0.7516907000, 0.7520027000,
1682 0.7523143000, 0.7526257000, 0.7529368000, 0.7532476000, 0.7535581000, 0.7538683000,
1683 0.7541783000, 0.7544879000, 0.7547973000, 0.7551064000, 0.7554152000, 0.7557237000,
1684 0.7560319000, 0.7563398000, 0.7566475000, 0.7569549000, 0.7572620000, 0.7575688000,
1685 0.7578753000, 0.7581816000, 0.7584875000, 0.7587932000, 0.7590986000, 0.7594038000,
1686 0.7597086000, 0.7600132000, 0.7603175000, 0.7606215000, 0.7609253000, 0.7612288000,
1687 0.7615320000, 0.7618349000, 0.7621376000, 0.7624399000, 0.7627420000, 0.7630439000,
1688 0.7633454000, 0.7636467000, 0.7639477000, 0.7642485000, 0.7645490000, 0.7648492000,
1689 0.7651491000, 0.7654488000, 0.7657482000, 0.7660473000, 0.7663462000, 0.7666448000,
1690 0.7669431000, 0.7672412000, 0.7675390000, 0.7678366000, 0.7681338000, 0.7684309000,
1691 0.7687276000, 0.7690241000, 0.7693203000, 0.7696163000, 0.7699120000, 0.7702074000,
1692 0.7705026000, 0.7707975000, 0.7710922000, 0.7713866000, 0.7716807000, 0.7719746000,
1693 0.7722682000, 0.7725616000, 0.7728547000, 0.7731475000, 0.7734401000, 0.7737325000,
1694 0.7740245000, 0.7743164000, 0.7746079000, 0.7748993000, 0.7751903000, 0.7754811000,
1695 0.7757717000, 0.7760620000, 0.7763521000, 0.7766419000, 0.7769314000, 0.7772207000,
1696 0.7775098000, 0.7777986000, 0.7780871000, 0.7783754000, 0.7786635000, 0.7789513000,
1697 0.7792389000, 0.7795262000, 0.7798132000, 0.7801000000, 0.7803866000, 0.7806729000,
1698 0.7809590000, 0.7812449000, 0.7815304000, 0.7818158000, 0.7821009000, 0.7823858000,
1699 0.7826704000, 0.7829547000, 0.7832389000, 0.7835228000, 0.7838064000, 0.7840898000,
1700 0.7843730000, 0.7846559000, 0.7849386000, 0.7852210000, 0.7855032000, 0.7857852000,
1701 0.7860669000, 0.7863484000, 0.7866297000, 0.7869107000, 0.7871915000, 0.7874720000,
1702 0.7877523000, 0.7880324000, 0.7883122000, 0.7885919000, 0.7888712000, 0.7891504000,
1703 0.7894293000, 0.7897079000, 0.7899864000, 0.7902646000, 0.7905425000, 0.7908203000,
1704 0.7910978000, 0.7913750000, 0.7916521000, 0.7919289000, 0.7922055000, 0.7924818000,
1705 0.7927579000, 0.7930338000, 0.7933095000, 0.7935849000, 0.7938601000, 0.7941351000,
1706 0.7944099000, 0.7946844000, 0.7949587000, 0.7952328000, 0.7955066000, 0.7957803000,
1707 0.7960537000, 0.7963268000, 0.7965998000, 0.7968725000, 0.7971450000, 0.7974173000,
1708 0.7976893000, 0.7979612000, 0.7982328000, 0.7985041000, 0.7987753000, 0.7990463000,
1709 0.7993170000, 0.7995875000, 0.7998578000, 0.8001278000, 0.8003977000, 0.8006673000,
1710 0.8009367000, 0.8012059000, 0.8014748000, 0.8017436000, 0.8020121000, 0.8022804000,
1711 0.8025485000, 0.8028164000, 0.8030841000, 0.8033515000, 0.8036187000, 0.8038857000,
1712 0.8041525000, 0.8044191000, 0.8046855000, 0.8049517000, 0.8052176000, 0.8054833000,
1713 0.8057488000, 0.8060141000, 0.8062792000, 0.8065441000, 0.8068088000, 0.8070732000,
1714 0.8073375000, 0.8076015000, 0.8078653000, 0.8081289000, 0.8083923000, 0.8086555000,
1715 0.8089185000, 0.8091813000, 0.8094438000, 0.8097062000, 0.8099683000, 0.8102303000,
1716 0.8104920000, 0.8107535000, 0.8110148000, 0.8112760000, 0.8115369000, 0.8117976000,
1717 0.8120580000, 0.8123183000, 0.8125784000, 0.8128383000, 0.8130980000, 0.8133574000,
1718 0.8136167000, 0.8138757000, 0.8141346000, 0.8143933000, 0.8146517000, 0.8149099000,
1719 0.8151680000, 0.8154258000, 0.8156835000, 0.8159409000, 0.8161982000, 0.8164552000,
1720 0.8167120000, 0.8169687000, 0.8172251000, 0.8174813000, 0.8177374000, 0.8179932000,
1721 0.8182488000, 0.8185043000, 0.8187595000, 0.8190146000, 0.8192694000, 0.8195241000,
1722 0.8197785000, 0.8200328000, 0.8202868000, 0.8205407000, 0.8207944000, 0.8210478000,
1723 0.8213011000, 0.8215542000, 0.8218071000, 0.8220598000, 0.8223122000, 0.8225645000,
1724 0.8228167000, 0.8230686000, 0.8233203000, 0.8235718000, 0.8238232000, 0.8240743000,
1725 0.8243252000, 0.8245760000, 0.8248266000, 0.8250769000, 0.8253271000, 0.8255771000,
1726 0.8258269000, 0.8260765000, 0.8263259000, 0.8265752000, 0.8268242000, 0.8270731000,
1727 0.8273217000, 0.8275702000, 0.8278185000, 0.8280666000, 0.8283145000, 0.8285622000,
1728 0.8288097000, 0.8290571000, 0.8293043000, 0.8295512000, 0.8297980000, 0.8300446000,
1729 0.8302910000, 0.8305373000, 0.8307833000, 0.8310292000, 0.8312748000, 0.8315203000,
1730 0.8317656000, 0.8320107000, 0.8322557000, 0.8325004000, 0.8327450000, 0.8329894000,
1731 0.8332336000, 0.8334776000, 0.8337214000, 0.8339651000, 0.8342085000, 0.8344518000,
1732 0.8346949000 };
1733 double width_a1 = a1_width[iii];
1734 a[0] = 1;
1735 a[1] = 0;
1736 b[0] = mass2 - sa;
1737 // b[1] = -mass*width*wid(mass2,mass,sa,sb,sc,r2,l);
1738 b[1] = -mass * width_a1;
1739 Com_Divide( a, b, prop );
1740}
1741void EvtD0ToKSKppimpi0::propagatorRBW_K1( double mass2, double mass, double width, double sa,
1742 double sb, double sc, double r2, int l,
1743 double prop[2] ) {
1744 double a[2], b[2];
1745 int iii = int( sa * 1000 / 2 ) - 300; //(4000-0)/2000=2
1746 if ( iii < 0 ) { iii = 0; }
1747 double K1270_width[2000] = {
1748 // 1700 bin, s = (0.6, 4)
1749 0.0000000014, 0.0000000052, 0.0000000127, 0.0000000252, 0.0000000442, 0.0000000709,
1750 0.0000001067, 0.0000001528, 0.0000002108, 0.0000002819, 0.0000003675, 0.0000004689,
1751 0.0000005876, 0.0000007249, 0.0000008822, 0.0000010609, 0.0000012624, 0.0000014882,
1752 0.0000017396, 0.0000020181, 0.0000023251, 0.0000026620, 0.0000030303, 0.0000034314,
1753 0.0000038668, 0.0000043379, 0.0000048463, 0.0000053933, 0.0000059805, 0.0000066094,
1754 0.0000072815, 0.0000079983, 0.0000087612, 0.0000095719, 0.0000104318, 0.0000113424,
1755 0.0000123054, 0.0000133223, 0.0000143946, 0.0000155239, 0.0000167117, 0.0000179597,
1756 0.0000192695, 0.0000206426, 0.0000220807, 0.0000235853, 0.0000251581, 0.0000268008,
1757 0.0000285149, 0.0000303022, 0.0000321643, 0.0000341028, 0.0000361195, 0.0000382160,
1758 0.0000403941, 0.0000426554, 0.0000450017, 0.0000474347, 0.0000499561, 0.0000525678,
1759 0.0000552715, 0.0000580689, 0.0000609619, 0.0000639522, 0.0000670417, 0.0000702323,
1760 0.0000735256, 0.0000769237, 0.0000804283, 0.0000840414, 0.0000877648, 0.0000916004,
1761 0.0000955503, 0.0000996163, 0.0001038003, 0.0001081045, 0.0001125306, 0.0001170808,
1762 0.0001217571, 0.0001265615, 0.0001314961, 0.0001365630, 0.0001417641, 0.0001471018,
1763 0.0001525781, 0.0001581951, 0.0001639551, 0.0001698602, 0.0001759127, 0.0001821148,
1764 0.0001884687, 0.0001949769, 0.0002016415, 0.0002084650, 0.0002154497, 0.0002225981,
1765 0.0002299124, 0.0002373952, 0.0002450491, 0.0002528764, 0.0002608798, 0.0002690618,
1766 0.0002774250, 0.0002859721, 0.0002947058, 0.0003036288, 0.0003127438, 0.0003220537,
1767 0.0003315613, 0.0003412695, 0.0003511812, 0.0003612993, 0.0003716269, 0.0003821671,
1768 0.0003929230, 0.0004038976, 0.0004150944, 0.0004265164, 0.0004381671, 0.0004500499,
1769 0.0004621681, 0.0004745254, 0.0004871253, 0.0004999714, 0.0005130674, 0.0005264172,
1770 0.0005400246, 0.0005538936, 0.0005680282, 0.0005824325, 0.0005971108, 0.0006120672,
1771 0.0006273063, 0.0006428325, 0.0006586504, 0.0006747647, 0.0006911803, 0.0007079021,
1772 0.0007249351, 0.0007422845, 0.0007599558, 0.0007779542, 0.0007962855, 0.0008149555,
1773 0.0008339699, 0.0008533351, 0.0008730571, 0.0008931425, 0.0009135979, 0.0009344302,
1774 0.0009556463, 0.0009772537, 0.0009992598, 0.0010216720, 0.0010444990, 0.0010677490,
1775 0.0010914300, 0.0011155510, 0.0011401210, 0.0011651510, 0.0011906480, 0.0012166250,
1776 0.0012430910, 0.0012700570, 0.0012975350, 0.0013255370, 0.0013540740, 0.0013831600,
1777 0.0014128090, 0.0014430330, 0.0014738470, 0.0015052670, 0.0015373090, 0.0015699870,
1778 0.0016033210, 0.0016373270, 0.0016720240, 0.0017074330, 0.0017435730, 0.0017804660,
1779 0.0018181340, 0.0018566010, 0.0018958920, 0.0019360320, 0.0019770490, 0.0020189710,
1780 0.0020618280, 0.0021056510, 0.0021504740, 0.0021963300, 0.0022432570, 0.0022912920,
1781 0.0023404760, 0.0023908500, 0.0024424590, 0.0024953480, 0.0025495680, 0.0026051670,
1782 0.0026622000, 0.0027207220, 0.0027807900, 0.0028424670, 0.0029058130, 0.0029708960,
1783 0.0030377810, 0.0031065390, 0.0031772420, 0.0032499620, 0.0033247750, 0.0034017560,
1784 0.0034809810, 0.0035625250, 0.0036464650, 0.0037328720, 0.0038218190, 0.0039133710,
1785 0.0040075910, 0.0041045360, 0.0042042540, 0.0043067850, 0.0044121600, 0.0045203990,
1786 0.0046315090, 0.0047454840, 0.0048623040, 0.0049819360, 0.0051043310, 0.0052294280,
1787 0.0053571500, 0.0054874080, 0.0056201030, 0.0057551230, 0.0058923500, 0.0060316590,
1788 0.0061729200, 0.0063159990, 0.0064607640, 0.0066070800, 0.0067548180, 0.0069038480,
1789 0.0070540500, 0.0072053060, 0.0073575050, 0.0075105450, 0.0076643270, 0.0078187630,
1790 0.0079737720, 0.0081292780, 0.0082852130, 0.0084415170, 0.0085981350, 0.0087550200,
1791 0.0089121270, 0.0090694200, 0.0092268660, 0.0093844380, 0.0095421120, 0.0096998670,
1792 0.0098576870, 0.0100155600, 0.0101734700, 0.0103314200, 0.0104894000, 0.0106474000,
1793 0.0108054200, 0.0109634700, 0.0111215500, 0.0112796500, 0.0114377900, 0.0115959700,
1794 0.0117541900, 0.0119124600, 0.0120707900, 0.0122292000, 0.0123876800, 0.0125462500,
1795 0.0127049200, 0.0128637100, 0.0130226100, 0.0131816400, 0.0133408200, 0.0135001600,
1796 0.0136596700, 0.0138193500, 0.0139792400, 0.0141393300, 0.0142996400, 0.0144601800,
1797 0.0146209700, 0.0147820200, 0.0149433400, 0.0151049500, 0.0152668600, 0.0154290800,
1798 0.0155916300, 0.0157545200, 0.0159177700, 0.0160813900, 0.0162453900, 0.0164097800,
1799 0.0165745900, 0.0167398200, 0.0169054800, 0.0170716000, 0.0172381900, 0.0174052600,
1800 0.0175728200, 0.0177408900, 0.0179094800, 0.0180786100, 0.0182483000, 0.0184185500,
1801 0.0185893800, 0.0187608000, 0.0189328400, 0.0191054900, 0.0192787900, 0.0194527400,
1802 0.0196273600, 0.0198026700, 0.0199786700, 0.0201553900, 0.0203328300, 0.0205110200,
1803 0.0206899700, 0.0208696900, 0.0210502100, 0.0212315300, 0.0214136700, 0.0215966500,
1804 0.0217804800, 0.0219651900, 0.0221507800, 0.0223372700, 0.0225246900, 0.0227130400,
1805 0.0229023400, 0.0230926100, 0.0232838700, 0.0234761300, 0.0236694200, 0.0238637400,
1806 0.0240591200, 0.0242555800, 0.0244531300, 0.0246517900, 0.0248515800, 0.0250525200,
1807 0.0252546300, 0.0254579200, 0.0256624200, 0.0258681400, 0.0260751100, 0.0262833500,
1808 0.0264928700, 0.0267036900, 0.0269158400, 0.0271293400, 0.0273442100, 0.0275604600,
1809 0.0277781300, 0.0279972300, 0.0282177800, 0.0284398100, 0.0286633400, 0.0288884000,
1810 0.0291149900, 0.0293431600, 0.0295729200, 0.0298042900, 0.0300373100, 0.0302719900,
1811 0.0305083600, 0.0307464500, 0.0309862700, 0.0312278600, 0.0314712500, 0.0317164500,
1812 0.0319635000, 0.0322124200, 0.0324632400, 0.0327159900, 0.0329706900, 0.0332273800,
1813 0.0334860800, 0.0337468200, 0.0340096300, 0.0342745400, 0.0345415800, 0.0348107800,
1814 0.0350821700, 0.0353557800, 0.0356316400, 0.0359097900, 0.0361902600, 0.0364730800,
1815 0.0367582700, 0.0370458900, 0.0373359500, 0.0376284900, 0.0379235500, 0.0382211500,
1816 0.0385213500, 0.0388241600, 0.0391296300, 0.0394377800, 0.0397486700, 0.0400623100,
1817 0.0403787600, 0.0406980400, 0.0410202000, 0.0413452600, 0.0416732700, 0.0420042700,
1818 0.0423382900, 0.0426753700, 0.0430155500, 0.0433588700, 0.0437053700, 0.0440550800,
1819 0.0444080500, 0.0447643100, 0.0451239000, 0.0454868700, 0.0458532400, 0.0462230700,
1820 0.0465963900, 0.0469732300, 0.0473536500, 0.0477376700, 0.0481253400, 0.0485166900,
1821 0.0489117600, 0.0493106000, 0.0497132400, 0.0501197100, 0.0505300600, 0.0509443200,
1822 0.0513625300, 0.0517847200, 0.0522109400, 0.0526412100, 0.0530755700, 0.0535140600,
1823 0.0539567000, 0.0544035400, 0.0548545900, 0.0553099000, 0.0557694900, 0.0562334000,
1824 0.0567016400, 0.0571742500, 0.0576512600, 0.0581326800, 0.0586185500, 0.0591088800,
1825 0.0596036900, 0.0601030200, 0.0606068600, 0.0611152500, 0.0616281900, 0.0621457100,
1826 0.0626678100, 0.0631945000, 0.0637258000, 0.0642617000, 0.0648022200, 0.0653473700,
1827 0.0658971300, 0.0664515200, 0.0670105300, 0.0675741500, 0.0681423900, 0.0687152300,
1828 0.0692926600, 0.0698746700, 0.0704612600, 0.0710524000, 0.0716480700, 0.0722482600,
1829 0.0728529400, 0.0734620900, 0.0740756900, 0.0746937100, 0.0753161100, 0.0759428600,
1830 0.0765739400, 0.0772093000, 0.0778489100, 0.0784927300, 0.0791407200, 0.0797928300,
1831 0.0804490200, 0.0811092400, 0.0817734600, 0.0824416000, 0.0831136300, 0.0837895000,
1832 0.0844691400, 0.0851525100, 0.0858395400, 0.0865301900, 0.0872243800, 0.0879220700,
1833 0.0886231900, 0.0893276700, 0.0900354700, 0.0907465100, 0.0914607200, 0.0921780600,
1834 0.0928984400, 0.0936218200, 0.0943481100, 0.0950772500, 0.0958091900, 0.1023867000,
1835 0.1055505000, 0.1081475000, 0.1104523000, 0.1125711000, 0.1145581000, 0.1164458000,
1836 0.1182550000, 0.1200004000, 0.1216926000, 0.1233395000, 0.1249474000, 0.1265212000,
1837 0.1280649000, 0.1295817000, 0.1310743000, 0.1325452000, 0.1339961000, 0.1354288000,
1838 0.1368448000, 0.1382454000, 0.1396316000, 0.1410044000, 0.1423647000, 0.1437133000,
1839 0.1450510000, 0.1463782000, 0.1476956000, 0.1490038000, 0.1503031000, 0.1515941000,
1840 0.1528770000, 0.1541523000, 0.1554203000, 0.1566812000, 0.1579354000, 0.1591832000,
1841 0.1604247000, 0.1616602000, 0.1628899000, 0.1641140000, 0.1653326000, 0.1665460000,
1842 0.1677543000, 0.1689576000, 0.1701561000, 0.1713499000, 0.1725391000, 0.1737239000,
1843 0.1749043000, 0.1760805000, 0.1772526000, 0.1784206000, 0.1795847000, 0.1807449000,
1844 0.1819014000, 0.1830541000, 0.1842031000, 0.1853487000, 0.1864907000, 0.1876292000,
1845 0.1887644000, 0.1898963000, 0.1910250000, 0.1921504000, 0.1932727000, 0.1943919000,
1846 0.1955080000, 0.1966211000, 0.1977313000, 0.1988386000, 0.1999430000, 0.2010446000,
1847 0.2021434000, 0.2032395000, 0.2043329000, 0.2054236000, 0.2065117000, 0.2075972000,
1848 0.2086801000, 0.2097605000, 0.2108384000, 0.2119139000, 0.2129870000, 0.2140576000,
1849 0.2151259000, 0.2161919000, 0.2172555000, 0.2183169000, 0.2193761000, 0.2204330000,
1850 0.2214878000, 0.2225404000, 0.2235909000, 0.2246392000, 0.2256855000, 0.2267297000,
1851 0.2277719000, 0.2288121000, 0.2298503000, 0.2308866000, 0.2319210000, 0.2329534000,
1852 0.2339840000, 0.2350126000, 0.2360395000, 0.2370646000, 0.2380878000, 0.2391093000,
1853 0.2401291000, 0.2411471000, 0.2421634000, 0.2431781000, 0.2441911000, 0.2452024000,
1854 0.2462122000, 0.2472203000, 0.2482269000, 0.2492319000, 0.2502354000, 0.2512374000,
1855 0.2522379000, 0.2532369000, 0.2542345000, 0.2552306000, 0.2562253000, 0.2572187000,
1856 0.2582106000, 0.2592012000, 0.2601905000, 0.2611784000, 0.2621651000, 0.2631504000,
1857 0.2641345000, 0.2651174000, 0.2660991000, 0.2670795000, 0.2680587000, 0.2690368000,
1858 0.2700137000, 0.2709895000, 0.2719642000, 0.2729378000, 0.2739103000, 0.2748817000,
1859 0.2758521000, 0.2768215000, 0.2777898000, 0.2787572000, 0.2797235000, 0.2806890000,
1860 0.2816534000, 0.2826170000, 0.2835796000, 0.2845414000, 0.2855022000, 0.2864623000,
1861 0.2874214000, 0.2883798000, 0.2893373000, 0.2902940000, 0.2912500000, 0.2922052000,
1862 0.2931597000, 0.2941134000, 0.2950664000, 0.2960187000, 0.2969704000, 0.2979213000,
1863 0.2988716000, 0.2998213000, 0.3007704000, 0.3017189000, 0.3026667000, 0.3036140000,
1864 0.3045608000, 0.3055070000, 0.3064527000, 0.3073978000, 0.3083425000, 0.3092867000,
1865 0.3102304000, 0.3111737000, 0.3121166000, 0.3130590000, 0.3140010000, 0.3149426000,
1866 0.3158839000, 0.3168248000, 0.3177653000, 0.3187056000, 0.3196455000, 0.3205851000,
1867 0.3215244000, 0.3224634000, 0.3234022000, 0.3243408000, 0.3252791000, 0.3262172000,
1868 0.3271551000, 0.3280928000, 0.3290303000, 0.3299677000, 0.3309050000, 0.3318421000,
1869 0.3327791000, 0.3337160000, 0.3346528000, 0.3355896000, 0.3365263000, 0.3374629000,
1870 0.3383995000, 0.3393362000, 0.3402728000, 0.3412094000, 0.3421460000, 0.3430827000,
1871 0.3440194000, 0.3449562000, 0.3458931000, 0.3468301000, 0.3477672000, 0.3487044000,
1872 0.3496418000, 0.3505793000, 0.3515170000, 0.3524548000, 0.3533929000, 0.3543311000,
1873 0.3552696000, 0.3562083000, 0.3571472000, 0.3580864000, 0.3590259000, 0.3599657000,
1874 0.3609057000, 0.3618461000, 0.3627868000, 0.3637279000, 0.3646693000, 0.3656111000,
1875 0.3665532000, 0.3674958000, 0.3684387000, 0.3693821000, 0.3703259000, 0.3712701000,
1876 0.3722149000, 0.3731600000, 0.3741057000, 0.3750519000, 0.3759986000, 0.3769458000,
1877 0.3778936000, 0.3788419000, 0.3797908000, 0.3807402000, 0.3816903000, 0.3826409000,
1878 0.3835922000, 0.3845441000, 0.3854966000, 0.3864498000, 0.3874037000, 0.3883583000,
1879 0.3893135000, 0.3902695000, 0.3912262000, 0.3921836000, 0.3931418000, 0.3941007000,
1880 0.3950604000, 0.3960209000, 0.3969822000, 0.3979443000, 0.3989072000, 0.3998710000,
1881 0.4008356000, 0.4018011000, 0.4027675000, 0.4037347000, 0.4047029000, 0.4056720000,
1882 0.4066420000, 0.4076129000, 0.4085848000, 0.4095576000, 0.4105315000, 0.4115063000,
1883 0.4124821000, 0.4134590000, 0.4144368000, 0.4154157000, 0.4163957000, 0.4173767000,
1884 0.4183588000, 0.4193420000, 0.4203264000, 0.4213118000, 0.4222983000, 0.4232860000,
1885 0.4242749000, 0.4252649000, 0.4262561000, 0.4272484000, 0.4282420000, 0.4292368000,
1886 0.4302328000, 0.4312301000, 0.4322286000, 0.4332284000, 0.4342295000, 0.4352318000,
1887 0.4362354000, 0.4372404000, 0.4382467000, 0.4392543000, 0.4402633000, 0.4412736000,
1888 0.4422853000, 0.4432984000, 0.4443129000, 0.4453288000, 0.4463461000, 0.4473648000,
1889 0.4483850000, 0.4494067000, 0.4504298000, 0.4514544000, 0.4524804000, 0.4535080000,
1890 0.4545371000, 0.4555678000, 0.4565999000, 0.4576337000, 0.4586689000, 0.4597058000,
1891 0.4607442000, 0.4617843000, 0.4628259000, 0.4638692000, 0.4649141000, 0.4659606000,
1892 0.4670088000, 0.4680586000, 0.4691101000, 0.4701633000, 0.4712183000, 0.4722749000,
1893 0.4733332000, 0.4743933000, 0.4754551000, 0.4765186000, 0.4775840000, 0.4786511000,
1894 0.4797199000, 0.4807906000, 0.4818631000, 0.4829374000, 0.4840135000, 0.4850915000,
1895 0.4861713000, 0.4872530000, 0.4883366000, 0.4894220000, 0.4905093000, 0.4915985000,
1896 0.4926897000, 0.4937827000, 0.4948777000, 0.4959747000, 0.4970736000, 0.4981744000,
1897 0.4992772000, 0.5003820000, 0.5014888000, 0.5025976000, 0.5037085000, 0.5048213000,
1898 0.5059362000, 0.5070531000, 0.5081720000, 0.5092931000, 0.5104162000, 0.5115413000,
1899 0.5126686000, 0.5137980000, 0.5149294000, 0.5160630000, 0.5171987000, 0.5183366000,
1900 0.5194766000, 0.5206187000, 0.5217630000, 0.5229095000, 0.5240581000, 0.5252090000,
1901 0.5263620000, 0.5275173000, 0.5286748000, 0.5298344000, 0.5309964000, 0.5321605000,
1902 0.5333269000, 0.5344956000, 0.5356665000, 0.5368397000, 0.5380152000, 0.5391930000,
1903 0.5403731000, 0.5415554000, 0.5427401000, 0.5439272000, 0.5451165000, 0.5463082000,
1904 0.5475022000, 0.5486986000, 0.5498973000, 0.5510984000, 0.5523019000, 0.5535077000,
1905 0.5547159000, 0.5559266000, 0.5571396000, 0.5583551000, 0.5595729000, 0.5607932000,
1906 0.5620159000, 0.5632411000, 0.5644686000, 0.5656987000, 0.5669312000, 0.5681661000,
1907 0.5694036000, 0.5706435000, 0.5718859000, 0.5731307000, 0.5743781000, 0.5756280000,
1908 0.5768804000, 0.5781352000, 0.5793926000, 0.5806526000, 0.5819150000, 0.5831800000,
1909 0.5844476000, 0.5857176000, 0.5869903000, 0.5882655000, 0.5895432000, 0.5908235000,
1910 0.5921064000, 0.5933919000, 0.5946799000, 0.5959706000, 0.5972638000, 0.5985596000,
1911 0.5998581000, 0.6011591000, 0.6024627000, 0.6037690000, 0.6050779000, 0.6063894000,
1912 0.6077036000, 0.6090203000, 0.6103397000, 0.6116618000, 0.6129865000, 0.6143139000,
1913 0.6156439000, 0.6169765000, 0.6183119000, 0.6196499000, 0.6209905000, 0.6223339000,
1914 0.6236799000, 0.6250286000, 0.6263800000, 0.6277340000, 0.6290908000, 0.6304503000,
1915 0.6318124000, 0.6331773000, 0.6345448000, 0.6359151000, 0.6372881000, 0.6386638000,
1916 0.6400422000, 0.6414233000, 0.6428071000, 0.6441937000, 0.6455830000, 0.6469750000,
1917 0.6483698000, 0.6497673000, 0.6511675000, 0.6525705000, 0.6539762000, 0.6553846000,
1918 0.6567958000, 0.6582098000, 0.6596264000, 0.6610459000, 0.6624681000, 0.6638930000,
1919 0.6653208000, 0.6667512000, 0.6681845000, 0.6696205000, 0.6710592000, 0.6725007000,
1920 0.6739450000, 0.6753921000, 0.6768419000, 0.6782945000, 0.6797499000, 0.6812080000,
1921 0.6826690000, 0.6841327000, 0.6855992000, 0.6870684000, 0.6885404000, 0.6900153000,
1922 0.6914929000, 0.6929732000, 0.6944564000, 0.6959423000, 0.6974311000, 0.6989226000,
1923 0.7004169000, 0.7019140000, 0.7034139000, 0.7049165000, 0.7064220000, 0.7079302000,
1924 0.7094412000, 0.7109550000, 0.7124716000, 0.7139910000, 0.7155132000, 0.7170382000,
1925 0.7185659000, 0.7200965000, 0.7216298000, 0.7231660000, 0.7247049000, 0.7262466000,
1926 0.7277911000, 0.7293383000, 0.7308884000, 0.7324413000, 0.7339969000, 0.7355553000,
1927 0.7371166000, 0.7386806000, 0.7402474000, 0.7418169000, 0.7433893000, 0.7449644000,
1928 0.7465424000, 0.7481231000, 0.7497066000, 0.7512929000, 0.7528819000, 0.7544737000,
1929 0.7560684000, 0.7576657000, 0.7592659000, 0.7608689000, 0.7624746000, 0.7640831000,
1930 0.7656943000, 0.7673084000, 0.7689252000, 0.7705448000, 0.7721671000, 0.7737922000,
1931 0.7754201000, 0.7770507000, 0.7786841000, 0.7803203000, 0.7819592000, 0.7836009000,
1932 0.7852454000, 0.7868926000, 0.7885425000, 0.7901952000, 0.7918507000, 0.7935089000,
1933 0.7951698000, 0.7968335000, 0.7985000000, 0.8001692000, 0.8018411000, 0.8035158000,
1934 0.8051932000, 0.8068733000, 0.8085562000, 0.8102418000, 0.8119301000, 0.8136212000,
1935 0.8153150000, 0.8170115000, 0.8187108000, 0.8204127000, 0.8221174000, 0.8238248000,
1936 0.8255349000, 0.8272477000, 0.8289633000, 0.8306815000, 0.8324025000, 0.8341261000,
1937 0.8358525000, 0.8375815000, 0.8393133000, 0.8410477000, 0.8427849000, 0.8445247000,
1938 0.8462672000, 0.8480124000, 0.8497603000, 0.8515108000, 0.8532641000, 0.8550200000,
1939 0.8567786000, 0.8585399000, 0.8603038000, 0.8620704000, 0.8638396000, 0.8656116000,
1940 0.8673861000, 0.8691634000, 0.8709433000, 0.8727258000, 0.8745110000, 0.8762988000,
1941 0.8780893000, 0.8798824000, 0.8816782000, 0.8834766000, 0.8852776000, 0.8870812000,
1942 0.8888875000, 0.8906964000, 0.8925079000, 0.8943221000, 0.8961388000, 0.8979582000,
1943 0.8997801000, 0.9016047000, 0.9034319000, 0.9052617000, 0.9070941000, 0.9089291000,
1944 0.9107666000, 0.9126068000, 0.9144495000, 0.9162949000, 0.9181428000, 0.9199933000,
1945 0.9218463000, 0.9237020000, 0.9255602000, 0.9274210000, 0.9292843000, 0.9311502000,
1946 0.9330186000, 0.9348896000, 0.9367632000, 0.9386393000, 0.9405179000, 0.9423991000,
1947 0.9442828000, 0.9461691000, 0.9480579000, 0.9499492000, 0.9518431000, 0.9537394000,
1948 0.9556383000, 0.9575397000, 0.9594436000, 0.9613501000, 0.9632590000, 0.9651704000,
1949 0.9670844000, 0.9690008000, 0.9709197000, 0.9728411000, 0.9747650000, 0.9766914000,
1950 0.9786203000, 0.9805516000, 0.9824854000, 0.9844217000, 0.9863605000, 0.9883017000,
1951 0.9902454000, 0.9921915000, 0.9941401000, 0.9960911000, 0.9980446000, 1.0000010000,
1952 1.0019590000, 1.0039200000, 1.0058830000, 1.0078490000, 1.0098170000, 1.0117870000,
1953 1.0137600000, 1.0157360000, 1.0177130000, 1.0196940000, 1.0216760000, 1.0236610000,
1954 1.0256490000, 1.0276380000, 1.0296310000, 1.0316250000, 1.0336220000, 1.0356210000,
1955 1.0376230000, 1.0396270000, 1.0416340000, 1.0436430000, 1.0456540000, 1.0476670000,
1956 1.0496830000, 1.0517020000, 1.0537220000, 1.0557450000, 1.0577710000, 1.0597980000,
1957 1.0618280000, 1.0638610000, 1.0658950000, 1.0679320000, 1.0699720000, 1.0720130000,
1958 1.0740570000, 1.0761040000, 1.0781520000, 1.0802030000, 1.0822560000, 1.0843120000,
1959 1.0863690000, 1.0884300000, 1.0904920000, 1.0925570000, 1.0946240000, 1.0966930000,
1960 1.0987640000, 1.1008380000, 1.1029140000, 1.1049920000, 1.1070730000, 1.1091560000,
1961 1.1112410000, 1.1133280000, 1.1154180000, 1.1175090000, 1.1196030000, 1.1217000000,
1962 1.1237980000, 1.1258990000, 1.1280020000, 1.1301070000, 1.1322140000, 1.1343240000,
1963 1.1364360000, 1.1385500000, 1.1406660000, 1.1427850000, 1.1449050000, 1.1470280000,
1964 1.1491530000, 1.1512800000, 1.1534100000, 1.1555410000, 1.1576750000, 1.1598110000,
1965 1.1619490000, 1.1640890000, 1.1662310000, 1.1683760000, 1.1705230000, 1.1726710000,
1966 1.1748220000, 1.1769750000, 1.1791310000, 1.1812880000, 1.1834480000, 1.1856090000,
1967 1.1877730000, 1.1899390000, 1.1921070000, 1.1942770000, 1.1964490000, 1.1986240000,
1968 1.2008000000, 1.2029780000, 1.2051590000, 1.2073420000, 1.2095270000, 1.2117130000,
1969 1.2139020000, 1.2160930000, 1.2182860000, 1.2204820000, 1.2226790000, 1.2248780000,
1970 1.2270790000, 1.2292830000, 1.2314880000, 1.2336960000, 1.2359050000, 1.2381170000,
1971 1.2403300000, 1.2425460000, 1.2447640000, 1.2469830000, 1.2492050000, 1.2514290000,
1972 1.2536540000, 1.2558820000, 1.2581120000, 1.2603430000, 1.2625770000, 1.2648130000,
1973 1.2670510000, 1.2692900000, 1.2715320000, 1.2737760000, 1.2760210000, 1.2782690000,
1974 1.2805180000, 1.2827700000, 1.2850230000, 1.2872790000, 1.2895360000, 1.2917950000,
1975 1.2940560000, 1.2963200000, 1.2985850000, 1.3008520000, 1.3031210000, 1.3053920000,
1976 1.3076650000, 1.3099390000, 1.3122160000, 1.3144940000, 1.3167750000, 1.3190570000,
1977 1.3213420000, 1.3236280000, 1.3259160000, 1.3282060000, 1.3304970000, 1.3327910000,
1978 1.3350870000, 1.3373840000, 1.3396830000, 1.3419850000, 1.3442880000, 1.3465920000,
1979 1.3488990000, 1.3512080000, 1.3535180000, 1.3558300000, 1.3581440000, 1.3604600000,
1980 1.3627780000, 1.3650970000, 1.3674190000, 1.3697420000, 1.3720670000, 1.3743940000,
1981 1.3767220000, 1.3790530000, 1.3813850000, 1.3837190000, 1.3860540000, 1.3883920000,
1982 1.3907310000, 1.3930720000, 1.3954150000, 1.3977600000, 1.4001060000, 1.4024540000,
1983 1.4048040000, 1.4071560000, 1.4095090000, 1.4118640000, 1.4142210000, 1.4165800000,
1984 1.4189400000, 1.4213020000, 1.4236660000, 1.4260310000, 1.4283990000, 1.4307670000,
1985 1.4331380000, 1.4355100000, 1.4378840000, 1.4402600000, 1.4426370000, 1.4450160000,
1986 1.4473970000, 1.4497800000, 1.4521640000, 1.4545490000, 1.4569370000, 1.4593260000,
1987 1.4617170000, 1.4641090000, 1.4665030000, 1.4688990000, 1.4712960000, 1.4736950000,
1988 1.4760960000, 1.4784980000, 1.4809020000, 1.4833070000, 1.4857140000, 1.4881230000,
1989 1.4905330000, 1.4929450000, 1.4953580000, 1.4977730000, 1.5001900000, 1.5026080000,
1990 1.5050280000, 1.5074490000, 1.5098720000, 1.5122970000, 1.5147230000, 1.5171500000,
1991 1.5195790000, 1.5220100000, 1.5244420000, 1.5268760000, 1.5293110000, 1.5317480000,
1992 1.5341870000, 1.5366270000, 1.5390680000, 1.5415110000, 1.5439550000, 1.5464010000,
1993 1.5488490000, 1.5512980000, 1.5537480000, 1.5562000000, 1.5586530000, 1.5611080000,
1994 1.5635650000, 1.5660230000, 1.5684820000, 1.5709430000, 1.5734050000, 1.5758680000,
1995 1.5783340000, 1.5808000000, 1.5832680000, 1.5857380000, 1.5882080000, 1.5906810000,
1996 1.5931540000, 1.5956290000, 1.5981060000, 1.6005840000, 1.6030630000, 1.6055440000,
1997 1.6080260000, 1.6105100000, 1.6129950000, 1.6154810000, 1.6179680000, 1.6204570000,
1998 1.6229480000, 1.6254400000, 1.6279330000, 1.6304270000, 1.6329230000, 1.6354200000,
1999 1.6379190000, 1.6404190000, 1.6429200000, 1.6454220000, 1.6479260000, 1.6504310000,
2000 1.6529380000, 1.6554450000, 1.6579540000, 1.6604650000, 1.6629760000, 1.6654890000,
2001 1.6680040000, 1.6705190000, 1.6730360000, 1.6755540000, 1.6780730000, 1.6805940000,
2002 1.6831160000, 1.6856390000, 1.6881630000, 1.6906890000, 1.6932160000, 1.6957440000,
2003 1.6982730000, 1.7008030000, 1.7033350000, 1.7058680000, 1.7084020000, 1.7109380000,
2004 1.7134740000, 1.7160120000, 1.7185510000, 1.7210910000, 1.7236330000, 1.7261750000,
2005 1.7287190000, 1.7312640000, 1.7338100000, 1.7363570000, 1.7389050000, 1.7414550000,
2006 1.7440050000, 1.7465570000, 1.7491100000, 1.7516640000, 1.7542190000, 1.7567760000,
2007 1.7593330000, 1.7618920000, 1.7644520000, 1.7670120000, 1.7695740000, 1.7721370000,
2008 1.7747010000, 1.7772660000, 1.7798330000, 1.7824000000, 1.7849680000, 1.7875380000,
2009 1.7901080000, 1.7926800000, 1.7952530000, 1.7978260000, 1.8004010000, 1.8029770000,
2010 1.8055540000, 1.8081320000, 1.8107100000, 1.8132900000, 1.8158710000, 1.8184530000,
2011 1.8210360000, 1.8236200000, 1.8262050000, 1.8287910000, 1.8313780000, 1.8339660000,
2012 1.8365550000, 1.8391450000, 1.8417360000, 1.8443270000, 1.8469200000, 1.8495140000,
2013 1.8521090000, 1.8547040000, 1.8573010000, 1.8598980000, 1.8624970000, 1.8650960000,
2014 1.8676970000, 1.8702980000, 1.8729000000, 1.8755030000, 1.8781070000, 1.8807120000,
2015 1.8833180000, 1.8859250000, 1.8885320000, 1.8911410000, 1.8937500000, 1.8963600000,
2016 1.8989710000, 1.9015830000, 1.9041960000, 1.9068100000, 1.9094240000, 1.9120400000,
2017 1.9146560000, 1.9172730000, 1.9198910000, 1.9225100000, 1.9251290000, 1.9277500000,
2018 1.9303710000, 1.9329930000, 1.9356160000, 1.9382400000, 1.9408640000, 1.9434890000,
2019 1.9461150000, 1.9487420000, 1.9513700000, 1.9539980000, 1.9566270000, 1.9592570000,
2020 1.9618880000, 1.9645190000, 1.9671520000, 1.9697850000, 1.9724180000, 1.9750530000,
2021 1.9776880000, 1.9803240000, 1.9829610000, 1.9855980000, 1.9882360000, 1.9908750000,
2022 1.9935140000, 1.9961550000, 1.9987960000, 2.0014370000, 2.0040800000, 2.0067230000,
2023 2.0093660000, 2.0120110000, 2.0146560000, 2.0173010000, 2.0199480000, 2.0225950000,
2024 2.0252430000, 2.0278910000, 2.0305400000, 2.0331900000, 2.0358400000, 2.0384910000,
2025 2.0411430000, 2.0437950000, 2.0464480000, 2.0491010000, 2.0517550000, 2.0544100000,
2026 2.0570650000, 2.0597210000, 2.0623780000, 2.0650350000, 2.0676930000, 2.0703510000,
2027 2.0730100000, 2.0756690000, 2.0783290000, 2.0809900000, 2.0836510000, 2.0863130000,
2028 2.0889750000, 2.0916380000, 2.0943010000, 2.0969650000, 2.0996290000, 2.1022940000,
2029 2.1049600000, 2.1076260000, 2.1102930000, 2.1129600000, 2.1156270000, 2.1182950000,
2030 2.1209640000, 2.1236330000, 2.1263030000, 2.1289730000, 2.1316440000, 2.1343150000,
2031 2.1369860000, 2.1396580000, 2.1423310000, 2.1450040000, 2.1476780000, 2.1503510000,
2032 2.1530260000, 2.1557010000 };
2033 double K1400_width[2000] = {
2034 // 1700 bin, s = (0.6, 4)
2035 0.0000000007, 0.0000000025, 0.0000000061, 0.0000000121, 0.0000000213, 0.0000000344,
2036 0.0000000519, 0.0000000747, 0.0000001034, 0.0000001389, 0.0000001818, 0.0000002329,
2037 0.0000002930, 0.0000003630, 0.0000004435, 0.0000005356, 0.0000006400, 0.0000007576,
2038 0.0000008892, 0.0000010359, 0.0000011985, 0.0000013780, 0.0000015753, 0.0000017914,
2039 0.0000020274, 0.0000022842, 0.0000025629, 0.0000028646, 0.0000031903, 0.0000035413,
2040 0.0000039185, 0.0000043232, 0.0000047566, 0.0000052199, 0.0000057143, 0.0000062411,
2041 0.0000068015, 0.0000073970, 0.0000080288, 0.0000086983, 0.0000094069, 0.0000101561,
2042 0.0000109474, 0.0000117821, 0.0000126619, 0.0000135884, 0.0000145631, 0.0000155876,
2043 0.0000166637, 0.0000177930, 0.0000189773, 0.0000202184, 0.0000215182, 0.0000228784,
2044 0.0000243010, 0.0000257879, 0.0000273412, 0.0000289629, 0.0000306551, 0.0000324199,
2045 0.0000342596, 0.0000361763, 0.0000381724, 0.0000402502, 0.0000424121, 0.0000446606,
2046 0.0000469982, 0.0000494275, 0.0000519511, 0.0000545718, 0.0000572923, 0.0000601155,
2047 0.0000630443, 0.0000660816, 0.0000692306, 0.0000724943, 0.0000758761, 0.0000793791,
2048 0.0000830068, 0.0000867626, 0.0000906501, 0.0000946730, 0.0000988349, 0.0001031398,
2049 0.0001075916, 0.0001121943, 0.0001169520, 0.0001218691, 0.0001269499, 0.0001321989,
2050 0.0001376208, 0.0001432202, 0.0001490021, 0.0001549715, 0.0001611335, 0.0001674933,
2051 0.0001740566, 0.0001808287, 0.0001878156, 0.0001950231, 0.0002024573, 0.0002101244,
2052 0.0002180310, 0.0002261837, 0.0002345892, 0.0002432546, 0.0002521872, 0.0002613945,
2053 0.0002708840, 0.0002806638, 0.0002907420, 0.0003011271, 0.0003118276, 0.0003228526,
2054 0.0003342113, 0.0003459132, 0.0003579682, 0.0003703863, 0.0003831781, 0.0003963544,
2055 0.0004099263, 0.0004239055, 0.0004383037, 0.0004531333, 0.0004684072, 0.0004841383,
2056 0.0005003405, 0.0005170277, 0.0005342145, 0.0005519161, 0.0005701480, 0.0005889265,
2057 0.0006082683, 0.0006281907, 0.0006487119, 0.0006698504, 0.0006916256, 0.0007140575,
2058 0.0007371670, 0.0007609758, 0.0007855062, 0.0008107815, 0.0008368260, 0.0008636648,
2059 0.0008913240, 0.0009198309, 0.0009492137, 0.0009795018, 0.0010107260, 0.0010429180,
2060 0.0010761110, 0.0011103390, 0.0011456390, 0.0011820480, 0.0012196060, 0.0012583530,
2061 0.0012983320, 0.0013395870, 0.0013821650, 0.0014261150, 0.0014714870, 0.0015183340,
2062 0.0015667120, 0.0016166800, 0.0016682970, 0.0017216280, 0.0017767390, 0.0018337000,
2063 0.0018925850, 0.0019534700, 0.0020164360, 0.0020815670, 0.0021489530, 0.0022186860,
2064 0.0022908640, 0.0023655910, 0.0024429750, 0.0025231290, 0.0026061740, 0.0026922350,
2065 0.0027814460, 0.0028739440, 0.0029698790, 0.0030694050, 0.0031726840, 0.0032798900,
2066 0.0033912030, 0.0035068140, 0.0036269240, 0.0037517470, 0.0038815040, 0.0040164320,
2067 0.0041567790, 0.0043028050, 0.0044547850, 0.0046130090, 0.0047777790, 0.0049494150,
2068 0.0051282500, 0.0053146360, 0.0055089370, 0.0057115350, 0.0059228300, 0.0061432340,
2069 0.0063731760, 0.0066130990, 0.0068634590, 0.0071247200, 0.0073973600, 0.0076818580,
2070 0.0079786990, 0.0082883670, 0.0086113370, 0.0089480750, 0.0092990290, 0.0096646210,
2071 0.0100452400, 0.0104412500, 0.0108529300, 0.0112805500, 0.0117242700, 0.0121842100,
2072 0.0126603900, 0.0131527500, 0.0136611300, 0.0141852900, 0.0147248800, 0.0152794700,
2073 0.0158485100, 0.0164314100, 0.0170274700, 0.0176359300, 0.0182559700, 0.0188867400,
2074 0.0195273500, 0.0201768900, 0.0208344400, 0.0214991000, 0.0221699600, 0.0228461600,
2075 0.0235268400, 0.0242112100, 0.0248985100, 0.0255880100, 0.0262790500, 0.0269710000,
2076 0.0276632900, 0.0283554000, 0.0290468400, 0.0297371800, 0.0304260300, 0.0311130400,
2077 0.0317978800, 0.0324802900, 0.0331600000, 0.0338368100, 0.0345105300, 0.0351809800,
2078 0.0358480300, 0.0365115500, 0.0371714500, 0.0378276300, 0.0384800200, 0.0391285600,
2079 0.0397732200, 0.0404139500, 0.0410507400, 0.0416835500, 0.0423124000, 0.0429372700,
2080 0.0435581700, 0.0441751100, 0.0447881100, 0.0453971800, 0.0460023400, 0.0466036300,
2081 0.0472010600, 0.0477946700, 0.0483845000, 0.0489705600, 0.0495529100, 0.0501315800,
2082 0.0507066000, 0.0512780200, 0.0518458700, 0.0524101800, 0.0529710200, 0.0535284000,
2083 0.0540823700, 0.0546329800, 0.0551802500, 0.0557242400, 0.0562649800, 0.0568025200,
2084 0.0573368800, 0.0578681100, 0.0583962500, 0.0589213300, 0.0594434000, 0.0599624900,
2085 0.0604786400, 0.0609918800, 0.0615022500, 0.0620097800, 0.0625145200, 0.0630164900,
2086 0.0635157300, 0.0640122700, 0.0645061500, 0.0649973900, 0.0654860300, 0.0659721100,
2087 0.0664556400, 0.0669366600, 0.0674152100, 0.0678913000, 0.0683649800, 0.0688362600,
2088 0.0693051800, 0.0697717600, 0.0702360300, 0.0706980200, 0.0711577500, 0.0716152500,
2089 0.0720705400, 0.0725236500, 0.0729746100, 0.0734234300, 0.0738701400, 0.0743147600,
2090 0.0747573300, 0.0751978500, 0.0756363500, 0.0760728500, 0.0765073800, 0.0769399600,
2091 0.0773706000, 0.0777993300, 0.0782261600, 0.0786511200, 0.0790742300, 0.0794955000,
2092 0.0799149600, 0.0803326100, 0.0807484900, 0.0811626100, 0.0815749800, 0.0819856300,
2093 0.0823945700, 0.0828018200, 0.0832073900, 0.0836113000, 0.0840135700, 0.0844142200,
2094 0.0848132500, 0.0852106900, 0.0856065500, 0.0860008400, 0.0863935900, 0.0867848000,
2095 0.0871744900, 0.0875626700, 0.0879493600, 0.0883345700, 0.0887183200, 0.0891006200,
2096 0.0894814800, 0.0898609100, 0.0902389300, 0.0906155600, 0.0909908000, 0.0913646700,
2097 0.0917371700, 0.0921083300, 0.0924781500, 0.0928466500, 0.0932138400, 0.0935797200,
2098 0.0939443200, 0.0943076400, 0.0946696900, 0.0950304800, 0.0953900400, 0.0957483600,
2099 0.0961054500, 0.0964613400, 0.0968160300, 0.0971695200, 0.0975218300, 0.0978729800,
2100 0.0982229600, 0.0985718000, 0.0989194900, 0.0992660600, 0.0996115100, 0.0999558400,
2101 0.1002991000, 0.1006412000, 0.1009823000, 0.1013223000, 0.1016612000, 0.1019991000,
2102 0.1023359000, 0.1026717000, 0.1030064000, 0.1033402000, 0.1036729000, 0.1040047000,
2103 0.1043354000, 0.1046652000, 0.1049939000, 0.1053218000, 0.1056486000, 0.1059745000,
2104 0.1062995000, 0.1066235000, 0.1069466000, 0.1072688000, 0.1075901000, 0.1079104000,
2105 0.1082299000, 0.1085485000, 0.1088662000, 0.1091831000, 0.1094990000, 0.1098142000,
2106 0.1101284000, 0.1104418000, 0.1107544000, 0.1110662000, 0.1113771000, 0.1116873000,
2107 0.1119966000, 0.1123051000, 0.1126128000, 0.1129198000, 0.1132259000, 0.1135313000,
2108 0.1138360000, 0.1141398000, 0.1144429000, 0.1147453000, 0.1150469000, 0.1153478000,
2109 0.1156480000, 0.1159474000, 0.1162461000, 0.1165442000, 0.1168415000, 0.1171381000,
2110 0.1174340000, 0.1177292000, 0.1180238000, 0.1183177000, 0.1186109000, 0.1189034000,
2111 0.1191953000, 0.1194865000, 0.1197771000, 0.1200670000, 0.1203563000, 0.1206449000,
2112 0.1209329000, 0.1212203000, 0.1215071000, 0.1217933000, 0.1220788000, 0.1223637000,
2113 0.1226480000, 0.1229318000, 0.1232149000, 0.1234974000, 0.1237794000, 0.1240607000,
2114 0.1243415000, 0.1246217000, 0.1249013000, 0.1251803000, 0.1254588000, 0.1257367000,
2115 0.1260140000, 0.1262908000, 0.1265670000, 0.1268427000, 0.1271178000, 0.1273923000,
2116 0.1276663000, 0.1279398000, 0.1282127000, 0.1284850000, 0.1287568000, 0.1290281000,
2117 0.1292988000, 0.1295690000, 0.1298387000, 0.1301078000, 0.1303764000, 0.1306445000,
2118 0.1309120000, 0.1311790000, 0.1314455000, 0.1317114000, 0.1319768000, 0.1322417000,
2119 0.1325060000, 0.1327699000, 0.1330332000, 0.1332960000, 0.1335582000, 0.1338199000,
2120 0.1340812000, 0.1343418000, 0.1346020000, 0.1348617000, 0.1351208000, 0.1355387000,
2121 0.1358629000, 0.1361711000, 0.1364708000, 0.1367648000, 0.1370546000, 0.1373411000,
2122 0.1376249000, 0.1379063000, 0.1381858000, 0.1384634000, 0.1387394000, 0.1390139000,
2123 0.1392870000, 0.1395588000, 0.1398294000, 0.1400988000, 0.1403671000, 0.1406344000,
2124 0.1409007000, 0.1411659000, 0.1414302000, 0.1416936000, 0.1419562000, 0.1422178000,
2125 0.1424786000, 0.1427386000, 0.1429977000, 0.1432561000, 0.1435136000, 0.1437704000,
2126 0.1440265000, 0.1442818000, 0.1445364000, 0.1447902000, 0.1450434000, 0.1452958000,
2127 0.1455475000, 0.1457986000, 0.1460490000, 0.1462986000, 0.1465477000, 0.1467960000,
2128 0.1470437000, 0.1472908000, 0.1475372000, 0.1477830000, 0.1480281000, 0.1482726000,
2129 0.1485165000, 0.1487597000, 0.1490024000, 0.1492444000, 0.1494858000, 0.1497266000,
2130 0.1499668000, 0.1502065000, 0.1504455000, 0.1506839000, 0.1509218000, 0.1511590000,
2131 0.1513957000, 0.1516318000, 0.1518674000, 0.1521023000, 0.1523367000, 0.1525706000,
2132 0.1528038000, 0.1530366000, 0.1532687000, 0.1535003000, 0.1537314000, 0.1539619000,
2133 0.1541918000, 0.1544213000, 0.1546501000, 0.1548785000, 0.1551063000, 0.1553336000,
2134 0.1555603000, 0.1557866000, 0.1560122000, 0.1562374000, 0.1564621000, 0.1566862000,
2135 0.1569098000, 0.1571329000, 0.1573555000, 0.1575776000, 0.1577992000, 0.1580203000,
2136 0.1582409000, 0.1584609000, 0.1586805000, 0.1588996000, 0.1591182000, 0.1593363000,
2137 0.1595539000, 0.1597710000, 0.1599876000, 0.1602038000, 0.1604194000, 0.1606346000,
2138 0.1608493000, 0.1610636000, 0.1612773000, 0.1614906000, 0.1617035000, 0.1619158000,
2139 0.1621277000, 0.1623391000, 0.1625501000, 0.1627606000, 0.1629706000, 0.1631802000,
2140 0.1633894000, 0.1635981000, 0.1638063000, 0.1640141000, 0.1642214000, 0.1644283000,
2141 0.1646347000, 0.1648407000, 0.1650463000, 0.1652514000, 0.1654561000, 0.1656604000,
2142 0.1658642000, 0.1660676000, 0.1662706000, 0.1664731000, 0.1666752000, 0.1668769000,
2143 0.1670782000, 0.1672790000, 0.1674794000, 0.1676794000, 0.1678790000, 0.1680782000,
2144 0.1682769000, 0.1684753000, 0.1686732000, 0.1688707000, 0.1690679000, 0.1692646000,
2145 0.1694609000, 0.1696568000, 0.1698523000, 0.1700474000, 0.1702421000, 0.1704364000,
2146 0.1706304000, 0.1708239000, 0.1710170000, 0.1712098000, 0.1714022000, 0.1715941000,
2147 0.1717857000, 0.1719769000, 0.1721677000, 0.1723582000, 0.1725482000, 0.1727379000,
2148 0.1729272000, 0.1731162000, 0.1733047000, 0.1734929000, 0.1736807000, 0.1738682000,
2149 0.1740552000, 0.1742420000, 0.1744283000, 0.1746143000, 0.1747999000, 0.1749851000,
2150 0.1751700000, 0.1753546000, 0.1755388000, 0.1757226000, 0.1759060000, 0.1760892000,
2151 0.1762719000, 0.1764543000, 0.1766364000, 0.1768181000, 0.1769995000, 0.1771805000,
2152 0.1773612000, 0.1775415000, 0.1777215000, 0.1779011000, 0.1780804000, 0.1782594000,
2153 0.1784380000, 0.1786163000, 0.1787943000, 0.1789719000, 0.1791492000, 0.1793262000,
2154 0.1795028000, 0.1796791000, 0.1798551000, 0.1800308000, 0.1802061000, 0.1803811000,
2155 0.1805558000, 0.1807301000, 0.1809042000, 0.1810779000, 0.1812513000, 0.1814244000,
2156 0.1815971000, 0.1817696000, 0.1819417000, 0.1821136000, 0.1822851000, 0.1824563000,
2157 0.1826272000, 0.1827978000, 0.1829681000, 0.1831380000, 0.1833077000, 0.1834771000,
2158 0.1836461000, 0.1838149000, 0.1839834000, 0.1841515000, 0.1843194000, 0.1844869000,
2159 0.1846542000, 0.1848212000, 0.1849879000, 0.1851543000, 0.1853203000, 0.1854862000,
2160 0.1856517000, 0.1858169000, 0.1859818000, 0.1861465000, 0.1863108000, 0.1864749000,
2161 0.1866387000, 0.1868022000, 0.1869654000, 0.1871284000, 0.1872910000, 0.1874534000,
2162 0.1876155000, 0.1877774000, 0.1879389000, 0.1881002000, 0.1882612000, 0.1884219000,
2163 0.1885824000, 0.1887426000, 0.1889025000, 0.1890621000, 0.1892215000, 0.1893806000,
2164 0.1895395000, 0.1896980000, 0.1898563000, 0.1900144000, 0.1901722000, 0.1903297000,
2165 0.1904869000, 0.1906439000, 0.1908007000, 0.1909571000, 0.1911133000, 0.1912693000,
2166 0.1914250000, 0.1915805000, 0.1917357000, 0.1918906000, 0.1920453000, 0.1921997000,
2167 0.1923539000, 0.1925078000, 0.1926615000, 0.1928149000, 0.1929681000, 0.1931211000,
2168 0.1932738000, 0.1934262000, 0.1935784000, 0.1937304000, 0.1938821000, 0.1940336000,
2169 0.1941848000, 0.1943358000, 0.1944865000, 0.1946370000, 0.1947873000, 0.1949374000,
2170 0.1950872000, 0.1952367000, 0.1953860000, 0.1955351000, 0.1956840000, 0.1958326000,
2171 0.1959810000, 0.1961292000, 0.1962771000, 0.1964248000, 0.1965723000, 0.1967196000,
2172 0.1968666000, 0.1970134000, 0.1971600000, 0.1973063000, 0.1974524000, 0.1975983000,
2173 0.1977440000, 0.1978894000, 0.1980347000, 0.1981797000, 0.1983245000, 0.1984690000,
2174 0.1986134000, 0.1987575000, 0.1989015000, 0.1990452000, 0.1991886000, 0.1993319000,
2175 0.1994750000, 0.1996178000, 0.1997604000, 0.1999029000, 0.2000451000, 0.2001871000,
2176 0.2003289000, 0.2004704000, 0.2006118000, 0.2007530000, 0.2008939000, 0.2010347000,
2177 0.2011752000, 0.2013156000, 0.2014557000, 0.2015956000, 0.2017353000, 0.2018749000,
2178 0.2020142000, 0.2021533000, 0.2022922000, 0.2024310000, 0.2025695000, 0.2027078000,
2179 0.2028459000, 0.2029839000, 0.2031216000, 0.2032591000, 0.2033965000, 0.2035336000,
2180 0.2036706000, 0.2038073000, 0.2039439000, 0.2040803000, 0.2042165000, 0.2043524000,
2181 0.2044883000, 0.2046239000, 0.2047593000, 0.2048945000, 0.2050296000, 0.2051644000,
2182 0.2052991000, 0.2054336000, 0.2055679000, 0.2057020000, 0.2058360000, 0.2059697000,
2183 0.2061033000, 0.2062367000, 0.2063699000, 0.2065029000, 0.2066358000, 0.2067684000,
2184 0.2069009000, 0.2070332000, 0.2071653000, 0.2072973000, 0.2074291000, 0.2075607000,
2185 0.2076921000, 0.2078233000, 0.2079544000, 0.2080853000, 0.2082161000, 0.2083466000,
2186 0.2084770000, 0.2086072000, 0.2087372000, 0.2088671000, 0.2089968000, 0.2091263000,
2187 0.2092557000, 0.2093849000, 0.2095139000, 0.2096428000, 0.2097715000, 0.2099000000,
2188 0.2100284000, 0.2101566000, 0.2102846000, 0.2104125000, 0.2105402000, 0.2106677000,
2189 0.2107951000, 0.2109223000, 0.2110494000, 0.2111763000, 0.2113030000, 0.2114296000,
2190 0.2115560000, 0.2116823000, 0.2118084000, 0.2119343000, 0.2120601000, 0.2121858000,
2191 0.2123112000, 0.2124366000, 0.2125617000, 0.2126868000, 0.2128116000, 0.2129363000,
2192 0.2130609000, 0.2131853000, 0.2133095000, 0.2134337000, 0.2135576000, 0.2136814000,
2193 0.2138051000, 0.2139286000, 0.2140519000, 0.2141751000, 0.2142982000, 0.2144211000,
2194 0.2145439000, 0.2146665000, 0.2147890000, 0.2149113000, 0.2150335000, 0.2151555000,
2195 0.2152774000, 0.2153992000, 0.2155208000, 0.2156423000, 0.2157636000, 0.2158848000,
2196 0.2160059000, 0.2161268000, 0.2162475000, 0.2163682000, 0.2164887000, 0.2166090000,
2197 0.2167292000, 0.2168493000, 0.2169693000, 0.2170891000, 0.2172088000, 0.2173283000,
2198 0.2174477000, 0.2175670000, 0.2176861000, 0.2178051000, 0.2179240000, 0.2180427000,
2199 0.2181613000, 0.2182798000, 0.2183982000, 0.2185164000, 0.2186345000, 0.2187524000,
2200 0.2188703000, 0.2189880000, 0.2191055000, 0.2192230000, 0.2193403000, 0.2194575000,
2201 0.2195745000, 0.2196915000, 0.2198083000, 0.2199250000, 0.2200416000, 0.2201580000,
2202 0.2202743000, 0.2203905000, 0.2205066000, 0.2206225000, 0.2207384000, 0.2208541000,
2203 0.2209697000, 0.2210851000, 0.2212005000, 0.2213157000, 0.2214308000, 0.2215458000,
2204 0.2216607000, 0.2217754000, 0.2218901000, 0.2220046000, 0.2221190000, 0.2222333000,
2205 0.2223475000, 0.2224615000, 0.2225755000, 0.2226893000, 0.2228030000, 0.2229166000,
2206 0.2230301000, 0.2231435000, 0.2232568000, 0.2233699000, 0.2234830000, 0.2235959000,
2207 0.2237088000, 0.2238215000, 0.2239341000, 0.2240466000, 0.2241590000, 0.2242713000,
2208 0.2243834000, 0.2244955000, 0.2246075000, 0.2247193000, 0.2248311000, 0.2249427000,
2209 0.2250543000, 0.2251657000, 0.2252770000, 0.2253883000, 0.2254994000, 0.2256104000,
2210 0.2257213000, 0.2258321000, 0.2259429000, 0.2260535000, 0.2261640000, 0.2262744000,
2211 0.2263847000, 0.2264949000, 0.2266050000, 0.2267151000, 0.2268250000, 0.2269348000,
2212 0.2270445000, 0.2271542000, 0.2272637000, 0.2273731000, 0.2274825000, 0.2275917000,
2213 0.2277009000, 0.2278099000, 0.2279189000, 0.2280277000, 0.2281365000, 0.2282452000,
2214 0.2283538000, 0.2284623000, 0.2285707000, 0.2286790000, 0.2287872000, 0.2288954000,
2215 0.2290034000, 0.2291114000, 0.2292192000, 0.2293270000, 0.2294347000, 0.2295423000,
2216 0.2296498000, 0.2297572000, 0.2298646000, 0.2299718000, 0.2300790000, 0.2301861000,
2217 0.2302931000, 0.2304000000, 0.2305068000, 0.2306136000, 0.2307202000, 0.2308268000,
2218 0.2309333000, 0.2310397000, 0.2311460000, 0.2312523000, 0.2313584000, 0.2314645000,
2219 0.2315705000, 0.2316765000, 0.2317823000, 0.2318881000, 0.2319938000, 0.2320994000,
2220 0.2322049000, 0.2323104000, 0.2324157000, 0.2325210000, 0.2326263000, 0.2327314000,
2221 0.2328365000, 0.2329415000, 0.2330464000, 0.2331513000, 0.2332560000, 0.2333607000,
2222 0.2334654000, 0.2335699000, 0.2336744000, 0.2337788000, 0.2338831000, 0.2339874000,
2223 0.2340916000, 0.2341957000, 0.2342998000, 0.2344038000, 0.2345077000, 0.2346115000,
2224 0.2347153000, 0.2348190000, 0.2349227000, 0.2350262000, 0.2351298000, 0.2352332000,
2225 0.2353366000, 0.2354399000, 0.2355431000, 0.2356463000, 0.2357494000, 0.2358525000,
2226 0.2359555000, 0.2360584000, 0.2361612000, 0.2362640000, 0.2363668000, 0.2364694000,
2227 0.2365721000, 0.2366746000, 0.2367771000, 0.2368795000, 0.2369819000, 0.2370842000,
2228 0.2371864000, 0.2372886000, 0.2373908000, 0.2374928000, 0.2375948000, 0.2376968000,
2229 0.2377987000, 0.2379005000, 0.2380023000, 0.2381041000, 0.2382057000, 0.2383074000,
2230 0.2384089000, 0.2385104000, 0.2386119000, 0.2387133000, 0.2388146000, 0.2389159000,
2231 0.2390172000, 0.2391184000, 0.2392195000, 0.2393206000, 0.2394216000, 0.2395226000,
2232 0.2396236000, 0.2397245000, 0.2398253000, 0.2399261000, 0.2400268000, 0.2401275000,
2233 0.2402281000, 0.2403287000, 0.2404293000, 0.2405298000, 0.2406302000, 0.2407306000,
2234 0.2408310000, 0.2409313000, 0.2410316000, 0.2411318000, 0.2412320000, 0.2413321000,
2235 0.2414322000, 0.2415323000, 0.2416323000, 0.2417323000, 0.2418322000, 0.2419321000,
2236 0.2420319000, 0.2421317000, 0.2422315000, 0.2423312000, 0.2424309000, 0.2425305000,
2237 0.2426301000, 0.2427297000, 0.2428292000, 0.2429287000, 0.2430281000, 0.2431275000,
2238 0.2432269000, 0.2433262000, 0.2434255000, 0.2435248000, 0.2436240000, 0.2437232000,
2239 0.2438224000, 0.2439215000, 0.2440206000, 0.2441196000, 0.2442187000, 0.2443176000,
2240 0.2444166000, 0.2445155000, 0.2446144000, 0.2447133000, 0.2448121000, 0.2449109000,
2241 0.2450097000, 0.2451084000, 0.2452071000, 0.2453058000, 0.2454044000, 0.2455031000,
2242 0.2456017000, 0.2457002000, 0.2457988000, 0.2458973000, 0.2459957000, 0.2460942000,
2243 0.2461926000, 0.2462910000, 0.2463894000, 0.2464878000, 0.2465861000, 0.2466844000,
2244 0.2467827000, 0.2468809000, 0.2469792000, 0.2470774000, 0.2471756000, 0.2472737000,
2245 0.2473719000, 0.2474700000, 0.2475681000, 0.2476662000, 0.2477642000, 0.2478623000,
2246 0.2479603000, 0.2480583000, 0.2481563000, 0.2482542000, 0.2483522000, 0.2484501000,
2247 0.2485480000, 0.2486459000, 0.2487437000, 0.2488416000, 0.2489394000, 0.2490372000,
2248 0.2491351000, 0.2492328000, 0.2493306000, 0.2494284000, 0.2495261000, 0.2496239000,
2249 0.2497216000, 0.2498193000, 0.2499170000, 0.2500146000, 0.2501123000, 0.2502100000,
2250 0.2503076000, 0.2504052000, 0.2505029000, 0.2506005000, 0.2506981000, 0.2507956000,
2251 0.2508932000, 0.2509908000, 0.2510883000, 0.2511859000, 0.2512834000, 0.2513810000,
2252 0.2514785000, 0.2515760000, 0.2516735000, 0.2517710000, 0.2518685000, 0.2519660000,
2253 0.2520635000, 0.2521610000, 0.2522584000, 0.2523559000, 0.2524534000, 0.2525508000,
2254 0.2526483000, 0.2527457000, 0.2528432000, 0.2529406000, 0.2530381000, 0.2531355000,
2255 0.2532330000, 0.2533304000, 0.2534278000, 0.2535252000, 0.2536227000, 0.2537201000,
2256 0.2538175000, 0.2539150000, 0.2540124000, 0.2541098000, 0.2542073000, 0.2543047000,
2257 0.2544021000, 0.2544996000, 0.2545970000, 0.2546944000, 0.2547919000, 0.2548893000,
2258 0.2549868000, 0.2550842000, 0.2551817000, 0.2552791000, 0.2553766000, 0.2554740000,
2259 0.2555715000, 0.2556690000, 0.2557665000, 0.2558639000, 0.2559614000, 0.2560589000,
2260 0.2561564000, 0.2562539000, 0.2563515000, 0.2564490000, 0.2565465000, 0.2566441000,
2261 0.2567416000, 0.2568392000, 0.2569367000, 0.2570343000, 0.2571319000, 0.2572295000,
2262 0.2573271000, 0.2574247000, 0.2575223000, 0.2576199000, 0.2577176000, 0.2578152000,
2263 0.2579129000, 0.2580106000, 0.2581082000, 0.2582059000, 0.2583036000, 0.2584014000,
2264 0.2584991000, 0.2585968000, 0.2586946000, 0.2587924000, 0.2588902000, 0.2589879000,
2265 0.2590858000, 0.2591836000, 0.2592814000, 0.2593793000, 0.2594771000, 0.2595750000,
2266 0.2596729000, 0.2597708000, 0.2598688000, 0.2599667000, 0.2600647000, 0.2601627000,
2267 0.2602606000, 0.2603587000, 0.2604567000, 0.2605547000, 0.2606528000, 0.2607509000,
2268 0.2608490000, 0.2609471000, 0.2610452000, 0.2611434000, 0.2612415000, 0.2613397000,
2269 0.2614379000, 0.2615361000, 0.2616344000, 0.2617326000, 0.2618309000, 0.2619292000,
2270 0.2620275000, 0.2621259000, 0.2622242000, 0.2623226000, 0.2624210000, 0.2625194000,
2271 0.2626179000, 0.2627164000, 0.2628148000, 0.2629133000, 0.2630119000, 0.2631104000,
2272 0.2632090000, 0.2633076000, 0.2634062000, 0.2635048000, 0.2636035000, 0.2637022000,
2273 0.2638009000, 0.2638996000, 0.2639984000, 0.2640971000, 0.2641959000, 0.2642947000,
2274 0.2643936000, 0.2644925000, 0.2645913000, 0.2646903000, 0.2647892000, 0.2648882000,
2275 0.2649871000, 0.2650861000, 0.2651852000, 0.2652842000, 0.2653833000, 0.2654824000,
2276 0.2655816000, 0.2656807000, 0.2657799000, 0.2658791000, 0.2659784000, 0.2660776000,
2277 0.2661769000, 0.2662762000, 0.2663755000, 0.2664749000, 0.2665743000, 0.2666737000,
2278 0.2667731000, 0.2668726000, 0.2669721000, 0.2670716000, 0.2671712000, 0.2672707000,
2279 0.2673703000, 0.2674700000, 0.2675696000, 0.2676693000, 0.2677690000, 0.2678687000,
2280 0.2679685000, 0.2680683000, 0.2681681000, 0.2682679000, 0.2683678000, 0.2684677000,
2281 0.2685676000, 0.2686675000, 0.2687675000, 0.2688675000, 0.2689676000, 0.2690676000,
2282 0.2691677000, 0.2692678000, 0.2693680000, 0.2694681000, 0.2695683000, 0.2696685000,
2283 0.2697688000, 0.2698691000, 0.2699694000, 0.2700697000, 0.2701701000, 0.2702705000,
2284 0.2703709000, 0.2704713000, 0.2705718000, 0.2706723000, 0.2707729000, 0.2708734000,
2285 0.2709740000, 0.2710746000, 0.2711753000, 0.2712759000, 0.2713766000, 0.2714774000,
2286 0.2715781000, 0.2716789000, 0.2717797000, 0.2718806000, 0.2719814000, 0.2720823000,
2287 0.2721833000, 0.2722842000, 0.2723852000, 0.2724862000, 0.2725873000, 0.2726883000,
2288 0.2727894000, 0.2728905000, 0.2729917000, 0.2730929000, 0.2731941000, 0.2732953000,
2289 0.2733966000, 0.2734979000, 0.2735992000, 0.2737005000, 0.2738019000, 0.2739033000,
2290 0.2740048000, 0.2741062000, 0.2742077000, 0.2743092000, 0.2744108000, 0.2745123000,
2291 0.2746139000, 0.2747156000, 0.2748172000, 0.2749189000, 0.2750206000, 0.2751223000,
2292 0.2752241000, 0.2753259000, 0.2754277000, 0.2755295000, 0.2756314000, 0.2757333000,
2293 0.2758352000, 0.2759372000, 0.2760392000, 0.2761412000, 0.2762432000, 0.2763453000,
2294 0.2764474000, 0.2765495000, 0.2766516000, 0.2767538000, 0.2768560000, 0.2769582000,
2295 0.2770604000, 0.2771627000, 0.2772650000, 0.2773673000, 0.2774697000, 0.2775720000,
2296 0.2776745000, 0.2777769000, 0.2778793000, 0.2779818000, 0.2780843000, 0.2781868000,
2297 0.2782894000, 0.2783920000, 0.2784946000, 0.2785972000, 0.2786999000, 0.2788026000,
2298 0.2789053000, 0.2790080000, 0.2791107000, 0.2792135000, 0.2793163000, 0.2794192000,
2299 0.2795220000, 0.2796249000, 0.2797278000, 0.2798307000, 0.2799337000, 0.2800366000,
2300 0.2801396000, 0.2802427000, 0.2803457000, 0.2804488000, 0.2805519000, 0.2806550000,
2301 0.2807581000, 0.2808613000, 0.2809645000, 0.2810677000, 0.2811709000, 0.2812742000,
2302 0.2813774000, 0.2814807000, 0.2815841000, 0.2816874000, 0.2817908000, 0.2818942000,
2303 0.2819976000, 0.2821010000, 0.2822044000, 0.2823079000, 0.2824114000, 0.2825149000,
2304 0.2826185000, 0.2827220000, 0.2828256000, 0.2829292000, 0.2830328000, 0.2831365000,
2305 0.2832401000, 0.2833438000, 0.2834475000, 0.2835512000, 0.2836550000, 0.2837588000,
2306 0.2838625000, 0.2839663000, 0.2840702000, 0.2841740000, 0.2842779000, 0.2843817000,
2307 0.2844856000, 0.2845896000, 0.2846935000, 0.2847975000, 0.2849014000, 0.2850054000,
2308 0.2851094000, 0.2852135000, 0.2853175000, 0.2854216000, 0.2855257000, 0.2856298000,
2309 0.2857339000, 0.2858380000, 0.2859422000, 0.2860464000, 0.2861505000, 0.2862548000,
2310 0.2863590000, 0.2864632000, 0.2865675000, 0.2866717000, 0.2867760000, 0.2868803000,
2311 0.2869847000, 0.2870890000, 0.2871933000, 0.2872977000, 0.2874021000, 0.2875065000,
2312 0.2876109000, 0.2877153000, 0.2878198000, 0.2879242000, 0.2880287000, 0.2881332000,
2313 0.2882377000, 0.2883422000, 0.2884468000, 0.2885513000, 0.2886559000, 0.2887604000,
2314 0.2888650000, 0.2889696000, 0.2890742000, 0.2891789000, 0.2892835000, 0.2893882000,
2315 0.2894928000, 0.2895975000, 0.2897022000, 0.2898069000, 0.2899116000, 0.2900163000,
2316 0.2901211000, 0.2902258000, 0.2903306000, 0.2904354000, 0.2905402000, 0.2906449000,
2317 0.2907498000, 0.2908546000, 0.2909594000, 0.2910642000, 0.2911691000, 0.2912740000,
2318 0.2913788000, 0.2914837000 };
2319 double width_K1;
2320 if ( mass < 1.35 ) { width_K1 = K1270_width[iii]; }
2321 else { width_K1 = K1400_width[iii]; }
2322 a[0] = 1;
2323 a[1] = 0;
2324 b[0] = mass2 - sa;
2325 // b[1] = -mass*width*wid(mass2,mass,sa,sb,sc,r2,l);
2326 b[1] = -mass * width_K1;
2327 Com_Divide( a, b, prop );
2328}
2329void EvtD0ToKSKppimpi0::propagatorGS( double mass2, double mass, double width, double sa,
2330 double sb, double sc, double r2, double prop[2] ) {
2331
2332 double GS1 = 0.636619783;
2333 double GS2 = 0.01860182466;
2334 double GS3 = 0.1591549458;
2335 double GS4 = 0.00620060822;
2336 double a[2], b[2];
2337 double tmp = sb - sc;
2338 double tmp1 = sa + tmp;
2339 double q2 = fabs( 0.25 * tmp1 * tmp1 / sa - sb );
2340 // if(q2<0) q2 = 1e-16;
2341
2342 double tmp2 = mass2 + tmp;
2343 double q02 = fabs( 0.25 * tmp2 * tmp2 / mass2 - sb );
2344 // if(q02<0) q02 = 1e-16;
2345
2346 double q = sqrt( q2 );
2347 double q0 = sqrt( q02 );
2348 double m = sqrt( sa );
2349 double q03 = q0 * q02;
2350 double tmp3 = log( mass + 2 * q0 ) + 1.2760418309; // log(mass_2Pion) = 1.2760418309;
2351
2352 double h = GS1 * q / m * ( log( m + 2 * q ) + 1.2760418309 );
2353 double h0 = GS1 * q0 / mass * tmp3;
2354 double dh = h0 * ( 0.125 / q02 - 0.5 / mass2 ) + GS3 / mass2;
2355 double d = GS2 / q02 * tmp3 + GS3 * mass / q0 - GS4 * mass / q03;
2356 double f = mass2 / q03 * ( q2 * ( h - h0 ) + ( mass2 - sa ) * q02 * dh );
2357
2358 a[0] = 1.0 + d * width / mass;
2359 a[1] = 0.0;
2360 b[0] = mass2 - sa + width * f;
2361 b[1] = -mass * width * widl1( mass2, mass, sa, sb, sc, r2 );
2362 Com_Divide( a, b, prop );
2363}
2364void EvtD0ToKSKppimpi0::rhoab( double sa, double sb, double sc, double res[2] ) {
2365 double tmp = sa + sb - sc;
2366 double q = 0.25 * tmp * tmp / sa - sb;
2367 if ( q >= 0 )
2368 {
2369 res[0] = 2.0 * sqrt( q / sa );
2370 res[1] = 0.0;
2371 }
2372 else
2373 {
2374 res[0] = 0.0;
2375 res[1] = 2.0 * sqrt( -q / sa );
2376 }
2377}
2378void EvtD0ToKSKppimpi0::rho4Pi( double sa, double res[2] ) {
2379 double temp = 1.0 - 0.3116765584 / sa; // 0.3116765584=0.13957*0.13957*16
2380 if ( temp >= 0 )
2381 {
2382 res[0] = sqrt( temp ) / ( 1.0 + exp( 9.8 - 3.5 * sa ) );
2383 res[1] = 0.0;
2384 }
2385 else
2386 {
2387 res[0] = 0.0;
2388 res[1] = sqrt( -temp ) / ( 1.0 + exp( 9.8 - 3.5 * sa ) );
2389 }
2390}
2391void EvtD0ToKSKppimpi0::propagatorsigma500( double sa, double sb, double sc, double prop[2] ) {
2392 double f = 0.5843 + 1.6663 * sa;
2393 const double M = 0.9264;
2394 const double mass2 = 0.85821696; // M*M
2395 const double mpi2d2 = 0.00973989245;
2396 double g1 = f * ( sa - mpi2d2 ) / ( mass2 - mpi2d2 ) * exp( ( mass2 - sa ) / 1.082 );
2397 double rho1s[2], rho1M[2], rho2s[2], rho2M[2], rho1[2], rho2[2];
2398 rhoab( sa, sb, sc, rho1s );
2399 rhoab( mass2, sb, sc, rho1M );
2400 rho4Pi( sa, rho2s );
2401 rho4Pi( mass2, rho2M );
2402 Com_Divide( rho1s, rho1M, rho1 );
2403 Com_Divide( rho2s, rho2M, rho2 );
2404 double a[2], b[2];
2405 a[0] = 1.0;
2406 a[1] = 0.0;
2407 b[0] = mass2 - sa + M * ( g1 * rho1[1] + 0.0024 * rho2[1] );
2408 b[1] = -M * ( g1 * rho1[0] + 0.0024 * rho2[0] );
2409 Com_Divide( a, b, prop );
2410}
2411void EvtD0ToKSKppimpi0::Flatte_rhoab( double sa, double sb, double sc, double rho[2] ) {
2412 double q = ( sa + sb - sc ) * ( sa + sb - sc ) / ( 4 * sa ) - sb;
2413 if ( q > 0 )
2414 {
2415 rho[0] = 2 * sqrt( q / sa );
2416 rho[1] = 0;
2417 }
2418 else if ( q < 0 )
2419 {
2420 rho[0] = 0;
2421 rho[1] = 2 * sqrt( -q / sa );
2422 }
2423}
2424void EvtD0ToKSKppimpi0::propagator980( double mass, double sx, double* sb, double* sc,
2425 double prop[2] ) {
2426 double unit[2] = { 1.0 };
2427 double ci[2] = { 0, 1 };
2428 double rho1[2];
2429 Flatte_rhoab( sx, sb[0], sc[0], rho1 );
2430 double rho2[2];
2431 Flatte_rhoab( sx, sb[1], sc[1], rho2 );
2432 double gK_f980 = 0.69465, gPi_f980 = 0.165;
2433 double tmp1[2] = { gK_f980, 0 };
2434 double tmp11[2];
2435 double tmp2[2] = { gPi_f980, 0 };
2436 double tmp22[2];
2437 Com_Multi( tmp1, rho1, tmp11 );
2438 Com_Multi( tmp2, rho2, tmp22 );
2439 double tmp3[2] = { tmp11[0] + tmp22[0], tmp11[1] + tmp22[1] };
2440 double tmp31[2];
2441 Com_Multi( tmp3, ci, tmp31 );
2442 double tmp4[2] = { mass * mass - sx - tmp31[0], -1.0 * tmp31[1] };
2443 Com_Divide( unit, tmp4, prop );
2444}
2445void EvtD0ToKSKppimpi0::propagatorFlatte( double mass, double width, double sa,
2446 double prop[2] ) {
2447 double q2_Pi, q2_Ka;
2448 double rhoPi[2], rhoKa[2];
2449 q2_Pi = 0.25 * sa - mPi * mPi;
2450 q2_Ka = 0.25 * sa - mKa * mKa;
2451 if ( q2_Pi > 0 )
2452 {
2453 rhoPi[0] = 2.0 * sqrt( q2_Pi / sa );
2454 rhoPi[1] = 0.0;
2455 }
2456 if ( q2_Pi <= 0 )
2457 {
2458 rhoPi[0] = 0.0;
2459 rhoPi[1] = 2.0 * sqrt( -q2_Pi / sa );
2460 }
2461 if ( q2_Ka > 0 )
2462 {
2463 rhoKa[0] = 2.0 * sqrt( q2_Ka / sa );
2464 rhoKa[1] = 0.0;
2465 }
2466 if ( q2_Ka <= 0 )
2467 {
2468 rhoKa[0] = 0.0;
2469 rhoKa[1] = 2.0 * sqrt( -q2_Ka / sa );
2470 }
2471 double a[2], b[2];
2472 a[0] = 1;
2473 a[1] = 0;
2474 b[0] = mass * mass - sa + 0.165 * rhoPi[1] + 0.69465 * rhoKa[1];
2475 b[1] = -( 0.165 * rhoPi[0] + 0.69465 * rhoKa[0] );
2476 Com_Divide( a, b, prop );
2477}
2478void EvtD0ToKSKppimpi0::KPiSLASS( double sa, double sb, double sc, double prop[2] ) {
2479 const double m1430 = 1.441;
2480 const double sa0 = 2.076481; // m1430*m1430;
2481 const double w1430 = 0.193;
2482 const double Lass1 = 0.25 / sa0;
2483 double tmp = sb - sc;
2484 double tmp1 = sa0 + tmp;
2485 double q0 = fabs( Lass1 * tmp1 * tmp1 - sb );
2486 // if(q0<0) q0 = 1e-16;
2487 double tmp2 = sa + tmp;
2488 double qs = fabs( 0.25 * tmp2 * tmp2 / sa - sb );
2489 double q = sqrt( qs );
2490 double width = w1430 * q * m1430 / sqrt( sa * q0 );
2491 double temp_R = atan( m1430 * width / ( sa0 - sa ) );
2492 if ( temp_R < 0 ) temp_R += math_pi;
2493 double deltaR = -109.7 * math_pi / 180.0 + temp_R;
2494 double temp_F =
2495 atan( 0.226 * q / ( 2.0 - 3.8194 * qs ) ); // 2.0*0.113 = 0.226; -33.8*0.113 = -3.8194
2496 if ( temp_F < 0 ) temp_F += math_pi;
2497 double deltaF = 0.1 * math_pi / 180.0 + temp_F;
2498 double deltaS = deltaR + 2.0 * deltaF;
2499 double t1 = 0.96 * sin( deltaF );
2500 double t2 = sin( deltaR );
2501 double CF[2], CS[2];
2502 CF[0] = cos( deltaF );
2503 CF[1] = sin( deltaF );
2504 CS[0] = cos( deltaS );
2505 CS[1] = sin( deltaS );
2506 prop[0] = t1 * CF[0] + t2 * CS[0];
2507 prop[1] = t1 * CF[1] + t2 * CS[1];
2508}
mass
double mass
Definition CosmicGenerator.cxx:128
f
TFile f("ana_bhabha660a_dqa_mcPat_zy_old.root")
g1
TF1 * g1
Definition DataBase/tau_mode.c:58
exp
EvtComplex exp(const EvtComplex &c)
Definition EvtComplex.hh:209
K1270_width
double K1270_width[]
Definition EvtDToKmPipPipPi0.cc:364
K1400_width
double K1400_width[]
Definition EvtDToKmPipPipPi0.cc:650
a1_width
double a1_width[]
Definition EvtDToKmPipPipPi0.cc:47
unit
*******INTEGER m_nBinMax INTEGER m_NdiMax !No of bins in histogram for cell exploration division $ !Last vertex $ !Last active cell $ !Last cell in buffer $ !No of sampling when dividing cell $ !No of function total $ !Flag for random ceel for $ !Flag for type of for WtMax $ !Flag which decides whether vertices are included in the sampling $ entire domain is hyp !Maximum effective eevents per saves r n generator level $ !Flag for chat level in !Latex Output unit
Definition FoamA.h:90
sin
double sin(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:185
cos
double cos(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:187
q
****INTEGER imax DOUBLE PRECISION m_pi *DOUBLE PRECISION m_amfin DOUBLE PRECISION m_Chfin DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_sinw2 DOUBLE PRECISION m_GFermi DOUBLE PRECISION m_MfinMin DOUBLE PRECISION m_ta2 INTEGER m_out INTEGER m_KeyFSR INTEGER m_KeyQCD *COMMON c_Semalib $ !copy of input $ !CMS energy $ !beam mass $ !final mass $ !beam charge $ !final charge $ !smallest final mass $ !Z mass $ !Z width $ !EW mixing angle $ !Gmu Fermi $ alphaQED at q
Definition KKsem.h:33
CS
TCrossPart * CS
Definition Mcgpj.cxx:51
EvtD0ToKSKppimpi0::clone
EvtDecayBase * clone()
Definition EvtD0ToKSKppimpi0.cc:36
EvtD0ToKSKppimpi0::getName
void getName(std::string &name)
Definition EvtD0ToKSKppimpi0.cc:34
EvtD0ToKSKppimpi0::decay
void decay(EvtParticle *p)
Definition EvtD0ToKSKppimpi0.cc:123
EvtDecayBase::checkSpinParent
void checkSpinParent(EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:479
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition EvtDecayBase.cc:261
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition EvtDecayBase.cc:465
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:68
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:449
EvtDecayProb::setProb
void setProb(double prob)
Definition EvtDecayProb.hh:40
EvtPDL::getStdHep
static int getStdHep(EvtId id)
Definition EvtPDL.hh:61
EvtParticle::getId
EvtId getId() const
Definition EvtParticle.cc:100
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:108
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:80
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1042
EvtVector4R::get
double get(int i) const
Definition EvtVector4R.hh:171
t
int t()
Definition t.c:1