| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // INCL++ intra-nuclear cascade model
27 // Alain Boudard, CEA-Saclay, France
28 // Joseph Cugnon, University of Liege, Belgium
29 // Jean-Christophe David, CEA-Saclay, France
30 // Pekka Kaitaniemi, CEA-Saclay, France, and H
31 // Sylvie Leray, CEA-Saclay, France
32 // Davide Mancusi, CEA-Saclay, France
33 //
34 #define INCLXX_IN_GEANT4_MODE 1
35
36 #include "globals.hh"
37
38 #include "G4INCLNKbToSpiChannel.hh"
39 #include "G4INCLKinematicsUtils.hh"
40 #include "G4INCLBinaryCollisionAvatar.hh"
41 #include "G4INCLRandom.hh"
42 #include "G4INCLGlobals.hh"
43 #include "G4INCLLogger.hh"
44 #include <algorithm>
45 #include "G4INCLPhaseSpaceGenerator.hh"
46
47 namespace G4INCL {
48
49 NKbToSpiChannel::NKbToSpiChannel(Particle *p
50 : particle1(p1), particle2(p2)
51 {}
52
53 NKbToSpiChannel::~NKbToSpiChannel(){}
54
55 void NKbToSpiChannel::fillFinalState(FinalSt
56
57 // p K0b -> S+ pi0 (2/3)
58 // p K0b -> S0 pi+ (2/3)
59 //
60 // p K- -> S+ pi- (1)
61 // p K- -> S0 pi0 (1/2)
62 // p K- -> S- pi+ (2/3)
63
64 Particle *nucleon;
65 Particle *kaon;
66
67 if(particle1->isNucleon()){
68 nucleon = particle1;
69 kaon = particle2;
70 }
71 else{
72 nucleon = particle2;
73 kaon = particle1;
74 }
75
76 ParticleType SigmaType, PionType;
77
78 const G4int iso = ParticleTable::getIsospi
79 const G4int iso_n = ParticleTable::getIsos
80 const G4double rdm = Random::shoot();
81
82 G4int WhichChannel = 99;
83
84 if(iso == 2 || iso == -2){
85 if(rdm < 0.5){
86 PionType = ParticleTable::getPionType(
87 SigmaType = SigmaZero;
88 WhichChannel = 0;
89 }
90 else{
91 SigmaType = ParticleTable::getSigmaTyp
92 PionType = PiZero;
93 WhichChannel = iso/2;
94 }
95 }
96 else{
97 if(rdm*13. < 6.){
98 PionType = ParticleTable::getPionType(
99 SigmaType = ParticleTable::getSigmaTyp
100 WhichChannel = iso_n;
101 }
102 else if(rdm*13. < 9.){
103 PionType = PiZero;
104 SigmaType = SigmaZero;
105 WhichChannel = 0;
106 }
107 else{
108 PionType = ParticleTable::getPionType(
109 SigmaType = ParticleTable::getSigmaTyp
110 WhichChannel = -iso_n;
111 }
112 }
113
114 ThreeVector mom_pion = KaonMomentum(kaon,n
115
116 kaon->setType(PionType);
117 nucleon->setType(SigmaType);
118
119 G4double norm = KinematicsUtils::momentumI
120
121 kaon->setMomentum(mom_pion*norm);
122 nucleon->setMomentum(-mom_pion*norm);
123
124 kaon->adjustEnergyFromMomentum();
125 nucleon->adjustEnergyFromMomentum();
126
127 fs->addModifiedParticle(nucleon);
128 fs->addModifiedParticle(kaon);
129
130 }
131
132 ThreeVector NKbToSpiChannel::KaonMomentum(Pa
133
134 const G4double pLab = KinematicsUtils::mom
135 G4double max = 0;
136
137 if(pLab < 435.) return Random::normVector(
138
139 if(WhichChannel == 0) max = 865;
140 else if(std::abs(WhichChannel) == 1) max =
141 else{
142 INCL_WARN("Bad channel chosen in Kb N To
143 return Random::normVector();
144 }
145
146 G4double cos_theta = 1.;
147 G4double sin_theta = 0.;
148 const G4double cos_phi = std::cos(Random::
149 const G4double sin_phi = std::sqrt(1-cos_p
150
151 const G4double x = kaon->getMomentum().get
152 const G4double y = kaon->getMomentum().get
153 const G4double z = kaon->getMomentum().get
154
155 const G4double r = std::sqrt(x*x+y*y+z*z);
156 const G4double rho = std::sqrt(x*x+y*y);
157
158 if(pLab >= max){
159 const G4double b = 12. * pLab/2375.; //
160 cos_theta = std::log(Random::shoot()*(st
161 sin_theta = std::sqrt(1-cos_theta*cos_th
162
163 }
164 else{
165 const G4double Legendre_coef[3][283][9]
166 {// K- p -> Sigma- pi+
167 {435,0.20268,0.89605,0.00315,-0.02043,
168 {440,0.2116,0.91169,0.00737,-0.02297,2
169 {445,0.21976,0.92725,0.01151,-0.02554,
170 {450,0.22691,0.94264,0.01546,-0.02817,
171 {455,0.23293,0.95777,0.0191,-0.03089,4
172 {460,0.23762,0.97259,0.02233,-0.03372,
173 {465,0.24084,0.98722,0.025,-0.03662,0.
174 {470,0.24244,1.00183,0.02695,-0.0395,0
175 {475,0.24243,1.0166,0.02804,-0.0423,0.
176 {480,0.2408,1.03171,0.0281,-0.04494,0.
177 {485,0.23759,1.04732,0.02704,-0.04728,
178 {490,0.23289,1.06364,0.02475,-0.04917,
179 {495,0.22693,1.08085,0.02112,-0.05045,
180 {500,0.21985,1.0991,0.01607,-0.05099,0
181 {505,0.21186,1.11838,0.00962,-0.05067,
182 {510,0.20311,1.13865,0.0018,-0.04938,0
183 {515,0.19378,1.15987,-0.00734,-0.04703
184 {520,0.18394,1.18199,-0.01776,-0.04353
185 {525,0.17362,1.20497,-0.02941,-0.03891
186 {530,0.16277,1.22874,-0.04223,-0.03321
187 {535,0.15135,1.25326,-0.05615,-0.02646
188 {540,0.13922,1.2785,-0.07107,-0.01873,
189 {545,0.12624,1.30447,-0.08684,-0.01017
190 {550,0.1123,1.3312,-0.10326,-0.00094,0
191 {555,0.09727,1.35871,-0.12015,0.00879,
192 {560,0.08117,1.38699,-0.13734,0.01891,
193 {565,0.064,1.41599,-0.15466,0.02935,-0
194 {570,0.0459,1.44562,-0.17197,0.04009,-
195 {575,0.02704,1.47583,-0.18911,0.05108,
196 {580,0.00769,1.50654,-0.20595,0.06228,
197 {585,-0.01183,1.53766,-0.22246,0.07364
198 {590,-0.03118,1.56911,-0.23858,0.08511
199 {595,-0.05002,1.60083,-0.25425,0.09665
200 {600,-0.06802,1.63271,-0.26943,0.10821
201 {605,-0.08492,1.66463,-0.28414,0.1197,
202 {610,-0.10045,1.69641,-0.29846,0.131,-
203 {615,-0.11445,1.72786,-0.31248,0.14199
204 {620,-0.12685,1.75882,-0.32626,0.15255
205 {625,-0.13761,1.78907,-0.3398,0.16263,
206 {630,-0.14665,1.81843,-0.353,0.17219,-
207 {635,-0.15404,1.84667,-0.36578,0.18122
208 {640,-0.15984,1.8736,-0.37806,0.18969,
209 {645,-0.1641,1.89903,-0.38974,0.19753,
210 {650,-0.1668,1.92283,-0.40077,0.20457,
211 {655,-0.16802,1.94484,-0.41108,0.2107,
212 {660,-0.16778,1.96491,-0.4206,0.21575,
213 {665,-0.16609,1.98295,-0.42932,0.21971
214 {670,-0.16293,1.99895,-0.4373,0.22266,
215 {675,-0.15835,2.01289,-0.4446,0.2247,-
216 {680,-0.15235,2.02475,-0.45128,0.22593
217 {685,-0.14502,2.03452,-0.45737,0.22638
218 {690,-0.1364,2.04219,-0.46282,0.22601,
219 {695,-0.12665,2.04773,-0.46761,0.22477
220 {700,-0.11594,2.05112,-0.47172,0.22263
221 {705,-0.10451,2.05233,-0.47513,0.21958
222 {710,-0.09257,2.05134,-0.47789,0.21563
223 {715,-0.08045,2.0481,-0.48004,0.21082,
224 {720,-0.06841,2.04258,-0.48164,0.20519
225 {725,-0.05677,2.03476,-0.4827,0.19876,
226 {730,-0.04577,2.02468,-0.48321,0.19159
227 {735,-0.03567,2.01241,-0.48312,0.18371
228 {740,-0.02666,1.99799,-0.4824,0.17518,
229 {745,-0.01887,1.98148,-0.48101,0.16603
230 {750,-0.01239,1.96293,-0.47891,0.15626
231 {755,-0.00733,1.9424,-0.47608,0.14587,
232 {760,-0.00369,1.91993,-0.47246,0.13483
233 {765,-0.0015,1.8956,-0.46801,0.12316,0
234 {770,-0.00083,1.86947,-0.46257,0.1109,
235 {775,-0.00169,1.84167,-0.45593,0.09812
236 {780,-0.0042,1.81233,-0.44791,0.0849,0
237 {785,-0.0085,1.78162,-0.43841,0.07136,
238 {790,-0.0149,1.74974,-0.4273,0.05758,0
239 {795,-0.02359,1.71687,-0.41448,0.04365
240 {800,-0.03492,1.68318,-0.3999,0.02967,
241 {805,-0.04918,1.64882,-0.38363,0.01571
242 {810,-0.06654,1.61392,-0.36573,0.00182
243 {815,-0.08712,1.57864,-0.34628,-0.0119
244 {820,-0.11012,1.54311,-0.32533,-0.0254
245 {825,-0.13493,1.50749,-0.30295,-0.0387
246 {830,-0.16027,1.47191,-0.2792,-0.05167
247 {835,-0.18509,1.43652,-0.25416,-0.0642
248 {840,-0.20852,1.40147,-0.22789,-0.0762
249 {845,-0.23013,1.3669,-0.20046,-0.08779
250 {850,-0.2499,1.33295,-0.17192,-0.09871
251 {855,-0.26818,1.29977,-0.14236,-0.1089
252 {860,-0.2854,1.2675,-0.11182,-0.11851,
253 {865,-0.30205,1.23628,-0.08039,-0.1272
254 {870,-0.31869,1.20621,-0.0481,-0.13519
255 {875,-0.33574,1.17732,-0.01498,-0.1423
256 {880,-0.35356,1.14966,0.01893,-0.14877
257 {885,-0.37245,1.12329,0.05362,-0.15447
258 {890,-0.3926,1.09818,0.08898,-0.15955,
259 {895,-0.41404,1.07429,0.12481,-0.16417
260 {900,-0.43679,1.05155,0.16084,-0.16846
261 {905,-0.46071,1.02994,0.1968,-0.17253,
262 {910,-0.4856,1.00944,0.23246,-0.17643,
263 {915,-0.51107,0.99007,0.26757,-0.18017
264 {920,-0.53674,0.97186,0.3019,-0.18373,
265 {925,-0.56217,0.95486,0.33525,-0.18707
266 {930,-0.58681,0.9392,0.36744,-0.19012,
267 {935,-0.61023,0.92501,0.39833,-0.1928,
268 {940,-0.63197,0.91241,0.42775,-0.19502
269 {945,-0.65175,0.90149,0.45562,-0.19679
270 {950,-0.66926,0.89231,0.4819,-0.19816,
271 {955,-0.68439,0.88494,0.50655,-0.19918
272 {960,-0.69707,0.87945,0.52952,-0.19989
273 {965,-0.70731,0.87589,0.55076,-0.20033
274 {970,-0.71516,0.87419,0.5703,-0.20045,
275 {975,-0.7207,0.87429,0.58813,-0.20021,
276 {980,-0.72394,0.87607,0.60419,-0.1995,
277 {985,-0.72503,0.87939,0.61843,-0.19819
278 {990,-0.72394,0.88412,0.63075,-0.19616
279 {995,-0.72067,0.89012,0.64109,-0.19327
280 {1000,-0.71516,0.89727,0.64936,-0.1894
281 {1005,-0.70731,0.90541,0.6555,-0.18445
282 {1010,-0.697,0.91431,0.65947,-0.17846,
283 {1015,-0.68411,0.92376,0.66132,-0.1714
284 {1020,-0.66835,0.93356,0.66107,-0.1634
285 {1025,-0.64976,0.94348,0.65877,-0.1545
286 {1030,-0.62823,0.95334,0.65446,-0.1445
287 {1035,-0.6038,0.96292,0.64817,-0.13375
288 {1040,-0.57654,0.97201,0.63996,-0.1219
289 {1045,-0.54672,0.9804,0.62985,-0.10933
290 {1050,-0.51475,0.98795,0.61792,-0.0957
291 {1055,-0.48114,0.99462,0.60434,-0.0812
292 {1060,-0.44627,1.00047,0.58928,-0.0657
293 {1065,-0.41097,1.00553,0.57295,-0.0492
294 {1070,-0.37567,1.00983,0.55554,-0.0315
295 {1075,-0.34091,1.01341,0.53725,-0.0127
296 {1080,-0.30701,1.0163,0.51825,0.00733,
297 {1085,-0.27413,1.01854,0.49874,0.0287,
298 {1090,-0.2425,1.02016,0.47892,0.0514,1
299 {1095,-0.21201,1.02123,0.45893,0.0754,
300 {1100,-0.18239,1.02178,0.43894,0.10065
301 {1105,-0.15367,1.02189,0.41909,0.12708
302 {1110,-0.1254,1.02159,0.39954,0.15464,
303 {1115,-0.09757,1.02096,0.38043,0.18328
304 {1120,-0.06982,1.02003,0.36193,0.21294
305 {1125,-0.04213,1.01887,0.34419,0.24357
306 {1130,-0.0142,1.01753,0.32732,0.27509,
307 {1135,0.01406,1.01605,0.31137,0.30741,
308 {1140,0.04282,1.01449,0.29633,0.3404,1
309 {1145,0.0719,1.0129,0.28221,0.37395,1.
310 {1150,0.10146,1.01132,0.26901,0.40794,
311 {1155,0.13127,1.00979,0.25673,0.44225,
312 {1160,0.1612,1.00837,0.24537,0.47676,1
313 {1165,0.19131,1.00711,0.23494,0.51136,
314 {1170,0.22156,1.00605,0.22546,0.54595,
315 {1175,0.25183,1.00528,0.21695,0.58048,
316 {1180,0.28191,1.00489,0.20947,0.61496,
317 {1185,0.31173,1.00498,0.20307,0.64937,
318 {1190,0.34118,1.00562,0.1978,0.68368,1
319 {1195,0.37021,1.0069,0.1937,0.71789,1.
320 {1200,0.39843,1.00892,0.19083,0.75199,
321 {1205,0.42584,1.01177,0.18923,0.78595,
322 {1210,0.45231,1.01553,0.18896,0.81978,
323 {1215,0.47765,1.02033,0.19014,0.85355,
324 {1220,0.50177,1.02629,0.19287,0.88735,
325 {1225,0.52457,1.03352,0.19729,0.92127,
326 {1230,0.54599,1.04214,0.2035,0.95541,1
327 {1235,0.56609,1.05228,0.21163,0.98984,
328 {1240,0.5847,1.06406,0.2218,1.02465,1.
329 {1245,0.60195,1.07759,0.23413,1.05995,
330 {1250,0.61786,1.09296,0.24867,1.09575,
331 {1255,0.63245,1.11008,0.26523,1.13191,
332 {1260,0.64572,1.1288,0.28356,1.16821,1
333 {1265,0.65767,1.14902,0.30341,1.20443,
334 {1270,0.6683,1.17065,0.32457,1.24045,1
335 {1275,0.67758,1.19375,0.34687,1.27626,
336 {1280,0.68536,1.21833,0.37014,1.31188,
337 {1285,0.69167,1.24444,0.39422,1.34729,
338 {1290,0.69645,1.27212,0.41894,1.38249,
339 {1295,0.69969,1.30139,0.44417,1.41745,
340 {1300,0.70139,1.33231,0.46975,1.45211,
341 {1305,0.70159,1.36489,0.49554,1.48643,
342 {1310,0.70037,1.39918,0.52138,1.52037,
343 {1315,0.69788,1.43522,0.54714,1.55387,
344 {1320,0.69426,1.47301,0.57267,1.5869,1
345 {1325,0.68978,1.51246,0.59789,1.61947,
346 {1330,0.68469,1.55333,0.62279,1.65156,
347 {1335,0.6793,1.59543,0.64739,1.68314,1
348 {1340,0.67393,1.63853,0.67168,1.71418,
349 {1345,0.66891,1.68242,0.69567,1.74465,
350 {1350,0.66455,1.72688,0.71937,1.77453,
351 {1355,0.66117,1.7717,0.74277,1.80377,1
352 {1360,0.65897,1.81664,0.76591,1.83233,
353 {1365,0.65815,1.86145,0.78889,1.86002,
354 {1370,0.65879,1.90585,0.81186,1.88663,
355 {1375,0.66083,1.94962,0.8348,1.91201,1
356 {1380,0.66406,1.9926,0.85756,1.93612,1
357 {1385,0.66812,2.03463,0.87996,1.95888,
358 {1390,0.67247,2.07556,0.90183,1.98025,
359 {1395,0.67646,2.11524,0.92297,2.00017,
360 {1400,0.67917,2.15351,0.94322,2.01858,
361 {1405,0.67999,2.19022,0.96239,2.03543,
362 {1410,0.67813,2.2252,0.98032,2.05066,1
363 {1415,0.67307,2.25832,0.99681,2.06423,
364 {1420,0.66456,2.28945,1.01176,2.07609,
365 {1425,0.65267,2.31863,1.02525,2.08636,
366 {1430,0.6378,2.34595,1.03744,2.09517,1
367 {1435,0.62062,2.37149,1.04848,2.10265,
368 {1440,0.60209,2.39534,1.05852,2.10894,
369 {1445,0.5832,2.41757,1.06772,2.11416,1
370 {1450,0.56496,2.43828,1.07622,2.11846,
371 {1455,0.54826,2.45754,1.08418,2.12197,
372 {1460,0.53374,2.47544,1.09175,2.12481,
373 {1465,0.52181,2.49206,1.09908,2.12713,
374 {1470,0.51254,2.50744,1.10619,2.12895,
375 {1475,0.50575,2.52157,1.11304,2.13027,
376 {1480,0.50085,2.53446,1.11958,2.13104,
377 {1485,0.49715,2.5461,1.12576,2.13126,1
378 {1490,0.49385,2.55652,1.13152,2.13089,
379 {1495,0.49019,2.56569,1.13682,2.12991,
380 {1500,0.48551,2.57364,1.14159,2.1283,1
381 {1505,0.47939,2.58035,1.14579,2.12603,
382 {1510,0.47148,2.58582,1.14937,2.12307,
383 {1515,0.46175,2.59007,1.15228,2.11941,
384 {1520,0.45033,2.59313,1.15453,2.11509,
385 {1525,0.43753,2.5951,1.15631,2.11032,1
386 {1530,0.42376,2.59609,1.15782,2.10533,
387 {1535,0.40952,2.59624,1.15926,2.10037,
388 {1540,0.39543,2.59565,1.16083,2.09565,
389 {1545,0.38215,2.59444,1.16275,2.09142,
390 {1550,0.37039,2.59273,1.1652,2.08789,1
391 {1555,0.36091,2.59062,1.16833,2.08526,
392 {1560,0.35456,2.58822,1.17224,2.08366,
393 {1565,0.35177,2.58562,1.17705,2.08327,
394 {1570,0.35337,2.58293,1.18288,2.08422,
395 {1575,0.35961,2.58025,1.18983,2.08669,
396 {1580,0.37063,2.57768,1.19803,2.09082,
397 {1585,0.3863,2.57533,1.20758,2.09678,1
398 {1590,0.40627,2.57329,1.2186,2.10471,1
399 {1595,0.43008,2.57166,1.2312,2.11477,1
400 {1600,0.45695,2.57056,1.2455,2.12712,1
401 {1605,0.48626,2.57007,1.2616,2.14191,1
402 {1610,0.51735,2.57029,1.27959,2.15926,
403 {1615,0.54962,2.57116,1.2992,2.17896,1
404 {1620,0.58245,2.57256,1.32004,2.20065,
405 {1625,0.61536,2.5744,1.34172,2.224,1.3
406 {1630,0.64786,2.57654,1.36383,2.24863,
407 {1635,0.67945,2.57889,1.38599,2.27422,
408 {1640,0.70969,2.58132,1.40779,2.30039,
409 {1645,0.73811,2.58373,1.42884,2.32681,
410 {1650,0.7643,2.58599,1.44873,2.35312,1
411 {1655,0.78784,2.58801,1.46709,2.37897,
412 {1660,0.80845,2.58964,1.4836,2.40404,1
413 {1665,0.8259,2.59072,1.49807,2.42802,1
414 {1670,0.84003,2.5911,1.51034,2.45063,1
415 {1675,0.85087,2.59061,1.52021,2.47156,
416 {1680,0.85853,2.58908,1.52752,2.49054,
417 {1685,0.86321,2.58636,1.53208,2.50725,
418 {1690,0.86521,2.58228,1.53373,2.52141,
419 {1695,0.86497,2.57669,1.53227,2.53272,
420 {1700,0.86261,2.56941,1.52753,2.54089,
421 {1705,0.85865,2.5603,1.51934,2.54562,1
422 {1710,0.85338,2.54921,1.50758,2.54666,
423 {1715,0.84709,2.53615,1.49237,2.54388,
424 {1720,0.83998,2.52115,1.4739,2.5372,1.
425 {1725,0.83223,2.50424,1.45236,2.52653,
426 {1730,0.82398,2.48546,1.42794,2.51179,
427 {1735,0.81532,2.46484,1.40083,2.49289,
428 {1740,0.80635,2.4424,1.37121,2.46973,1
429 {1745,0.79717,2.4182,1.33927,2.44227,1
430 {1750,0.78787,2.39228,1.30516,2.41059,
431 {1755,0.77864,2.36474,1.26903,2.37487,
432 {1760,0.76947,2.33565,1.231,2.33528,1.
433 {1765,0.76063,2.3051,1.19124,2.29197,1
434 {1770,0.75229,2.27315,1.14986,2.24514,
435 {1775,0.74462,2.23989,1.10702,2.19493,
436 {1780,0.73776,2.20541,1.06285,2.14153,
437 {1785,0.73178,2.16977,1.0175,2.0851,1.
438 {1790,0.72674,2.13307,0.9711,2.02582,1
439 {1795,0.7226,2.09537,0.9238,1.96385,1.
440 {1800,0.7193,2.05676,0.87573,1.89936,0
441 {1805,0.71674,2.01732,0.82704,1.83253,
442 {1810,0.71485,1.97712,0.77787,1.76351,
443 {1815,0.7135,1.93625,0.72835,1.69249,0
444 {1820,0.71257,1.89479,0.67863,1.61963,
445 {1825,0.71198,1.85282,0.62885,1.5451,0
446 {1830,0.71165,1.81041,0.57915,1.46907,
447 {1835,0.7115,1.76765,0.52966,1.39172,0
448 {1840,0.71149,1.72461,0.48053,1.3132,0
449 {1845,0.71158,1.68139,0.43186,1.23376,
450 {// K- p -> Sigma0 pi0
451 {435,0.01082,1.10989,0.01663,-0.00335,
452 {440,-0.01888,1.07693,-0.00805,-0.0120
453 {445,-0.04853,1.04408,-0.03276,-0.0207
454 {450,-0.07803,1.01152,-0.0575,-0.02944
455 {455,-0.10733,0.97938,-0.0823,-0.03809
456 {460,-0.13633,0.9478,-0.10716,-0.0467,
457 {465,-0.16491,0.91681,-0.13197,-0.0552
458 {470,-0.19294,0.88642,-0.15661,-0.0635
459 {475,-0.2203,0.85665,-0.18094,-0.07154
460 {480,-0.24683,0.82754,-0.20479,-0.0792
461 {485,-0.27226,0.79925,-0.22775,-0.0863
462 {490,-0.29631,0.77198,-0.24937,-0.0926
463 {495,-0.31868,0.74592,-0.26919,-0.0980
464 {500,-0.33912,0.72125,-0.28682,-0.1022
465 {505,-0.35756,0.69815,-0.30212,-0.1048
466 {510,-0.37395,0.67677,-0.31501,-0.1054
467 {515,-0.38825,0.65729,-0.32543,-0.1037
468 {520,-0.40045,0.63982,-0.33324,-0.0996
469 {525,-0.41054,0.62436,-0.33823,-0.0934
470 {530,-0.41855,0.61093,-0.34015,-0.0858
471 {535,-0.42447,0.59951,-0.33876,-0.0771
472 {540,-0.42842,0.59011,-0.33407,-0.0679
473 {545,-0.43073,0.5827,-0.32663,-0.05813
474 {550,-0.43174,0.57727,-0.31706,-0.0478
475 {555,-0.4318,0.57379,-0.30599,-0.03697
476 {560,-0.43126,0.57219,-0.29405,-0.0257
477 {565,-0.43044,0.57227,-0.28185,-0.0143
478 {570,-0.4298,0.57397,-0.26991,-0.00328
479 {575,-0.42982,0.57729,-0.25874,0.0072,
480 {580,-0.4308,0.58214,-0.24872,0.01687,
481 {585,-0.4328,0.58835,-0.2401,0.02571,0
482 {590,-0.43612,0.59583,-0.23321,0.03376
483 {595,-0.44115,0.6045,-0.22839,0.04108,
484 {600,-0.4482,0.61432,-0.22597,0.04772,
485 {605,-0.45692,0.62534,-0.22614,0.0535,
486 {610,-0.46673,0.63752,-0.22906,0.05801
487 {615,-0.47702,0.65085,-0.23489,0.06082
488 {620,-0.48723,0.66531,-0.24375,0.06153
489 {625,-0.4968,0.68072,-0.25539,0.06011,
490 {630,-0.50526,0.69682,-0.26923,0.05678
491 {635,-0.51223,0.71341,-0.28474,0.05178
492 {640,-0.51759,0.73044,-0.30159,0.04536
493 {645,-0.52116,0.74789,-0.31929,0.0379,
494 {650,-0.52267,0.76582,-0.33727,0.02987
495 {655,-0.52186,0.78423,-0.35493,0.02172
496 {660,-0.51846,0.80318,-0.3717,0.01394,
497 {665,-0.51218,0.82275,-0.387,0.007,0,0
498 {670,-0.50273,0.84317,-0.40028,0.00144
499 {675,-0.48982,0.86469,-0.41096,-0.0022
500 {680,-0.47326,0.88743,-0.41875,-0.0034
501 {685,-0.45292,0.91119,-0.42381,-0.0022
502 {690,-0.42867,0.93557,-0.42646,0.00161
503 {695,-0.40062,0.96016,-0.42703,0.00789
504 {700,-0.36903,0.98453,-0.42585,0.01652
505 {705,-0.33405,1.00813,-0.42325,0.02726
506 {710,-0.29565,1.03016,-0.41947,0.03953
507 {715,-0.25378,1.0498,-0.41477,0.05274,
508 {720,-0.20851,1.06626,-0.40938,0.06618
509 {725,-0.16053,1.07916,-0.40351,0.07896
510 {730,-0.11059,1.08856,-0.39716,0.09047
511 {735,-0.05947,1.09453,-0.3903,0.10015,
512 {740,-0.00792,1.09717,-0.38288,0.10747
513 {745,0.0433,1.09656,-0.37485,0.11203,0
514 {750,0.09344,1.09277,-0.36603,0.11403,
515 {755,0.14178,1.0859,-0.35625,0.11383,0
516 {760,0.18761,1.07615,-0.34549,0.1117,0
517 {765,0.23037,1.06383,-0.33392,0.10794,
518 {770,0.26959,1.04936,-0.32217,0.1028,0
519 {775,0.30481,1.03313,-0.311,0.09648,0,
520 {780,0.33569,1.0155,-0.3011,0.08904,0,
521 {785,0.36209,0.9969,-0.29309,0.08029,0
522 {790,0.38386,0.97771,-0.28758,0.07006,
523 {795,0.40088,0.95833,-0.28517,0.05815,
524 {800,0.41307,0.93883,-0.28627,0.04427,
525 {805,0.42048,0.91889,-0.29099,0.02801,
526 {810,0.42339,0.89844,-0.29914,0.0094,0
527 {815,0.42217,0.87751,-0.31045,-0.01136
528 {820,0.4172,0.85614,-0.32463,-0.03409,
529 {825,0.40886,0.83439,-0.34137,-0.05857
530 {830,0.39754,0.8123,-0.36041,-0.08462,
531 {835,0.38361,0.7899,-0.38143,-0.11203,
532 {840,0.36745,0.76725,-0.40416,-0.14061
533 {845,0.34944,0.74439,-0.4283,-0.17016,
534 {850,0.32996,0.72136,-0.45357,-0.20047
535 {855,0.30939,0.6982,-0.47967,-0.23135,
536 {860,0.2881,0.67496,-0.50631,-0.26261,
537 {865,0.26649,0.65168,-0.53321,-0.29404
538 {0,0,0,0,0,0,0,0,0},
539 {0,0,0,0,0,0,0,0,0},
540 {0,0,0,0,0,0,0,0,0},
541 {0,0,0,0,0,0,0,0,0},
542 {0,0,0,0,0,0,0,0,0},
543 {0,0,0,0,0,0,0,0,0},
544 {0,0,0,0,0,0,0,0,0},
545 {0,0,0,0,0,0,0,0,0},
546 {0,0,0,0,0,0,0,0,0},
547 {0,0,0,0,0,0,0,0,0},
548 {0,0,0,0,0,0,0,0,0},
549 {0,0,0,0,0,0,0,0,0},
550 {0,0,0,0,0,0,0,0,0},
551 {0,0,0,0,0,0,0,0,0},
552 {0,0,0,0,0,0,0,0,0},
553 {0,0,0,0,0,0,0,0,0},
554 {0,0,0,0,0,0,0,0,0},
555 {0,0,0,0,0,0,0,0,0},
556 {0,0,0,0,0,0,0,0,0},
557 {0,0,0,0,0,0,0,0,0},
558 {0,0,0,0,0,0,0,0,0},
559 {0,0,0,0,0,0,0,0,0},
560 {0,0,0,0,0,0,0,0,0},
561 {0,0,0,0,0,0,0,0,0},
562 {0,0,0,0,0,0,0,0,0},
563 {0,0,0,0,0,0,0,0,0},
564 {0,0,0,0,0,0,0,0,0},
565 {0,0,0,0,0,0,0,0,0},
566 {0,0,0,0,0,0,0,0,0},
567 {0,0,0,0,0,0,0,0,0},
568 {0,0,0,0,0,0,0,0,0},
569 {0,0,0,0,0,0,0,0,0},
570 {0,0,0,0,0,0,0,0,0},
571 {0,0,0,0,0,0,0,0,0},
572 {0,0,0,0,0,0,0,0,0},
573 {0,0,0,0,0,0,0,0,0},
574 {0,0,0,0,0,0,0,0,0},
575 {0,0,0,0,0,0,0,0,0},
576 {0,0,0,0,0,0,0,0,0},
577 {0,0,0,0,0,0,0,0,0},
578 {0,0,0,0,0,0,0,0,0},
579 {0,0,0,0,0,0,0,0,0},
580 {0,0,0,0,0,0,0,0,0},
581 {0,0,0,0,0,0,0,0,0},
582 {0,0,0,0,0,0,0,0,0},
583 {0,0,0,0,0,0,0,0,0},
584 {0,0,0,0,0,0,0,0,0},
585 {0,0,0,0,0,0,0,0,0},
586 {0,0,0,0,0,0,0,0,0},
587 {0,0,0,0,0,0,0,0,0},
588 {0,0,0,0,0,0,0,0,0},
589 {0,0,0,0,0,0,0,0,0},
590 {0,0,0,0,0,0,0,0,0},
591 {0,0,0,0,0,0,0,0,0},
592 {0,0,0,0,0,0,0,0,0},
593 {0,0,0,0,0,0,0,0,0},
594 {0,0,0,0,0,0,0,0,0},
595 {0,0,0,0,0,0,0,0,0},
596 {0,0,0,0,0,0,0,0,0},
597 {0,0,0,0,0,0,0,0,0},
598 {0,0,0,0,0,0,0,0,0},
599 {0,0,0,0,0,0,0,0,0},
600 {0,0,0,0,0,0,0,0,0},
601 {0,0,0,0,0,0,0,0,0},
602 {0,0,0,0,0,0,0,0,0},
603 {0,0,0,0,0,0,0,0,0},
604 {0,0,0,0,0,0,0,0,0},
605 {0,0,0,0,0,0,0,0,0},
606 {0,0,0,0,0,0,0,0,0},
607 {0,0,0,0,0,0,0,0,0},
608 {0,0,0,0,0,0,0,0,0},
609 {0,0,0,0,0,0,0,0,0},
610 {0,0,0,0,0,0,0,0,0},
611 {0,0,0,0,0,0,0,0,0},
612 {0,0,0,0,0,0,0,0,0},
613 {0,0,0,0,0,0,0,0,0},
614 {0,0,0,0,0,0,0,0,0},
615 {0,0,0,0,0,0,0,0,0},
616 {0,0,0,0,0,0,0,0,0},
617 {0,0,0,0,0,0,0,0,0},
618 {0,0,0,0,0,0,0,0,0},
619 {0,0,0,0,0,0,0,0,0},
620 {0,0,0,0,0,0,0,0,0},
621 {0,0,0,0,0,0,0,0,0},
622 {0,0,0,0,0,0,0,0,0},
623 {0,0,0,0,0,0,0,0,0},
624 {0,0,0,0,0,0,0,0,0},
625 {0,0,0,0,0,0,0,0,0},
626 {0,0,0,0,0,0,0,0,0},
627 {0,0,0,0,0,0,0,0,0},
628 {0,0,0,0,0,0,0,0,0},
629 {0,0,0,0,0,0,0,0,0},
630 {0,0,0,0,0,0,0,0,0},
631 {0,0,0,0,0,0,0,0,0},
632 {0,0,0,0,0,0,0,0,0},
633 {0,0,0,0,0,0,0,0,0}},
634 {// K- p -> Sigma+ pi-
635 {435,-0.41506,0.51794,-0.39061,-0.1725
636 {440,-0.41233,0.48944,-0.38429,-0.1651
637 {445,-0.4096,0.46112,-0.37795,-0.15764
638 {450,-0.40685,0.4332,-0.37154,-0.15026
639 {455,-0.40408,0.40591,-0.36504,-0.1430
640 {460,-0.40131,0.37946,-0.35841,-0.1359
641 {465,-0.39874,0.35403,-0.35162,-0.1290
642 {470,-0.39658,0.32979,-0.34467,-0.1224
643 {475,-0.39505,0.30691,-0.33753,-0.1160
644 {480,-0.3944,0.28554,-0.3302,-0.11002,
645 {485,-0.39493,0.26577,-0.32275,-0.1042
646 {490,-0.39696,0.24765,-0.31529,-0.0988
647 {495,-0.4008,0.23126,-0.30789,-0.09372
648 {500,-0.40672,0.21665,-0.30065,-0.0889
649 {505,-0.41477,0.20379,-0.2936,-0.08456
650 {510,-0.42493,0.19264,-0.2868,-0.0806,
651 {515,-0.4372,0.18316,-0.28027,-0.07711
652 {520,-0.45145,0.17527,-0.27403,-0.0741
653 {525,-0.46723,0.16882,-0.268,-0.07163,
654 {530,-0.48408,0.16365,-0.26212,-0.0696
655 {535,-0.50152,0.15962,-0.25632,-0.0680
656 {540,-0.5191,0.1566,-0.25056,-0.06694,
657 {545,-0.53644,0.15456,-0.24489,-0.0662
658 {550,-0.55315,0.15349,-0.23935,-0.0660
659 {555,-0.56883,0.15338,-0.23399,-0.0662
660 {560,-0.58327,0.15419,-0.22886,-0.0669
661 {565,-0.59663,0.15588,-0.2239,-0.06794
662 {570,-0.60911,0.15838,-0.21909,-0.0693
663 {575,-0.62091,0.16164,-0.21439,-0.0711
664 {580,-0.63222,0.1656,-0.20978,-0.07333
665 {585,-0.64322,0.17023,-0.20527,-0.0758
666 {590,-0.65408,0.17551,-0.20089,-0.0786
667 {595,-0.66499,0.18138,-0.19664,-0.0818
668 {600,-0.6761,0.18784,-0.19254,-0.08532
669 {605,-0.68744,0.19486,-0.18869,-0.0891
670 {610,-0.69895,0.2025,-0.18519,-0.09324
671 {615,-0.71054,0.21077,-0.18219,-0.0976
672 {620,-0.72215,0.21971,-0.17979,-0.1024
673 {625,-0.73381,0.22936,-0.17813,-0.1075
674 {630,-0.74564,0.23978,-0.17729,-0.1129
675 {635,-0.75775,0.25101,-0.17738,-0.1187
676 {640,-0.77025,0.26312,-0.17851,-0.1249
677 {645,-0.78311,0.27609,-0.18077,-0.1316
678 {650,-0.79619,0.28987,-0.18424,-0.1386
679 {655,-0.80933,0.30439,-0.18901,-0.1462
680 {660,-0.8224,0.3196,-0.19518,-0.15426,
681 {665,-0.83513,0.3354,-0.2028,-0.16277,
682 {670,-0.84713,0.35164,-0.21197,-0.1718
683 {675,-0.85802,0.36815,-0.22273,-0.1813
684 {680,-0.86734,0.38479,-0.23516,-0.1914
685 {685,-0.87414,0.4013,-0.24915,-0.2022,
686 {690,-0.87701,0.41738,-0.26452,-0.2135
687 {695,-0.87457,0.43271,-0.28104,-0.2254
688 {700,-0.86542,0.447,-0.29851,-0.238,-0
689 {705,-0.8484,0.46,-0.31674,-0.25122,-0
690 {710,-0.82293,0.47163,-0.33558,-0.2650
691 {715,-0.7885,0.48183,-0.35485,-0.2795,
692 {720,-0.74459,0.49053,-0.37441,-0.2944
693 {725,-0.69112,0.49769,-0.39405,-0.31,-
694 {730,-0.62899,0.50329,-0.41345,-0.3260
695 {735,-0.55919,0.50731,-0.43227,-0.3424
696 {740,-0.48275,0.50974,-0.45017,-0.3594
697 {745,-0.40083,0.51061,-0.46688,-0.3767
698 {750,-0.31522,0.51019,-0.48236,-0.3943
699 {755,-0.22787,0.50878,-0.49662,-0.4122
700 {760,-0.14072,0.5067,-0.5097,-0.43031,
701 {765,-0.05569,0.50425,-0.52159,-0.4484
702 {770,0.0256,0.50172,-0.53224,-0.46666,
703 {775,0.10169,0.49939,-0.54158,-0.48485
704 {780,0.17164,0.49748,-0.54958,-0.503,-
705 {785,0.23528,0.49609,-0.55635,-0.52109
706 {790,0.29246,0.49532,-0.56201,-0.53909
707 {795,0.34305,0.49529,-0.56666,-0.55697
708 {800,0.38717,0.49602,-0.57044,-0.5747,
709 {805,0.42531,0.49745,-0.57352,-0.59223
710 {810,0.45797,0.49951,-0.57609,-0.60951
711 {815,0.48565,0.50213,-0.57833,-0.6265,
712 {820,0.50886,0.50523,-0.58041,-0.64314
713 {825,0.5281,0.50875,-0.58253,-0.65939,
714 {830,0.54388,0.51262,-0.58485,-0.6752,
715 {835,0.55669,0.51675,-0.58757,-0.69052
716 {840,0.56705,0.5211,-0.59085,-0.70531,
717 {845,0.57545,0.52557,-0.59488,-0.71951
718 {850,0.5824,0.53011,-0.59985,-0.73309,
719 {855,0.5884,0.53464,-0.60592,-0.74598,
720 {860,0.59395,0.53909,-0.61329,-0.75815
721 {865,0.59954,0.54338,-0.62212,-0.76955
722 {870,0.60541,0.54745,-0.63256,-0.78013
723 {875,0.61159,0.55119,-0.64468,-0.78984
724 {880,0.61811,0.55452,-0.65859,-0.79866
725 {885,0.62501,0.55734,-0.67437,-0.80652
726 {890,0.63226,0.55962,-0.69201,-0.81339
727 {895,0.63991,0.56144,-0.71124,-0.81921
728 {900,0.6482,0.56292,-0.73165,-0.82391,
729 {905,0.65734,0.56418,-0.7528,-0.82746,
730 {910,0.66707,0.56537,-0.77433,-0.82982
731 {915,0.67677,0.56667,-0.79588,-0.83094
732 {920,0.68581,0.56824,-0.81712,-0.83079
733 {925,0.6938,0.57027,-0.83765,-0.82936,
734 {930,0.70068,0.57295,-0.85705,-0.8266,
735 {935,0.70637,0.57649,-0.8749,-0.8225,-
736 {940,0.71082,0.58109,-0.89077,-0.81703
737 {945,0.71395,0.58698,-0.90423,-0.81016
738 {950,0.71565,0.59441,-0.91489,-0.80182
739 {955,0.71583,0.60361,-0.92234,-0.79199
740 {960,0.71439,0.61483,-0.92615,-0.78061
741 {965,0.71134,0.62829,-0.92599,-0.76765
742 {970,0.70709,0.64411,-0.92166,-0.7531,
743 {975,0.7021,0.66238,-0.91307,-0.73694,
744 {980,0.69672,0.68301,-0.90027,-0.7192,
745 {985,0.69117,0.70589,-0.8834,-0.69991,
746 {990,0.68572,0.73088,-0.8626,-0.6791,-
747 {995,0.6806,0.75785,-0.838,-0.6568,0.0
748 {1000,0.67606,0.78666,-0.80974,-0.6330
749 {1005,0.67236,0.81718,-0.77798,-0.6078
750 {1010,0.66984,0.84916,-0.74299,-0.5813
751 {1015,0.66874,0.88239,-0.70509,-0.5534
752 {1020,0.66925,0.91668,-0.66453,-0.5242
753 {1025,0.67157,0.95184,-0.62161,-0.4938
754 {1030,0.67589,0.98766,-0.57661,-0.4622
755 {1035,0.68241,1.02396,-0.52979,-0.4295
756 {1040,0.69132,1.06053,-0.48143,-0.3958
757 {1045,0.70281,1.09718,-0.43183,-0.3610
758 {1050,0.71699,1.13373,-0.38124,-0.3254
759 {1055,0.73356,1.17007,-0.32997,-0.2888
760 {1060,0.75214,1.20612,-0.2783,-0.25144
761 {1065,0.77234,1.24178,-0.22653,-0.2133
762 {1070,0.79378,1.27697,-0.17493,-0.1744
763 {1075,0.81607,1.31159,-0.1238,-0.13504
764 {1080,0.83883,1.34556,-0.07343,-0.0950
765 {1085,0.86167,1.37879,-0.02411,-0.0545
766 {1090,0.88422,1.41119,0.02389,-0.01366
767 {1095,0.90616,1.44276,0.07039,0.02762,
768 {1100,0.92718,1.47349,0.11521,0.06924,
769 {1105,0.94696,1.50338,0.15818,0.11115,
770 {1110,0.9652,1.53241,0.19913,0.15331,0
771 {1115,0.9816,1.5606,0.23789,0.19569,0.
772 {1120,0.99584,1.58794,0.27429,0.23824,
773 {1125,1.00761,1.61442,0.30816,0.28092,
774 {1130,1.01668,1.64004,0.33938,0.32369,
775 {1135,1.02313,1.66484,0.36806,0.36651,
776 {1140,1.02709,1.68883,0.39433,0.40934,
777 {1145,1.02872,1.71205,0.41837,0.45217,
778 {1150,1.02817,1.73453,0.4403,0.49494,0
779 {1155,1.02559,1.75629,0.4603,0.53762,0
780 {1160,1.02111,1.77737,0.47851,0.58018,
781 {1165,1.01489,1.7978,0.49508,0.62259,0
782 {1170,1.00709,1.81759,0.51017,0.66481,
783 {1175,0.99786,1.83673,0.52401,0.70683,
784 {1180,0.98737,1.8552,0.53683,0.74862,0
785 {1185,0.9758,1.87296,0.54887,0.79018,0
786 {1190,0.96331,1.89,0.56034,0.8315,0.74
787 {1195,0.95008,1.90628,0.57149,0.87255,
788 {1200,0.93626,1.92179,0.58256,0.91332,
789 {1205,0.92203,1.9365,0.59376,0.95379,0
790 {1210,0.90757,1.95038,0.60531,0.99396,
791 {1215,0.89304,1.96341,0.61726,1.03379,
792 {1220,0.87864,1.97556,0.62967,1.07325,
793 {1225,0.86453,1.9868,0.64256,1.11232,0
794 {1230,0.85091,1.9971,0.65596,1.15095,0
795 {1235,0.83795,2.00644,0.66993,1.18913,
796 {1240,0.82582,2.01479,0.68448,1.22682,
797 {1245,0.81473,2.02212,0.69966,1.26398,
798 {1250,0.80486,2.0284,0.71551,1.3006,0.
799 {1255,0.7966,2.03353,0.73205,1.33662,0
800 {1260,0.79032,2.03741,0.74932,1.372,1.
801 {1265,0.78642,2.03993,0.76735,1.40669,
802 {1270,0.78492,2.04109,0.78615,1.44068,
803 {1275,0.78537,2.04098,0.80574,1.47397,
804 {1280,0.78729,2.03971,0.8261,1.50657,1
805 {1285,0.79021,2.03737,0.84724,1.53849,
806 {1290,0.79375,2.03407,0.86916,1.56972,
807 {1295,0.79786,2.02989,0.89186,1.60028,
808 {1300,0.80263,2.02488,0.91536,1.63014,
809 {1305,0.80811,2.01914,0.93965,1.65929,
810 {1310,0.81435,2.01272,0.96476,1.68774,
811 {1315,0.82144,2.0057,0.99069,1.71547,1
812 {1320,0.82936,1.99815,1.01744,1.74248,
813 {1325,0.83779,1.99015,1.04496,1.76877,
814 {1330,0.84665,1.9818,1.07325,1.79437,1
815 {1335,0.85598,1.97322,1.10228,1.81929,
816 {1340,0.86584,1.96452,1.13207,1.84353,
817 {1345,0.87626,1.9558,1.16259,1.86713,1
818 {1350,0.88729,1.94717,1.19385,1.89008,
819 {1355,0.89897,1.93876,1.22584,1.9124,1
820 {1360,0.91133,1.93063,1.25852,1.9341,1
821 {1365,0.92428,1.92278,1.29175,1.95515,
822 {1370,0.93773,1.91515,1.32538,1.97551,
823 {1375,0.95159,1.90768,1.35924,1.99516,
824 {1380,0.96583,1.90033,1.39318,2.01405,
825 {1385,0.98042,1.89304,1.42703,2.03215,
826 {1390,0.9953,1.88573,1.46066,2.04944,1
827 {1395,1.01046,1.87835,1.49391,2.06587,
828 {1400,1.02585,1.87084,1.52662,2.08142,
829 {1405,1.04143,1.86314,1.55865,2.09605,
830 {1410,1.05718,1.8552,1.58983,2.10973,1
831 {1415,1.07304,1.84694,1.62001,2.12243,
832 {1420,1.08898,1.83832,1.64908,2.13412,
833 {1425,1.10484,1.82934,1.67699,2.14484,
834 {1430,1.12048,1.81999,1.70374,2.15463,
835 {1435,1.13574,1.81027,1.72932,2.16352,
836 {1440,1.15046,1.80018,1.75373,2.17157,
837 {1445,1.16449,1.78973,1.77696,2.17881,
838 {1450,1.17767,1.77892,1.799,2.18528,1.
839 {1455,1.18985,1.76774,1.81985,2.19103,
840 {1460,1.20086,1.75621,1.8395,2.1961,1.
841 {1465,1.21057,1.74432,1.85795,2.20053,
842 {1470,1.219,1.7322,1.87524,2.20436,1.6
843 {1475,1.22629,1.72006,1.89148,2.20761,
844 {1480,1.2326,1.7081,1.90674,2.21031,1.
845 {1485,1.23807,1.69653,1.92113,2.2125,1
846 {1490,1.24286,1.68556,1.93474,2.21421,
847 {1495,1.24712,1.67541,1.94765,2.21547,
848 {1500,1.251,1.66628,1.95997,2.21629,1.
849 {1505,1.25464,1.65837,1.97178,2.21673,
850 {1510,1.25822,1.65191,1.98318,2.2168,1
851 {1515,1.26186,1.64709,1.99426,2.21654,
852 {1520,1.26574,1.64404,2.00509,2.21596,
853 {1525,1.27001,1.6427,2.01568,2.21509,1
854 {1530,1.27484,1.64297,2.02605,2.21392,
855 {1535,1.28041,1.64477,2.03621,2.21246,
856 {1540,1.28686,1.64799,2.04616,2.21072,
857 {1545,1.29439,1.65256,2.05591,2.2087,1
858 {1550,1.3031,1.65836,2.06547,2.20642,1
859 {1555,1.3129,1.66525,2.07475,2.20382,1
860 {1560,1.32356,1.67306,2.08365,2.20085,
861 {1565,1.33485,1.68162,2.09205,2.19746,
862 {1570,1.34654,1.69074,2.09985,2.19357,
863 {1575,1.35842,1.70026,2.10691,2.18914,
864 {1580,1.37025,1.71001,2.11314,2.1841,1
865 {1585,1.38181,1.71981,2.11843,2.17841,
866 {1590,1.39288,1.72949,2.12265,2.17199,
867 {1595,1.40323,1.73887,2.1257,2.16479,1
868 {1600,1.41263,1.74778,2.12746,2.15675,
869 {1605,1.42086,1.75605,2.12782,2.14782,
870 {1610,1.42772,1.76352,2.12669,2.13793,
871 {1615,1.43328,1.77017,2.12406,2.12708,
872 {1620,1.43768,1.77599,2.11996,2.11527,
873 {1625,1.4411,1.78103,2.11445,2.10249,1
874 {1630,1.44369,1.7853,2.10756,2.08876,1
875 {1635,1.44561,1.78882,2.09933,2.07406,
876 {1640,1.44703,1.79161,2.0898,2.05839,1
877 {1645,1.4481,1.7937,2.07902,2.04178,1.
878 {1650,1.449,1.79509,2.06703,2.0242,1.8
879 {1655,1.44987,1.79582,2.05386,2.00566,
880 {1660,1.4508,1.79591,2.03956,1.98617,1
881 {1665,1.45177,1.79536,2.02422,1.96574,
882 {1670,1.45273,1.79421,2.0079,1.94436,1
883 {1675,1.45366,1.79245,1.99068,1.92205,
884 {1680,1.45453,1.79012,1.97262,1.89882,
885 {1685,1.4553,1.78723,1.95379,1.87466,1
886 {1690,1.45594,1.78379,1.93428,1.84959,
887 {1695,1.45642,1.77982,1.91414,1.82361,
888 {1700,1.4567,1.77534,1.89345,1.79673,1
889 {1705,1.45676,1.77036,1.87229,1.76895,
890 {1710,1.45653,1.76489,1.85069,1.74029,
891 {1715,1.45585,1.75888,1.8286,1.71072,1
892 {1720,1.45454,1.75226,1.80593,1.68023,
893 {1725,1.45241,1.74498,1.78259,1.6488,1
894 {1730,1.44926,1.73696,1.75849,1.61643,
895 {1735,1.44492,1.72815,1.73354,1.58309,
896 {1740,1.43919,1.71847,1.70766,1.54877,
897 {1745,1.43192,1.70788,1.68076,1.51346,
898 {1750,1.42321,1.69643,1.6529,1.4772,1.
899 {1755,1.41325,1.68419,1.62417,1.44003,
900 {1760,1.40226,1.67124,1.59466,1.40201,
901 {1765,1.39045,1.65768,1.56447,1.3632,1
902 {1770,1.37802,1.64358,1.53368,1.32364,
903 {1775,1.36518,1.62903,1.5024,1.28339,1
904 {1780,1.35214,1.61411,1.47071,1.24251,
905 {1785,1.3391,1.59891,1.43872,1.20104,1
906 {1790,1.32627,1.58351,1.4065,1.15905,1
907 {1795,1.31385,1.56799,1.37416,1.11658,
908 {1800,1.30206,1.55243,1.3418,1.07368,1
909 {1805,1.29111,1.53692,1.3095,1.03042,1
910 {1810,1.28119,1.52155,1.27735,0.98684,
911 {1815,1.27252,1.50639,1.24546,0.943,1.
912 {1820,1.2653,1.49154,1.21391,0.89895,1
913 {1825,1.25974,1.47706,1.18281,0.85475,
914 {1830,1.25605,1.46305,1.15224,0.81044,
915 {1835,1.25443,1.4496,1.12229,0.76609,1
916 {1840,1.25509,1.43678,1.09306,0.72174,
917 {1845,1.25803,1.42461,1.06458,0.67744,
918
919
920
921 const G4int coef_ener = G4int((pLab-Lege
922 const G4double sup_ener = pLab/5. - coef
923
924 // assert(pLab >= Legendre_coef[WhichChannel+1
925
926 // Legendre coefficient normalized
927 const G4double A0 = 1.;
928 const G4double A1 = (1-sup_ener)*Legendr
929 const G4double A2 = (1-sup_ener)*Legendr
930 const G4double A3 = (1-sup_ener)*Legendr
931 const G4double A4 = (1-sup_ener)*Legendr
932 const G4double A5 = (1-sup_ener)*Legendr
933 const G4double A6 = (1-sup_ener)*Legendr
934 const G4double A7 = (1-sup_ener)*Legendr
935 const G4double A8 = (1-sup_ener)*Legendr
936
937 // Theoritical max if all Ai > 0 (often
938 const G4double A = std::fabs(A0) + std::
939
940 G4bool success = false;
941 G4int maxloop = 0;
942
943 while(!success && maxloop < 1000){
944
945 cos_theta = Random::shoot()*2-1.; // n
946
947 // Legendre Polynomial
948 G4double P0 = A0;
949 G4double P1 = A1*cos_theta;
950 G4double P2 = A2/2.*(3*std::pow(cos_th
951 G4double P3 = A3/2.*(5*std::pow(cos_th
952 G4double P4 = A4/8.*(35*std::pow(cos_t
953 G4double P5 = A5/8.*(63*std::pow(cos_t
954 G4double P6 = A6/16.*(231*std::pow(cos
955 G4double P7 = A7/16.*(429*std::pow(cos
956 G4double P8 = A8/128.*(6435*std::pow(c
957
958 G4double P = (P0 + P1 + P2 + P3 + P4 +
959
960 if(Random::shoot()*A < P) success = tr
961 maxloop +=1 ;
962 if(maxloop==1000) cos_theta = std::log
963 }
964 sin_theta = std::sqrt(1-cos_theta*cos_th
965 }
966
967 if(rho == 0) return ThreeVector(sin_theta*
968 // Rotation in the direction of the incide
969 const G4double px = x/r*cos_theta - y/rho*
970 const G4double py = y/r*cos_theta + x/rho*
971 const G4double pz = z/r*cos_theta - rho/r*
972
973 return ThreeVector(px,py,pz);
974 }
975
976 }
977