| 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 "G4INCLNKbToLpiChannel.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 NKbToLpiChannel::NKbToLpiChannel(Particle *p1, Particle *p2)
50 : particle1(p1), particle2(p2)
51 {}
52
53 NKbToLpiChannel::~NKbToLpiChannel(){}
54
55 void NKbToLpiChannel::fillFinalState(FinalState *fs) {
56
57 Particle *nucleon;
58 Particle *kaon;
59
60 if(particle1->isNucleon()){
61 nucleon = particle1;
62 kaon = particle2;
63 }
64 else{
65 nucleon = particle2;
66 kaon = particle1;
67 }
68 const G4int iso = ParticleTable::getIsospin(nucleon->getType()) + ParticleTable::getIsospin(kaon->getType());
69
70
71 ThreeVector mom_pion = KaonMomentum(kaon,nucleon); // based on K- p -> Lambda pi0
72
73 kaon->setType(ParticleTable::getPionType(iso));
74 nucleon->setType(Lambda);
75
76 G4double norm = KinematicsUtils::momentumInCM(kaon,nucleon);
77
78 kaon->setMomentum(mom_pion*norm);
79 nucleon->setMomentum(-mom_pion*norm);
80
81 kaon->adjustEnergyFromMomentum();
82 nucleon->adjustEnergyFromMomentum();
83
84 fs->addModifiedParticle(nucleon);
85 fs->addModifiedParticle(kaon);
86
87 }
88
89 ThreeVector NKbToLpiChannel::KaonMomentum(Particle const * const kaon, Particle const * const nucleon){
90
91 const G4double pLab = KinematicsUtils::momentumInLab(kaon,nucleon);
92
93 if(pLab < 435.) return Random::normVector(); // isotropic
94
95 G4double cos_theta = 1.;
96 G4double sin_theta = 0.;
97 const G4double cos_phi = std::cos(Random::shoot()*Math::twoPi);
98 const G4double sin_phi = std::sqrt(1-cos_phi*cos_phi);
99
100 const G4double x = kaon->getMomentum().getX();
101 const G4double y = kaon->getMomentum().getY();
102 const G4double z = kaon->getMomentum().getZ();
103
104 const G4double r = std::sqrt(x*x+y*y+z*z);
105 const G4double rho = std::sqrt(x*x+y*y);
106
107 if(pLab >= 1845.){
108 const G4double b = 12. * pLab/2375.; // correspond to the forward slope description at 2375 MeV/c in K- p elastic
109 cos_theta = std::log(Random::shoot()*(std::exp(b)-std::exp(-b))+std::exp(-b))/b;
110 sin_theta = std::sqrt(1-cos_theta*cos_theta);
111
112 }
113 else{
114 const G4double Legendre_coef[283][8] = {
115 {435,1.29381,0.87668,0.24531,0.00844,-0.0067,4e-05,0},
116 {440,1.28561,0.87491,0.24024,0.00529,-0.00609,3e-05,0},
117 {445,1.27743,0.87315,0.23522,0.00212,-0.00547,3e-05,0},
118 {450,1.26928,0.87139,0.2303,-0.00105,-0.00485,2e-05,0},
119 {455,1.2612,0.86965,0.22553,-0.00424,-0.0042,1e-05,0},
120 {460,1.25317,0.86792,0.22099,-0.00745,-0.00354,0,0},
121 {465,1.24512,0.86623,0.21678,-0.01069,-0.00284,0,0},
122 {470,1.23692,0.86461,0.21304,-0.01397,-0.00211,-1e-05,0},
123 {475,1.22848,0.86311,0.20991,-0.01729,-0.00132,-2e-05,0},
124 {480,1.21967,0.86172,0.20749,-0.02066,-0.00048,-3e-05,0},
125 {485,1.21035,0.86044,0.20577,-0.02407,0.00044,-3e-05,0},
126 {490,1.2004,0.85923,0.20473,-0.02748,0.00143,-4e-05,0},
127 {495,1.18967,0.85804,0.20434,-0.03087,0.00251,-5e-05,0},
128 {500,1.17802,0.85684,0.20455,-0.03423,0.00369,-6e-05,0},
129 {505,1.16533,0.8556,0.20535,-0.03753,0.00497,-7e-05,0},
130 {510,1.15145,0.85427,0.20669,-0.04075,0.00637,-8e-05,0},
131 {515,1.13624,0.85282,0.20855,-0.04386,0.0079,-9e-05,0},
132 {520,1.11972,0.85121,0.21086,-0.0468,0.00955,-1e-04,0},
133 {525,1.10208,0.84953,0.21379,-0.04947,0.01133,-0.00011,0},
134 {530,1.0835,0.84789,0.21759,-0.05178,0.01323,-0.00012,0},
135 {535,1.06417,0.84639,0.2225,-0.05361,0.01525,-0.00013,0},
136 {540,1.04423,0.8451,0.22862,-0.05493,0.01738,-0.00014,0},
137 {545,1.02374,0.844,0.23576,-0.05583,0.0196,-0.00015,0},
138 {550,1.00274,0.84308,0.2437,-0.05642,0.02189,-0.00016,0},
139 {555,0.98125,0.8423,0.25219,-0.05681,0.02423,-0.00017,0},
140 {560,0.95932,0.84165,0.26101,-0.0571,0.02659,-0.00017,0},
141 {565,0.93696,0.84108,0.26988,-0.05743,0.02895,-0.00018,0},
142 {570,0.91408,0.84053,0.27837,-0.05801,0.03126,-0.00018,0},
143 {575,0.89056,0.83994,0.28602,-0.05907,0.03345,-0.00018,0},
144 {580,0.86629,0.8393,0.29249,-0.06069,0.03543,-0.00017,0},
145 {585,0.84122,0.83869,0.29765,-0.06277,0.0371,-0.00016,0},
146 {590,0.81531,0.83818,0.30154,-0.06516,0.03842,-0.00015,0},
147 {595,0.78853,0.83781,0.30425,-0.06771,0.03943,-0.00012,0},
148 {600,0.76086,0.83762,0.30587,-0.07025,0.04015,-9e-05,1e-05},
149 {605,0.73242,0.83777,0.30661,-0.07259,0.0406,-6e-05,1e-05},
150 {610,0.7034,0.8385,0.30673,-0.07453,0.04081,-1e-05,1e-05},
151 {615,0.67402,0.84003,0.3065,-0.07588,0.0408,5e-05,1e-05},
152 {620,0.6445,0.8426,0.30618,-0.07642,0.0406,0.00012,1e-05},
153 {625,0.61503,0.84637,0.30581,-0.07599,0.04021,2e-04,2e-05},
154 {630,0.5858,0.85147,0.30525,-0.07444,0.03962,3e-04,2e-05},
155 {635,0.557,0.85801,0.30443,-0.07169,0.03881,0.00041,2e-05},
156 {640,0.5289,0.86611,0.30346,-0.06786,0.03772,0.00053,3e-05},
157 {645,0.50176,0.87586,0.30245,-0.06312,0.03626,0.00068,3e-05},
158 {650,0.47582,0.8874,0.30153,-0.05762,0.03438,0.00084,3e-05},
159 {655,0.45135,0.90083,0.30082,-0.05148,0.03201,0.00101,4e-05},
160 {660,0.42865,0.91629,0.30054,-0.0447,0.0292,0.0012,4e-05},
161 {665,0.40794,0.93393,0.30079,-0.03707,0.02606,0.00141,4e-05},
162 {670,0.38927,0.95387,0.30149,-0.0283,0.02276,0.00163,5e-05},
163 {675,0.37268,0.97623,0.30251,-0.01806,0.01946,0.00187,5e-05},
164 {680,0.35818,1.0011,0.30373,-0.00607,0.01635,0.0021,5e-05},
165 {685,0.34574,1.02858,0.30527,0.00796,0.0136,0.00234,5e-05},
166 {690,0.33527,1.05871,0.30741,0.02426,0.01138,0.00257,5e-05},
167 {695,0.32658,1.09138,0.31026,0.04292,0.00982,0.0028,5e-05},
168 {700,0.31939,1.12634,0.31382,0.06382,0.00895,0.003,5e-05},
169 {705,0.31344,1.16334,0.31815,0.08676,0.00879,0.00317,4e-05},
170 {710,0.3086,1.20215,0.32345,0.11138,0.00935,0.00328,4e-05},
171 {715,0.30471,1.24255,0.32994,0.13731,0.01061,0.00334,3e-05},
172 {720,0.30161,1.28419,0.33773,0.16421,0.01265,0.00332,2e-05},
173 {725,0.29909,1.32659,0.34682,0.19187,0.01564,0.0032,0},
174 {730,0.29695,1.36945,0.35725,0.22005,0.01956,0.00296,-2e-05},
175 {735,0.29502,1.41255,0.36915,0.24854,0.0243,0.00258,-4e-05},
176 {740,0.29311,1.45571,0.38262,0.27709,0.02974,0.00202,-7e-05},
177 {745,0.29109,1.49871,0.39781,0.3055,0.03575,0.00127,-0.00011},
178 {750,0.28898,1.54138,0.41503,0.33366,0.0422,0.00029,-0.00015},
179 {755,0.28683,1.58345,0.43457,0.36147,0.04894,-0.00094,-0.00019},
180 {760,0.28451,1.6243,0.45625,0.38871,0.05591,-0.00244,-0.00025},
181 {765,0.28175,1.66337,0.47979,0.41518,0.06308,-0.00425,-0.00031},
182 {770,0.27835,1.70025,0.50482,0.44067,0.07042,-0.00639,-0.00037},
183 {775,0.27417,1.73456,0.53095,0.46502,0.07785,-0.00889,-0.00044},
184 {780,0.26902,1.76599,0.55773,0.48808,0.08531,-0.01176,-0.00052},
185 {785,0.26276,1.79435,0.58464,0.50969,0.09264,-0.01503,-6e-04},
186 {790,0.2552,1.81947,0.61114,0.52972,0.09971,-0.0187,-0.00068},
187 {795,0.24617,1.84117,0.6367,0.54802,0.10638,-0.02279,-0.00077},
188 {800,0.23554,1.85933,0.66076,0.56451,0.11244,-0.0273,-0.00086},
189 {805,0.22322,1.87397,0.68275,0.57917,0.11765,-0.03222,-0.00095},
190 {810,0.20914,1.8851,0.70209,0.592,0.12173,-0.03752,-0.00104},
191 {815,0.1932,1.89272,0.7182,0.60299,0.12443,-0.04319,-0.00113},
192 {820,0.17533,1.89686,0.73052,0.61214,0.1255,-0.04922,-0.00121},
193 {825,0.15565,1.89756,0.73887,0.61947,0.12491,-0.05557,-0.00128},
194 {830,0.13439,1.89493,0.74336,0.62504,0.12283,-0.06219,-0.00134},
195 {835,0.11178,1.88905,0.74408,0.62889,0.11942,-0.06905,-0.00139},
196 {840,0.08807,1.88001,0.74115,0.63106,0.11484,-0.07609,-0.00143},
197 {845,0.0635,1.86791,0.73467,0.63161,0.10925,-0.08327,-0.00144},
198 {850,0.03831,1.85283,0.72475,0.63058,0.10282,-0.09055,-0.00144},
199 {855,0.01278,1.83488,0.7116,0.62802,0.09574,-0.09784,-0.00142},
200 {860,-0.01274,1.81419,0.69565,0.624,0.0883,-0.105,-0.00136},
201 {865,-0.03789,1.79087,0.67736,0.61858,0.08079,-0.11186,-0.00127},
202 {870,-0.06232,1.76507,0.65719,0.61182,0.07351,-0.11824,-0.00115},
203 {875,-0.08566,1.73691,0.63561,0.60379,0.06676,-0.12399,-0.00099},
204 {880,-0.10755,1.70652,0.61308,0.59454,0.06082,-0.12894,-0.00077},
205 {885,-0.12763,1.67402,0.59005,0.58414,0.05598,-0.13293,-0.00051},
206 {890,-0.1458,1.63947,0.56679,0.57245,0.05242,-0.13584,-0.00019},
207 {895,-0.16195,1.60305,0.54339,0.5594,0.05012,-0.13773,0.00019},
208 {900,-0.17602,1.56493,0.51984,0.54502,0.04902,-0.13871,0.00064},
209 {905,-0.18798,1.52525,0.4961,0.52929,0.049,-0.13888,0.00117},
210 {910,-0.19778,1.48414,0.47216,0.51223,0.04998,-0.13838,0.00177},
211 {915,-0.20538,1.44176,0.44798,0.49384,0.05185,-0.13732,0.00247},
212 {920,-0.21074,1.39824,0.42355,0.47412,0.05451,-0.13581,0.00325},
213 {925,-0.21388,1.35383,0.39886,0.45319,0.05787,-0.13396,0.00414},
214 {930,-0.21484,1.30892,0.37401,0.4313,0.06178,-0.13183,0.00512},
215 {935,-0.21369,1.2639,0.34905,0.40872,0.06613,-0.1295,0.0062},
216 {940,-0.2105,1.21913,0.32406,0.38569,0.07078,-0.12704,0.00739},
217 {945,-0.20536,1.17491,0.29905,0.36232,0.07556,-0.12458,0.00869},
218 {950,-0.19836,1.13145,0.27401,0.33866,0.08027,-0.12224,0.01009},
219 {955,-0.18958,1.08898,0.24895,0.31474,0.08473,-0.12018,0.01158},
220 {960,-0.17912,1.04775,0.22386,0.2906,0.08874,-0.11852,0.01318},
221 {965,-0.16707,1.00794,0.19872,0.26627,0.09217,-0.11741,0.01487},
222 {970,-0.15345,0.96965,0.17358,0.24171,0.09504,-0.11702,0.01662},
223 {975,-0.1383,0.93292,0.14845,0.21687,0.09744,-0.1175,0.01841},
224 {980,-0.12165,0.89782,0.12337,0.19178,0.09944,-0.11883,0.0202},
225 {985,-0.10352,0.86439,0.09835,0.16649,0.10108,-0.12091,0.02197},
226 {990,-0.08393,0.83266,0.07345,0.14106,0.1024,-0.12362,0.02367},
227 {995,-0.06289,0.80271,0.04868,0.11554,0.10346,-0.12687,0.02527},
228 {1000,-0.04044,0.77456,0.02407,0.09,0.10431,-0.13055,0.02675},
229 {1005,-0.0166,0.74824,-3e-04,0.06448,0.10496,-0.13453,0.02801},
230 {1010,0.00852,0.72368,-0.02409,0.03902,0.1053,-0.13858,0.02872},
231 {1015,0.03479,0.70081,-0.0472,0.01366,0.10528,-0.14262,0.02878},
232 {1020,0.06212,0.67958,-0.06962,-0.01158,0.10492,-0.14664,0.02821},
233 {1025,0.0904,0.65992,-0.09133,-0.03664,0.10421,-0.15064,0.02703},
234 {1030,0.11951,0.64179,-0.11234,-0.0615,0.10315,-0.15461,0.02526},
235 {1035,0.14935,0.62513,-0.13261,-0.08611,0.10174,-0.15854,0.02291},
236 {1040,0.1798,0.60989,-0.15215,-0.11044,0.09999,-0.16241,0.02},
237 {1045,0.21078,0.59602,-0.17094,-0.13445,0.0979,-0.16622,0.01656},
238 {1050,0.24215,0.58346,-0.18897,-0.15811,0.09547,-0.16996,0.01259},
239 {1055,0.27383,0.57216,-0.20623,-0.18137,0.0927,-0.17363,0.00812},
240 {1060,0.30569,0.56207,-0.2227,-0.20421,0.08959,-0.1772,0.00316},
241 {1065,0.33764,0.55313,-0.23838,-0.22657,0.08615,-0.18068,-0.00227},
242 {1070,0.36956,0.54529,-0.25326,-0.24843,0.08237,-0.18405,-0.00815},
243 {1075,0.40134,0.53849,-0.26732,-0.26974,0.07826,-0.1873,-0.01447},
244 {1080,0.43288,0.53269,-0.28055,-0.29047,0.07382,-0.19043,-0.02119},
245 {1085,0.46407,0.52783,-0.29295,-0.31059,0.06905,-0.19342,-0.02832},
246 {1090,0.49481,0.52386,-0.30452,-0.33004,0.06396,-0.19627,-0.03582},
247 {1095,0.52503,0.52073,-0.31538,-0.34881,0.05856,-0.1989,-0.04359},
248 {1100,0.55465,0.51841,-0.3257,-0.36683,0.05286,-0.20125,-0.05154},
249 {1105,0.58361,0.51686,-0.3356,-0.38409,0.0469,-0.20327,-0.05957},
250 {1110,0.61184,0.51604,-0.34525,-0.40052,0.04068,-0.20488,-0.06757},
251 {1115,0.63927,0.5159,-0.35478,-0.4161,0.03422,-0.20602,-0.07545},
252 {1120,0.66583,0.51641,-0.36435,-0.43079,0.02755,-0.20662,-0.08311},
253 {1125,0.69146,0.51753,-0.37411,-0.44454,0.02068,-0.20662,-0.09043},
254 {1130,0.71609,0.51922,-0.38418,-0.45732,0.01364,-0.20594,-0.09736},
255 {1135,0.73967,0.5214,-0.3946,-0.46913,0.00649,-0.20451,-0.10389},
256 {1140,0.76215,0.52403,-0.40542,-0.47994,-0.00069,-0.20222,-0.11009},
257 {1145,0.78348,0.52703,-0.41665,-0.48977,-0.00782,-0.19898,-0.11598},
258 {1150,0.80362,0.53033,-0.42833,-0.49859,-0.01484,-0.19469,-0.1216},
259 {1155,0.82251,0.53388,-0.44048,-0.50642,-0.02166,-0.18927,-0.12701},
260 {1160,0.8401,0.53761,-0.45313,-0.51323,-0.02821,-0.18261,-0.13223},
261 {1165,0.85635,0.54144,-0.46631,-0.51903,-0.03442,-0.17462,-0.13732},
262 {1170,0.87122,0.54533,-0.48004,-0.52382,-0.04023,-0.16525,-0.14231},
263 {1175,0.88472,0.54924,-0.49427,-0.52765,-0.0457,-0.15454,-0.1472},
264 {1180,0.89689,0.55314,-0.50896,-0.5306,-0.0509,-0.1426,-0.15202},
265 {1185,0.90774,0.55701,-0.52406,-0.53273,-0.05591,-0.12951,-0.15675},
266 {1190,0.91731,0.56082,-0.53952,-0.53411,-0.06079,-0.11537,-0.16142},
267 {1195,0.92561,0.56455,-0.55528,-0.53481,-0.06563,-0.10027,-0.16601},
268 {1200,0.93269,0.56816,-0.57129,-0.53488,-0.07049,-0.08429,-0.17056},
269 {1205,0.93856,0.57163,-0.58751,-0.53441,-0.07546,-0.06754,-0.17505},
270 {1210,0.94325,0.57494,-0.60389,-0.53345,-0.0806,-0.0501,-0.17949},
271 {1215,0.94679,0.57806,-0.62037,-0.53208,-0.086,-0.03207,-0.1839},
272 {1220,0.94921,0.58096,-0.6369,-0.53036,-0.09172,-0.01354,-0.18827},
273 {1225,0.95052,0.58361,-0.65343,-0.52836,-0.09784,0.00541,-0.19262},
274 {1230,0.95077,0.58599,-0.66992,-0.52614,-0.10444,0.02468,-0.19694},
275 {1235,0.94997,0.58808,-0.68631,-0.52378,-0.11159,0.04417,-0.20126},
276 {1240,0.94816,0.58984,-0.70255,-0.52133,-0.11936,0.06379,-0.20557},
277 {1245,0.94535,0.59125,-0.71859,-0.51888,-0.12783,0.08346,-0.20987},
278 {1250,0.94157,0.59229,-0.73437,-0.51648,-0.13706,0.10305,-0.21416},
279 {1255,0.93683,0.59293,-0.74976,-0.51419,-0.1471,0.12231,-0.21828},
280 {1260,0.93113,0.59315,-0.76463,-0.51207,-0.15797,0.14096,-0.22206},
281 {1265,0.92447,0.59294,-0.77883,-0.51019,-0.1697,0.15872,-0.22534},
282 {1270,0.91689,0.59232,-0.79218,-0.50845,-0.18223,0.17547,-0.22804},
283 {1275,0.90848,0.5914,-0.80445,-0.5066,-0.19545,0.19131,-0.23018},
284 {1280,0.89933,0.59026,-0.81539,-0.50438,-0.20919,0.20634,-0.2318},
285 {1285,0.88952,0.58902,-0.82477,-0.50151,-0.22331,0.22066,-0.23293},
286 {1290,0.87914,0.58776,-0.83234,-0.49778,-0.23768,0.23437,-0.23361},
287 {1295,0.86827,0.58655,-0.83797,-0.49307,-0.25217,0.24756,-0.23388},
288 {1300,0.85696,0.58546,-0.84148,-0.48733,-0.26668,0.26028,-0.23375},
289 {1305,0.8453,0.58454,-0.84274,-0.4805,-0.2811,0.27262,-0.23327},
290 {1310,0.83334,0.58385,-0.84159,-0.47251,-0.29532,0.28466,-0.23245},
291 {1315,0.82116,0.58345,-0.8379,-0.46331,-0.30922,0.29645,-0.23134},
292 {1320,0.80881,0.58341,-0.8315,-0.45284,-0.32269,0.30807,-0.22995},
293 {1325,0.79636,0.58376,-0.82231,-0.44106,-0.33567,0.31957,-0.22833},
294 {1330,0.78383,0.58452,-0.81043,-0.42806,-0.34819,0.33087,-0.22653},
295 {1335,0.77121,0.58565,-0.79603,-0.41395,-0.36036,0.34185,-0.22459},
296 {1340,0.75852,0.58715,-0.77928,-0.39886,-0.37227,0.3524,-0.22258},
297 {1345,0.74577,0.589,-0.76033,-0.38289,-0.38399,0.36242,-0.22055},
298 {1350,0.73295,0.59118,-0.73935,-0.36617,-0.39563,0.37178,-0.21857},
299 {1355,0.72008,0.59368,-0.71651,-0.3488,-0.40726,0.38039,-0.21667},
300 {1360,0.70715,0.59646,-0.69194,-0.3309,-0.41893,0.38818,-0.2149},
301 {1365,0.6941,0.59943,-0.66566,-0.31253,-0.43043,0.39531,-0.21314},
302 {1370,0.68084,0.60248,-0.63764,-0.29373,-0.44151,0.40198,-0.21126},
303 {1375,0.66737,0.60553,-0.60795,-0.27456,-0.45201,0.40833,-0.20918},
304 {1380,0.65371,0.60856,-0.57672,-0.25501,-0.46189,0.41438,-0.2069},
305 {1385,0.63991,0.61158,-0.5441,-0.2351,-0.47112,0.42015,-0.20442},
306 {1390,0.62603,0.61456,-0.51023,-0.21484,-0.47966,0.42567,-0.20173},
307 {1395,0.6121,0.6175,-0.47525,-0.19424,-0.48749,0.43095,-0.19884},
308 {1400,0.59816,0.62039,-0.43932,-0.17332,-0.49457,0.43601,-0.19575},
309 {1405,0.58428,0.62321,-0.40257,-0.15208,-0.50088,0.44087,-0.19246},
310 {1410,0.57048,0.62596,-0.36515,-0.13053,-0.50637,0.44556,-0.18897},
311 {1415,0.55682,0.62862,-0.32721,-0.1087,-0.51102,0.4501,-0.18528},
312 {1420,0.54335,0.63118,-0.28889,-0.08658,-0.5148,0.4545,-0.18139},
313 {1425,0.5301,0.63363,-0.25033,-0.06419,-0.51767,0.45878,-0.1773},
314 {1430,0.51712,0.63597,-0.21167,-0.04154,-0.5196,0.46298,-0.173},
315 {1435,0.50446,0.63817,-0.17307,-0.01864,-0.52057,0.4671,-0.16851},
316 {1440,0.49217,0.64024,-0.13467,0.0045,-0.52053,0.47116,-0.16382},
317 {1445,0.48028,0.64215,-0.09661,0.02786,-0.51946,0.4752,-0.15892},
318 {1450,0.46885,0.64391,-0.05904,0.05144,-0.51733,0.47922,-0.15383},
319 {1455,0.45792,0.64548,-0.02209,0.07522,-0.51411,0.48325,-0.14854},
320 {1460,0.44754,0.64688,0.01408,0.09919,-0.50975,0.48732,-0.14305},
321 {1465,0.43774,0.64808,0.04933,0.12335,-0.50424,0.49143,-0.13736},
322 {1470,0.42855,0.64906,0.08361,0.14764,-0.49758,0.49559,-0.13154},
323 {1475,0.41997,0.6498,0.11691,0.17201,-0.4898,0.49979,-0.12566},
324 {1480,0.41199,0.65025,0.14924,0.19638,-0.48094,0.50402,-0.11983},
325 {1485,0.4046,0.65039,0.1806,0.2207,-0.47103,0.50827,-0.11413},
326 {1490,0.39781,0.65018,0.211,0.2449,-0.4601,0.51253,-0.10866},
327 {1495,0.39161,0.6496,0.24043,0.26892,-0.44819,0.51679,-0.10352},
328 {1500,0.38599,0.64862,0.26889,0.29269,-0.43534,0.52104,-0.09879},
329 {1505,0.38096,0.64719,0.2964,0.31615,-0.42158,0.52526,-0.09457},
330 {1510,0.37651,0.6453,0.32295,0.33923,-0.40693,0.52944,-0.09095},
331 {1515,0.37264,0.6429,0.34854,0.36186,-0.39144,0.53359,-0.08803},
332 {1520,0.36933,0.63998,0.37315,0.38401,-0.37512,0.53773,-0.08583},
333 {1525,0.36653,0.63654,0.39667,0.40564,-0.35794,0.54203,-0.08422},
334 {1530,0.3642,0.63257,0.41901,0.42676,-0.33984,0.54668,-0.08304},
335 {1535,0.36228,0.62809,0.44005,0.44733,-0.3208,0.55185,-0.08216},
336 {1540,0.36072,0.62309,0.45968,0.46736,-0.30076,0.55773,-0.08141},
337 {1545,0.35948,0.61757,0.47781,0.48683,-0.27969,0.5645,-0.08065},
338 {1550,0.35851,0.61154,0.49434,0.50572,-0.25756,0.57231,-0.07974},
339 {1555,0.35781,0.60502,0.50929,0.52405,-0.23439,0.58116,-0.07861},
340 {1560,0.35738,0.59804,0.52273,0.54186,-0.21023,0.59098,-0.07724},
341 {1565,0.35725,0.59063,0.53473,0.55915,-0.18514,0.60169,-0.07558},
342 {1570,0.35742,0.58282,0.54538,0.57597,-0.15918,0.61322,-0.07362},
343 {1575,0.35792,0.57463,0.55473,0.59234,-0.13239,0.62549,-0.07132},
344 {1580,0.35874,0.5661,0.56287,0.60828,-0.10484,0.63843,-0.06864},
345 {1585,0.35991,0.55725,0.56986,0.62383,-0.07657,0.65196,-0.06557},
346 {1590,0.36144,0.54811,0.57578,0.63899,-0.04764,0.66601,-0.06206},
347 {1595,0.36335,0.53871,0.5807,0.65382,-0.0181,0.6805,-0.05809},
348 {1600,0.36564,0.52909,0.5847,0.66832,0.01199,0.69536,-0.05363},
349 {1605,0.36833,0.51926,0.58783,0.68253,0.04258,0.7105,-0.04864},
350 {1610,0.37143,0.50926,0.59019,0.69647,0.07361,0.72586,-0.0431},
351 {1615,0.37496,0.49911,0.59183,0.71017,0.10503,0.74137,-0.03697},
352 {1620,0.37894,0.48885,0.59284,0.72365,0.13679,0.75693,-0.03023},
353 {1625,0.38336,0.47851,0.59329,0.73695,0.16883,0.77249,-0.02284},
354 {1630,0.38826,0.46811,0.59324,0.75008,0.2011,0.78797,-0.01477},
355 {1635,0.39363,0.45767,0.59277,0.76308,0.23354,0.80328,-0.00599},
356 {1640,0.39951,0.44724,0.59195,0.77597,0.26611,0.81836,0.00352},
357 {1645,0.40588,0.43684,0.59086,0.78878,0.29874,0.83312,0.01381},
358 {1650,0.41278,0.4265,0.58956,0.80153,0.33138,0.84751,0.0249},
359 {1655,0.42022,0.41624,0.58813,0.81425,0.36397,0.86143,0.03681},
360 {1660,0.4282,0.40614,0.58665,0.82695,0.39646,0.87483,0.04954},
361 {1665,0.43672,0.39632,0.58517,0.83966,0.42876,0.88768,0.06298},
362 {1670,0.44579,0.38689,0.58378,0.85238,0.46079,0.89994,0.07703},
363 {1675,0.45542,0.37799,0.58253,0.86511,0.49245,0.91158,0.09161},
364 {1680,0.46559,0.36972,0.5815,0.87787,0.52367,0.92256,0.1066},
365 {1685,0.47631,0.36223,0.58075,0.89067,0.55437,0.93285,0.12193},
366 {1690,0.48759,0.35562,0.58036,0.90351,0.58446,0.94242,0.13748},
367 {1695,0.49943,0.35002,0.58039,0.91642,0.61386,0.95122,0.15316},
368 {1700,0.51182,0.34555,0.58092,0.92938,0.64248,0.95923,0.16887},
369 {1705,0.52478,0.34234,0.58201,0.94243,0.67025,0.96642,0.18453},
370 {1710,0.5383,0.34049,0.58374,0.95556,0.69713,0.97277,0.20007},
371 {1715,0.55239,0.34002,0.58616,0.96883,0.72336,0.97841,0.21569},
372 {1720,0.56707,0.34092,0.58936,0.98229,0.74924,0.98349,0.23159},
373 {1725,0.58235,0.34319,0.59342,0.99599,0.77504,0.98814,0.24802},
374 {1730,0.59825,0.34683,0.59839,1.00998,0.80106,0.99252,0.26518},
375 {1735,0.61478,0.35185,0.60435,1.02432,0.82759,0.99678,0.28332},
376 {1740,0.63195,0.35823,0.61138,1.03905,0.85493,1.00106,0.30266},
377 {1745,0.64977,0.36598,0.61954,1.05422,0.88333,1.00549,0.32339},
378 {1750,0.66822,0.37504,0.6288,1.06985,0.91282,1.01013,0.34556},
379 {1755,0.68725,0.38533,0.6391,1.08592,0.94333,1.01496,0.36913},
380 {1760,0.70685,0.39679,0.65037,1.10245,0.9748,1.01998,0.39406},
381 {1765,0.72695,0.40934,0.66257,1.11943,1.00714,1.02519,0.42031},
382 {1770,0.74754,0.42291,0.67563,1.13684,1.0403,1.03058,0.44784},
383 {1775,0.76857,0.43743,0.68947,1.15471,1.07419,1.03615,0.47661},
384 {1780,0.79,0.45282,0.70406,1.17301,1.10875,1.04191,0.50658},
385 {1785,0.81179,0.46902,0.71931,1.19174,1.14391,1.04784,0.53771},
386 {1790,0.83392,0.48594,0.73518,1.21092,1.17959,1.05394,0.56997},
387 {1795,0.85634,0.50352,0.75159,1.23052,1.21573,1.06022,0.60332},
388 {1800,0.87901,0.52169,0.76849,1.25055,1.25225,1.06666,0.63771},
389 {1805,0.9019,0.54038,0.78582,1.27101,1.28908,1.07326,0.6731},
390 {1810,0.92496,0.5595,0.80351,1.2919,1.32615,1.08003,0.70946},
391 {1815,0.94817,0.57899,0.82151,1.31321,1.36339,1.08695,0.74675},
392 {1820,0.97149,0.59878,0.83975,1.33493,1.40073,1.09403,0.78493},
393 {1825,0.99487,0.61879,0.85816,1.35707,1.4381,1.10125,0.82396},
394 {1830,1.01828,0.63895,0.8767,1.37963,1.47543,1.10863,0.86381},
395 {1835,1.04168,0.65919,0.89529,1.4026,1.51264,1.11616,0.90442},
396 {1840,1.06504,0.67943,0.91388,1.42598,1.54966,1.12382,0.94576},
397 {1845,1.08832,0.69962,0.93241,1.44973,1.58645,1.13162,0.98775}};
398
399 const G4int coef_ener = G4int((pLab-Legendre_coef[0][0])/5);
400 const G4double sup_ener = pLab/5. - coef_ener -Legendre_coef[0][0]/5;
401
402 // assert(pLab >= Legendre_coef[coef_ener][0] && pLab < Legendre_coef[coef_ener+1][0]);
403
404 // Legendre coefficient normalized
405 const G4double A0 = 1.;
406 const G4double A1 = (1-sup_ener)*Legendre_coef[coef_ener][1] + sup_ener*Legendre_coef[coef_ener+1][1];
407 const G4double A2 = (1-sup_ener)*Legendre_coef[coef_ener][2] + sup_ener*Legendre_coef[coef_ener+1][2];
408 const G4double A3 = (1-sup_ener)*Legendre_coef[coef_ener][3] + sup_ener*Legendre_coef[coef_ener+1][3];
409 const G4double A4 = (1-sup_ener)*Legendre_coef[coef_ener][4] + sup_ener*Legendre_coef[coef_ener+1][4];
410 const G4double A5 = (1-sup_ener)*Legendre_coef[coef_ener][5] + sup_ener*Legendre_coef[coef_ener+1][5];
411 const G4double A6 = (1-sup_ener)*Legendre_coef[coef_ener][6] + sup_ener*Legendre_coef[coef_ener+1][6];
412 const G4double A7 = (1-sup_ener)*Legendre_coef[coef_ener][7] + sup_ener*Legendre_coef[coef_ener+1][7];
413
414 // Theoritical max if all Ai > 0 (often the case)
415 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);
416
417 G4bool success = false;
418 G4int maxloop = 0;
419
420 while(!success && maxloop < 1000){
421
422 cos_theta = Random::shoot()*2-1.; // not optimized
423
424 // Legendre Polynomial
425 G4double P0 = A0;
426 G4double P1 = A1*cos_theta;
427 G4double P2 = A2/2.*(3*std::pow(cos_theta,2)-1);
428 G4double P3 = A3/2.*(5*std::pow(cos_theta,3)-3*cos_theta);
429 G4double P4 = A4/8.*(35*std::pow(cos_theta,4)-30*std::pow(cos_theta,2)+3);
430 G4double P5 = A5/8.*(63*std::pow(cos_theta,5)-70*std::pow(cos_theta,3)+15*cos_theta);
431 G4double P6 = A6/16.*(231*std::pow(cos_theta,6)-315*std::pow(cos_theta,4)+105*std::pow(cos_theta,2)-5);
432 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);
433
434 G4double P = (P0 + P1 + P2 + P3 + P4 + P5 + P6 + P7)/2; // /2 for the normalisation
435
436 if(Random::shoot()*A < P) success = true;
437 maxloop +=1 ;
438 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
439 }
440 sin_theta = std::sqrt(1-cos_theta*cos_theta);
441 }
442
443 if(rho == 0) return ThreeVector(sin_theta*cos_phi,sin_theta*sin_phi,cos_theta);
444 // Rotation in the direction of the incident kaon
445 const G4double px = x/r*cos_theta - y/rho*sin_theta*cos_phi + z/r*x/rho*sin_theta*sin_phi;
446 const G4double py = y/r*cos_theta + x/rho*sin_theta*cos_phi + z/r*y/rho*sin_theta*sin_phi;
447 const G4double pz = z/r*cos_theta - rho/r*sin_theta*sin_phi;
448
449 return ThreeVector(px,py,pz);
450 }
451
452 }
453