1 | ------------------------------------------------------------------- |
---|
2 | |
---|
3 | ========================================================== |
---|
4 | Geant4 - an Object-Oriented Toolkit for Physics Simulation |
---|
5 | ========================================================== |
---|
6 | $Id: History,v 1.32 2010/10/29 06:54:23 gunter Exp $ |
---|
7 | --------------------------------------------------------------------- |
---|
8 | |
---|
9 | History file for the Liege cascade INCL Model |
---|
10 | --------------------------------------------- |
---|
11 | |
---|
12 | This file should be used to summarize modifications introduced in the |
---|
13 | code and to keep track of all tags. |
---|
14 | |
---|
15 | --------------------------------------------------------------- |
---|
16 | * Please list in reverse chronological order (last date on top) |
---|
17 | --------------------------------------------------------------- |
---|
18 | |
---|
19 | 29 October 2010 - Gunter Folger (hadr-incl-V09-03-05) |
---|
20 | - Fixed several more compilation warnings for conversion of double to int |
---|
21 | in G4Incl.cc, G4InclAblaCascadeInterface.cc, and G4InclAblaLightIonInterface.cc. |
---|
22 | |
---|
23 | 28 October 2010 - Gabriele Cosmo (hadr-incl-V09-03-04) |
---|
24 | ----------------------------------------------------- |
---|
25 | - Fixed compilation warning in G4InclDataDefs.hh for conversion of double |
---|
26 | to integer in initialization. |
---|
27 | |
---|
28 | 25 October 2010 - Pekka Kaitaniemi (hadr-incl-V09-03-03) |
---|
29 | -------------------------------------------------------- |
---|
30 | - Refactored INCL input data structure to its own class G4InclInput |
---|
31 | o Made data members private |
---|
32 | o Access to input data only through accessor methods |
---|
33 | o INCL now uses integer A and Z internally as well as in the interface |
---|
34 | |
---|
35 | 14 September 2010 - Pekka Kaitaniemi (hadr-incl-V09-03-02) |
---|
36 | ---------------------------------------------------------- |
---|
37 | - Migrated to integer A and Z when using G4Nucleus |
---|
38 | - Bugfixes: |
---|
39 | o Make sure we use Geant4 random number generators everywhere |
---|
40 | o Silenced some unnecessary warnings |
---|
41 | o Make sure we use INCL for all light ion projectiles up to Carbon-12 |
---|
42 | |
---|
43 | 16 August 2010 - Pekka Kaitaniemi |
---|
44 | --------------------------------- |
---|
45 | - Moved the INCL input initialization to the INCL input class |
---|
46 | - Implemented inverse-kinematics treatment for ion-hydrogen collisions |
---|
47 | |
---|
48 | 01 June 2010 - Pekka Kaitaniemi |
---|
49 | ------------------------------- |
---|
50 | - Improved the Fermi break-up configuration |
---|
51 | |
---|
52 | 26 February 2010 - Pekka Kaitaniemi |
---|
53 | ----------------------------------- |
---|
54 | - Merged INCL light ion projectile support to the integration branch |
---|
55 | o Support infrastructure for light ions up to Carbon |
---|
56 | o Use Geant4 Fermi break-up for projectile spectators and light cascade |
---|
57 | remnants |
---|
58 | |
---|
59 | 14 June 2010 - Gabriele Cosmo (hadr-incl-V09-03-01) |
---|
60 | ------------------------------------------------------ |
---|
61 | - Added missing virtual destructor to G4AblaFissionBase, |
---|
62 | G4InclAblaVirtualData and G4InclRandomNumbers. |
---|
63 | |
---|
64 | 27 April 2010 - Pekka Kaitaniemi (hadr-incl-V09-03-00) |
---|
65 | ------------------------------------------------------ |
---|
66 | - Fix: Nuclear deformation and shell effect data reader now constructs |
---|
67 | the corresponding data tables correctly. |
---|
68 | - Added some internal logging facilities to INCL (disabled by default) |
---|
69 | |
---|
70 | 09 December 2009 - Pekka Kaitaniemi (hadr-incl-V09-02-08) |
---|
71 | --------------------------------------------------------- |
---|
72 | - Fix: Added a safeguard against division by zero in INCL nucleon |
---|
73 | transmission probability calculation. |
---|
74 | |
---|
75 | 08 December 2009 - Pekka Kaitaniemi (hadr-incl-V09-02-07) |
---|
76 | --------------------------------------------------------- |
---|
77 | - Fix: Added a safeguard against division by zero in INCL |
---|
78 | initialization. |
---|
79 | |
---|
80 | 08 December 2009 - Pekka Kaitaniemi (hadr-incl-V09-02-06) |
---|
81 | --------------------------------------------------------- |
---|
82 | - Fix: Added a safeguard to avoid (hopefully rare) negative energies |
---|
83 | of INCL output particles. This would result in an attempt to take |
---|
84 | square root of a negative number. |
---|
85 | |
---|
86 | 04 December 2009 - Pekka Kaitaniemi (hadr-incl-V09-02-05) |
---|
87 | --------------------------------------------------------- |
---|
88 | - Fix: Corrected the G4InclAblaLightIonInterface error handling. In |
---|
89 | case INCL is unable to produce a valid cascade it will return the |
---|
90 | original projectile particle back to Geant4. Unfortunately in case |
---|
91 | of the light ion interface it returned a particle with NULL pointer |
---|
92 | to particle definition. |
---|
93 | - Added safeguards to both G4InclAblaCascadeInterface and |
---|
94 | G4InclAblaLightIonInterface to prevent returning of particles with |
---|
95 | NULL definition. |
---|
96 | |
---|
97 | 03 December 2009 - Pekka Kaitaniemi (hadr-incl-V09-02-04) |
---|
98 | --------------------------------------------------------- |
---|
99 | - Fix: Added a safeguard against attempts to compute the logarithm of zero in |
---|
100 | INCL |
---|
101 | |
---|
102 | 02 December 2009 - Pekka Kaitaniemi (hadr-incl-V09-02-03) |
---|
103 | --------------------------------------------------------- |
---|
104 | - Fix: Corrected datatype problem in INCL initialization |
---|
105 | |
---|
106 | 22 November 2009 - Pekka Kaitaniemi (hadr-incl-V09-02-02) |
---|
107 | --------------------------------------------------------- |
---|
108 | - Fix: Activated Geant4 random number generator instead of the built-in one. |
---|
109 | |
---|
110 | 18 November 2009 - Pekka Kaitaniemi (hadr-incl-V09-02-01) |
---|
111 | --------------------------------------------------------- |
---|
112 | - Added safeguard division by zero (or negative) energy in ABLA |
---|
113 | fission fragment handling |
---|
114 | - Minor fix in the INCL particle reflection time calculation |
---|
115 | |
---|
116 | 17 October 2009 - Pekka Kaitaniemi (hadr-incl-V09-02-00) |
---|
117 | -------------------------------------------------------- |
---|
118 | - Bugfixes to INCL nuclear potential handling: |
---|
119 | The interpolation function produced discontinuous results. This was |
---|
120 | due to an array index off-by-one bug. This bugfix corrects the |
---|
121 | impact parameter distribution. |
---|
122 | - Fixed several variable initialization issues in INCL. |
---|
123 | - Fixed severalFORTRAN to C++ translation issues in INCL. |
---|
124 | |
---|
125 | 06 November 2008 - Pekka Kaitaniemi (hadr-incl-V09-01-05) |
---|
126 | --------------------------------------------------------- |
---|
127 | - Fixed GCC 4.3.2 "src/G4Incl.cc:1712: warning: array subscript is above array |
---|
128 | bounds" problem. The problem occured because array size was used as an |
---|
129 | index. |
---|
130 | |
---|
131 | 06 November 2008 - Gabriele Cosmo (hadr-incl-V09-01-04) |
---|
132 | ------------------------------------------------------- |
---|
133 | - Fix for compilation errors on Windows-VC++; include <time.h> in |
---|
134 | G4AblaFission and G4AblaFissionSimfis18 translation units to make them |
---|
135 | self-consistent. |
---|
136 | |
---|
137 | 24 October 2008 - Dennis Wright (hadr-incl-V09-01-03) |
---|
138 | ----------------------------------------------------- |
---|
139 | - use GetNuclearMass() instead of GetAtomicMass() in G4AblaEvaporation.cc |
---|
140 | |
---|
141 | 10 October 2008 - Pekka Kaitaniemi |
---|
142 | ------------------------------------ |
---|
143 | - Added ability to use Geant4 Fermi break-up |
---|
144 | for light cascade remnants |
---|
145 | |
---|
146 | 12 September 2008 - Pekka Kaitaniemi (hadr-incl-V09-01-02) |
---|
147 | ------------------------------------ |
---|
148 | - Bug fixes in INCL cascade: |
---|
149 | o Fixed target nucleus initialization for targets with A >= 0.19 and |
---|
150 | A <= 28. |
---|
151 | o Fixed target nucleus radius calculation for targets 6 <= A < 19. |
---|
152 | o Fixed conditional statements to avoid uninitialized variables. |
---|
153 | o Several array index handling (off-by-one) bugs fixed. |
---|
154 | o Several possible double/integer datatype problems fixed. |
---|
155 | |
---|
156 | 25 June 2008 - Pekka Kaitaniemi (hadr-incl-V09-01-01) |
---|
157 | ------------------------------ |
---|
158 | - Portability fixes and code clean-up: |
---|
159 | o Removed all assert() calls. |
---|
160 | o Removed all isnan() and isinf() calls. |
---|
161 | o All math functions now use std:: namespace. |
---|
162 | o Silenced the fission models. |
---|
163 | |
---|
164 | 26 May 2008 - Pekka Kaitaniemi (hadr-incl-V09-01-00) |
---|
165 | ------------------------------ |
---|
166 | - Bug fixes: |
---|
167 | o Fixed serious data reading bug affecting the reading of data |
---|
168 | files: flalpha.dat, frldm.dat and vgsld.dat. |
---|
169 | - Fission reorganization: |
---|
170 | o Two fission models are provided. By default GSI SimFis3 fission is |
---|
171 | used. Optional SimFis18 is provided. |
---|
172 | |
---|
173 | 22 April 2008 - Pekka Kaitaniemi |
---|
174 | ----------------------------------- |
---|
175 | - INCL/ABLA changes: |
---|
176 | o Fix: Partial fix to the Alpha evaporation energies. |
---|
177 | o Fix: Neutron energies. |
---|
178 | o Fix: Moved remnant direction calculation block and added safeguards. |
---|
179 | o ABLA bugfix: switched the transformation matrix indices. |
---|
180 | |
---|
181 | 27 February 2008 - Pekka Kaitaniemi |
---|
182 | ----------------------------------- |
---|
183 | - INCL/ABLA changes since Geant4 9.1: |
---|
184 | o Initialized variables to zero and improved coding style. |
---|
185 | o Variable "bet" needs to be part of G4Fiss. |
---|
186 | o Fixed a variable definition. Variable "homega" is part of "G4Fiss". |
---|
187 | o Fixed an index off-by-one bug in G4Abla. |
---|
188 | o Fix an infinite loop. |
---|
189 | o Move the stored event number to the class definition. |
---|
190 | o Added a possible (commented) fix for ABLA output off-by-one problem. |
---|
191 | o Cleaned up G4Abla::qrot comments. |
---|
192 | - Changes to the test suite since Geant4 9.1: |
---|
193 | o Added more tests. |
---|
194 | o New test system. |
---|
195 | o Compare FORTRAN and C++ results. |
---|
196 | o Preliminary GUI for running plotting scripts. |
---|
197 | |
---|
198 | 04 December 2007 - Gabriele Cosmo (hadr-incl-V09-00-09) |
---|
199 | --------------------------------- |
---|
200 | - Again commented out call to isnan() in G4Incl.cc for allow |
---|
201 | porting to platforms other than Linux.... |
---|
202 | |
---|
203 | 03 December 2007 - Pekka Kaitaniemi |
---|
204 | ----------------------------------- |
---|
205 | - Added bugfixes to INCL nucleon array handling. Several off-by-one |
---|
206 | bugs were fixed. These fixes help with INCL neutron energy spectra. |
---|
207 | - 4-momentum conservation problem was noticed to be mainly ABLA |
---|
208 | related. |
---|
209 | o INCL alone violates energy conservation only by less than 20 MeV. |
---|
210 | o ABLA has problems with baryon number conservation. Some produced |
---|
211 | nucleons are not given back to INCL and Geant4. This is probably the |
---|
212 | main reason for 4-momentum conservation violations. |
---|
213 | - Added 4-momentum and baryon number conservation checking code to the |
---|
214 | interfaces (behind flag DEBUGINCL). |
---|
215 | |
---|
216 | 15 November 2007 - Pekka Kaitaniemi |
---|
217 | ----------------------------------- |
---|
218 | - First fix to secondary particle Pt, Px and Phi distributions. |
---|
219 | o Momentum x-component was transformed incorrectly from spherical to |
---|
220 | cartesian coordinates. |
---|
221 | |
---|
222 | 08 November 2007 - Pekka Kaitaniemi |
---|
223 | ----------------------------------- |
---|
224 | - Started using G4UniformRand as a random number generator. |
---|
225 | o It is much more sophisticated generator than the original INCL one |
---|
226 | was. |
---|
227 | o This should eliminate infinite loop in test runs on 64-bit |
---|
228 | plattforms. The old generator probably was not 64-bit clean. |
---|
229 | |
---|
230 | 31 October 2007 - Gabriele Cosmo (hadr-incl-V09-00-05) |
---|
231 | -------------------------------- |
---|
232 | - Again commented out calls to isnan() and isinf() in G4Abla.cc for allow |
---|
233 | porting to platforms other than Linux. |
---|
234 | |
---|
235 | 30 October 2007 - Pekka Kaitaniemi |
---|
236 | ---------------------------------- |
---|
237 | - Abla now runs. Integrated energy spectrums for protons and neutrons |
---|
238 | seem to be OK. Isotope production is still problematic. There seems |
---|
239 | to be too much fission. |
---|
240 | |
---|
241 | 24 October 2007 - Pekka Kaitaniemi (hadr-incl-V09-00-04) |
---|
242 | ---------------------------------- |
---|
243 | - Modified INCL to filter out unphysical events containing particles |
---|
244 | that have more energy than the energy of the projectile. This should |
---|
245 | fix problem we had with outcoming neutrons (and protons) that have |
---|
246 | more energy than the projectile particle. |
---|
247 | - Added G4InclAblaDataFile class that reads data used by INCL/ABLA |
---|
248 | from datafiles located in $G4ABLADATA. If this environment variable |
---|
249 | is not correctly set INCL/ABLA returns G4Exception. |
---|
250 | - Problematic hardcoded data has been removed. |
---|
251 | |
---|
252 | 16 October 2007 - Pekka Kaitaniemi |
---|
253 | ------------------------------------------------------ |
---|
254 | - Silenced the rather verbose output. Now nothing should be printed on |
---|
255 | terminal by default. |
---|
256 | - Fixed INCL internal random number generator bug that sometimes |
---|
257 | (rarely) produced random numbers that were greater than 1.0. |
---|
258 | - Random number seeds are now taken from CLHEP::HepRandom. |
---|
259 | - Cleaned up some commented code. |
---|
260 | - Found new NaN bugs in ABLA. |
---|
261 | |
---|
262 | 13 October 2007 - Pekka Kaitaniemi |
---|
263 | ------------------------------------------------------ |
---|
264 | - Interface cleanup. G4InclCascadeInterface and |
---|
265 | G4InclLightIonInterface should be fully operational. |
---|
266 | - INCL was reverted back to using simple random number generator for |
---|
267 | test result compatibility (to produce exactly same results as the |
---|
268 | FORTRAN version does). |
---|
269 | - Fixes for ABLA evaporation/fission: |
---|
270 | o Fission barrier problem fixed. Now fission code is called during |
---|
271 | the run. |
---|
272 | o Several NaN bugs fixed. |
---|
273 | o Random seed initialization bug fixed. |
---|
274 | o Stand-alone INCL/ABLA program (test/standaloneinclabla) mostly |
---|
275 | works (produces results that are sometimes in good agreement with |
---|
276 | FORTRAN version and sometimes in not so good agreement). With |
---|
277 | heavy targets and lots of energy fission is overestimated. |
---|
278 | o INCL+ABLA interface not yet fully operational. G4 simulation run |
---|
279 | crashes. |
---|
280 | |
---|
281 | 11 October 2007 - Gabriele Cosmo (hadr-incl-V09-00-01) |
---|
282 | ------------------------------------------------------ |
---|
283 | - Fixes for porting: |
---|
284 | o Commented out usage of isnan() and isinf() in assertions. |
---|
285 | o Replaced usage of M_PI with CLHEP::pi. |
---|
286 | o Use std namespace for standard mathematical functions. |
---|
287 | o Casted to double arguments used in std::pow() function. |
---|
288 | o Commented out initialisation of data in G4InclAblaHardcodedData.cc |
---|
289 | for Windows; temporary solution to allow for compilation, pending |
---|
290 | solution of loading of data from file ! |
---|
291 | |
---|
292 | 23 May 2007 - Aatos Heikkinen (hadr-incl-V09-00-00) |
---|
293 | --------------------------------------------------- |
---|
294 | - First version of the Liege cascade INCL Model. |
---|
295 | |
---|
296 | |
---|