| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
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 Helsinki Institute of Physics, Finland
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 *p1, Particle *p2)
50 : particle1(p1), particle2(p2)
51 {}
52
53 NKbToSpiChannel::~NKbToSpiChannel(){}
54
55 void NKbToSpiChannel::fillFinalState(FinalState *fs) {
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::getIsospin(nucleon->getType()) + ParticleTable::getIsospin(kaon->getType());
79 const G4int iso_n = ParticleTable::getIsospin(nucleon->getType());
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(iso);
87 SigmaType = SigmaZero;
88 WhichChannel = 0;
89 }
90 else{
91 SigmaType = ParticleTable::getSigmaType(iso);
92 PionType = PiZero;
93 WhichChannel = iso/2;
94 }
95 }
96 else{
97 if(rdm*13. < 6.){
98 PionType = ParticleTable::getPionType(-2.*iso_n);
99 SigmaType = ParticleTable::getSigmaType(2.*iso_n);
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(2.*iso_n);
109 SigmaType = ParticleTable::getSigmaType(-2.*iso_n);
110 WhichChannel = -iso_n;
111 }
112 }
113
114 ThreeVector mom_pion = KaonMomentum(kaon,nucleon,WhichChannel); // based on K- p collision
115
116 kaon->setType(PionType);
117 nucleon->setType(SigmaType);
118
119 G4double norm = KinematicsUtils::momentumInCM(kaon,nucleon);
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(Particle const * const kaon, Particle const * const nucleon, G4int WhichChannel){
133
134 const G4double pLab = KinematicsUtils::momentumInLab(kaon,nucleon);
135 G4double max = 0;
136
137 if(pLab < 435.) return Random::normVector(); // isotropic
138
139 if(WhichChannel == 0) max = 865;
140 else if(std::abs(WhichChannel) == 1) max = 1845;
141 else{
142 INCL_WARN("Bad channel chosen in Kb N To Sigma pi angular distribution determination. Return an isotropic distribution.\n");
143 return Random::normVector();
144 }
145
146 G4double cos_theta = 1.;
147 G4double sin_theta = 0.;
148 const G4double cos_phi = std::cos(Random::shoot()*Math::twoPi);
149 const G4double sin_phi = std::sqrt(1-cos_phi*cos_phi);
150
151 const G4double x = kaon->getMomentum().getX();
152 const G4double y = kaon->getMomentum().getY();
153 const G4double z = kaon->getMomentum().getZ();
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.; // correspond to the forward slope description at 2375 MeV/c
160 cos_theta = std::log(Random::shoot()*(std::exp(b)-std::exp(-b))+std::exp(-b))/b;
161 sin_theta = std::sqrt(1-cos_theta*cos_theta);
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,0.00013,1e-05,0,0},
168 {440,0.2116,0.91169,0.00737,-0.02297,2e-04,0,0,0},
169 {445,0.21976,0.92725,0.01151,-0.02554,0.00027,-1e-05,0,0},
170 {450,0.22691,0.94264,0.01546,-0.02817,0.00033,-2e-05,0,0},
171 {455,0.23293,0.95777,0.0191,-0.03089,4e-04,-3e-05,0,0},
172 {460,0.23762,0.97259,0.02233,-0.03372,0.00046,-4e-05,0,0},
173 {465,0.24084,0.98722,0.025,-0.03662,0.00053,-5e-05,0,0},
174 {470,0.24244,1.00183,0.02695,-0.0395,0.00059,-6e-05,0,0},
175 {475,0.24243,1.0166,0.02804,-0.0423,0.00065,-7e-05,0,0},
176 {480,0.2408,1.03171,0.0281,-0.04494,0.00071,-8e-05,0,0},
177 {485,0.23759,1.04732,0.02704,-0.04728,0.00076,-8e-05,0,0},
178 {490,0.23289,1.06364,0.02475,-0.04917,0.00081,-9e-05,0,0},
179 {495,0.22693,1.08085,0.02112,-0.05045,0.00086,-1e-04,0,0},
180 {500,0.21985,1.0991,0.01607,-0.05099,0.00089,-0.00011,0,0},
181 {505,0.21186,1.11838,0.00962,-0.05067,0.00092,-0.00012,0,0},
182 {510,0.20311,1.13865,0.0018,-0.04938,0.00094,-0.00012,0,0},
183 {515,0.19378,1.15987,-0.00734,-0.04703,0.00095,-0.00013,0,0},
184 {520,0.18394,1.18199,-0.01776,-0.04353,0.00094,-0.00013,0,0},
185 {525,0.17362,1.20497,-0.02941,-0.03891,0.00092,-0.00014,0,0},
186 {530,0.16277,1.22874,-0.04223,-0.03321,0.00088,-0.00014,1e-05,0},
187 {535,0.15135,1.25326,-0.05615,-0.02646,0.00082,-0.00014,1e-05,0},
188 {540,0.13922,1.2785,-0.07107,-0.01873,0.00073,-0.00014,1e-05,0},
189 {545,0.12624,1.30447,-0.08684,-0.01017,0.00061,-0.00013,1e-05,0},
190 {550,0.1123,1.3312,-0.10326,-0.00094,0.00047,-0.00012,1e-05,0},
191 {555,0.09727,1.35871,-0.12015,0.00879,0.00029,-0.00011,1e-05,0},
192 {560,0.08117,1.38699,-0.13734,0.01891,7e-05,-9e-05,1e-05,0},
193 {565,0.064,1.41599,-0.15466,0.02935,-0.00019,-7e-05,2e-05,0},
194 {570,0.0459,1.44562,-0.17197,0.04009,-5e-04,-5e-05,2e-05,0},
195 {575,0.02704,1.47583,-0.18911,0.05108,-0.00085,-2e-05,2e-05,0},
196 {580,0.00769,1.50654,-0.20595,0.06228,-0.00126,2e-05,2e-05,0},
197 {585,-0.01183,1.53766,-0.22246,0.07364,-0.00172,6e-05,2e-05,0},
198 {590,-0.03118,1.56911,-0.23858,0.08511,-0.00224,0.00011,2e-05,0},
199 {595,-0.05002,1.60083,-0.25425,0.09665,-0.00282,0.00016,3e-05,0},
200 {600,-0.06802,1.63271,-0.26943,0.10821,-0.00346,0.00023,3e-05,0},
201 {605,-0.08492,1.66463,-0.28414,0.1197,-0.00417,3e-04,3e-05,0},
202 {610,-0.10045,1.69641,-0.29846,0.131,-0.00494,0.00038,3e-05,0},
203 {615,-0.11445,1.72786,-0.31248,0.14199,-0.00578,0.00047,3e-05,0},
204 {620,-0.12685,1.75882,-0.32626,0.15255,-0.00669,0.00057,3e-05,0},
205 {625,-0.13761,1.78907,-0.3398,0.16263,-0.00766,0.00068,3e-05,0},
206 {630,-0.14665,1.81843,-0.353,0.17219,-0.00869,8e-04,3e-05,0},
207 {635,-0.15404,1.84667,-0.36578,0.18122,-0.00977,0.00093,3e-05,0},
208 {640,-0.15984,1.8736,-0.37806,0.18969,-0.0109,0.00107,3e-05,0},
209 {645,-0.1641,1.89903,-0.38974,0.19753,-0.01207,0.00121,3e-05,0},
210 {650,-0.1668,1.92283,-0.40077,0.20457,-0.01326,0.00137,3e-05,0},
211 {655,-0.16802,1.94484,-0.41108,0.2107,-0.01447,0.00153,3e-05,0},
212 {660,-0.16778,1.96491,-0.4206,0.21575,-0.01567,0.0017,2e-05,0},
213 {665,-0.16609,1.98295,-0.42932,0.21971,-0.01684,0.00187,2e-05,1e-05},
214 {670,-0.16293,1.99895,-0.4373,0.22266,-0.01797,0.00204,1e-05,1e-05},
215 {675,-0.15835,2.01289,-0.4446,0.2247,-0.01902,0.00221,0,1e-05},
216 {680,-0.15235,2.02475,-0.45128,0.22593,-0.01997,0.00238,-1e-05,1e-05},
217 {685,-0.14502,2.03452,-0.45737,0.22638,-0.02077,0.00255,-2e-05,1e-05},
218 {690,-0.1364,2.04219,-0.46282,0.22601,-0.02139,0.0027,-3e-05,1e-05},
219 {695,-0.12665,2.04773,-0.46761,0.22477,-0.02178,0.00284,-4e-05,1e-05},
220 {700,-0.11594,2.05112,-0.47172,0.22263,-0.0219,0.00296,-6e-05,1e-05},
221 {705,-0.10451,2.05233,-0.47513,0.21958,-0.0217,0.00305,-8e-05,1e-05},
222 {710,-0.09257,2.05134,-0.47789,0.21563,-0.02111,0.00311,-1e-04,1e-05},
223 {715,-0.08045,2.0481,-0.48004,0.21082,-0.02008,0.00314,-0.00012,1e-05},
224 {720,-0.06841,2.04258,-0.48164,0.20519,-0.01853,0.00312,-0.00015,1e-05},
225 {725,-0.05677,2.03476,-0.4827,0.19876,-0.01641,0.00305,-0.00018,1e-05},
226 {730,-0.04577,2.02468,-0.48321,0.19159,-0.01364,0.00292,-0.00021,1e-05},
227 {735,-0.03567,2.01241,-0.48312,0.18371,-0.01013,0.00271,-0.00024,1e-05},
228 {740,-0.02666,1.99799,-0.4824,0.17518,-0.0058,0.00243,-0.00028,1e-05},
229 {745,-0.01887,1.98148,-0.48101,0.16603,-0.00058,0.00205,-0.00031,1e-05},
230 {750,-0.01239,1.96293,-0.47891,0.15626,0.00561,0.00157,-0.00035,1e-05},
231 {755,-0.00733,1.9424,-0.47608,0.14587,0.01287,0.00098,-0.00039,1e-05},
232 {760,-0.00369,1.91993,-0.47246,0.13483,0.02128,0.00026,-0.00043,0},
233 {765,-0.0015,1.8956,-0.46801,0.12316,0.03091,-0.00059,-0.00048,0},
234 {770,-0.00083,1.86947,-0.46257,0.1109,0.04184,-0.0016,-0.00052,0},
235 {775,-0.00169,1.84167,-0.45593,0.09812,0.05414,-0.00277,-0.00056,0},
236 {780,-0.0042,1.81233,-0.44791,0.0849,0.06785,-0.0041,-0.00061,-1e-05},
237 {785,-0.0085,1.78162,-0.43841,0.07136,0.08299,-0.00562,-0.00065,-1e-05},
238 {790,-0.0149,1.74974,-0.4273,0.05758,0.09958,-0.00733,-0.00068,-2e-05},
239 {795,-0.02359,1.71687,-0.41448,0.04365,0.11765,-0.00925,-0.00072,-2e-05},
240 {800,-0.03492,1.68318,-0.3999,0.02967,0.13718,-0.01137,-0.00074,-3e-05},
241 {805,-0.04918,1.64882,-0.38363,0.01571,0.15809,-0.0137,-0.00076,-4e-05},
242 {810,-0.06654,1.61392,-0.36573,0.00182,0.18032,-0.01625,-0.00077,-4e-05},
243 {815,-0.08712,1.57864,-0.34628,-0.01193,0.20379,-0.01901,-0.00077,-5e-05},
244 {820,-0.11012,1.54311,-0.32533,-0.02547,0.22841,-0.02199,-0.00075,-6e-05},
245 {825,-0.13493,1.50749,-0.30295,-0.03874,0.25413,-0.02519,-0.00072,-7e-05},
246 {830,-0.16027,1.47191,-0.2792,-0.05167,0.28087,-0.02861,-0.00068,-8e-05},
247 {835,-0.18509,1.43652,-0.25416,-0.0642,0.30854,-0.03226,-0.00061,-9e-05},
248 {840,-0.20852,1.40147,-0.22789,-0.07626,0.33709,-0.03613,-0.00053,-1e-04},
249 {845,-0.23013,1.3669,-0.20046,-0.08779,0.36642,-0.04024,-0.00042,-0.00012},
250 {850,-0.2499,1.33295,-0.17192,-0.09871,0.39647,-0.04458,-3e-04,-0.00013},
251 {855,-0.26818,1.29977,-0.14236,-0.10898,0.42717,-0.04916,-0.00014,-0.00015},
252 {860,-0.2854,1.2675,-0.11182,-0.11851,0.45844,-0.05397,4e-05,-0.00016},
253 {865,-0.30205,1.23628,-0.08039,-0.12725,0.4902,-0.05903,0.00025,-0.00018},
254 {870,-0.31869,1.20621,-0.0481,-0.13519,0.52231,-0.06438,0.00049,-2e-04},
255 {875,-0.33574,1.17732,-0.01498,-0.14236,0.5546,-0.07012,0.00076,-0.00021},
256 {880,-0.35356,1.14966,0.01893,-0.14877,0.58687,-0.07634,0.00107,-0.00023},
257 {885,-0.37245,1.12329,0.05362,-0.15447,0.61893,-0.08312,0.00141,-0.00025},
258 {890,-0.3926,1.09818,0.08898,-0.15955,0.65061,-0.09049,0.00179,-0.00027},
259 {895,-0.41404,1.07429,0.12481,-0.16417,0.68173,-0.09828,0.00221,-0.00029},
260 {900,-0.43679,1.05155,0.16084,-0.16846,0.7122,-0.10616,0.00266,-0.00031},
261 {905,-0.46071,1.02994,0.1968,-0.17253,0.74193,-0.11383,0.00314,-0.00033},
262 {910,-0.4856,1.00944,0.23246,-0.17643,0.77089,-0.12104,0.00366,-0.00034},
263 {915,-0.51107,0.99007,0.26757,-0.18017,0.79911,-0.12764,0.00421,-0.00035},
264 {920,-0.53674,0.97186,0.3019,-0.18373,0.8266,-0.13347,0.00478,-0.00036},
265 {925,-0.56217,0.95486,0.33525,-0.18707,0.85342,-0.13846,0.00537,-0.00036},
266 {930,-0.58681,0.9392,0.36744,-0.19012,0.87966,-0.14268,0.00598,-0.00036},
267 {935,-0.61023,0.92501,0.39833,-0.1928,0.9054,-0.14621,0.00659,-0.00035},
268 {940,-0.63197,0.91241,0.42775,-0.19502,0.93074,-0.14913,0.0072,-0.00034},
269 {945,-0.65175,0.90149,0.45562,-0.19679,0.95586,-0.15152,0.0078,-0.00031},
270 {950,-0.66926,0.89231,0.4819,-0.19816,0.98094,-0.15347,0.00838,-0.00027},
271 {955,-0.68439,0.88494,0.50655,-0.19918,1.00619,-0.15508,0.00891,-0.00022},
272 {960,-0.69707,0.87945,0.52952,-0.19989,1.03181,-0.15642,0.00939,-0.00016},
273 {965,-0.70731,0.87589,0.55076,-0.20033,1.05793,-0.15758,0.0098,-8e-05},
274 {970,-0.71516,0.87419,0.5703,-0.20045,1.08454,-0.15852,0.01012,2e-05},
275 {975,-0.7207,0.87429,0.58813,-0.20021,1.11152,-0.1592,0.01032,0.00013},
276 {980,-0.72394,0.87607,0.60419,-0.1995,1.13868,-0.1595,0.01035,0.00027},
277 {985,-0.72503,0.87939,0.61843,-0.19819,1.16573,-0.15928,0.01018,0.00043},
278 {990,-0.72394,0.88412,0.63075,-0.19616,1.1924,-0.15842,0.00978,0.00061},
279 {995,-0.72067,0.89012,0.64109,-0.19327,1.21844,-0.15676,0.0091,0.00081},
280 {1000,-0.71516,0.89727,0.64936,-0.1894,1.24357,-0.15418,0.0081,0.00105},
281 {1005,-0.70731,0.90541,0.6555,-0.18445,1.26753,-0.15055,0.00676,0.00131},
282 {1010,-0.697,0.91431,0.65947,-0.17846,1.29013,-0.1459,0.00501,0.00161},
283 {1015,-0.68411,0.92376,0.66132,-0.17145,1.31127,-0.14023,0.00282,0.00193},
284 {1020,-0.66835,0.93356,0.66107,-0.16347,1.3309,-0.13354,0.00014,0.00228},
285 {1025,-0.64976,0.94348,0.65877,-0.15451,1.34895,-0.12582,-0.0031,0.00266},
286 {1030,-0.62823,0.95334,0.65446,-0.14459,1.36535,-0.11708,-0.00694,0.00307},
287 {1035,-0.6038,0.96292,0.64817,-0.13375,1.38003,-0.10731,-0.01143,0.00351},
288 {1040,-0.57654,0.97201,0.63996,-0.12199,1.39292,-0.0965,-0.01662,0.00397},
289 {1045,-0.54672,0.9804,0.62985,-0.10933,1.40396,-0.08465,-0.02256,0.00446},
290 {1050,-0.51475,0.98795,0.61792,-0.09577,1.41311,-0.07177,-0.02927,0.00498},
291 {1055,-0.48114,0.99462,0.60434,-0.08128,1.42048,-0.05788,-0.03657,0.00552},
292 {1060,-0.44627,1.00047,0.58928,-0.06579,1.42621,-0.04301,-0.04426,0.00608},
293 {1065,-0.41097,1.00553,0.57295,-0.04925,1.43042,-0.02719,-0.05212,0.00665},
294 {1070,-0.37567,1.00983,0.55554,-0.03158,1.43327,-0.01044,-0.05994,0.00724},
295 {1075,-0.34091,1.01341,0.53725,-0.01274,1.43488,0.00719,-0.06751,0.00782},
296 {1080,-0.30701,1.0163,0.51825,0.00733,1.4354,0.02567,-0.07463,0.00841},
297 {1085,-0.27413,1.01854,0.49874,0.0287,1.43496,0.04499,-0.08107,0.009},
298 {1090,-0.2425,1.02016,0.47892,0.0514,1.4337,0.06511,-0.08666,0.00957},
299 {1095,-0.21201,1.02123,0.45893,0.0754,1.43171,0.08607,-0.09129,0.01013},
300 {1100,-0.18239,1.02178,0.43894,0.10065,1.4291,0.10791,-0.09489,0.01066},
301 {1105,-0.15367,1.02189,0.41909,0.12708,1.42598,0.13065,-0.09738,0.01114},
302 {1110,-0.1254,1.02159,0.39954,0.15464,1.42243,0.15435,-0.09869,0.01156},
303 {1115,-0.09757,1.02096,0.38043,0.18328,1.41856,0.17903,-0.09875,0.01192},
304 {1120,-0.06982,1.02003,0.36193,0.21294,1.41446,0.20473,-0.09748,0.0122},
305 {1125,-0.04213,1.01887,0.34419,0.24357,1.41025,0.23151,-0.0948,0.01238},
306 {1130,-0.0142,1.01753,0.32732,0.27509,1.40599,0.25935,-0.09067,0.01246},
307 {1135,0.01406,1.01605,0.31137,0.30741,1.40167,0.28815,-0.08514,0.01242},
308 {1140,0.04282,1.01449,0.29633,0.3404,1.39724,0.31776,-0.07828,0.01222},
309 {1145,0.0719,1.0129,0.28221,0.37395,1.39267,0.34804,-0.07016,0.01185},
310 {1150,0.10146,1.01132,0.26901,0.40794,1.3879,0.37885,-0.06087,0.01128},
311 {1155,0.13127,1.00979,0.25673,0.44225,1.38291,0.41003,-0.05049,0.0105},
312 {1160,0.1612,1.00837,0.24537,0.47676,1.37764,0.44145,-0.03908,0.00947},
313 {1165,0.19131,1.00711,0.23494,0.51136,1.37206,0.47296,-0.02672,0.00818},
314 {1170,0.22156,1.00605,0.22546,0.54595,1.36612,0.5044,-0.01351,0.0066},
315 {1175,0.25183,1.00528,0.21695,0.58048,1.35981,0.53559,0.00041,0.00471},
316 {1180,0.28191,1.00489,0.20947,0.61496,1.35311,0.56633,0.01488,0.00247},
317 {1185,0.31173,1.00498,0.20307,0.64937,1.34598,0.59642,0.02974,-0.00015},
318 {1190,0.34118,1.00562,0.1978,0.68368,1.33842,0.62566,0.04483,-0.00318},
319 {1195,0.37021,1.0069,0.1937,0.71789,1.33041,0.65385,0.05998,-0.00664},
320 {1200,0.39843,1.00892,0.19083,0.75199,1.32192,0.6808,0.07505,-0.01057},
321 {1205,0.42584,1.01177,0.18923,0.78595,1.31293,0.7063,0.08986,-0.015},
322 {1210,0.45231,1.01553,0.18896,0.81978,1.30344,0.7302,0.10428,-0.01995},
323 {1215,0.47765,1.02033,0.19014,0.85355,1.29349,0.7525,0.11823,-0.02545},
324 {1220,0.50177,1.02629,0.19287,0.88735,1.28312,0.77325,0.13164,-0.03151},
325 {1225,0.52457,1.03352,0.19729,0.92127,1.27241,0.79248,0.14447,-0.03814},
326 {1230,0.54599,1.04214,0.2035,0.95541,1.26139,0.81024,0.15665,-0.04537},
327 {1235,0.56609,1.05228,0.21163,0.98984,1.25012,0.82657,0.16812,-0.05321},
328 {1240,0.5847,1.06406,0.2218,1.02465,1.23865,0.84152,0.17883,-0.06168},
329 {1245,0.60195,1.07759,0.23413,1.05995,1.22704,0.85513,0.18872,-0.07079},
330 {1250,0.61786,1.09296,0.24867,1.09575,1.21531,0.86741,0.19769,-0.08061},
331 {1255,0.63245,1.11008,0.26523,1.13191,1.20336,0.87829,0.20553,-0.09138},
332 {1260,0.64572,1.1288,0.28356,1.16821,1.19107,0.88763,0.21202,-0.10339},
333 {1265,0.65767,1.14902,0.30341,1.20443,1.1783,0.89534,0.2169,-0.11694},
334 {1270,0.6683,1.17065,0.32457,1.24045,1.16511,0.90145,0.22018,-0.13214},
335 {1275,0.67758,1.19375,0.34687,1.27626,1.15175,0.90623,0.22212,-0.14893},
336 {1280,0.68536,1.21833,0.37014,1.31188,1.1385,0.90994,0.22304,-0.16725},
337 {1285,0.69167,1.24444,0.39422,1.34729,1.12563,0.91284,0.22321,-0.187},
338 {1290,0.69645,1.27212,0.41894,1.38249,1.11339,0.91518,0.22293,-0.20808},
339 {1295,0.69969,1.30139,0.44417,1.41745,1.10188,0.91705,0.22235,-0.23029},
340 {1300,0.70139,1.33231,0.46975,1.45211,1.09119,0.91853,0.22163,-0.25337},
341 {1305,0.70159,1.36489,0.49554,1.48643,1.0814,0.91967,0.22091,-0.27708},
342 {1310,0.70037,1.39918,0.52138,1.52037,1.07258,0.92056,0.22033,-0.30119},
343 {1315,0.69788,1.43522,0.54714,1.55387,1.06482,0.92125,0.22005,-0.32545},
344 {1320,0.69426,1.47301,0.57267,1.5869,1.05819,0.92182,0.22018,-0.34964},
345 {1325,0.68978,1.51246,0.59789,1.61947,1.05272,0.92232,0.22064,-0.37368},
346 {1330,0.68469,1.55333,0.62279,1.65156,1.04837,0.92281,0.22135,-0.39738},
347 {1335,0.6793,1.59543,0.64739,1.68314,1.04513,0.92337,0.22227,-0.4205},
348 {1340,0.67393,1.63853,0.67168,1.71418,1.04298,0.92406,0.22335,-0.44283},
349 {1345,0.66891,1.68242,0.69567,1.74465,1.04188,0.92496,0.22454,-0.46412},
350 {1350,0.66455,1.72688,0.71937,1.77453,1.04181,0.92613,0.22578,-0.48416},
351 {1355,0.66117,1.7717,0.74277,1.80377,1.04274,0.92764,0.22704,-0.50271},
352 {1360,0.65897,1.81664,0.76591,1.83233,1.04468,0.92957,0.2283,-0.51949},
353 {1365,0.65815,1.86145,0.78889,1.86002,1.0477,0.93203,0.2298,-0.53402},
354 {1370,0.65879,1.90585,0.81186,1.88663,1.05191,0.93512,0.23178,-0.54579},
355 {1375,0.66083,1.94962,0.8348,1.91201,1.05725,0.93889,0.23444,-0.55462},
356 {1380,0.66406,1.9926,0.85756,1.93612,1.06349,0.94336,0.23782,-0.56076},
357 {1385,0.66812,2.03463,0.87996,1.95888,1.07037,0.94852,0.24197,-0.56449},
358 {1390,0.67247,2.07556,0.90183,1.98025,1.07762,0.95437,0.24693,-0.5661},
359 {1395,0.67646,2.11524,0.92297,2.00017,1.08501,0.9609,0.25275,-0.56588},
360 {1400,0.67917,2.15351,0.94322,2.01858,1.09227,0.96812,0.25948,-0.56411},
361 {1405,0.67999,2.19022,0.96239,2.03543,1.09914,0.97601,0.26715,-0.56106},
362 {1410,0.67813,2.2252,0.98032,2.05066,1.10537,0.98459,0.27582,-0.55704},
363 {1415,0.67307,2.25832,0.99681,2.06423,1.11071,0.99384,0.28552,-0.55231},
364 {1420,0.66456,2.28945,1.01176,2.07609,1.11496,1.00374,0.29628,-0.54712},
365 {1425,0.65267,2.31863,1.02525,2.08636,1.11811,1.01425,0.30798,-0.54148},
366 {1430,0.6378,2.34595,1.03744,2.09517,1.12022,1.02527,0.32046,-0.53535},
367 {1435,0.62062,2.37149,1.04848,2.10265,1.12135,1.03672,0.33357,-0.52871},
368 {1440,0.60209,2.39534,1.05852,2.10894,1.12156,1.04854,0.34718,-0.52152},
369 {1445,0.5832,2.41757,1.06772,2.11416,1.1209,1.06063,0.36112,-0.51375},
370 {1450,0.56496,2.43828,1.07622,2.11846,1.11943,1.07293,0.37526,-0.50535},
371 {1455,0.54826,2.45754,1.08418,2.12197,1.1172,1.08535,0.38943,-0.49631},
372 {1460,0.53374,2.47544,1.09175,2.12481,1.11428,1.09782,0.40349,-0.48658},
373 {1465,0.52181,2.49206,1.09908,2.12713,1.11072,1.11025,0.41731,-0.47614},
374 {1470,0.51254,2.50744,1.10619,2.12895,1.10663,1.12251,0.4309,-0.46499},
375 {1475,0.50575,2.52157,1.11304,2.13027,1.10217,1.13444,0.44441,-0.4532},
376 {1480,0.50085,2.53446,1.11958,2.13104,1.09748,1.14586,0.45799,-0.44082},
377 {1485,0.49715,2.5461,1.12576,2.13126,1.09271,1.15662,0.47178,-0.4279},
378 {1490,0.49385,2.55652,1.13152,2.13089,1.08803,1.16654,0.48593,-0.4145},
379 {1495,0.49019,2.56569,1.13682,2.12991,1.08357,1.17545,0.50059,-0.40067},
380 {1500,0.48551,2.57364,1.14159,2.1283,1.07949,1.1832,0.51591,-0.38647},
381 {1505,0.47939,2.58035,1.14579,2.12603,1.07594,1.18959,0.53203,-0.37195},
382 {1510,0.47148,2.58582,1.14937,2.12307,1.07308,1.19448,0.5491,-0.35716},
383 {1515,0.46175,2.59007,1.15228,2.11941,1.07105,1.1977,0.56728,-0.34216},
384 {1520,0.45033,2.59313,1.15453,2.11509,1.06998,1.19914,0.58667,-0.32702},
385 {1525,0.43753,2.5951,1.15631,2.11032,1.06994,1.19888,0.60734,-0.31178},
386 {1530,0.42376,2.59609,1.15782,2.10533,1.071,1.19702,0.62932,-0.29653},
387 {1535,0.40952,2.59624,1.15926,2.10037,1.07322,1.19365,0.65266,-0.28133},
388 {1540,0.39543,2.59565,1.16083,2.09565,1.07667,1.18887,0.6774,-0.26624},
389 {1545,0.38215,2.59444,1.16275,2.09142,1.08141,1.18276,0.7036,-0.25134},
390 {1550,0.37039,2.59273,1.1652,2.08789,1.08751,1.17544,0.73127,-0.2367},
391 {1555,0.36091,2.59062,1.16833,2.08526,1.09503,1.16716,0.76032,-0.22237},
392 {1560,0.35456,2.58822,1.17224,2.08366,1.10406,1.15822,0.79059,-0.20842},
393 {1565,0.35177,2.58562,1.17705,2.08327,1.11465,1.14894,0.82195,-0.19492},
394 {1570,0.35337,2.58293,1.18288,2.08422,1.1269,1.13961,0.85423,-0.18192},
395 {1575,0.35961,2.58025,1.18983,2.08669,1.14085,1.13056,0.88728,-0.1695},
396 {1580,0.37063,2.57768,1.19803,2.09082,1.1566,1.12209,0.92096,-0.1577},
397 {1585,0.3863,2.57533,1.20758,2.09678,1.17421,1.1145,0.95511,-0.1466},
398 {1590,0.40627,2.57329,1.2186,2.10471,1.19376,1.10811,0.98957,-0.13626},
399 {1595,0.43008,2.57166,1.2312,2.11477,1.21531,1.10322,1.0242,-0.12673},
400 {1600,0.45695,2.57056,1.2455,2.12712,1.23895,1.10014,1.05885,-0.11809},
401 {1605,0.48626,2.57007,1.2616,2.14191,1.26474,1.09918,1.09336,-0.1104},
402 {1610,0.51735,2.57029,1.27959,2.15926,1.29271,1.1006,1.12758,-0.1037},
403 {1615,0.54962,2.57116,1.2992,2.17896,1.3226,1.10434,1.16131,-0.09794},
404 {1620,0.58245,2.57256,1.32004,2.20065,1.35402,1.11019,1.19434,-0.09302},
405 {1625,0.61536,2.5744,1.34172,2.224,1.38657,1.11794,1.22646,-0.08883},
406 {1630,0.64786,2.57654,1.36383,2.24863,1.41985,1.12739,1.25746,-0.08529},
407 {1635,0.67945,2.57889,1.38599,2.27422,1.45346,1.13832,1.28713,-0.08227},
408 {1640,0.70969,2.58132,1.40779,2.30039,1.48701,1.15051,1.31525,-0.07969},
409 {1645,0.73811,2.58373,1.42884,2.32681,1.52009,1.16377,1.34161,-0.07744},
410 {1650,0.7643,2.58599,1.44873,2.35312,1.55232,1.17789,1.36601,-0.07541},
411 {1655,0.78784,2.58801,1.46709,2.37897,1.58329,1.19264,1.38822,-0.07352},
412 {1660,0.80845,2.58964,1.4836,2.40404,1.61266,1.20785,1.40808,-0.07166},
413 {1665,0.8259,2.59072,1.49807,2.42802,1.64013,1.22332,1.42543,-0.06977},
414 {1670,0.84003,2.5911,1.51034,2.45063,1.6654,1.2389,1.44012,-0.06776},
415 {1675,0.85087,2.59061,1.52021,2.47156,1.68818,1.25439,1.45202,-0.06558},
416 {1680,0.85853,2.58908,1.52752,2.49054,1.7082,1.26964,1.46098,-0.06314},
417 {1685,0.86321,2.58636,1.53208,2.50725,1.72514,1.28447,1.46686,-0.06037},
418 {1690,0.86521,2.58228,1.53373,2.52141,1.73872,1.29869,1.4695,-0.05719},
419 {1695,0.86497,2.57669,1.53227,2.53272,1.74866,1.31215,1.46878,-0.05354},
420 {1700,0.86261,2.56941,1.52753,2.54089,1.75465,1.32465,1.46453,-0.04935},
421 {1705,0.85865,2.5603,1.51934,2.54562,1.7564,1.33604,1.45663,-0.04453},
422 {1710,0.85338,2.54921,1.50758,2.54666,1.75368,1.34611,1.44499,-0.03906},
423 {1715,0.84709,2.53615,1.49237,2.54388,1.74639,1.35461,1.42984,-0.03313},
424 {1720,0.83998,2.52115,1.4739,2.5372,1.73451,1.36127,1.41148,-0.02693},
425 {1725,0.83223,2.50424,1.45236,2.52653,1.718,1.36581,1.3902,-0.02071},
426 {1730,0.82398,2.48546,1.42794,2.51179,1.69681,1.36795,1.36629,-0.01468},
427 {1735,0.81532,2.46484,1.40083,2.49289,1.67092,1.36742,1.34005,-0.00906},
428 {1740,0.80635,2.4424,1.37121,2.46973,1.6403,1.36394,1.31177,-0.00407},
429 {1745,0.79717,2.4182,1.33927,2.44227,1.60491,1.35727,1.28174,8e-05},
430 {1750,0.78787,2.39228,1.30516,2.41059,1.56494,1.34739,1.25012,0.00337},
431 {1755,0.77864,2.36474,1.26903,2.37487,1.52061,1.33437,1.21704,0.00581},
432 {1760,0.76947,2.33565,1.231,2.33528,1.47217,1.31829,1.18262,0.00745},
433 {1765,0.76063,2.3051,1.19124,2.29197,1.41987,1.29923,1.14697,0.00832},
434 {1770,0.75229,2.27315,1.14986,2.24514,1.36393,1.27726,1.11021,0.00845},
435 {1775,0.74462,2.23989,1.10702,2.19493,1.3046,1.25247,1.07248,0.00788},
436 {1780,0.73776,2.20541,1.06285,2.14153,1.24212,1.22492,1.03388,0.00663},
437 {1785,0.73178,2.16977,1.0175,2.0851,1.17673,1.1947,0.99455,0.00475},
438 {1790,0.72674,2.13307,0.9711,2.02582,1.10867,1.16189,0.95459,0.00226},
439 {1795,0.7226,2.09537,0.9238,1.96385,1.03819,1.12655,0.91414,-0.00079},
440 {1800,0.7193,2.05676,0.87573,1.89936,0.96551,1.08878,0.8733,-0.00439},
441 {1805,0.71674,2.01732,0.82704,1.83253,0.89088,1.04864,0.83221,-0.00849},
442 {1810,0.71485,1.97712,0.77787,1.76351,0.81455,1.00621,0.79099,-0.01305},
443 {1815,0.7135,1.93625,0.72835,1.69249,0.73675,0.96158,0.74975,-0.01806},
444 {1820,0.71257,1.89479,0.67863,1.61963,0.65772,0.91481,0.70861,-0.02347},
445 {1825,0.71198,1.85282,0.62885,1.5451,0.5777,0.86598,0.66769,-0.02924},
446 {1830,0.71165,1.81041,0.57915,1.46907,0.49694,0.81518,0.62713,-0.03536},
447 {1835,0.7115,1.76765,0.52966,1.39172,0.41567,0.76248,0.58703,-0.04177},
448 {1840,0.71149,1.72461,0.48053,1.3132,0.33413,0.70795,0.54752,-0.04846},
449 {1845,0.71158,1.68139,0.43186,1.23376,0.25256,0.65182,0.50865,-0.05537}},
450 {// K- p -> Sigma0 pi0
451 {435,0.01082,1.10989,0.01663,-0.00335,0,0,0,0},
452 {440,-0.01888,1.07693,-0.00805,-0.01206,0,0,0,0},
453 {445,-0.04853,1.04408,-0.03276,-0.02075,0,0,0,0},
454 {450,-0.07803,1.01152,-0.0575,-0.02944,0,0,0,0},
455 {455,-0.10733,0.97938,-0.0823,-0.03809,0,0,0,0},
456 {460,-0.13633,0.9478,-0.10716,-0.0467,0,0,0,0},
457 {465,-0.16491,0.91681,-0.13197,-0.0552,0,0,0,0},
458 {470,-0.19294,0.88642,-0.15661,-0.06351,0,0,0,0},
459 {475,-0.2203,0.85665,-0.18094,-0.07154,0,0,0,0},
460 {480,-0.24683,0.82754,-0.20479,-0.07921,0,0,0,0},
461 {485,-0.27226,0.79925,-0.22775,-0.08632,0,0,0,0},
462 {490,-0.29631,0.77198,-0.24937,-0.09267,0,0,0,0},
463 {495,-0.31868,0.74592,-0.26919,-0.09804,0,0,0,0},
464 {500,-0.33912,0.72125,-0.28682,-0.10221,0,0,0,0},
465 {505,-0.35756,0.69815,-0.30212,-0.10482,0,0,0,0},
466 {510,-0.37395,0.67677,-0.31501,-0.10547,0,0,0,0},
467 {515,-0.38825,0.65729,-0.32543,-0.10379,0,0,0,0},
468 {520,-0.40045,0.63982,-0.33324,-0.09964,0,0,0,0},
469 {525,-0.41054,0.62436,-0.33823,-0.09347,0,0,0,0},
470 {530,-0.41855,0.61093,-0.34015,-0.0858,0,0,0,0},
471 {535,-0.42447,0.59951,-0.33876,-0.07716,0,0,0,0},
472 {540,-0.42842,0.59011,-0.33407,-0.06792,0,0,0,0},
473 {545,-0.43073,0.5827,-0.32663,-0.05813,0,0,0,0},
474 {550,-0.43174,0.57727,-0.31706,-0.04781,0,0,0,0},
475 {555,-0.4318,0.57379,-0.30599,-0.03697,0,0,0,0},
476 {560,-0.43126,0.57219,-0.29405,-0.0257,0,0,0,0},
477 {565,-0.43044,0.57227,-0.28185,-0.01436,0,0,0,0},
478 {570,-0.4298,0.57397,-0.26991,-0.00328,0,0,0,0},
479 {575,-0.42982,0.57729,-0.25874,0.0072,0,0,0,0},
480 {580,-0.4308,0.58214,-0.24872,0.01687,0,0,0,0},
481 {585,-0.4328,0.58835,-0.2401,0.02571,0,0,0,0},
482 {590,-0.43612,0.59583,-0.23321,0.03376,0,0,0,0},
483 {595,-0.44115,0.6045,-0.22839,0.04108,0,0,0,0},
484 {600,-0.4482,0.61432,-0.22597,0.04772,0,0,0,0},
485 {605,-0.45692,0.62534,-0.22614,0.0535,0,0,0,0},
486 {610,-0.46673,0.63752,-0.22906,0.05801,0,0,0,0},
487 {615,-0.47702,0.65085,-0.23489,0.06082,0,0,0,0},
488 {620,-0.48723,0.66531,-0.24375,0.06153,0,0,0,0},
489 {625,-0.4968,0.68072,-0.25539,0.06011,0,0,0,0},
490 {630,-0.50526,0.69682,-0.26923,0.05678,0,0,0,0},
491 {635,-0.51223,0.71341,-0.28474,0.05178,0,0,0,0},
492 {640,-0.51759,0.73044,-0.30159,0.04536,0,0,0,0},
493 {645,-0.52116,0.74789,-0.31929,0.0379,0,0,0,0},
494 {650,-0.52267,0.76582,-0.33727,0.02987,0,0,0,0},
495 {655,-0.52186,0.78423,-0.35493,0.02172,0,0,0,0},
496 {660,-0.51846,0.80318,-0.3717,0.01394,0,0,0,0},
497 {665,-0.51218,0.82275,-0.387,0.007,0,0,0,0},
498 {670,-0.50273,0.84317,-0.40028,0.00144,0,0,0,0},
499 {675,-0.48982,0.86469,-0.41096,-0.0022,0,0,0,0},
500 {680,-0.47326,0.88743,-0.41875,-0.00347,0,0,0,0},
501 {685,-0.45292,0.91119,-0.42381,-0.0022,0,0,0,0},
502 {690,-0.42867,0.93557,-0.42646,0.00161,0,0,0,0},
503 {695,-0.40062,0.96016,-0.42703,0.00789,0,0,0,0},
504 {700,-0.36903,0.98453,-0.42585,0.01652,0,0,0,0},
505 {705,-0.33405,1.00813,-0.42325,0.02726,0,0,0,0},
506 {710,-0.29565,1.03016,-0.41947,0.03953,0,0,0,0},
507 {715,-0.25378,1.0498,-0.41477,0.05274,0,0,0,0},
508 {720,-0.20851,1.06626,-0.40938,0.06618,0,0,0,0},
509 {725,-0.16053,1.07916,-0.40351,0.07896,0,0,0,0},
510 {730,-0.11059,1.08856,-0.39716,0.09047,0,0,0,0},
511 {735,-0.05947,1.09453,-0.3903,0.10015,0,0,0,0},
512 {740,-0.00792,1.09717,-0.38288,0.10747,0,0,0,0},
513 {745,0.0433,1.09656,-0.37485,0.11203,0,0,0,0},
514 {750,0.09344,1.09277,-0.36603,0.11403,0,0,0,0},
515 {755,0.14178,1.0859,-0.35625,0.11383,0,0,0,0},
516 {760,0.18761,1.07615,-0.34549,0.1117,0,0,0,0},
517 {765,0.23037,1.06383,-0.33392,0.10794,0,0,0,0},
518 {770,0.26959,1.04936,-0.32217,0.1028,0,0,0,0},
519 {775,0.30481,1.03313,-0.311,0.09648,0,0,0,0},
520 {780,0.33569,1.0155,-0.3011,0.08904,0,0,0,0},
521 {785,0.36209,0.9969,-0.29309,0.08029,0,0,0,0},
522 {790,0.38386,0.97771,-0.28758,0.07006,0,0,0,0},
523 {795,0.40088,0.95833,-0.28517,0.05815,0,0,0,0},
524 {800,0.41307,0.93883,-0.28627,0.04427,0,0,0,0},
525 {805,0.42048,0.91889,-0.29099,0.02801,0,0,0,0},
526 {810,0.42339,0.89844,-0.29914,0.0094,0,0,0,0},
527 {815,0.42217,0.87751,-0.31045,-0.01136,0,0,0,0},
528 {820,0.4172,0.85614,-0.32463,-0.03409,0,0,0,0},
529 {825,0.40886,0.83439,-0.34137,-0.05857,0,0,0,0},
530 {830,0.39754,0.8123,-0.36041,-0.08462,0,0,0,0},
531 {835,0.38361,0.7899,-0.38143,-0.11203,0,0,0,0},
532 {840,0.36745,0.76725,-0.40416,-0.14061,0,0,0,0},
533 {845,0.34944,0.74439,-0.4283,-0.17016,0,0,0,0},
534 {850,0.32996,0.72136,-0.45357,-0.20047,0,0,0,0},
535 {855,0.30939,0.6982,-0.47967,-0.23135,0,0,0,0},
536 {860,0.2881,0.67496,-0.50631,-0.26261,0,0,0,0},
537 {865,0.26649,0.65168,-0.53321,-0.29404,0,0,0,0},
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.17259,-0.00101,-0.00057,0,1e-05},
636 {440,-0.41233,0.48944,-0.38429,-0.1651,-0.00059,-0.00043,0,1e-05},
637 {445,-0.4096,0.46112,-0.37795,-0.15764,-0.00017,-0.00029,0,0},
638 {450,-0.40685,0.4332,-0.37154,-0.15026,0.00025,-0.00016,0,0},
639 {455,-0.40408,0.40591,-0.36504,-0.14301,0.00067,-2e-05,0,0},
640 {460,-0.40131,0.37946,-0.35841,-0.13594,0.00109,0.00012,0,0},
641 {465,-0.39874,0.35403,-0.35162,-0.12907,0.00152,0.00025,0,0},
642 {470,-0.39658,0.32979,-0.34467,-0.12244,0.00194,0.00039,0,0},
643 {475,-0.39505,0.30691,-0.33753,-0.11608,0.00237,0.00053,0,0},
644 {480,-0.3944,0.28554,-0.3302,-0.11002,0.00279,0.00067,0,-1e-05},
645 {485,-0.39493,0.26577,-0.32275,-0.10426,0.00322,0.00081,0,-1e-05},
646 {490,-0.39696,0.24765,-0.31529,-0.09883,0.00364,0.00095,0,-1e-05},
647 {495,-0.4008,0.23126,-0.30789,-0.09372,0.00407,0.00109,0,-1e-05},
648 {500,-0.40672,0.21665,-0.30065,-0.08895,0.00449,0.00124,0,-1e-05},
649 {505,-0.41477,0.20379,-0.2936,-0.08456,0.00491,0.00138,0,-1e-05},
650 {510,-0.42493,0.19264,-0.2868,-0.0806,0.00533,0.00152,0,-1e-05},
651 {515,-0.4372,0.18316,-0.28027,-0.07711,0.00575,0.00167,0,-2e-05},
652 {520,-0.45145,0.17527,-0.27403,-0.07413,0.00616,0.00181,0,-2e-05},
653 {525,-0.46723,0.16882,-0.268,-0.07163,0.00656,0.00195,0,-2e-05},
654 {530,-0.48408,0.16365,-0.26212,-0.06961,0.00696,0.0021,0,-2e-05},
655 {535,-0.50152,0.15962,-0.25632,-0.06805,0.00734,0.00224,0,-2e-05},
656 {540,-0.5191,0.1566,-0.25056,-0.06694,0.00772,0.00238,0,-2e-05},
657 {545,-0.53644,0.15456,-0.24489,-0.06628,0.00807,0.00251,0,-3e-05},
658 {550,-0.55315,0.15349,-0.23935,-0.06605,0.00841,0.00264,0,-3e-05},
659 {555,-0.56883,0.15338,-0.23399,-0.06626,0.00872,0.00277,0,-3e-05},
660 {560,-0.58327,0.15419,-0.22886,-0.0669,0.00901,0.0029,0,-3e-05},
661 {565,-0.59663,0.15588,-0.2239,-0.06794,0.00927,0.00301,0,-3e-05},
662 {570,-0.60911,0.15838,-0.21909,-0.06938,0.00949,0.00312,0,-3e-05},
663 {575,-0.62091,0.16164,-0.21439,-0.07118,0.00967,0.00322,0,-4e-05},
664 {580,-0.63222,0.1656,-0.20978,-0.07333,0.00981,0.00331,0,-4e-05},
665 {585,-0.64322,0.17023,-0.20527,-0.07583,0.0099,0.00339,0,-4e-05},
666 {590,-0.65408,0.17551,-0.20089,-0.07867,0.00993,0.00345,0,-4e-05},
667 {595,-0.66499,0.18138,-0.19664,-0.08184,0.0099,0.0035,0,-4e-05},
668 {600,-0.6761,0.18784,-0.19254,-0.08532,0.0098,0.00353,0,-4e-05},
669 {605,-0.68744,0.19486,-0.18869,-0.08913,0.00962,0.00354,0,-4e-05},
670 {610,-0.69895,0.2025,-0.18519,-0.09324,0.00936,0.00353,0,-4e-05},
671 {615,-0.71054,0.21077,-0.18219,-0.09768,0.009,0.00349,0,-4e-05},
672 {620,-0.72215,0.21971,-0.17979,-0.10242,0.00855,0.00343,0,-4e-05},
673 {625,-0.73381,0.22936,-0.17813,-0.1075,0.00798,0.00333,0,-4e-05},
674 {630,-0.74564,0.23978,-0.17729,-0.11293,0.0073,0.00321,0,-4e-05},
675 {635,-0.75775,0.25101,-0.17738,-0.11874,0.0065,0.00305,0,-4e-05},
676 {640,-0.77025,0.26312,-0.17851,-0.12496,0.00555,0.00285,1e-05,-4e-05},
677 {645,-0.78311,0.27609,-0.18077,-0.1316,0.00446,0.00261,1e-05,-4e-05},
678 {650,-0.79619,0.28987,-0.18424,-0.13869,0.00322,0.00232,1e-05,-3e-05},
679 {655,-0.80933,0.30439,-0.18901,-0.14623,0.00181,0.00199,1e-05,-3e-05},
680 {660,-0.8224,0.3196,-0.19518,-0.15426,0.00022,0.0016,2e-05,-3e-05},
681 {665,-0.83513,0.3354,-0.2028,-0.16277,-0.00156,0.00116,2e-05,-3e-05},
682 {670,-0.84713,0.35164,-0.21197,-0.1718,-0.00353,0.00066,2e-05,-2e-05},
683 {675,-0.85802,0.36815,-0.22273,-0.18137,-0.00571,1e-04,3e-05,-2e-05},
684 {680,-0.86734,0.38479,-0.23516,-0.19149,-0.00811,-0.00053,3e-05,-1e-05},
685 {685,-0.87414,0.4013,-0.24915,-0.2022,-0.01074,-0.00123,4e-05,0},
686 {690,-0.87701,0.41738,-0.26452,-0.2135,-0.0136,-0.00201,4e-05,0},
687 {695,-0.87457,0.43271,-0.28104,-0.22543,-0.01672,-0.00286,4e-05,1e-05},
688 {700,-0.86542,0.447,-0.29851,-0.238,-0.02008,-0.00379,4e-05,2e-05},
689 {705,-0.8484,0.46,-0.31674,-0.25122,-0.02371,-0.00481,5e-05,3e-05},
690 {710,-0.82293,0.47163,-0.33558,-0.26507,-0.02761,-0.00591,5e-05,4e-05},
691 {715,-0.7885,0.48183,-0.35485,-0.2795,-0.03179,-0.00711,4e-05,5e-05},
692 {720,-0.74459,0.49053,-0.37441,-0.29449,-0.03625,-0.0084,4e-05,7e-05},
693 {725,-0.69112,0.49769,-0.39405,-0.31,-0.041,-0.00979,3e-05,8e-05},
694 {730,-0.62899,0.50329,-0.41345,-0.32601,-0.04604,-0.01127,2e-05,9e-05},
695 {735,-0.55919,0.50731,-0.43227,-0.34248,-0.05136,-0.01286,1e-05,0.00011},
696 {740,-0.48275,0.50974,-0.45017,-0.3594,-0.05697,-0.01455,-1e-05,0.00013},
697 {745,-0.40083,0.51061,-0.46688,-0.37671,-0.06287,-0.01634,-3e-05,0.00015},
698 {750,-0.31522,0.51019,-0.48236,-0.39436,-0.06905,-0.01823,-6e-05,0.00017},
699 {755,-0.22787,0.50878,-0.49662,-0.41225,-0.0755,-0.02023,-1e-04,0.00019},
700 {760,-0.14072,0.5067,-0.5097,-0.43031,-0.08221,-0.02233,-0.00014,0.00021},
701 {765,-0.05569,0.50425,-0.52159,-0.44847,-0.08917,-0.02453,-0.00019,0.00023},
702 {770,0.0256,0.50172,-0.53224,-0.46666,-0.09637,-0.02683,-0.00024,0.00026},
703 {775,0.10169,0.49939,-0.54158,-0.48485,-0.10377,-0.02923,-0.00031,0.00028},
704 {780,0.17164,0.49748,-0.54958,-0.503,-0.11136,-0.03171,-0.00038,0.00031},
705 {785,0.23528,0.49609,-0.55635,-0.52109,-0.11908,-0.03426,-0.00045,0.00034},
706 {790,0.29246,0.49532,-0.56201,-0.53909,-0.1269,-0.03689,-0.00053,0.00037},
707 {795,0.34305,0.49529,-0.56666,-0.55697,-0.13479,-0.03958,-0.00061,4e-04},
708 {800,0.38717,0.49602,-0.57044,-0.5747,-0.14269,-0.04231,-0.00069,0.00043},
709 {805,0.42531,0.49745,-0.57352,-0.59223,-0.15055,-0.04508,-0.00077,0.00046},
710 {810,0.45797,0.49951,-0.57609,-0.60951,-0.15829,-0.04787,-0.00084,0.00049},
711 {815,0.48565,0.50213,-0.57833,-0.6265,-0.16587,-0.05066,-0.00089,0.00052},
712 {820,0.50886,0.50523,-0.58041,-0.64314,-0.17322,-0.05344,-0.00093,0.00056},
713 {825,0.5281,0.50875,-0.58253,-0.65939,-0.18027,-0.05618,-0.00095,0.00059},
714 {830,0.54388,0.51262,-0.58485,-0.6752,-0.18697,-0.05889,-0.00095,0.00062},
715 {835,0.55669,0.51675,-0.58757,-0.69052,-0.19326,-0.06154,-0.00091,0.00065},
716 {840,0.56705,0.5211,-0.59085,-0.70531,-0.19907,-0.06411,-0.00084,0.00069},
717 {845,0.57545,0.52557,-0.59488,-0.71951,-0.20435,-0.06659,-0.00074,0.00072},
718 {850,0.5824,0.53011,-0.59985,-0.73309,-0.20902,-0.06896,-0.00059,0.00075},
719 {855,0.5884,0.53464,-0.60592,-0.74598,-0.21304,-0.07122,-4e-04,0.00078},
720 {860,0.59395,0.53909,-0.61329,-0.75815,-0.21633,-0.07334,-0.00016,0.00081},
721 {865,0.59954,0.54338,-0.62212,-0.76955,-0.21885,-0.0753,0.00013,0.00084},
722 {870,0.60541,0.54745,-0.63256,-0.78013,-0.22056,-0.0771,0.00048,0.00087},
723 {875,0.61159,0.55119,-0.64468,-0.78984,-0.22145,-0.07869,0.00089,9e-04},
724 {880,0.61811,0.55452,-0.65859,-0.79866,-0.22151,-0.08006,0.00137,0.00093},
725 {885,0.62501,0.55734,-0.67437,-0.80652,-0.22076,-0.08119,0.00192,0.00095},
726 {890,0.63226,0.55962,-0.69201,-0.81339,-0.21917,-0.08204,0.00254,0.00097},
727 {895,0.63991,0.56144,-0.71124,-0.81921,-0.21677,-0.08265,0.00322,0.00099},
728 {900,0.6482,0.56292,-0.73165,-0.82391,-0.21351,-0.08304,0.00396,0.00101},
729 {905,0.65734,0.56418,-0.7528,-0.82746,-0.20938,-0.08324,0.00472,0.00102},
730 {910,0.66707,0.56537,-0.77433,-0.82982,-0.20437,-0.08329,0.0055,0.00103},
731 {915,0.67677,0.56667,-0.79588,-0.83094,-0.19849,-0.08322,0.00626,0.00103},
732 {920,0.68581,0.56824,-0.81712,-0.83079,-0.19176,-0.08306,0.00696,0.00102},
733 {925,0.6938,0.57027,-0.83765,-0.82936,-0.18417,-0.08286,0.00755,0.00101},
734 {930,0.70068,0.57295,-0.85705,-0.8266,-0.17573,-0.08266,0.00799,0.00099},
735 {935,0.70637,0.57649,-0.8749,-0.8225,-0.16642,-0.08249,0.00819,0.00097},
736 {940,0.71082,0.58109,-0.89077,-0.81703,-0.15624,-0.08242,0.00809,0.00093},
737 {945,0.71395,0.58698,-0.90423,-0.81016,-0.14519,-0.08246,0.0076,0.00088},
738 {950,0.71565,0.59441,-0.91489,-0.80182,-0.13326,-0.08261,0.00663,0.00082},
739 {955,0.71583,0.60361,-0.92234,-0.79199,-0.12046,-0.08291,0.00506,0.00075},
740 {960,0.71439,0.61483,-0.92615,-0.78061,-0.10679,-0.08336,0.0028,0.00067},
741 {965,0.71134,0.62829,-0.92599,-0.76765,-0.09225,-0.08398,-0.00027,0.00056},
742 {970,0.70709,0.64411,-0.92166,-0.7531,-0.07689,-0.08479,-0.00425,0.00045},
743 {975,0.7021,0.66238,-0.91307,-0.73694,-0.06076,-0.0858,-0.00928,0.00031},
744 {980,0.69672,0.68301,-0.90027,-0.7192,-0.04391,-0.08702,-0.01547,0.00016},
745 {985,0.69117,0.70589,-0.8834,-0.69991,-0.02642,-0.08844,-0.02292,-1e-05},
746 {990,0.68572,0.73088,-0.8626,-0.6791,-0.00834,-0.09004,-0.03174,-0.00021},
747 {995,0.6806,0.75785,-0.838,-0.6568,0.01026,-0.09181,-0.04202,-0.00043},
748 {1000,0.67606,0.78666,-0.80974,-0.63304,0.02931,-0.09375,-0.05388,-0.00067},
749 {1005,0.67236,0.81718,-0.77798,-0.60786,0.04877,-0.09583,-0.06739,-0.00094},
750 {1010,0.66984,0.84916,-0.74299,-0.5813,0.06857,-0.09799,-0.08243,-0.00124},
751 {1015,0.66874,0.88239,-0.70509,-0.5534,0.08867,-0.10021,-0.09858,-0.00156},
752 {1020,0.66925,0.91668,-0.66453,-0.52422,0.10903,-0.10243,-0.11538,-0.00191},
753 {1025,0.67157,0.95184,-0.62161,-0.49382,0.12959,-0.10462,-0.13236,-0.00229},
754 {1030,0.67589,0.98766,-0.57661,-0.46226,0.15032,-0.10673,-0.14906,-0.00271},
755 {1035,0.68241,1.02396,-0.52979,-0.42958,0.17118,-0.10873,-0.165,-0.00316},
756 {1040,0.69132,1.06053,-0.48143,-0.39584,0.1921,-0.11057,-0.17973,-0.00364},
757 {1045,0.70281,1.09718,-0.43183,-0.36109,0.21306,-0.11223,-0.19278,-0.00415},
758 {1050,0.71699,1.13373,-0.38124,-0.3254,0.234,-0.11366,-0.2038,-0.0047},
759 {1055,0.73356,1.17007,-0.32997,-0.28883,0.25489,-0.11484,-0.21289,-0.00529},
760 {1060,0.75214,1.20612,-0.2783,-0.25144,0.27572,-0.11578,-0.22029,-0.00591},
761 {1065,0.77234,1.24178,-0.22653,-0.2133,0.29644,-0.11646,-0.22622,-0.00657},
762 {1070,0.79378,1.27697,-0.17493,-0.17448,0.31703,-0.11688,-0.23091,-0.00727},
763 {1075,0.81607,1.31159,-0.1238,-0.13504,0.33746,-0.11701,-0.23458,-0.008},
764 {1080,0.83883,1.34556,-0.07343,-0.09504,0.3577,-0.11687,-0.23746,-0.00877},
765 {1085,0.86167,1.37879,-0.02411,-0.05456,0.37772,-0.11643,-0.23978,-0.00959},
766 {1090,0.88422,1.41119,0.02389,-0.01366,0.39751,-0.11569,-0.24169,-0.01044},
767 {1095,0.90616,1.44276,0.07039,0.02762,0.41704,-0.11467,-0.24307,-0.01132},
768 {1100,0.92718,1.47349,0.11521,0.06924,0.43632,-0.11339,-0.24371,-0.01225},
769 {1105,0.94696,1.50338,0.15818,0.11115,0.45534,-0.11187,-0.24342,-0.01321},
770 {1110,0.9652,1.53241,0.19913,0.15331,0.47409,-0.11012,-0.242,-0.0142},
771 {1115,0.9816,1.5606,0.23789,0.19569,0.49258,-0.10817,-0.23926,-0.01523},
772 {1120,0.99584,1.58794,0.27429,0.23824,0.51078,-0.10603,-0.23498,-0.01629},
773 {1125,1.00761,1.61442,0.30816,0.28092,0.52871,-0.10373,-0.22898,-0.01738},
774 {1130,1.01668,1.64004,0.33938,0.32369,0.54636,-0.10127,-0.22118,-0.0185},
775 {1135,1.02313,1.66484,0.36806,0.36651,0.56375,-0.09867,-0.21198,-0.01964},
776 {1140,1.02709,1.68883,0.39433,0.40934,0.58091,-0.09593,-0.20188,-0.02081},
777 {1145,1.02872,1.71205,0.41837,0.45217,0.59787,-0.09304,-0.1914,-0.02199},
778 {1150,1.02817,1.73453,0.4403,0.49494,0.61465,-0.09001,-0.18106,-0.02319},
779 {1155,1.02559,1.75629,0.4603,0.53762,0.63128,-0.08684,-0.17138,-0.0244},
780 {1160,1.02111,1.77737,0.47851,0.58018,0.64779,-0.08352,-0.16287,-0.02561},
781 {1165,1.01489,1.7978,0.49508,0.62259,0.6642,-0.08006,-0.15604,-0.02682},
782 {1170,1.00709,1.81759,0.51017,0.66481,0.68055,-0.07647,-0.15123,-0.02803},
783 {1175,0.99786,1.83673,0.52401,0.70683,0.69689,-0.07271,-0.14805,-0.02923},
784 {1180,0.98737,1.8552,0.53683,0.74862,0.71327,-0.06878,-0.14595,-0.03041},
785 {1185,0.9758,1.87296,0.54887,0.79018,0.72977,-0.06466,-0.14436,-0.03156},
786 {1190,0.96331,1.89,0.56034,0.8315,0.74642,-0.06033,-0.14272,-0.03268},
787 {1195,0.95008,1.90628,0.57149,0.87255,0.7633,-0.05577,-0.14046,-0.03375},
788 {1200,0.93626,1.92179,0.58256,0.91332,0.78046,-0.05097,-0.13703,-0.03476},
789 {1205,0.92203,1.9365,0.59376,0.95379,0.79796,-0.04591,-0.13185,-0.03572},
790 {1210,0.90757,1.95038,0.60531,0.99396,0.81585,-0.04057,-0.12448,-0.0366},
791 {1215,0.89304,1.96341,0.61726,1.03379,0.83412,-0.03495,-0.1149,-0.0374},
792 {1220,0.87864,1.97556,0.62967,1.07325,0.85278,-0.02905,-0.10319,-0.0381},
793 {1225,0.86453,1.9868,0.64256,1.11232,0.87181,-0.02286,-0.08944,-0.03869},
794 {1230,0.85091,1.9971,0.65596,1.15095,0.8912,-0.01639,-0.07374,-0.03915},
795 {1235,0.83795,2.00644,0.66993,1.18913,0.91096,-0.00962,-0.05619,-0.03947},
796 {1240,0.82582,2.01479,0.68448,1.22682,0.93106,-0.00255,-0.03687,-0.03963},
797 {1245,0.81473,2.02212,0.69966,1.26398,0.9515,0.00481,-0.01587,-0.03962},
798 {1250,0.80486,2.0284,0.71551,1.3006,0.97227,0.01247,0.00658,-0.03943},
799 {1255,0.7966,2.03353,0.73205,1.33662,0.99334,0.0204,0.02967,-0.03908},
800 {1260,0.79032,2.03741,0.74932,1.372,1.0147,0.02857,0.05245,-0.0386},
801 {1265,0.78642,2.03993,0.76735,1.40669,1.03632,0.03694,0.07397,-0.038},
802 {1270,0.78492,2.04109,0.78615,1.44068,1.05815,0.04549,0.09368,-0.03732},
803 {1275,0.78537,2.04098,0.80574,1.47397,1.08018,0.05421,0.11146,-0.03658},
804 {1280,0.78729,2.03971,0.8261,1.50657,1.10237,0.0631,0.12727,-0.0358},
805 {1285,0.79021,2.03737,0.84724,1.53849,1.12468,0.07215,0.14103,-0.035},
806 {1290,0.79375,2.03407,0.86916,1.56972,1.14706,0.08134,0.15264,-0.03421},
807 {1295,0.79786,2.02989,0.89186,1.60028,1.16947,0.09069,0.16178,-0.03347},
808 {1300,0.80263,2.02488,0.91536,1.63014,1.19182,0.10018,0.16808,-0.03285},
809 {1305,0.80811,2.01914,0.93965,1.65929,1.21405,0.10981,0.17119,-0.03239},
810 {1310,0.81435,2.01272,0.96476,1.68774,1.23609,0.11959,0.17074,-0.03214},
811 {1315,0.82144,2.0057,0.99069,1.71547,1.25787,0.12949,0.16635,-0.03216},
812 {1320,0.82936,1.99815,1.01744,1.74248,1.27932,0.13954,0.15783,-0.03249},
813 {1325,0.83779,1.99015,1.04496,1.76877,1.30041,0.14973,0.14596,-0.03314},
814 {1330,0.84665,1.9818,1.07325,1.79437,1.32111,0.1601,0.13175,-0.03405},
815 {1335,0.85598,1.97322,1.10228,1.81929,1.3414,0.17067,0.11614,-0.03517},
816 {1340,0.86584,1.96452,1.13207,1.84353,1.36125,0.18147,0.10006,-0.03645},
817 {1345,0.87626,1.9558,1.16259,1.86713,1.38065,0.19253,0.08445,-0.03781},
818 {1350,0.88729,1.94717,1.19385,1.89008,1.39956,0.20387,0.07024,-0.0392},
819 {1355,0.89897,1.93876,1.22584,1.9124,1.41797,0.21552,0.05839,-0.04056},
820 {1360,0.91133,1.93063,1.25852,1.9341,1.43584,0.22751,0.04964,-0.04184},
821 {1365,0.92428,1.92278,1.29175,1.95515,1.45313,0.23983,0.04407,-0.043},
822 {1370,0.93773,1.91515,1.32538,1.97551,1.46982,0.25248,0.04159,-0.04401},
823 {1375,0.95159,1.90768,1.35924,1.99516,1.48588,0.26547,0.0418,-0.04485},
824 {1380,0.96583,1.90033,1.39318,2.01405,1.50132,0.27876,0.04395,-0.04548},
825 {1385,0.98042,1.89304,1.42703,2.03215,1.51615,0.29237,0.04726,-0.0459},
826 {1390,0.9953,1.88573,1.46066,2.04944,1.53039,0.30627,0.05095,-0.04606},
827 {1395,1.01046,1.87835,1.49391,2.06587,1.54404,0.32046,0.05424,-0.04595},
828 {1400,1.02585,1.87084,1.52662,2.08142,1.55712,0.33492,0.05634,-0.04555},
829 {1405,1.04143,1.86314,1.55865,2.09605,1.56963,0.34966,0.05648,-0.04482},
830 {1410,1.05718,1.8552,1.58983,2.10973,1.58159,0.36465,0.05388,-0.04374},
831 {1415,1.07304,1.84694,1.62001,2.12243,1.59301,0.37988,0.04775,-0.04229},
832 {1420,1.08898,1.83832,1.64908,2.13412,1.6039,0.39535,0.03758,-0.04044},
833 {1425,1.10484,1.82934,1.67699,2.14484,1.6143,0.41106,0.02383,-0.03819},
834 {1430,1.12048,1.81999,1.70374,2.15463,1.62425,0.427,0.00724,-0.03552},
835 {1435,1.13574,1.81027,1.72932,2.16352,1.63379,0.44316,-0.01145,-0.03242},
836 {1440,1.15046,1.80018,1.75373,2.17157,1.64296,0.45955,-0.03152,-0.02889},
837 {1445,1.16449,1.78973,1.77696,2.17881,1.65179,0.47617,-0.05224,-0.0249},
838 {1450,1.17767,1.77892,1.799,2.18528,1.66033,0.493,-0.07288,-0.02046},
839 {1455,1.18985,1.76774,1.81985,2.19103,1.66861,0.51005,-0.09269,-0.01555},
840 {1460,1.20086,1.75621,1.8395,2.1961,1.67668,0.52732,-0.11096,-0.01015},
841 {1465,1.21057,1.74432,1.85795,2.20053,1.68457,0.54479,-0.12698,-0.00427},
842 {1470,1.219,1.7322,1.87524,2.20436,1.69232,0.56245,-0.14044,0.00207},
843 {1475,1.22629,1.72006,1.89148,2.20761,1.69992,0.58024,-0.15137,0.00883},
844 {1480,1.2326,1.7081,1.90674,2.21031,1.70741,0.59813,-0.15979,0.01594},
845 {1485,1.23807,1.69653,1.92113,2.2125,1.71479,0.61607,-0.16572,0.02336},
846 {1490,1.24286,1.68556,1.93474,2.21421,1.72208,0.63402,-0.16918,0.03103},
847 {1495,1.24712,1.67541,1.94765,2.21547,1.7293,0.65192,-0.17019,0.03891},
848 {1500,1.251,1.66628,1.95997,2.21629,1.73646,0.66975,-0.16879,0.04693},
849 {1505,1.25464,1.65837,1.97178,2.21673,1.74357,0.68744,-0.16499,0.05505},
850 {1510,1.25822,1.65191,1.98318,2.2168,1.75066,0.70496,-0.15881,0.06321},
851 {1515,1.26186,1.64709,1.99426,2.21654,1.75773,0.72227,-0.15028,0.07137},
852 {1520,1.26574,1.64404,2.00509,2.21596,1.7648,0.73931,-0.13954,0.07946},
853 {1525,1.27001,1.6427,2.01568,2.21509,1.77186,0.75605,-0.12701,0.08744},
854 {1530,1.27484,1.64297,2.02605,2.21392,1.77888,0.77244,-0.11314,0.09525},
855 {1535,1.28041,1.64477,2.03621,2.21246,1.78586,0.78844,-0.09836,0.10284},
856 {1540,1.28686,1.64799,2.04616,2.21072,1.79279,0.80401,-0.08314,0.11014},
857 {1545,1.29439,1.65256,2.05591,2.2087,1.79966,0.8191,-0.06791,0.11712},
858 {1550,1.3031,1.65836,2.06547,2.20642,1.80644,0.83367,-0.05309,0.12371},
859 {1555,1.3129,1.66525,2.07475,2.20382,1.81312,0.84769,-0.03882,0.1299},
860 {1560,1.32356,1.67306,2.08365,2.20085,1.81969,0.86115,-0.02508,0.1357},
861 {1565,1.33485,1.68162,2.09205,2.19746,1.82612,0.87402,-0.0119,0.14112},
862 {1570,1.34654,1.69074,2.09985,2.19357,1.83239,0.88627,0.00073,0.14615},
863 {1575,1.35842,1.70026,2.10691,2.18914,1.83849,0.89789,0.01281,0.1508},
864 {1580,1.37025,1.71001,2.11314,2.1841,1.8444,0.90885,0.02433,0.15507},
865 {1585,1.38181,1.71981,2.11843,2.17841,1.8501,0.91914,0.03529,0.15897},
866 {1590,1.39288,1.72949,2.12265,2.17199,1.85557,0.92872,0.04568,0.1625},
867 {1595,1.40323,1.73887,2.1257,2.16479,1.86078,0.93759,0.0555,0.16566},
868 {1600,1.41263,1.74778,2.12746,2.15675,1.86573,0.9457,0.06474,0.16846},
869 {1605,1.42086,1.75605,2.12782,2.14782,1.8704,0.95306,0.07341,0.1709},
870 {1610,1.42772,1.76352,2.12669,2.13793,1.87475,0.95962,0.08151,0.17299},
871 {1615,1.43328,1.77017,2.12406,2.12708,1.87879,0.96539,0.08916,0.17472},
872 {1620,1.43768,1.77599,2.11996,2.11527,1.88249,0.97035,0.09653,0.17611},
873 {1625,1.4411,1.78103,2.11445,2.10249,1.88584,0.97449,0.10377,0.17715},
874 {1630,1.44369,1.7853,2.10756,2.08876,1.88883,0.97781,0.11105,0.17786},
875 {1635,1.44561,1.78882,2.09933,2.07406,1.89145,0.98031,0.11854,0.17823},
876 {1640,1.44703,1.79161,2.0898,2.05839,1.89369,0.98196,0.1264,0.17827},
877 {1645,1.4481,1.7937,2.07902,2.04178,1.89552,0.98277,0.1348,0.17798},
878 {1650,1.449,1.79509,2.06703,2.0242,1.89693,0.98272,0.14389,0.17737},
879 {1655,1.44987,1.79582,2.05386,2.00566,1.89793,0.98181,0.15384,0.17644},
880 {1660,1.4508,1.79591,2.03956,1.98617,1.89846,0.98003,0.16463,0.17518},
881 {1665,1.45177,1.79536,2.02422,1.96574,1.89851,0.97737,0.17604,0.17354},
882 {1670,1.45273,1.79421,2.0079,1.94436,1.89801,0.9738,0.18781,0.1715},
883 {1675,1.45366,1.79245,1.99068,1.92205,1.89694,0.96931,0.1997,0.16901},
884 {1680,1.45453,1.79012,1.97262,1.89882,1.89525,0.9639,0.21145,0.16602},
885 {1685,1.4553,1.78723,1.95379,1.87466,1.8929,0.95754,0.22282,0.16251},
886 {1690,1.45594,1.78379,1.93428,1.84959,1.88984,0.95022,0.23355,0.15843},
887 {1695,1.45642,1.77982,1.91414,1.82361,1.88604,0.94193,0.24341,0.15375},
888 {1700,1.4567,1.77534,1.89345,1.79673,1.88144,0.93265,0.25213,0.14842},
889 {1705,1.45676,1.77036,1.87229,1.76895,1.87602,0.92237,0.25948,0.1424},
890 {1710,1.45653,1.76489,1.85069,1.74029,1.86972,0.91108,0.26518,0.13566},
891 {1715,1.45585,1.75888,1.8286,1.71072,1.86253,0.89877,0.26888,0.12821},
892 {1720,1.45454,1.75226,1.80593,1.68023,1.85444,0.88545,0.27022,0.12005},
893 {1725,1.45241,1.74498,1.78259,1.6488,1.84541,0.87112,0.26884,0.11119},
894 {1730,1.44926,1.73696,1.75849,1.61643,1.83543,0.8558,0.26436,0.10163},
895 {1735,1.44492,1.72815,1.73354,1.58309,1.82449,0.83947,0.25641,0.0914},
896 {1740,1.43919,1.71847,1.70766,1.54877,1.81256,0.82215,0.24464,0.08048},
897 {1745,1.43192,1.70788,1.68076,1.51346,1.79963,0.80384,0.22872,0.06891},
898 {1750,1.42321,1.69643,1.6529,1.4772,1.78573,0.78458,0.20876,0.05669},
899 {1755,1.41325,1.68419,1.62417,1.44003,1.77092,0.76442,0.18503,0.04387},
900 {1760,1.40226,1.67124,1.59466,1.40201,1.75523,0.74341,0.15781,0.03049},
901 {1765,1.39045,1.65768,1.56447,1.3632,1.73872,0.72159,0.12737,0.01657},
902 {1770,1.37802,1.64358,1.53368,1.32364,1.72144,0.69903,0.094,0.00214},
903 {1775,1.36518,1.62903,1.5024,1.28339,1.70344,0.67576,0.05798,-0.01274},
904 {1780,1.35214,1.61411,1.47071,1.24251,1.68477,0.65185,0.01957,-0.02806},
905 {1785,1.3391,1.59891,1.43872,1.20104,1.66547,0.62732,-0.02093,-0.04378},
906 {1790,1.32627,1.58351,1.4065,1.15905,1.6456,0.60225,-0.06326,-0.05985},
907 {1795,1.31385,1.56799,1.37416,1.11658,1.6252,0.57668,-0.10713,-0.07626},
908 {1800,1.30206,1.55243,1.3418,1.07368,1.60433,0.55065,-0.15227,-0.09296},
909 {1805,1.29111,1.53692,1.3095,1.03042,1.58303,0.52422,-0.1984,-0.10991},
910 {1810,1.28119,1.52155,1.27735,0.98684,1.56136,0.49744,-0.24524,-0.12709},
911 {1815,1.27252,1.50639,1.24546,0.943,1.53936,0.47035,-0.29251,-0.14447},
912 {1820,1.2653,1.49154,1.21391,0.89895,1.51708,0.44301,-0.33993,-0.16199},
913 {1825,1.25974,1.47706,1.18281,0.85475,1.49457,0.41547,-0.38723,-0.17964},
914 {1830,1.25605,1.46305,1.15224,0.81044,1.47189,0.38777,-0.43413,-0.19737},
915 {1835,1.25443,1.4496,1.12229,0.76609,1.44907,0.35997,-0.48035,-0.21516},
916 {1840,1.25509,1.43678,1.09306,0.72174,1.42618,0.33212,-0.5256,-0.23297},
917 {1845,1.25803,1.42461,1.06458,0.67744,1.40325,0.30425,-0.56974,-0.25076}}};
918
919
920
921 const G4int coef_ener = G4int((pLab-Legendre_coef[WhichChannel+1][0][0])/5);
922 const G4double sup_ener = pLab/5. - coef_ener - Legendre_coef[WhichChannel+1][0][0]/5;
923
924 // assert(pLab >= Legendre_coef[WhichChannel+1][coef_ener][0] && pLab < Legendre_coef[WhichChannel+1][coef_ener+1][0]);
925
926 // Legendre coefficient normalized
927 const G4double A0 = 1.;
928 const G4double A1 = (1-sup_ener)*Legendre_coef[WhichChannel+1][coef_ener][1] + sup_ener*Legendre_coef[WhichChannel+1][coef_ener+1][1];
929 const G4double A2 = (1-sup_ener)*Legendre_coef[WhichChannel+1][coef_ener][2] + sup_ener*Legendre_coef[WhichChannel+1][coef_ener+1][2];
930 const G4double A3 = (1-sup_ener)*Legendre_coef[WhichChannel+1][coef_ener][3] + sup_ener*Legendre_coef[WhichChannel+1][coef_ener+1][3];
931 const G4double A4 = (1-sup_ener)*Legendre_coef[WhichChannel+1][coef_ener][4] + sup_ener*Legendre_coef[WhichChannel+1][coef_ener+1][4];
932 const G4double A5 = (1-sup_ener)*Legendre_coef[WhichChannel+1][coef_ener][5] + sup_ener*Legendre_coef[WhichChannel+1][coef_ener+1][5];
933 const G4double A6 = (1-sup_ener)*Legendre_coef[WhichChannel+1][coef_ener][6] + sup_ener*Legendre_coef[WhichChannel+1][coef_ener+1][6];
934 const G4double A7 = (1-sup_ener)*Legendre_coef[WhichChannel+1][coef_ener][7] + sup_ener*Legendre_coef[WhichChannel+1][coef_ener+1][7];
935 const G4double A8 = (1-sup_ener)*Legendre_coef[WhichChannel+1][coef_ener][8] + sup_ener*Legendre_coef[WhichChannel+1][coef_ener+1][8];
936
937 // Theoritical max if all Ai > 0 (often the case)
938 const G4double A = std::fabs(A0) + std::fabs(A1) + std::fabs(A2) + std::fabs(A3) + std::fabs(A4) + std::fabs(A5) + std::fabs(A6) + std::fabs(A7) + std::fabs(A8);
939
940 G4bool success = false;
941 G4int maxloop = 0;
942
943 while(!success && maxloop < 1000){
944
945 cos_theta = Random::shoot()*2-1.; // not optimized
946
947 // Legendre Polynomial
948 G4double P0 = A0;
949 G4double P1 = A1*cos_theta;
950 G4double P2 = A2/2.*(3*std::pow(cos_theta,2)-1);
951 G4double P3 = A3/2.*(5*std::pow(cos_theta,3)-3*cos_theta);
952 G4double P4 = A4/8.*(35*std::pow(cos_theta,4)-30*std::pow(cos_theta,2)+3);
953 G4double P5 = A5/8.*(63*std::pow(cos_theta,5)-70*std::pow(cos_theta,3)+15*cos_theta);
954 G4double P6 = A6/16.*(231*std::pow(cos_theta,6)-315*std::pow(cos_theta,4)+105*std::pow(cos_theta,2)-5);
955 G4double P7 = A7/16.*(429*std::pow(cos_theta,7)-693*std::pow(cos_theta,5)+315*std::pow(cos_theta,3)-35*cos_theta);
956 G4double P8 = A8/128.*(6435*std::pow(cos_theta,8)-12012*std::pow(cos_theta,6)+6930*std::pow(cos_theta,4)-1260*std::pow(cos_theta,2)+35);
957
958 G4double P = (P0 + P1 + P2 + P3 + P4 + P5 + P6 + P7 + P8)/2; // /2 for the normalisation
959
960 if(Random::shoot()*A < P) success = true;
961 maxloop +=1 ;
962 if(maxloop==1000) cos_theta = std::log(Random::shoot()*(std::exp(10.)-std::exp(-10.))+std::exp(-10.))/10.; // if no success in 1E4 shoot, probably angulard distribution piked very forward
963 }
964 sin_theta = std::sqrt(1-cos_theta*cos_theta);
965 }
966
967 if(rho == 0) return ThreeVector(sin_theta*cos_phi,sin_theta*sin_phi,cos_theta);
968 // Rotation in the direction of the incident kaon
969 const G4double px = x/r*cos_theta - y/rho*sin_theta*cos_phi + z/r*x/rho*sin_theta*sin_phi;
970 const G4double py = y/r*cos_theta + x/rho*sin_theta*cos_phi + z/r*y/rho*sin_theta*sin_phi;
971 const G4double pz = z/r*cos_theta - rho/r*sin_theta*sin_phi;
972
973 return ThreeVector(px,py,pz);
974 }
975
976 }
977