source: trunk/source/processes/hadronic/models/pre_equilibrium/History @ 1253

Last change on this file since 1253 was 1228, checked in by garnier, 15 years ago

update geant4.9.3 tag
File size: 10.8 KB
Line 
1-------------------------------------------------------------------
2
3     ==========================================================
4     Geant4 - an Object-Oriented Toolkit for Physics Simulation
5     ==========================================================
6
7          History file for hadronic/models/de_excitation
8          ----------------------------------------------
9
10This file should be used to summarize modifications introduced in the
11code and to keep track of all tags.
12
13     ---------------------------------------------------------------
14     * Please list in reverse chronological order (last date on top)
15     ---------------------------------------------------------------
16
1721-November 2009 V.Ivanchenko   hadr-pre-V09-02-07
18---------------------------------------------------
19G4PreCompoundTransitions - return back transition probability
20
2119-November 2009 V.Ivanchenko   hadr-pre-V09-02-06
22---------------------------------------------------
23G4PreCompoundModel       - (JMQ) CEM transition probabilities are set
24                                 as a default
25G4PreCompoundTransitions - (JMQ) CEM transition probabilities have been
26                                 renormalized (transitions with Delta N=+2
27                                 are increased by factor 5)
28
2912-November 2009 G.Cosmo        hadr-pre-V09-02-05
30---------------------------------------------------
31G4PrecompoundEmission - fixed compilation error on Windows.
32
3312-November 2009 G.Folger       hadr-pre-V09-02-04
34---------------------------------------------------
35G4PrecompoundEmission - modifiy AngularDistribution() to avoid floating point
36                        exceptions.
37                      - remove calls to log in rho() replaced by use of new
38                        logfactorial()
39                               
4013-February 2009 V.Ivanchenko   hadr-pre-V09-02-03
41---------------------------------------------------
42G4PreCompoundXXX - changed the shape of probabilities (return back to 9.2) 
43                   for d, t, He3, alpha (JMQ)
44
4512-February 2009 V.Ivanchenko   hadr-pre-V09-02-02
46---------------------------------------------------
47G4PreCompoundXXX - changed the shape of probabilities
48                   for d, t, He3, alpha near the Coulomb barrier (JMQ)
49
5011-February 2009 V.Ivanchenko   hadr-pre-V09-02-01
51---------------------------------------------------
52G4PreCompoundXXX - set default Opt3 back, add decrease Coulomb barrier
53                   for d, t, a, he3 (JMQ)
54
5510-February 2009 V.Ivanchenko   hadr-pre-V09-02-00
56---------------------------------------------------
57Some clean up of comments
58G4PreCompoundIon - fixed probability of light ion emmision (JMQ)
59G4PreCompoundXXX - by default Opt1 is used for d, t, a, he3,
60                   Opt3 for n, p (JMQ)
61
6209-December 2008 A.Howard   hadr-pre-V09-01-15
63---------------------------------------------------
64Added protection for close to zero excitation energy in
65G4PreCompoundModel.cc (according to JMQ, MAC) to not try to de-excite.
66
6709-December 2008 A.Howard   hadr-pre-V09-01-14
68---------------------------------------------------
69Added protection for close to zero excitation energy in
70G4PreCompoundTransitions.cc (returns 0.0), prevents FPE later on.
71
7227-November 2008 A.Howard   hadr-pre-V09-01-13
73---------------------------------------------------
74Added data member initialisation to G4VPreCompoundTransitions.hh
75
7619-November 2008 A.Howard   hadr-pre-V09-01-12
77---------------------------------------------------
78JMQ fix to G4PreCompoundNeutron.cc for Zirconium (if (nu < 0.)nu=-nu).
79
8030-September 2008 A.Howard   hadr-pre-V09-01-11
81---------------------------------------------------
82Trivial protection against negative probabilities for incident protons on
83targets with A < Carbon.
84
8522-September 2008 A.Howard   hadr-pre-V09-01-10
86---------------------------------------------------
87JMQ's latest developments - which are an extension of hadr-pre-V09-01-08 with
88cross-section options (equivalent to the de-excitation modifications) and other
89options (SICB, Never Go Back,...).  In addition Gunter's fix to the factorials
90is included from hadr-pre-V09-01-09.
91Added new files: G4PreCompoundAlpha.cc G4PreCompoundDeuteron.cc
92G4PreCompoundHe3.cc G4PreCompoundNeutron.cc G4PreCompoundProton.cc
93G4PreCompoundTriton.cc (Classes previously existed.  Now they have more
94involved cross-section calculations).
95
9619-September 2008 A.Howard   hadr-pre-V09-01-09
97---------------------------------------------------
98Including Gunter's fix (see 11-August below) on top of ref-07
99(hadr-pre-V09-00-04), i.e. without JMQ's developments.
100
10111-August 2008 G.Folger     
102-----------------------------------------------
103Rewrite algorithm in G4PreCompoundEmission::rho() to avoid frequent
104floating point overflow when using Precompund in combination
105with Binary Cascade.
106
107
10824 July 2008 J. M. Quesada   hadr-pre-V09-01-08
109---------------------------------------------------
110Minor fixes and Coulomb barrier just for Wellisch's proton cross section
111(OPT=2)
112
113
11423 July 2008 V.Ivanchenko    hadr-pre-V09-01-07
115---------------------------------------------------
116Fixed mistake in tagging
117
11823 July 2008 V.Ivanchenko    hadr-pre-V09-01-06
119---------------------------------------------------
120Return back Coulomb barrier initialisation (J.M.Quesada)
121
12227 June 2008 V.Ivanchenko    hadr-pre-V09-01-05
123---------------------------------------------------
124Fixed G4PreCompoundNucleon and G4PreCompoundIon for
125Coulomb barrier, added protections for zero cross sections,
126OPT=2 is used (J.M.Quesada)
127
12805 June 2008 J. M. Quesada   hadr-pre-V09-01-04
129---------------------------------------------------
130Bug fixed in OPT=1 (Chatterjee) charged particle cross sections.
131Unphysical values at very low emission energies have been corrected (set
132to 0); OPT=1 ( OPT=2 in previous tag)
133
134
13515-May 2008 J. M.Quesada     hadr-pre-V09-01-03
136-----------------------------------------------
137The retrieval of transition probabilities in G4PreCompoundModel.cc
138(for protection against unphysically crossing the landa_+=landa_- condition)
139was misplaced.  Now it has been placed inside the preequilibrium transitions
140loop.
141
1428-May 2008 J. M.Quesada     hadr-pre-V09-01-02
143-----------------------------------------------
144Protection against non physical situation has been set: Equilibrium exciton number
145Neq (when reached, equilibrium regime starts) should correspond to equal transition
146probabilities "back" and "forth". Nevertheless,  for heavy target (Neq is big even
147for low incident energies) after first emission this condition (equal trans. prob.)
148is reached far before reaching corresponding Neq of the residual. Unless this "jump"
149to equilibrium is forced in this case, preequilibrium will be
150spuriously prolonged (as a side effect, with huge CPU consumption). 
151G4PreCompoundModel and G4PreCompoundTransitions classes have been modified.
152
153
1541-May 2008 J. M.Quesada     hadr-pre-V09-01-01
155-----------------------------------------------
156- First trial with cvs. 00 skipped by mistake.
157New inverse cross sections:
158        OPT=1 Chetterjee's parameterization to reaction cross sections from optical potential global fittings.
159        OPT=2 as OPT=1 but for protons the Wellisch's parameterization for protons is used  (DEFAULT)
160        OPT=3 Kalbach's modifications of Chatterjee's parameterization of cross sections
161        OPT=4 as OPT=3 but for protons the Wellisch's parameterization for protons is used
162
163Coulomb barrier has been suppressed as it enters through inverse cross sections. Also methods related to former Dostrovski's cross sections (alpha , beta and C parameters) have been suppressed. In all particle header files:
164G4PreCompoundNucleon.hh ,G4PreCompoundNeutron.hh,G4PreCompoundProton.hh, G4PreCompoundIon.hh, G4PreCompoundDeuteron.hh, G4PreCompoundTriton.hh ,G4PreCompoundHe3.hh , G4PreCompoundAlpha.hh
165
166
167Soft cutoff transition from equilibrium to equilibrium has been suppressed (unphysical and not necessary anymore).In G4PreCompoundModel.
168
169New transition probabilities :
170        OPT=1  Gudima's formulae
171        OPT=2  Machner parameterization of matrix elements for transitions (DEFAULT)
172 
173Level density is set to A/10 (Gudima's prescription).
174
175Several bugs fixed (average xs in nucleon-nucleon xs in G4PreCompundTransitions.cc and  emission probability formula in G4PreCompoundIon.cc ) . 
176
17725-October 2007 A.Howard     hadr-pre-V09-00-04
178-----------------------------------------------
179- Removed compilation warning on unused parameter in G4PreCompoundIon.cc
180
181
18211-October 2007 A.Howard     hadr-pre-V09-00-03
183-----------------------------------------------
184- Corrected for integer arithmetic bug in G4PreCompoundModel.cc
185  (ratio of charged to neutral excitons)
186
187
18801-October 2007 A.Howard     hadr-pre-V09-00-02
189-----------------------------------------------
190- Protected Rj (Q_beta) factor against divide by zeroes (should never happen?)
191  The following files have been modified:
192  G4PreCompoundAlpha.hh, G4PreCompoundDeuteron.hh, G4PreCompoundHe3.hh,
193  G4PreCompoundTriton.hh, G4PreCompoundNeutron.hh, G4PreCompoundProton.hh
194 
195
19623-August 2007 A.Howard     hadr-pre-V09-00-01
197----------------------------------------------
198- Development tag incorporating further suggestions/corrections from
199  Jose-Manuel Quesada particularly with respect to light ion production.
200  The following files have been modified:
201  G4PreCompoundIon.cc, G4PreCompoundDeuteron.hh, G4PreCompoundHe3.hh,
202  G4PreCompoundTriton.hh, G4PreCompoundAlpha.hh, G4PreCompoundIon.hh,
203  G4VPreCompoundIon.hh (unused?)
204
205
20623-July 2007 A.Howard       hadr-pre-V09-00-00
207----------------------------------------------
208- Development tag incorporating the suggestions/corrections of Jose-Manuel
209  Quesada.  These should reflect the literature and units expected for the
210  exciton model.  The following files have been modified:
211  G4PreCompoundNucleon.hh, G4VPreCompoundNucleon.hh, G4PreCompoundAlpha.hh,
212  G4PreCompoundHe3.hh, G4PreCompoundDeuteron.hh, G4PreCompoundNeutron.hh,
213  G4PreCompoundProton.hh, G4PreCompoundTriton.hh
214  G4VPreCompoundIon.cc, G4PreCompoundIon.cc, G4PreCompoundModel.cc,
215  G4VPreCompoundFragment.cc, G4PreCompoundTransitions.cc,
216  G4PreCompoundNucleon.cc, G4VPreCompoundNucleon.cc
217  The main changes are:
218  - introducing an Rj factor into the probability calculation,
219  - correcting the units for the Pauli Blocking Factor - missing a g factor,
220  - decreasing the number of charged particles if the number excitons is
221    decreased (previously this number was only increased with increasing
222    excitons).
223
224
22504-May 2007 G.Folger       hadr-pre-V08-02-00
226----------------------------------------------
227-  G4PreCompoundTransitions.cc: Protect against divide by 0
228
22919 Aug 2006 Dennis Wright (hadr-pre-V08-01-00)
230----------------------------------------------
231- exciton_model
232   G4VPreCompoundFragment.hh, .icc : remove meaningless const from
233   methods GetA, GetZ, GetRestA, GetRestZ, GetCoulombBarrier, GetBindingEnergy,
234   GetMaximalKineticEnergy, GetEnergyThreshold, GetEmissionProbability,
235   GetNuclearMass, GetRestNuclearMass and GetReducedMass
236
Note: See TracBrowser for help on using the repository browser.