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