| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // INCL++ intra-nuclear cascade model 26 // INCL++ intra-nuclear cascade model
27 // Alain Boudard, CEA-Saclay, France << 27 // Pekka Kaitaniemi, CEA and Helsinki Institute of Physics
28 // Joseph Cugnon, University of Liege, Belgium << 28 // Davide Mancusi, CEA
29 // Jean-Christophe David, CEA-Saclay, France << 29 // Alain Boudard, CEA
30 // Pekka Kaitaniemi, CEA-Saclay, France, and H << 30 // Sylvie Leray, CEA
31 // Sylvie Leray, CEA-Saclay, France << 31 // Joseph Cugnon, University of Liege
32 // Davide Mancusi, CEA-Saclay, France <<
33 // 32 //
34 #define INCLXX_IN_GEANT4_MODE 1 33 #define INCLXX_IN_GEANT4_MODE 1
35 34
36 #include "globals.hh" 35 #include "globals.hh"
37 36
38 /** \file G4INCLEventInfo.hh 37 /** \file G4INCLEventInfo.hh
39 * \brief Simple container for output of event 38 * \brief Simple container for output of event results.
40 * 39 *
41 * Contains the results of an INCL cascade. 40 * Contains the results of an INCL cascade.
42 * 41 *
43 * \date 21 January 2011 42 * \date 21 January 2011
44 * \author Davide Mancusi 43 * \author Davide Mancusi
45 */ 44 */
46 45
47 #ifndef G4INCLEVENTINFO_HH_HH << 46 #ifndef G4INCLEVENTINFO_HH
48 #define G4INCLEVENTINFO_HH_HH 1 << 47 #define G4INCLEVENTINFO_HH 1
49 48
50 #include "G4INCLParticleType.hh" 49 #include "G4INCLParticleType.hh"
51 #ifdef INCL_ROOT_USE 50 #ifdef INCL_ROOT_USE
52 #include <Rtypes.h> 51 #include <Rtypes.h>
53 #endif 52 #endif
54 #include <string> 53 #include <string>
55 #include <vector> 54 #include <vector>
56 #include <algorithm> 55 #include <algorithm>
57 56
58 namespace G4INCL { 57 namespace G4INCL {
59 #ifndef INCL_ROOT_USE 58 #ifndef INCL_ROOT_USE
60 typedef G4int Int_t; 59 typedef G4int Int_t;
61 typedef short Short_t; 60 typedef short Short_t;
62 typedef G4float Float_t; 61 typedef G4float Float_t;
63 typedef G4double Double_t; 62 typedef G4double Double_t;
64 typedef G4bool Bool_t; 63 typedef G4bool Bool_t;
65 #endif 64 #endif
66 65
67 struct EventInfo { 66 struct EventInfo {
68 EventInfo() : 67 EventInfo() :
69 nParticles(0), << 68 projectileType(UnknownParticle),
70 event(0), << 69 At(0), Zt(0), Ap(0), Zp(0),
71 eventBias((Float_t)0.0), << 70 Ep(0.),
72 nRemnants(0), << 71 impactParameter(0.0), nCollisions(0), stoppingTime(0.0),
73 projectileType(0), << 72 EBalance(0.0), pLongBalance(0.0), pTransBalance(0.0),
74 At(0), << 73 nCascadeParticles(0), nRemnants(0), nParticles(0),
75 Zt(0), << 74 transparent(true),
76 St(0), <<
77 Ap(0), <<
78 Zp(0), <<
79 Sp(0), <<
80 Ep((Float_t)0.0), <<
81 impactParameter((Float_t)0.0), <<
82 nCollisions(0), <<
83 stoppingTime((Float_t)0.0), <<
84 EBalance((Float_t)0.0), <<
85 firstEBalance((Float_t)0.0), <<
86 pLongBalance((Float_t)0.0), <<
87 pTransBalance((Float_t)0.0), <<
88 nCascadeParticles(0), <<
89 transparent(false), <<
90 annihilationP(false), <<
91 annihilationN(false), <<
92 forcedCompoundNucleus(false), 75 forcedCompoundNucleus(false),
93 nucleonAbsorption(false), << 76 nucleonAbsorption(false), pionAbsorption(false), nDecays(0),
94 pionAbsorption(false), << 77 nBlockedCollisions(0), nBlockedDecays(0),
95 nDecays(0), << 78 effectiveImpactParameter(0.0),
96 nSrcCollisions(0), <<
97 nSrcPairs(0), <<
98 nBlockedCollisions(0), <<
99 nBlockedDecays(0), <<
100 effectiveImpactParameter((Float_t)0.0) <<
101 deltasInside(false), 79 deltasInside(false),
102 sigmasInside(false), <<
103 kaonsInside(false), <<
104 antikaonsInside(false), <<
105 lambdasInside(false), <<
106 forcedDeltasInside(false), 80 forcedDeltasInside(false),
107 forcedDeltasOutside(false), 81 forcedDeltasOutside(false),
108 forcedPionResonancesOutside(false), <<
109 absorbedStrangeParticle(false), <<
110 forcedSigmaOutside(false), <<
111 forcedStrangeInside(false), <<
112 emitLambda(0), <<
113 emitKaon(false), <<
114 clusterDecay(false), 82 clusterDecay(false),
115 firstCollisionTime((Float_t)0.0), << 83 firstCollisionTime(0.),
116 firstCollisionXSec((Float_t)0.0), << 84 firstCollisionXSec(0.),
117 firstCollisionSpectatorPosition((Float <<
118 firstCollisionSpectatorMomentum((Float <<
119 firstCollisionIsElastic(false), <<
120 nReflectionAvatars(0), 85 nReflectionAvatars(0),
121 nCollisionAvatars(0), 86 nCollisionAvatars(0),
122 nDecayAvatars(0), 87 nDecayAvatars(0),
123 nUnmergedSpectators(0), << 88 nUnmergedSpectators(0)
124 nEnergyViolationInteraction(0) <<
125 <<
126 { 89 {
>> 90 std::fill_n(ARem, maxSizeRemnants, 0);
>> 91 std::fill_n(ZRem, maxSizeRemnants, 0);
>> 92 std::fill_n(EStarRem, maxSizeRemnants, ((Float_t)0.));
>> 93 std::fill_n(JRem, maxSizeRemnants, ((Float_t)0.));
>> 94 std::fill_n(EKinRem, maxSizeRemnants, ((Float_t)0.));
>> 95 std::fill_n(pxRem, maxSizeRemnants, ((Float_t)0.));
>> 96 std::fill_n(pyRem, maxSizeRemnants, ((Float_t)0.));
>> 97 std::fill_n(pzRem, maxSizeRemnants, ((Float_t)0.));
>> 98 std::fill_n(thetaRem, maxSizeRemnants, ((Float_t)0.));
>> 99 std::fill_n(phiRem, maxSizeRemnants, ((Float_t)0.));
>> 100 std::fill_n(jxRem, maxSizeRemnants, ((Float_t)0.));
>> 101 std::fill_n(jyRem, maxSizeRemnants, ((Float_t)0.));
>> 102 std::fill_n(jzRem, maxSizeRemnants, ((Float_t)0.));
>> 103
127 std::fill_n(A, maxSizeParticles, 0); 104 std::fill_n(A, maxSizeParticles, 0);
128 std::fill_n(Z, maxSizeParticles, 0); 105 std::fill_n(Z, maxSizeParticles, 0);
129 std::fill_n(S, maxSizeParticles, 0); << 106 std::fill_n(emissionTime, maxSizeParticles, ((Float_t)0.));
130 std::fill_n(J, maxSizeParticles, 0); << 107 std::fill_n(EKin, maxSizeParticles, ((Float_t)0.));
131 std::fill_n(PDGCode, maxSizeParticles, << 108 std::fill_n(px, maxSizeParticles, ((Float_t)0.));
132 std::fill_n(ParticleBias, maxSizeParti << 109 std::fill_n(py, maxSizeParticles, ((Float_t)0.));
133 std::fill_n(EKin, maxSizeParticles, (F << 110 std::fill_n(pz, maxSizeParticles, ((Float_t)0.));
134 std::fill_n(px, maxSizeParticles, (Flo << 111 std::fill_n(theta, maxSizeParticles, ((Float_t)0.));
135 std::fill_n(py, maxSizeParticles, (Flo << 112 std::fill_n(phi, maxSizeParticles, ((Float_t)0.));
136 std::fill_n(pz, maxSizeParticles, (Flo <<
137 std::fill_n(theta, maxSizeParticles, ( <<
138 std::fill_n(phi, maxSizeParticles, (Fl <<
139 std::fill_n(origin, maxSizeParticles, 113 std::fill_n(origin, maxSizeParticles, 0);
140 std::fill_n(parentResonancePDGCode, ma << 114 };
141 std::fill_n(parentResonanceID, maxSize <<
142 std::fill_n(emissionTime, maxSizeParti <<
143 std::fill_n(ARem, maxSizeRemnants, 0); <<
144 std::fill_n(ZRem, maxSizeRemnants, 0); <<
145 std::fill_n(SRem, maxSizeRemnants, 0); <<
146 std::fill_n(EStarRem, maxSizeRemnants, <<
147 std::fill_n(JRem, maxSizeRemnants, (Fl <<
148 std::fill_n(EKinRem, maxSizeRemnants, <<
149 std::fill_n(pxRem, maxSizeRemnants, (F <<
150 std::fill_n(pyRem, maxSizeRemnants, (F <<
151 std::fill_n(pzRem, maxSizeRemnants, (F <<
152 std::fill_n(thetaRem, maxSizeRemnants, <<
153 std::fill_n(phiRem, maxSizeRemnants, ( <<
154 std::fill_n(jxRem, maxSizeRemnants, (F <<
155 std::fill_n(jyRem, maxSizeRemnants, (F <<
156 std::fill_n(jzRem, maxSizeRemnants, (F <<
157 std::fill_n(EKinPrime, maxSizeParticle <<
158 std::fill_n(pzPrime, maxSizeParticles, <<
159 std::fill_n(thetaPrime, maxSizeParticl <<
160 } <<
161 115
162 /** \brief Number of the event */ 116 /** \brief Number of the event */
163 static G4ThreadLocal Int_t eventNumber; << 117 static Int_t eventNumber;
>> 118 /** \brief Protjectile particle type */
>> 119 ParticleType projectileType;
164 120
165 /** \brief Maximum array size for remnan <<
166 static const Short_t maxSizeRemnants = 1 <<
167 <<
168 /** \brief Maximum array size for emitte <<
169 static const Short_t maxSizeParticles = <<
170 <<
171 /** \brief Number of particles in the fi <<
172 Short_t nParticles; <<
173 /** \brief Sequential number of the even <<
174 Int_t event; <<
175 /** \brief Particle mass number */ <<
176 Short_t A[maxSizeParticles]; <<
177 /** \brief Particle charge number */ <<
178 Short_t Z[maxSizeParticles]; <<
179 /** \brief Particle strangeness number * <<
180 Short_t S[maxSizeParticles]; <<
181 /** \brief Particle angular momemtum */ <<
182 Short_t J[maxSizeParticles]; <<
183 /** \brief PDG numbering of the particle <<
184 Int_t PDGCode[maxSizeParticles]; <<
185 /** \brief Event bias */ <<
186 Float_t eventBias; <<
187 /** \brief Particle weight due to the bi <<
188 Float_t ParticleBias[maxSizeParticles]; <<
189 /** \brief Particle kinetic energy [MeV] <<
190 Float_t EKin[maxSizeParticles]; <<
191 /** \brief Particle momentum, x componen <<
192 Float_t px[maxSizeParticles]; <<
193 /** \brief Particle momentum, y componen <<
194 Float_t py[maxSizeParticles]; <<
195 /** \brief Particle momentum, z componen <<
196 Float_t pz[maxSizeParticles]; <<
197 /** \brief Particle momentum polar angle <<
198 Float_t theta[maxSizeParticles]; <<
199 /** \brief Particle momentum azimuthal a <<
200 Float_t phi[maxSizeParticles]; <<
201 /** \brief Origin of the particle <<
202 * <<
203 * Should be -1 for cascade particles, o <<
204 * de-excitation particles. */ <<
205 Short_t origin[maxSizeParticles]; <<
206 /** \brief Particle's parent resonance P <<
207 Int_t parentResonancePDGCode[maxSizePart <<
208 /** \brief Particle's parent resonance u <<
209 Int_t parentResonanceID[maxSizeParticles <<
210 /** \brief History of the particle <<
211 * <<
212 * Condensed information about the de-ex <<
213 * cascade particles, it is just an empt <<
214 * from the de-excitation of a cascade r <<
215 * characters. Each character represents <<
216 * the de-excitation process. The curren <<
217 * values and their meanings are the fol <<
218 * <<
219 * e: evaporation product <<
220 * E: evaporation residue <<
221 * m: multifragmentation <<
222 * a: light partner in asymmetric fissio <<
223 * A: heavy partner in asymmetric fissio <<
224 * f: light partner in fission <<
225 * F: heavy partner in fission <<
226 * s: saddle-to-scission emission <<
227 * n: non-statistical emission (decay) * <<
228 std::vector<std::string> history; <<
229 /** \brief Number of remnants */ <<
230 Short_t nRemnants; <<
231 /** \brief Projectile particle type */ <<
232 Int_t projectileType; <<
233 /** \brief Mass number of the target nuc 121 /** \brief Mass number of the target nucleus */
234 Short_t At; 122 Short_t At;
235 /** \brief Charge number of the target n 123 /** \brief Charge number of the target nucleus */
236 Short_t Zt; 124 Short_t Zt;
237 /** \brief Strangeness number of the tar << 125
238 Short_t St; <<
239 /** \brief Mass number of the projectile 126 /** \brief Mass number of the projectile nucleus */
240 Short_t Ap; 127 Short_t Ap;
241 /** \brief Charge number of the projecti 128 /** \brief Charge number of the projectile nucleus */
242 Short_t Zp; 129 Short_t Zp;
243 /** \brief Strangeness number of the pro <<
244 Short_t Sp; <<
245 /** \brief Projectile kinetic energy giv 130 /** \brief Projectile kinetic energy given as input */
246 Float_t Ep; 131 Float_t Ep;
>> 132
247 /** \brief Impact parameter [fm] */ 133 /** \brief Impact parameter [fm] */
248 Float_t impactParameter; 134 Float_t impactParameter;
249 /** \brief Number of accepted two-body c 135 /** \brief Number of accepted two-body collisions */
250 Int_t nCollisions; 136 Int_t nCollisions;
251 /** \brief Cascade stopping time [fm/c] 137 /** \brief Cascade stopping time [fm/c] */
252 Float_t stoppingTime; 138 Float_t stoppingTime;
>> 139
253 /** \brief Energy-conservation balance [ 140 /** \brief Energy-conservation balance [MeV] */
254 Float_t EBalance; 141 Float_t EBalance;
255 /** \brief First value for the energy-co <<
256 Float_t firstEBalance; <<
257 /** \brief Longitudinal momentum-conserv 142 /** \brief Longitudinal momentum-conservation balance [MeV/c] */
258 Float_t pLongBalance; 143 Float_t pLongBalance;
259 /** \brief Transverse momentum-conservat 144 /** \brief Transverse momentum-conservation balance [MeV/c] */
260 Float_t pTransBalance; 145 Float_t pTransBalance;
>> 146
261 /** \brief Number of cascade particles * 147 /** \brief Number of cascade particles */
262 Short_t nCascadeParticles; 148 Short_t nCascadeParticles;
>> 149 /** \brief Number of remnants */
>> 150 Int_t nRemnants;
>> 151 /** \brief Total number of emitted particles */
>> 152 Int_t nParticles;
>> 153
263 /** \brief True if the event is transpar 154 /** \brief True if the event is transparent */
264 Bool_t transparent; 155 Bool_t transparent;
265 /** \brief True if annihilation at rest <<
266 Bool_t annihilationP; <<
267 /** \brief True if annihilation at rest <<
268 Bool_t annihilationN; <<
269 /** \brief True if the event is a forced 156 /** \brief True if the event is a forced CN */
270 Bool_t forcedCompoundNucleus; 157 Bool_t forcedCompoundNucleus;
271 /** \brief True if the event is a nucleo << 158 /** \brief True if the event is absorption */
272 Bool_t nucleonAbsorption; 159 Bool_t nucleonAbsorption;
273 /** \brief True if the event is a pion a << 160 /** \brief True if the event is absorption */
274 Bool_t pionAbsorption; 161 Bool_t pionAbsorption;
275 /** \brief Number of accepted Delta deca 162 /** \brief Number of accepted Delta decays */
276 Int_t nDecays; 163 Int_t nDecays;
277 /** \brief Number of accepted SRC collis <<
278 Int_t nSrcCollisions; <<
279 /** \brief Number of src pairs */ <<
280 Int_t nSrcPairs; <<
281 /** \brief Number of two-body collisions 164 /** \brief Number of two-body collisions blocked by Pauli or CDPP */
282 Int_t nBlockedCollisions; 165 Int_t nBlockedCollisions;
283 /** \brief Number of decays blocked by P 166 /** \brief Number of decays blocked by Pauli or CDPP */
284 Int_t nBlockedDecays; 167 Int_t nBlockedDecays;
>> 168 /** \brief Number of reflection avatars */
285 /** \brief Effective (Coulomb-distorted) 169 /** \brief Effective (Coulomb-distorted) impact parameter [fm] */
286 Float_t effectiveImpactParameter; 170 Float_t effectiveImpactParameter;
287 /** \brief Event involved deltas in the << 171
>> 172 /// \brief Event involved deltas in the nucleus at the end of the cascade
288 Bool_t deltasInside; 173 Bool_t deltasInside;
289 /** \brief Event involved sigmas in the << 174 /// \brief Event involved forced delta decays inside the nucleus
290 Bool_t sigmasInside; <<
291 /** \brief Event involved kaons in the n <<
292 Bool_t kaonsInside; <<
293 /** \brief Event involved antikaons in t <<
294 Bool_t antikaonsInside; <<
295 /** \brief Event involved lambdas in the <<
296 Bool_t lambdasInside; <<
297 /** \brief Event involved forced delta d <<
298 Bool_t forcedDeltasInside; 175 Bool_t forcedDeltasInside;
299 /** \brief Event involved forced delta d << 176 /// \brief Event involved forced delta decays outside the nucleus
300 Bool_t forcedDeltasOutside; 177 Bool_t forcedDeltasOutside;
301 /** \brief Event involved forced eta/ome << 178 /// \brief Event involved cluster decay
302 Bool_t forcedPionResonancesOutside; <<
303 /** \brief Event involved forced strange <<
304 Bool_t absorbedStrangeParticle; <<
305 /** \brief Event involved forced Sigma Z <<
306 Bool_t forcedSigmaOutside; <<
307 /** \brief Event involved forced antiKao <<
308 Bool_t forcedStrangeInside; <<
309 /** \brief Number of forced Lambda emit <<
310 Int_t emitLambda; <<
311 /** \brief Event involved forced Kaon em <<
312 Bool_t emitKaon; <<
313 /** \brief Event involved cluster decay <<
314 Bool_t clusterDecay; 179 Bool_t clusterDecay;
>> 180
315 /** \brief Time of the first collision [ 181 /** \brief Time of the first collision [fm/c] */
316 Float_t firstCollisionTime; 182 Float_t firstCollisionTime;
317 /** \brief Cross section of the first co 183 /** \brief Cross section of the first collision (mb) */
318 Float_t firstCollisionXSec; 184 Float_t firstCollisionXSec;
319 /** \brief Position of the spectator on << 185
320 Float_t firstCollisionSpectatorPosition; <<
321 /** \brief Momentum of the spectator on <<
322 Float_t firstCollisionSpectatorMomentum; <<
323 /** \brief True if the first collision w <<
324 Bool_t firstCollisionIsElastic; <<
325 /** \brief Number of reflection avatars <<
326 Int_t nReflectionAvatars; 186 Int_t nReflectionAvatars;
327 /** \brief Number of collision avatars * 187 /** \brief Number of collision avatars */
328 Int_t nCollisionAvatars; 188 Int_t nCollisionAvatars;
329 /** \brief Number of decay avatars */ 189 /** \brief Number of decay avatars */
330 Int_t nDecayAvatars; 190 Int_t nDecayAvatars;
331 /** \brief Number of dynamical spectator << 191
>> 192 /// \brief Number of dynamical spectators that were merged back into the projectile remnant
332 Int_t nUnmergedSpectators; 193 Int_t nUnmergedSpectators;
333 /** \brief Number of attempted collision << 194
334 Int_t nEnergyViolationInteraction; << 195 /** \brief Maximum array size for remnants */
335 /** \brief Emission time [fm/c] */ << 196 static const Short_t maxSizeRemnants = 10;
336 Float_t emissionTime[maxSizeParticles]; <<
337 /** \brief Remnant mass number */ 197 /** \brief Remnant mass number */
338 Short_t ARem[maxSizeRemnants]; 198 Short_t ARem[maxSizeRemnants];
339 /** \brief Remnant charge number */ 199 /** \brief Remnant charge number */
340 Short_t ZRem[maxSizeRemnants]; 200 Short_t ZRem[maxSizeRemnants];
341 /** \brief Remnant strangeness number */ <<
342 Short_t SRem[maxSizeRemnants]; <<
343 /** \brief Remnant excitation energy [Me 201 /** \brief Remnant excitation energy [MeV] */
344 Float_t EStarRem[maxSizeRemnants]; 202 Float_t EStarRem[maxSizeRemnants];
345 /** \brief Remnant spin [\f$\hbar\f$] */ 203 /** \brief Remnant spin [\f$\hbar\f$] */
346 Float_t JRem[maxSizeRemnants]; 204 Float_t JRem[maxSizeRemnants];
347 /** \brief Remnant kinetic energy [MeV] 205 /** \brief Remnant kinetic energy [MeV] */
348 Float_t EKinRem[maxSizeRemnants]; 206 Float_t EKinRem[maxSizeRemnants];
349 /** \brief Remnant momentum, x component 207 /** \brief Remnant momentum, x component [MeV/c] */
350 Float_t pxRem[maxSizeRemnants]; 208 Float_t pxRem[maxSizeRemnants];
351 /** \brief Remnant momentum, y component 209 /** \brief Remnant momentum, y component [MeV/c] */
352 Float_t pyRem[maxSizeRemnants]; 210 Float_t pyRem[maxSizeRemnants];
353 /** \brief Remnant momentum, z component 211 /** \brief Remnant momentum, z component [MeV/c] */
354 Float_t pzRem[maxSizeRemnants]; 212 Float_t pzRem[maxSizeRemnants];
355 /** \brief Remnant momentum polar angle 213 /** \brief Remnant momentum polar angle [radians] */
356 Float_t thetaRem[maxSizeRemnants]; 214 Float_t thetaRem[maxSizeRemnants];
357 /** \brief Remnant momentum azimuthal an 215 /** \brief Remnant momentum azimuthal angle [radians] */
358 Float_t phiRem[maxSizeRemnants]; 216 Float_t phiRem[maxSizeRemnants];
359 /** \brief Remnant angular momentum, x c << 217 /** \brief Remnant angular momentum, x component [hbar] */
360 Float_t jxRem[maxSizeRemnants]; 218 Float_t jxRem[maxSizeRemnants];
361 /** \brief Remnant angular momentum, y c << 219 /** \brief Remnant angular momentum, y component [hbar] */
362 Float_t jyRem[maxSizeRemnants]; 220 Float_t jyRem[maxSizeRemnants];
363 /** \brief Remnant angular momentum, z c << 221 /** \brief Remnant angular momentum, z component [hbar] */
364 Float_t jzRem[maxSizeRemnants]; 222 Float_t jzRem[maxSizeRemnants];
>> 223
>> 224 /** \brief Maximum array size for emitted particles */
>> 225 static const Short_t maxSizeParticles = 1000;
>> 226 /** \brief Particle mass number */
>> 227 Short_t A[maxSizeParticles];
>> 228 /** \brief Particle charge number */
>> 229 Short_t Z[maxSizeParticles];
>> 230 /** \brief Emission time [fm/c] */
>> 231 Float_t emissionTime[maxSizeParticles];
>> 232 /** \brief Particle kinetic energy [MeV] */
>> 233 Float_t EKin[maxSizeParticles];
>> 234 /** \brief Particle momentum, x component [MeV/c] */
>> 235 Float_t px[maxSizeParticles];
>> 236 /** \brief Particle momentum, y component [MeV/c] */
>> 237 Float_t py[maxSizeParticles];
>> 238 /** \brief Particle momentum, z component [MeV/c] */
>> 239 Float_t pz[maxSizeParticles];
>> 240 /** \brief Particle momentum polar angle [radians] */
>> 241 Float_t theta[maxSizeParticles];
>> 242 /** \brief Particle momentum azimuthal angle [radians] */
>> 243 Float_t phi[maxSizeParticles];
>> 244 /** \brief Origin of the particle
>> 245 *
>> 246 * Should be -1 for cascade particles, or the number of the remnant for
>> 247 * de-excitation particles.
>> 248 *
>> 249 */
>> 250 Short_t origin[maxSizeParticles];
>> 251 /** \brief History of the particle
>> 252 *
>> 253 * Condensed information about the de-excitation chain of a particle. For
>> 254 * cascade particles, it is just an empty string. For particles arising
>> 255 * from the de-excitation of a cascade remnant, it is a string of
>> 256 * characters. Each character represents one or more identical steps in
>> 257 * the de-excitation process. The currently defined possible character
>> 258 * values and their meanings are the following:
>> 259 *
>> 260 * e: evaporation product
>> 261 * E: evaporation residue
>> 262 * m: multifragmentation
>> 263 * a: light partner in asymmetric fission or IMF emission
>> 264 * A: heavy partner in asymmetric fission or IMF emission
>> 265 * f: light partner in fission
>> 266 * F: heavy partner in fission
>> 267 * s: saddle-to-scission emission
>> 268 * n: non-statistical emission (decay)
>> 269 */
>> 270 std::vector<std::string> history;
>> 271
>> 272 #ifdef INCL_INVERSE_KINEMATICS
365 /** \brief Particle kinetic energy, in i 273 /** \brief Particle kinetic energy, in inverse kinematics [MeV] */
366 Float_t EKinPrime[maxSizeParticles]; 274 Float_t EKinPrime[maxSizeParticles];
367 /** \brief Particle momentum, z componen 275 /** \brief Particle momentum, z component, in inverse kinematics [MeV/c] */
368 Float_t pzPrime[maxSizeParticles]; 276 Float_t pzPrime[maxSizeParticles];
369 /** \brief Particle momentum polar angle 277 /** \brief Particle momentum polar angle, in inverse kinematics [radians] */
370 Float_t thetaPrime[maxSizeParticles]; 278 Float_t thetaPrime[maxSizeParticles];
>> 279 #endif // INCL_INVERSE_KINEMATICS
371 280
372 /** \brief Reset the EventInfo members * 281 /** \brief Reset the EventInfo members */
373 void reset() { 282 void reset() {
374 nParticles = 0; <<
375 event = 0; <<
376 eventBias = (Float_t)0.0; <<
377 history.clear(); <<
378 nRemnants = 0; <<
379 projectileType = 0; <<
380 At = 0; <<
381 Zt = 0; <<
382 St = 0; <<
383 Ap = 0; 283 Ap = 0;
384 Zp = 0; 284 Zp = 0;
385 Sp = 0; << 285 At = 0;
386 Ep = (Float_t)0.0; << 286 Zt = 0;
387 impactParameter = (Float_t)0.0; << 287 impactParameter = 0.0;
>> 288 effectiveImpactParameter = 0.0;
>> 289 stoppingTime = 0.0;
>> 290 EBalance = 0.0;
>> 291 pLongBalance = 0.0;
>> 292 pTransBalance = 0.0;
388 nCollisions = 0; 293 nCollisions = 0;
389 stoppingTime = (Float_t)0.0; << 294 nBlockedCollisions = 0;
390 EBalance = (Float_t)0.0; << 295 nDecays = 0;
391 firstEBalance = (Float_t)0.0; << 296 nBlockedDecays= 0;
392 pLongBalance = (Float_t)0.0; << 297 nDecays = 0;
393 pTransBalance = (Float_t)0.0; <<
394 nCascadeParticles = 0; 298 nCascadeParticles = 0;
395 transparent = false; << 299 nRemnants = 0;
396 annihilationP = false; << 300 nParticles = 0;
397 annihilationN = false; << 301 transparent = true;
398 forcedCompoundNucleus = false; 302 forcedCompoundNucleus = false;
399 nucleonAbsorption = false; << 303 nucleonAbsorption = false;
400 pionAbsorption = false; << 304 pionAbsorption = false;
401 nDecays = 0; <<
402 nSrcCollisions = 0; <<
403 nSrcPairs = 0; <<
404 nBlockedCollisions = 0; <<
405 nBlockedDecays = 0; <<
406 effectiveImpactParameter = (Float_t)0. <<
407 deltasInside = false; <<
408 sigmasInside = false; <<
409 kaonsInside = false; <<
410 antikaonsInside = false; <<
411 lambdasInside = false; <<
412 forcedDeltasInside = false; 305 forcedDeltasInside = false;
413 forcedDeltasOutside = false; 306 forcedDeltasOutside = false;
414 forcedPionResonancesOutside = false; << 307 deltasInside = false;
415 absorbedStrangeParticle = false; <<
416 forcedSigmaOutside = false; <<
417 forcedStrangeInside = false; <<
418 emitLambda = 0; <<
419 emitKaon = false; <<
420 clusterDecay = false; 308 clusterDecay = false;
421 firstCollisionTime = (Float_t)0.0; <<
422 firstCollisionXSec = (Float_t)0.0; <<
423 firstCollisionSpectatorPosition = (Flo <<
424 firstCollisionSpectatorMomentum = (Flo <<
425 firstCollisionIsElastic = false; <<
426 nReflectionAvatars = 0; <<
427 nCollisionAvatars = 0; <<
428 nDecayAvatars = 0; <<
429 nUnmergedSpectators = 0; 309 nUnmergedSpectators = 0;
430 nEnergyViolationInteraction = 0; <<
431 <<
432 } 310 }
433 311
434 /// \brief Move a remnant to the particl << 312 #ifdef INCL_INVERSE_KINEMATICS
435 void remnantToParticle(const G4int remna <<
436 <<
437 /// \brief Fill the variables describing <<
438 void fillInverseKinematics(const Double_ 313 void fillInverseKinematics(const Double_t gamma);
>> 314 #endif // INCL_INVERSE_KINEMATICS
439 }; 315 };
440 } 316 }
441 317
442 #endif /* G4INCLEVENTINFO_HH_HH */ << 318 #endif /* G4INCLEVENTINFO_HH */
443 319