| 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_HH
48 #define G4INCLEVENTINFO_HH_HH 1 47 #define G4INCLEVENTINFO_HH_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 nParticles(0),
70 event(0), <<
71 eventBias((Float_t)0.0), <<
72 nRemnants(0), 69 nRemnants(0),
73 projectileType(0), 70 projectileType(0),
74 At(0), 71 At(0),
75 Zt(0), 72 Zt(0),
76 St(0), <<
77 Ap(0), 73 Ap(0),
78 Zp(0), 74 Zp(0),
79 Sp(0), <<
80 Ep((Float_t)0.0), 75 Ep((Float_t)0.0),
81 impactParameter((Float_t)0.0), 76 impactParameter((Float_t)0.0),
82 nCollisions(0), 77 nCollisions(0),
83 stoppingTime((Float_t)0.0), 78 stoppingTime((Float_t)0.0),
84 EBalance((Float_t)0.0), 79 EBalance((Float_t)0.0),
85 firstEBalance((Float_t)0.0), <<
86 pLongBalance((Float_t)0.0), 80 pLongBalance((Float_t)0.0),
87 pTransBalance((Float_t)0.0), 81 pTransBalance((Float_t)0.0),
88 nCascadeParticles(0), 82 nCascadeParticles(0),
89 transparent(false), 83 transparent(false),
90 annihilationP(false), <<
91 annihilationN(false), <<
92 forcedCompoundNucleus(false), 84 forcedCompoundNucleus(false),
93 nucleonAbsorption(false), 85 nucleonAbsorption(false),
94 pionAbsorption(false), 86 pionAbsorption(false),
95 nDecays(0), 87 nDecays(0),
96 nSrcCollisions(0), <<
97 nSrcPairs(0), <<
98 nBlockedCollisions(0), 88 nBlockedCollisions(0),
99 nBlockedDecays(0), 89 nBlockedDecays(0),
100 effectiveImpactParameter((Float_t)0.0) 90 effectiveImpactParameter((Float_t)0.0),
101 deltasInside(false), 91 deltasInside(false),
102 sigmasInside(false), <<
103 kaonsInside(false), <<
104 antikaonsInside(false), <<
105 lambdasInside(false), <<
106 forcedDeltasInside(false), 92 forcedDeltasInside(false),
107 forcedDeltasOutside(false), 93 forcedDeltasOutside(false),
108 forcedPionResonancesOutside(false), <<
109 absorbedStrangeParticle(false), <<
110 forcedSigmaOutside(false), <<
111 forcedStrangeInside(false), <<
112 emitLambda(0), <<
113 emitKaon(false), <<
114 clusterDecay(false), 94 clusterDecay(false),
115 firstCollisionTime((Float_t)0.0), 95 firstCollisionTime((Float_t)0.0),
116 firstCollisionXSec((Float_t)0.0), 96 firstCollisionXSec((Float_t)0.0),
117 firstCollisionSpectatorPosition((Float 97 firstCollisionSpectatorPosition((Float_t)0.0),
118 firstCollisionSpectatorMomentum((Float 98 firstCollisionSpectatorMomentum((Float_t)0.0),
119 firstCollisionIsElastic(false), 99 firstCollisionIsElastic(false),
120 nReflectionAvatars(0), 100 nReflectionAvatars(0),
121 nCollisionAvatars(0), 101 nCollisionAvatars(0),
122 nDecayAvatars(0), 102 nDecayAvatars(0),
123 nUnmergedSpectators(0), 103 nUnmergedSpectators(0),
124 nEnergyViolationInteraction(0) 104 nEnergyViolationInteraction(0)
>> 105 #ifdef INCL_INVERSE_KINEMATICS
125 106
>> 107 #endif
126 { 108 {
127 std::fill_n(A, maxSizeParticles, 0); 109 std::fill_n(A, maxSizeParticles, 0);
128 std::fill_n(Z, maxSizeParticles, 0); 110 std::fill_n(Z, maxSizeParticles, 0);
129 std::fill_n(S, maxSizeParticles, 0); <<
130 std::fill_n(J, maxSizeParticles, 0); <<
131 std::fill_n(PDGCode, maxSizeParticles, <<
132 std::fill_n(ParticleBias, maxSizeParti <<
133 std::fill_n(EKin, maxSizeParticles, (F 111 std::fill_n(EKin, maxSizeParticles, (Float_t)0.0);
134 std::fill_n(px, maxSizeParticles, (Flo 112 std::fill_n(px, maxSizeParticles, (Float_t)0.0);
135 std::fill_n(py, maxSizeParticles, (Flo 113 std::fill_n(py, maxSizeParticles, (Float_t)0.0);
136 std::fill_n(pz, maxSizeParticles, (Flo 114 std::fill_n(pz, maxSizeParticles, (Float_t)0.0);
137 std::fill_n(theta, maxSizeParticles, ( 115 std::fill_n(theta, maxSizeParticles, (Float_t)0.0);
138 std::fill_n(phi, maxSizeParticles, (Fl 116 std::fill_n(phi, maxSizeParticles, (Float_t)0.0);
139 std::fill_n(origin, maxSizeParticles, 117 std::fill_n(origin, maxSizeParticles, 0);
140 std::fill_n(parentResonancePDGCode, ma <<
141 std::fill_n(parentResonanceID, maxSize <<
142 std::fill_n(emissionTime, maxSizeParti 118 std::fill_n(emissionTime, maxSizeParticles, (Float_t)0.0);
143 std::fill_n(ARem, maxSizeRemnants, 0); 119 std::fill_n(ARem, maxSizeRemnants, 0);
144 std::fill_n(ZRem, maxSizeRemnants, 0); 120 std::fill_n(ZRem, maxSizeRemnants, 0);
145 std::fill_n(SRem, maxSizeRemnants, 0); <<
146 std::fill_n(EStarRem, maxSizeRemnants, 121 std::fill_n(EStarRem, maxSizeRemnants, (Float_t)0.0);
147 std::fill_n(JRem, maxSizeRemnants, (Fl 122 std::fill_n(JRem, maxSizeRemnants, (Float_t)0.0);
148 std::fill_n(EKinRem, maxSizeRemnants, 123 std::fill_n(EKinRem, maxSizeRemnants, (Float_t)0.0);
149 std::fill_n(pxRem, maxSizeRemnants, (F 124 std::fill_n(pxRem, maxSizeRemnants, (Float_t)0.0);
150 std::fill_n(pyRem, maxSizeRemnants, (F 125 std::fill_n(pyRem, maxSizeRemnants, (Float_t)0.0);
151 std::fill_n(pzRem, maxSizeRemnants, (F 126 std::fill_n(pzRem, maxSizeRemnants, (Float_t)0.0);
152 std::fill_n(thetaRem, maxSizeRemnants, 127 std::fill_n(thetaRem, maxSizeRemnants, (Float_t)0.0);
153 std::fill_n(phiRem, maxSizeRemnants, ( 128 std::fill_n(phiRem, maxSizeRemnants, (Float_t)0.0);
154 std::fill_n(jxRem, maxSizeRemnants, (F 129 std::fill_n(jxRem, maxSizeRemnants, (Float_t)0.0);
155 std::fill_n(jyRem, maxSizeRemnants, (F 130 std::fill_n(jyRem, maxSizeRemnants, (Float_t)0.0);
156 std::fill_n(jzRem, maxSizeRemnants, (F 131 std::fill_n(jzRem, maxSizeRemnants, (Float_t)0.0);
>> 132 #ifdef INCL_INVERSE_KINEMATICS
157 std::fill_n(EKinPrime, maxSizeParticle 133 std::fill_n(EKinPrime, maxSizeParticles, (Float_t)0.0);
158 std::fill_n(pzPrime, maxSizeParticles, 134 std::fill_n(pzPrime, maxSizeParticles, (Float_t)0.0);
159 std::fill_n(thetaPrime, maxSizeParticl 135 std::fill_n(thetaPrime, maxSizeParticles, (Float_t)0.0);
>> 136 #endif
160 } 137 }
161 138
162 /** \brief Number of the event */ 139 /** \brief Number of the event */
163 static G4ThreadLocal Int_t eventNumber; 140 static G4ThreadLocal Int_t eventNumber;
164 141
165 /** \brief Maximum array size for remnan 142 /** \brief Maximum array size for remnants */
166 static const Short_t maxSizeRemnants = 1 143 static const Short_t maxSizeRemnants = 10;
167 144
168 /** \brief Maximum array size for emitte 145 /** \brief Maximum array size for emitted particles */
169 static const Short_t maxSizeParticles = 146 static const Short_t maxSizeParticles = 1000;
170 147
171 /** \brief Number of particles in the fi 148 /** \brief Number of particles in the final state */
172 Short_t nParticles; 149 Short_t nParticles;
173 /** \brief Sequential number of the even <<
174 Int_t event; <<
175 /** \brief Particle mass number */ 150 /** \brief Particle mass number */
176 Short_t A[maxSizeParticles]; 151 Short_t A[maxSizeParticles];
177 /** \brief Particle charge number */ 152 /** \brief Particle charge number */
178 Short_t Z[maxSizeParticles]; 153 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] 154 /** \brief Particle kinetic energy [MeV] */
190 Float_t EKin[maxSizeParticles]; 155 Float_t EKin[maxSizeParticles];
191 /** \brief Particle momentum, x componen 156 /** \brief Particle momentum, x component [MeV/c] */
192 Float_t px[maxSizeParticles]; 157 Float_t px[maxSizeParticles];
193 /** \brief Particle momentum, y componen 158 /** \brief Particle momentum, y component [MeV/c] */
194 Float_t py[maxSizeParticles]; 159 Float_t py[maxSizeParticles];
195 /** \brief Particle momentum, z componen 160 /** \brief Particle momentum, z component [MeV/c] */
196 Float_t pz[maxSizeParticles]; 161 Float_t pz[maxSizeParticles];
197 /** \brief Particle momentum polar angle 162 /** \brief Particle momentum polar angle [radians] */
198 Float_t theta[maxSizeParticles]; 163 Float_t theta[maxSizeParticles];
199 /** \brief Particle momentum azimuthal a 164 /** \brief Particle momentum azimuthal angle [radians] */
200 Float_t phi[maxSizeParticles]; 165 Float_t phi[maxSizeParticles];
201 /** \brief Origin of the particle 166 /** \brief Origin of the particle
202 * 167 *
203 * Should be -1 for cascade particles, o 168 * Should be -1 for cascade particles, or the number of the remnant for
204 * de-excitation particles. */ 169 * de-excitation particles. */
205 Short_t origin[maxSizeParticles]; 170 Short_t origin[maxSizeParticles];
206 /** \brief Particle's parent resonance P << 171 /** \brief Emission time [fm/c] */
207 Int_t parentResonancePDGCode[maxSizePart << 172 Float_t emissionTime[maxSizeParticles];
208 /** \brief Particle's parent resonance u <<
209 Int_t parentResonanceID[maxSizeParticles <<
210 /** \brief History of the particle 173 /** \brief History of the particle
211 * 174 *
212 * Condensed information about the de-ex 175 * Condensed information about the de-excitation chain of a particle. For
213 * cascade particles, it is just an empt 176 * cascade particles, it is just an empty string. For particles arising
214 * from the de-excitation of a cascade r 177 * from the de-excitation of a cascade remnant, it is a string of
215 * characters. Each character represents 178 * characters. Each character represents one or more identical steps in
216 * the de-excitation process. The curren 179 * the de-excitation process. The currently defined possible character
217 * values and their meanings are the fol 180 * values and their meanings are the following:
218 * 181 *
219 * e: evaporation product 182 * e: evaporation product
220 * E: evaporation residue 183 * E: evaporation residue
221 * m: multifragmentation 184 * m: multifragmentation
222 * a: light partner in asymmetric fissio 185 * a: light partner in asymmetric fission or IMF emission
223 * A: heavy partner in asymmetric fissio 186 * A: heavy partner in asymmetric fission or IMF emission
224 * f: light partner in fission 187 * f: light partner in fission
225 * F: heavy partner in fission 188 * F: heavy partner in fission
226 * s: saddle-to-scission emission 189 * s: saddle-to-scission emission
227 * n: non-statistical emission (decay) * 190 * n: non-statistical emission (decay) */
228 std::vector<std::string> history; 191 std::vector<std::string> history;
229 /** \brief Number of remnants */ 192 /** \brief Number of remnants */
230 Short_t nRemnants; << 193 Int_t nRemnants;
>> 194 /** \brief Remnant mass number */
>> 195 Short_t ARem[maxSizeRemnants];
>> 196 /** \brief Remnant charge number */
>> 197 Short_t ZRem[maxSizeRemnants];
>> 198 /** \brief Remnant excitation energy [MeV] */
>> 199 Float_t EStarRem[maxSizeRemnants];
>> 200 /** \brief Remnant spin [\f$\hbar\f$] */
>> 201 Float_t JRem[maxSizeRemnants];
>> 202 /** \brief Remnant kinetic energy [MeV] */
>> 203 Float_t EKinRem[maxSizeRemnants];
>> 204 /** \brief Remnant momentum, x component [MeV/c] */
>> 205 Float_t pxRem[maxSizeRemnants];
>> 206 /** \brief Remnant momentum, y component [MeV/c] */
>> 207 Float_t pyRem[maxSizeRemnants];
>> 208 /** \brief Remnant momentum, z component [MeV/c] */
>> 209 Float_t pzRem[maxSizeRemnants];
>> 210 /** \brief Remnant momentum polar angle [radians] */
>> 211 Float_t thetaRem[maxSizeRemnants];
>> 212 /** \brief Remnant momentum azimuthal angle [radians] */
>> 213 Float_t phiRem[maxSizeRemnants];
>> 214 /** \brief Remnant angular momentum, x component [\f$\hbar\f$] */
>> 215 Float_t jxRem[maxSizeRemnants];
>> 216 /** \brief Remnant angular momentum, y component [\f$\hbar\f$] */
>> 217 Float_t jyRem[maxSizeRemnants];
>> 218 /** \brief Remnant angular momentum, z component [\f$\hbar\f$] */
>> 219 Float_t jzRem[maxSizeRemnants];
231 /** \brief Projectile particle type */ 220 /** \brief Projectile particle type */
232 Int_t projectileType; 221 Int_t projectileType;
233 /** \brief Mass number of the target nuc 222 /** \brief Mass number of the target nucleus */
234 Short_t At; 223 Short_t At;
235 /** \brief Charge number of the target n 224 /** \brief Charge number of the target nucleus */
236 Short_t Zt; 225 Short_t Zt;
237 /** \brief Strangeness number of the tar <<
238 Short_t St; <<
239 /** \brief Mass number of the projectile 226 /** \brief Mass number of the projectile nucleus */
240 Short_t Ap; 227 Short_t Ap;
241 /** \brief Charge number of the projecti 228 /** \brief Charge number of the projectile nucleus */
242 Short_t Zp; 229 Short_t Zp;
243 /** \brief Strangeness number of the pro <<
244 Short_t Sp; <<
245 /** \brief Projectile kinetic energy giv 230 /** \brief Projectile kinetic energy given as input */
246 Float_t Ep; 231 Float_t Ep;
247 /** \brief Impact parameter [fm] */ 232 /** \brief Impact parameter [fm] */
248 Float_t impactParameter; 233 Float_t impactParameter;
249 /** \brief Number of accepted two-body c 234 /** \brief Number of accepted two-body collisions */
250 Int_t nCollisions; 235 Int_t nCollisions;
251 /** \brief Cascade stopping time [fm/c] 236 /** \brief Cascade stopping time [fm/c] */
252 Float_t stoppingTime; 237 Float_t stoppingTime;
253 /** \brief Energy-conservation balance [ 238 /** \brief Energy-conservation balance [MeV] */
254 Float_t EBalance; 239 Float_t EBalance;
255 /** \brief First value for the energy-co <<
256 Float_t firstEBalance; <<
257 /** \brief Longitudinal momentum-conserv 240 /** \brief Longitudinal momentum-conservation balance [MeV/c] */
258 Float_t pLongBalance; 241 Float_t pLongBalance;
259 /** \brief Transverse momentum-conservat 242 /** \brief Transverse momentum-conservation balance [MeV/c] */
260 Float_t pTransBalance; 243 Float_t pTransBalance;
261 /** \brief Number of cascade particles * 244 /** \brief Number of cascade particles */
262 Short_t nCascadeParticles; 245 Short_t nCascadeParticles;
263 /** \brief True if the event is transpar 246 /** \brief True if the event is transparent */
264 Bool_t transparent; 247 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 248 /** \brief True if the event is a forced CN */
270 Bool_t forcedCompoundNucleus; 249 Bool_t forcedCompoundNucleus;
271 /** \brief True if the event is a nucleo 250 /** \brief True if the event is a nucleon absorption */
272 Bool_t nucleonAbsorption; 251 Bool_t nucleonAbsorption;
273 /** \brief True if the event is a pion a 252 /** \brief True if the event is a pion absorption */
274 Bool_t pionAbsorption; 253 Bool_t pionAbsorption;
275 /** \brief Number of accepted Delta deca 254 /** \brief Number of accepted Delta decays */
276 Int_t nDecays; 255 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 256 /** \brief Number of two-body collisions blocked by Pauli or CDPP */
282 Int_t nBlockedCollisions; 257 Int_t nBlockedCollisions;
283 /** \brief Number of decays blocked by P 258 /** \brief Number of decays blocked by Pauli or CDPP */
284 Int_t nBlockedDecays; 259 Int_t nBlockedDecays;
285 /** \brief Effective (Coulomb-distorted) 260 /** \brief Effective (Coulomb-distorted) impact parameter [fm] */
286 Float_t effectiveImpactParameter; 261 Float_t effectiveImpactParameter;
287 /** \brief Event involved deltas in the 262 /** \brief Event involved deltas in the nucleus at the end of the cascade */
288 Bool_t deltasInside; 263 Bool_t deltasInside;
289 /** \brief Event involved sigmas in the <<
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 264 /** \brief Event involved forced delta decays inside the nucleus */
298 Bool_t forcedDeltasInside; 265 Bool_t forcedDeltasInside;
299 /** \brief Event involved forced delta d 266 /** \brief Event involved forced delta decays outside the nucleus */
300 Bool_t forcedDeltasOutside; 267 Bool_t forcedDeltasOutside;
301 /** \brief Event involved forced eta/ome <<
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 268 /** \brief Event involved cluster decay */
314 Bool_t clusterDecay; 269 Bool_t clusterDecay;
315 /** \brief Time of the first collision [ 270 /** \brief Time of the first collision [fm/c] */
316 Float_t firstCollisionTime; 271 Float_t firstCollisionTime;
317 /** \brief Cross section of the first co 272 /** \brief Cross section of the first collision (mb) */
318 Float_t firstCollisionXSec; 273 Float_t firstCollisionXSec;
319 /** \brief Position of the spectator on 274 /** \brief Position of the spectator on the first collision (fm) */
320 Float_t firstCollisionSpectatorPosition; 275 Float_t firstCollisionSpectatorPosition;
321 /** \brief Momentum of the spectator on 276 /** \brief Momentum of the spectator on the first collision (fm) */
322 Float_t firstCollisionSpectatorMomentum; 277 Float_t firstCollisionSpectatorMomentum;
323 /** \brief True if the first collision w 278 /** \brief True if the first collision was elastic */
324 Bool_t firstCollisionIsElastic; 279 Bool_t firstCollisionIsElastic;
325 /** \brief Number of reflection avatars 280 /** \brief Number of reflection avatars */
326 Int_t nReflectionAvatars; 281 Int_t nReflectionAvatars;
327 /** \brief Number of collision avatars * 282 /** \brief Number of collision avatars */
328 Int_t nCollisionAvatars; 283 Int_t nCollisionAvatars;
329 /** \brief Number of decay avatars */ 284 /** \brief Number of decay avatars */
330 Int_t nDecayAvatars; 285 Int_t nDecayAvatars;
331 /** \brief Number of dynamical spectator 286 /** \brief Number of dynamical spectators that were merged back into the projectile remnant */
332 Int_t nUnmergedSpectators; 287 Int_t nUnmergedSpectators;
333 /** \brief Number of attempted collision 288 /** \brief Number of attempted collisions/decays for which the energy-conservation algorithm failed to find a solution. */
334 Int_t nEnergyViolationInteraction; 289 Int_t nEnergyViolationInteraction;
335 /** \brief Emission time [fm/c] */ << 290
336 Float_t emissionTime[maxSizeParticles]; << 291 #ifdef INCL_INVERSE_KINEMATICS
337 /** \brief Remnant mass number */ <<
338 Short_t ARem[maxSizeRemnants]; <<
339 /** \brief Remnant charge number */ <<
340 Short_t ZRem[maxSizeRemnants]; <<
341 /** \brief Remnant strangeness number */ <<
342 Short_t SRem[maxSizeRemnants]; <<
343 /** \brief Remnant excitation energy [Me <<
344 Float_t EStarRem[maxSizeRemnants]; <<
345 /** \brief Remnant spin [\f$\hbar\f$] */ <<
346 Float_t JRem[maxSizeRemnants]; <<
347 /** \brief Remnant kinetic energy [MeV] <<
348 Float_t EKinRem[maxSizeRemnants]; <<
349 /** \brief Remnant momentum, x component <<
350 Float_t pxRem[maxSizeRemnants]; <<
351 /** \brief Remnant momentum, y component <<
352 Float_t pyRem[maxSizeRemnants]; <<
353 /** \brief Remnant momentum, z component <<
354 Float_t pzRem[maxSizeRemnants]; <<
355 /** \brief Remnant momentum polar angle <<
356 Float_t thetaRem[maxSizeRemnants]; <<
357 /** \brief Remnant momentum azimuthal an <<
358 Float_t phiRem[maxSizeRemnants]; <<
359 /** \brief Remnant angular momentum, x c <<
360 Float_t jxRem[maxSizeRemnants]; <<
361 /** \brief Remnant angular momentum, y c <<
362 Float_t jyRem[maxSizeRemnants]; <<
363 /** \brief Remnant angular momentum, z c <<
364 Float_t jzRem[maxSizeRemnants]; <<
365 /** \brief Particle kinetic energy, in i 292 /** \brief Particle kinetic energy, in inverse kinematics [MeV] */
366 Float_t EKinPrime[maxSizeParticles]; 293 Float_t EKinPrime[maxSizeParticles];
367 /** \brief Particle momentum, z componen 294 /** \brief Particle momentum, z component, in inverse kinematics [MeV/c] */
368 Float_t pzPrime[maxSizeParticles]; 295 Float_t pzPrime[maxSizeParticles];
369 /** \brief Particle momentum polar angle 296 /** \brief Particle momentum polar angle, in inverse kinematics [radians] */
370 Float_t thetaPrime[maxSizeParticles]; 297 Float_t thetaPrime[maxSizeParticles];
>> 298 #endif // INCL_INVERSE_KINEMATICS
371 299
372 /** \brief Reset the EventInfo members * 300 /** \brief Reset the EventInfo members */
373 void reset() { 301 void reset() {
374 nParticles = 0; 302 nParticles = 0;
375 event = 0; <<
376 eventBias = (Float_t)0.0; <<
377 history.clear(); 303 history.clear();
378 nRemnants = 0; 304 nRemnants = 0;
379 projectileType = 0; 305 projectileType = 0;
380 At = 0; 306 At = 0;
381 Zt = 0; 307 Zt = 0;
382 St = 0; <<
383 Ap = 0; 308 Ap = 0;
384 Zp = 0; 309 Zp = 0;
385 Sp = 0; <<
386 Ep = (Float_t)0.0; 310 Ep = (Float_t)0.0;
387 impactParameter = (Float_t)0.0; 311 impactParameter = (Float_t)0.0;
388 nCollisions = 0; 312 nCollisions = 0;
389 stoppingTime = (Float_t)0.0; 313 stoppingTime = (Float_t)0.0;
390 EBalance = (Float_t)0.0; 314 EBalance = (Float_t)0.0;
391 firstEBalance = (Float_t)0.0; <<
392 pLongBalance = (Float_t)0.0; 315 pLongBalance = (Float_t)0.0;
393 pTransBalance = (Float_t)0.0; 316 pTransBalance = (Float_t)0.0;
394 nCascadeParticles = 0; 317 nCascadeParticles = 0;
395 transparent = false; 318 transparent = false;
396 annihilationP = false; <<
397 annihilationN = false; <<
398 forcedCompoundNucleus = false; 319 forcedCompoundNucleus = false;
399 nucleonAbsorption = false; 320 nucleonAbsorption = false;
400 pionAbsorption = false; 321 pionAbsorption = false;
401 nDecays = 0; 322 nDecays = 0;
402 nSrcCollisions = 0; <<
403 nSrcPairs = 0; <<
404 nBlockedCollisions = 0; 323 nBlockedCollisions = 0;
405 nBlockedDecays = 0; 324 nBlockedDecays = 0;
406 effectiveImpactParameter = (Float_t)0. 325 effectiveImpactParameter = (Float_t)0.0;
407 deltasInside = false; 326 deltasInside = false;
408 sigmasInside = false; <<
409 kaonsInside = false; <<
410 antikaonsInside = false; <<
411 lambdasInside = false; <<
412 forcedDeltasInside = false; 327 forcedDeltasInside = false;
413 forcedDeltasOutside = false; 328 forcedDeltasOutside = false;
414 forcedPionResonancesOutside = false; <<
415 absorbedStrangeParticle = false; <<
416 forcedSigmaOutside = false; <<
417 forcedStrangeInside = false; <<
418 emitLambda = 0; <<
419 emitKaon = false; <<
420 clusterDecay = false; 329 clusterDecay = false;
421 firstCollisionTime = (Float_t)0.0; 330 firstCollisionTime = (Float_t)0.0;
422 firstCollisionXSec = (Float_t)0.0; 331 firstCollisionXSec = (Float_t)0.0;
423 firstCollisionSpectatorPosition = (Flo 332 firstCollisionSpectatorPosition = (Float_t)0.0;
424 firstCollisionSpectatorMomentum = (Flo 333 firstCollisionSpectatorMomentum = (Float_t)0.0;
425 firstCollisionIsElastic = false; 334 firstCollisionIsElastic = false;
426 nReflectionAvatars = 0; 335 nReflectionAvatars = 0;
427 nCollisionAvatars = 0; 336 nCollisionAvatars = 0;
428 nDecayAvatars = 0; 337 nDecayAvatars = 0;
429 nUnmergedSpectators = 0; 338 nUnmergedSpectators = 0;
430 nEnergyViolationInteraction = 0; 339 nEnergyViolationInteraction = 0;
>> 340 #ifdef INCL_INVERSE_KINEMATICS
431 341
>> 342 #endif
432 } 343 }
433 344
434 /// \brief Move a remnant to the particl 345 /// \brief Move a remnant to the particle array
435 void remnantToParticle(const G4int remna 346 void remnantToParticle(const G4int remnantIndex);
436 347
>> 348 #ifdef INCL_INVERSE_KINEMATICS
437 /// \brief Fill the variables describing 349 /// \brief Fill the variables describing the reaction in inverse kinematics
438 void fillInverseKinematics(const Double_ 350 void fillInverseKinematics(const Double_t gamma);
>> 351 #endif // INCL_INVERSE_KINEMATICS
439 }; 352 };
440 } 353 }
441 354
442 #endif /* G4INCLEVENTINFO_HH_HH */ 355 #endif /* G4INCLEVENTINFO_HH_HH */
443 356