| 1 | ------------------------------------------------------------------- |
|---|
| 2 | |
|---|
| 3 | ================================================== |
|---|
| 4 | Geant4 - an Object-Oriented Toolkit for Simulation |
|---|
| 5 | ================================================== |
|---|
| 6 | |
|---|
| 7 | History file for Low Energy Parameterized Models |
|---|
| 8 | ------------------------------------------------ |
|---|
| 9 | |
|---|
| 10 | This file should be used to summarize modifications introduced in the |
|---|
| 11 | code and to keep track of all tags. |
|---|
| 12 | |
|---|
| 13 | --------------------------------------------------------------- |
|---|
| 14 | * Please list in reverse chronological order (last date on top) |
|---|
| 15 | --------------------------------------------------------------- |
|---|
| 16 | |
|---|
| 17 | 27 November 2008 - V.Ivanchenko (hadr-cohe-V09-01-08) |
|---|
| 18 | ---------------------------------------------------- |
|---|
| 19 | G4ChargeExchangeProcess - added high and low energy protections into |
|---|
| 20 | cross section |
|---|
| 21 | |
|---|
| 22 | 20 November 2008 - V.Ivanchenko (hadr-cohe-V09-01-07) |
|---|
| 23 | ---------------------------------------------------- |
|---|
| 24 | G4ChargeExchangeProcess and G4ChargeExchange model prepaed to be used |
|---|
| 25 | in Physics Lists; energy threshold is set to 1.0 MeV |
|---|
| 26 | |
|---|
| 27 | 22 August 2008 - V.Ivanchenko (hadr-cohe-V09-01-06) |
|---|
| 28 | ---------------------------------------------------- |
|---|
| 29 | - set process sub-types according to enumerator for G4UHadronElasticProcess |
|---|
| 30 | and for G4ChargeExchangeProcess |
|---|
| 31 | |
|---|
| 32 | 05 August 2008 - V.Ivanchenko (hadr-cohe-V09-01-05) |
|---|
| 33 | ---------------------------------------------------- |
|---|
| 34 | G4HadronElastic - added protection inside LElastic model against rare |
|---|
| 35 | unphysical values of sampled t |
|---|
| 36 | |
|---|
| 37 | 08 Jul 2008 - D.H. Wright (hadr-cohe-V09-01-04) |
|---|
| 38 | ------------------------------------------------ |
|---|
| 39 | - set process sub-type to 11 for G4UHadronElasticProcess and 16 for |
|---|
| 40 | G4ChargeExchangeProcess |
|---|
| 41 | |
|---|
| 42 | 19 May 2008 - V.Ivanchenko (hadr-cohe-V09-01-03) |
|---|
| 43 | ------------------------------------------------ |
|---|
| 44 | G4HadronElastic - added counts of abnornal sampling results and cout |
|---|
| 45 | at the destructor of the class |
|---|
| 46 | G4UHadronElasticProcess and G4ChargeExchangeProcess - added call to |
|---|
| 47 | BuildPhysicsTable of the base class G4HadronicProcess |
|---|
| 48 | |
|---|
| 49 | 10 Apr 2008 - V.Ivanchenko (hadr-cohe-V09-01-02) |
|---|
| 50 | ------------------------------------------------ |
|---|
| 51 | G4HadronElastic - added extra protection against precision lost |
|---|
| 52 | when projectile almostly lost kinetic energy; |
|---|
| 53 | set lowestEnergyLimit=1.e-6*eV (was 0.0) |
|---|
| 54 | |
|---|
| 55 | 29 Jan 2008 - V.Ivanchenko (hadr-cohe-V09-01-01) |
|---|
| 56 | ------------------------------------------------- |
|---|
| 57 | G4UHadronElasticProcess - call to G4QElasticCrossSection with option |
|---|
| 58 | "false", so parameters for sampling of |
|---|
| 59 | scattering are precalculated |
|---|
| 60 | G4HadronElastic - added extra protection - if sampled t out of the |
|---|
| 61 | physics range then no scattering |
|---|
| 62 | |
|---|
| 63 | 16 Jan 2008 - V.Ivanchenko (hadr-cohe-V09-01-00) |
|---|
| 64 | ------------------------------------------------ |
|---|
| 65 | G4DiffuseElastic - remove dependence on other elastic models |
|---|
| 66 | G4ElasticHadrNucleusHE - fixed non-reproducibility, removed commented |
|---|
| 67 | part of old code; add more formatting |
|---|
| 68 | G4HadronElastic - return back threshold to G4ElasticHadrNucleusHE |
|---|
| 69 | model 1 GeV |
|---|
| 70 | |
|---|
| 71 | 10 Dec 2007 - V.Ivanchenko (hadr-cohe-V09-00-07) |
|---|
| 72 | ------------------------------------------------ |
|---|
| 73 | G4HadronElastic - set threshold to G4ElasticHadrNucleusHE model |
|---|
| 74 | DBL_MAX |
|---|
| 75 | |
|---|
| 76 | 15 Nov 2007 - V.Ivanchenko (hadr-cohe-V09-00-06) |
|---|
| 77 | ------------------------------------------------ |
|---|
| 78 | G4UHadronElasticProcess - introduced "ifdef VERBOSE" for verbosity |
|---|
| 79 | checks inside methods called at each step |
|---|
| 80 | |
|---|
| 81 | 14 Nov 2007 - V.Ivanchenko (hadr-cohe-V09-00-04) |
|---|
| 82 | ------------------------------------------------ |
|---|
| 83 | G4DiffuseElastic - add control on verbosity at initialisation |
|---|
| 84 | |
|---|
| 85 | 13 Nov 2007 - V.Ivanchenko (hadr-cohe-V09-00-03) |
|---|
| 86 | ------------------------------------------------ |
|---|
| 87 | G4ElasticHadrNucleusHE - (N.Starkov) updated: |
|---|
| 88 | - the difference in scattering off protons and |
|---|
| 89 | neutrons is included; |
|---|
| 90 | - optized computation of the distribution function |
|---|
| 91 | (important for heavy targets); |
|---|
| 92 | - fixed model for scattering off hydrogen; |
|---|
| 93 | - tunned model parameters for better description |
|---|
| 94 | of total and elastic cross sections |
|---|
| 95 | |
|---|
| 96 | 06 Nov 2007 - V.Ivanchenko (hadr-cohe-V09-00-02) |
|---|
| 97 | ------------------------------------------------ |
|---|
| 98 | G4DiffuseElastic - (V.Grichine) initialisation of the model for |
|---|
| 99 | each new target nucleus on-fly |
|---|
| 100 | |
|---|
| 101 | 05 Nov 2007 - V.Ivanchenko (hadr-cohe-V09-00-01) |
|---|
| 102 | ------------------------------------------------ |
|---|
| 103 | |
|---|
| 104 | G4DiffuseElastic - (V.Grichine) initialisation of the model in |
|---|
| 105 | constructor with particle definition for the |
|---|
| 106 | vector of available elements; tuning |
|---|
| 107 | at high energies. |
|---|
| 108 | |
|---|
| 109 | 17 Oct 2007 - V.Ivanchenko (hadr-cohe-V09-00-00) |
|---|
| 110 | ------------------------------------------------ |
|---|
| 111 | G4LEnp - removed erroneous code for identity |
|---|
| 112 | exchange of particles ( F.W. Jones) |
|---|
| 113 | |
|---|
| 114 | G4DiffuseElastic - added Coulomb elastic scattering amplitude |
|---|
| 115 | using Wentzel model (V. Grichine) |
|---|
| 116 | |
|---|
| 117 | 14 June 2007 - V.Ivanchenko (hadr-cohe-V08-03-07) |
|---|
| 118 | ------------------------------------------------ |
|---|
| 119 | G4UHadronElasticProcess - use "true" flag colling CHIPS to compute cross |
|---|
| 120 | section |
|---|
| 121 | |
|---|
| 122 | 14 June 2007 - V.Ivanchenko (hadr-cohe-V08-03-06) |
|---|
| 123 | ----------------------------------------------- |
|---|
| 124 | G4ElasticHadrNucleusHE - fixed compilation nuclear parameters for heavy |
|---|
| 125 | nucleus |
|---|
| 126 | G4DiffuseElastic - new class |
|---|
| 127 | |
|---|
| 128 | 31 May 2007 - V.Ivanchenko (hadr-cohe-V08-03-05) |
|---|
| 129 | ----------------------------------------------- |
|---|
| 130 | G4ElasticHadrNucleusHE - fixed compilation problem at Windows |
|---|
| 131 | use only double values inside std::pow |
|---|
| 132 | |
|---|
| 133 | 31 May 2007 - V.Ivanchenko (hadr-cohe-V08-03-04) |
|---|
| 134 | ----------------------------------------------- |
|---|
| 135 | G4ElasticHadrNucleusHE - fixed compilation problem at SUN |
|---|
| 136 | move static const out of switch |
|---|
| 137 | |
|---|
| 138 | 25 May 2007 - D.H. Wright (hadr-cohe-V08-03-03) |
|---|
| 139 | ----------------------------------------------- |
|---|
| 140 | add model names to constructors in G4LEnp, G4LEpp, G4ChargeExchange.cc, |
|---|
| 141 | G4HadronElastic, and G4ElasticHadrNucleusHE |
|---|
| 142 | |
|---|
| 143 | 24 May 2007 - V.Ivanchenko (hadr-cohe-V08-03-02) |
|---|
| 144 | ----------------------------------------------- |
|---|
| 145 | G4ElasticHadrNucleusHE - optimised computation of tables |
|---|
| 146 | G4HadronElastic - use HElastic only for mesons |
|---|
| 147 | |
|---|
| 148 | 15 May 2007 - V.Ivanchenko (hadr-cohe-V08-03-01) |
|---|
| 149 | ----------------------------------------------- |
|---|
| 150 | G4HadronElastic - use QElasticCrossSection for p, n incident for any energy |
|---|
| 151 | and any target; use HElastic for the rest |
|---|
| 152 | |
|---|
| 153 | 15 May 2007 - V.Ivanchenko (hadr-cohe-V08-03-00) |
|---|
| 154 | ----------------------------------------------- |
|---|
| 155 | G4ElasticHadrNucleusHE - general cleanup and redesign of private methods |
|---|
| 156 | |
|---|
| 157 | 5 May 2007 - V.Ivanchenko (hadr-cohe-V08-02-06) |
|---|
| 158 | ----------------------------------------------- |
|---|
| 159 | G4HadronElastic - do not use HE model for hydrogen target to avoid NaN; |
|---|
| 160 | use QElasticCrossSection for p, n incident for any energy |
|---|
| 161 | only for p, He targets |
|---|
| 162 | |
|---|
| 163 | 4 May 2007 - V.Ivanchenko (hadr-cohe-V08-02-05) |
|---|
| 164 | ----------------------------------------------- |
|---|
| 165 | G4HadronElastic - do not use HE model for hydrogen target to avoid NaN; |
|---|
| 166 | use QElasticCrossSection for p, n incident for any energy |
|---|
| 167 | |
|---|
| 168 | 2 May 2007 - V.Ivanchenko (hadr-cohe-V08-02-04) |
|---|
| 169 | ----------------------------------------------- |
|---|
| 170 | G4ElasticHadrNucleusHE - scale sampled t as ptot^2 |
|---|
| 171 | G4HadronElastic - use QElasticCrossSection for p, n < GeV and |
|---|
| 172 | G4ElasticHadrNucleusHE model above |
|---|
| 173 | still use QElasticCrossSEction for H and He targets |
|---|
| 174 | |
|---|
| 175 | 2 Apr 2007 - V.Ivanchenko (hadr-cohe-V08-02-03) |
|---|
| 176 | ----------------------------------------------- |
|---|
| 177 | G4ElasticHadrNucleusHE - use iterator search on bin instead of |
|---|
| 178 | computing log |
|---|
| 179 | G4HadronElastic - use the most probable isotope for given Z to call |
|---|
| 180 | to HE model (CPU/memory more efficient) |
|---|
| 181 | |
|---|
| 182 | 30 Jan 2007 - V.Ivanchenko (hadr-cohe-V08-02-02) |
|---|
| 183 | ------------------------------------------------ |
|---|
| 184 | G4UHadronElasticProcess - fix bug in number of neutrons in CHIPS interface |
|---|
| 185 | G4HadronElastic, G4ChargeExchange - more accuratelly compute Z and A |
|---|
| 186 | |
|---|
| 187 | 30 Jan 2007 - V.Ivanchenko (hadr-cohe-V08-02-01) |
|---|
| 188 | ------------------------------------------------ |
|---|
| 189 | G4UHadronElasticProcess, G4ChargeExchangeProcess - fix problem |
|---|
| 190 | for the case, when isotope vector is not defined |
|---|
| 191 | |
|---|
| 192 | 23 Jan 2007 - V.Ivanchenko (hadr-cohe-V08-02-00) |
|---|
| 193 | ------------------------------------------------ |
|---|
| 194 | G4UHadronElasticProcess - migrade to interface to cross section with Z and A |
|---|
| 195 | G4ChargeExchangeProcess - migrade to interface to cross section with Z and A |
|---|
| 196 | do not use QElastic total cross sections |
|---|
| 197 | |
|---|
| 198 | 23 Nov 2006 - V.Ivanchenko (hadr-cohe-V08-01-12) |
|---|
| 199 | ------------------------------------------------ |
|---|
| 200 | G4ElasticHadrNucleusHE - remove unused variables; |
|---|
| 201 | - set ONQ0=3 (instead of 5) as it was in the |
|---|
| 202 | tag 05 in order to have better CPU performance |
|---|
| 203 | - cleanup comments |
|---|
| 204 | |
|---|
| 205 | 22 Nov 2006 - V.Ivanchenko (hadr-cohe-V08-01-11) |
|---|
| 206 | ------------------------------------------------ |
|---|
| 207 | G4HadronElastic: set DBL_MAX low limit for the HE model |
|---|
| 208 | |
|---|
| 209 | 21 Nov 2006 - V.Ivanchenko (hadr-cohe-V08-01-10) |
|---|
| 210 | ------------------------------------------------ |
|---|
| 211 | Cleanup of G4ElasticHadrNucleusHE |
|---|
| 212 | |
|---|
| 213 | 18 Nov 2006 - V.Ivanchenko (hadr-cohe-V08-01-09) |
|---|
| 214 | ------------------------------------------------ |
|---|
| 215 | N.Starkov: Fix problem of |
|---|
| 216 | |
|---|
| 217 | 17 Nov 2006 - V.Ivanchenko (hadr-cohe-V08-01-08) |
|---|
| 218 | ------------------------------------------------ |
|---|
| 219 | N.Starkov: Fix problem of initialization of tables, |
|---|
| 220 | Fix problem of CPU overhead |
|---|
| 221 | |
|---|
| 222 | 17 Nov 2006 - V.Ivanchenko (hadr-cohe-V08-01-07) |
|---|
| 223 | ------------------------------------------------ |
|---|
| 224 | N.Starkov: Code optimization is perfomed |
|---|
| 225 | |
|---|
| 226 | 16 Nov 2006 - V.Ivanchenko (hadr-cohe-V08-01-06) |
|---|
| 227 | ------------------------------------------------ |
|---|
| 228 | G4HadronElastic - set threshold 0.4 GeV for HE model (VI) |
|---|
| 229 | G4ElasticHadrNucleusHE - add corrections allowing to work down |
|---|
| 230 | to 0.4 GeV (N.Starkov) |
|---|
| 231 | |
|---|
| 232 | 24 Oct 2006 - V.Ivanchenko (hadr-cohe-V08-01-05) |
|---|
| 233 | ------------------------------------------------ |
|---|
| 234 | Set lowestEnergy=0 for neutral particles, =eV for changed (VI) |
|---|
| 235 | G4HadronElastic - set threshold 1 GeV for G4ElasticHadrNucleusHE model (VI) |
|---|
| 236 | G4ElasticHadrNucleusHE - complete review and update of the class, |
|---|
| 237 | data tables are computed on-flight, |
|---|
| 238 | G4ELASTICDATA are not used anymore, |
|---|
| 239 | the model is applicable for all nuclei, except |
|---|
| 240 | deuteron, triton, He3, and A>238 (N.Starkov) |
|---|
| 241 | |
|---|
| 242 | 26 Sept 2006 - V.Ivanchenko (hadr-cohe-V08-01-04) |
|---|
| 243 | ------------------------------------------------- |
|---|
| 244 | Introduce extra protection - below lowestEnergy (10-6 eV) not active |
|---|
| 245 | |
|---|
| 246 | 10 August 2006 - V.Ivanchenko (hadr-cohe-V08-01-03) |
|---|
| 247 | --------------------------------------------------- |
|---|
| 248 | Introduce NaN check inside G4HadronElastic, G4ChargeExchange, |
|---|
| 249 | and G4UHadronElasticProcess; |
|---|
| 250 | |
|---|
| 251 | 24 July 2006 - V.Ivanchenko (hadr-cohe-V08-01-02) |
|---|
| 252 | ----------------------------------------------- |
|---|
| 253 | Remove HP data from G4UHadronElasticProcess (it is set in physics lists) |
|---|
| 254 | Add low energy threshold for QElastic 20 MeV inside G4HadronElastic |
|---|
| 255 | |
|---|
| 256 | 13 July 2006 - V.Ivanchenko (hadr-cohe-V08-01-01) |
|---|
| 257 | ----------------------------------------------- |
|---|
| 258 | Fix bug: remove 1 keV limit from G4UHadronElasticProcess |
|---|
| 259 | |
|---|
| 260 | 07 July 2006 - V.Ivanchenko (hadr-cohe-V08-01-00) |
|---|
| 261 | ----------------------------------------------- |
|---|
| 262 | Fix bug: the rotation of final state need to be done in PostStep of |
|---|
| 263 | the process, not on the model level (V.Ivanchenko) |
|---|
| 264 | |
|---|
| 265 | 21 June 2006 - V.Ivanchenko (hadr-cohe-V08-00-09) |
|---|
| 266 | ----------------------------------------------- |
|---|
| 267 | Fix - remove nan values n scattering on deuteron (V.Ivanchenko) |
|---|
| 268 | |
|---|
| 269 | 17 June 2006 - V.Ivanchenko (ghad-cohe-V08-00-08) |
|---|
| 270 | ----------------------------------------------- |
|---|
| 271 | Use sampling of final state for n, p scattering on deuteron |
|---|
| 272 | and alpha from LElastic (V.Ivanchenko) |
|---|
| 273 | |
|---|
| 274 | 7 June 2006 - V.Ivanchenko (ghad-cohe-V08-00-07) |
|---|
| 275 | ----------------------------------------------- |
|---|
| 276 | Fix complilation warnings at WINDOWS (V.Ivanchenko) |
|---|
| 277 | |
|---|
| 278 | 2 June 2006 - V.Ivanchenko (ghad-cohe-V08-00-06) |
|---|
| 279 | ----------------------------------------------- |
|---|
| 280 | Fix complilation warnings at SUN (V.Ivanchenko) |
|---|
| 281 | |
|---|
| 282 | 1 June 2006 - V.Ivanchenko (ghad-cohe-V08-00-05) |
|---|
| 283 | ----------------------------------------------- |
|---|
| 284 | Add regime of HElastic which is not using data files (N.Starkov) |
|---|
| 285 | Cleanup (V.Ivanchenko) |
|---|
| 286 | |
|---|
| 287 | 30 May 2006 - V.Ivanchenko (ghad-cohe-V08-00-04) |
|---|
| 288 | ----------------------------------------------- |
|---|
| 289 | Add ChargeExchangeProcess, cleanup G4HadronElastic |
|---|
| 290 | |
|---|
| 291 | 2 May 2006 - V.Ivanchenko (ghad-cohe-V08-00-03) |
|---|
| 292 | ----------------------------------------------- |
|---|
| 293 | Rename G4LElasticB -> G4HadronElastic |
|---|
| 294 | |
|---|
| 295 | 25 April 2006 - V.Ivanchenko (ghad-cohe-V08-00-02) |
|---|
| 296 | ------------------------------------------------- |
|---|
| 297 | G4UHadronElasticProcess and G4LElasticB: use np scattering from CHIPS, |
|---|
| 298 | fix run time warnings |
|---|
| 299 | Rename G4LElasticB by G4HadronElastic |
|---|
| 300 | |
|---|
| 301 | |
|---|
| 302 | 24 April 2006 - V.Ivanchenko (ghad-cohe-V08-00-01) |
|---|
| 303 | ------------------------------------------------- |
|---|
| 304 | G4UHadronElasticProcess - new process allows to use G4Isotopes composition |
|---|
| 305 | of elements defined by user and uses cross sections from |
|---|
| 306 | CHIPS and HP models when appropriate |
|---|
| 307 | |
|---|
| 308 | 13 April 2006 - V.Ivanchenko (ghad-cohe-V08-00-00) |
|---|
| 309 | ------------------------------------------------- |
|---|
| 310 | G4LElasticB - moved from low_energy; bug fixes - remove charge |
|---|
| 311 | exchange option, remove low limit on primary momentum; |
|---|
| 312 | add limit on secondary kinetic energy; add only s-wave |
|---|
| 313 | regime for low momentum scattering |
|---|
| 314 | GNUmakfile - introduce dependence on chips |
|---|
| 315 | |
|---|
| 316 | 28 Nov 2005 - G. Cosmo (hadr-cohe-V07-01-05) |
|---|
| 317 | --------------------------------------------- |
|---|
| 318 | Migrated to <sstream> from G4ElasticHadrNucleusHE. |
|---|
| 319 | |
|---|
| 320 | 25 Nov 2005 - D. Wright (hadr-cohe-V07-01-04) |
|---|
| 321 | --------------------------------------------- |
|---|
| 322 | From Solaris testing: line 83 of G4DiffElasticHadrScattering.hh: |
|---|
| 323 | multiple declaration of variable c1. Fixed. |
|---|
| 324 | |
|---|
| 325 | 25 Nov 2005 - D. Wright (hadr-cohe-V07-01-03) |
|---|
| 326 | --------------------------------------------- |
|---|
| 327 | Add detailed comments of N. Starkov from previous tags to History |
|---|
| 328 | file: |
|---|
| 329 | G4DiffElasticHadrNucleus.hh, .cc : |
|---|
| 330 | add high energy elastic scattering on protons |
|---|
| 331 | add new method HadronProtonDiffCrSec |
|---|
| 332 | |
|---|
| 333 | G4ElasticHadrNucleusHE.hh, .cc : |
|---|
| 334 | add high energy elastic scattering on protons |
|---|
| 335 | |
|---|
| 336 | G4HadronValues.hh, .cc : |
|---|
| 337 | use PDG value instead of static particle pointers |
|---|
| 338 | |
|---|
| 339 | From V. Ivantchenko: fixed memory leaks |
|---|
| 340 | 4-momentum balance fixed, tested OK |
|---|
| 341 | fixed compiler warnings |
|---|
| 342 | |
|---|
| 343 | 25 Nov 2005 - V.Ivanchenko (ghad-cohe-V07-01-03) |
|---|
| 344 | ------------------------------------------------ |
|---|
| 345 | Integration of recent development of the subpackage |
|---|
| 346 | Remove of the G4Rutherford for the time being |
|---|
| 347 | |
|---|
