Changeset 1340 for trunk/source/processes/hadronic/models/cascade/History

Timestamp:

Nov 5, 2010, 3:45:55 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

• trunk/source/processes/hadronic/models/cascade/History (modified) (2 diffs)

trunk/source/processes/hadronic/models/cascade/History

@@ -1 @@
-$Id: History,v 1.91 2010/06/11 17:07:23 mkelsey Exp $
+$Id: History,v 1.203 2010/10/22 20:41:05 mkelsey Exp $
 -------------------------------------------------------------------
 
@@ -16 +16 @@
      ---------------------------------------------------------------
 
-11 June 2010 Michael Keley (hadr-casc-V09-03-43)
-------------------------------------------------
+22 October 2010 Michael Kelsey (hadr-casc-V09-03-84)
+----------------------------------------------------
+- G4WatcherGun.cc: Change first numeric argument to G4NuclWatcher to int.
+
+NOTE:  Only this change is include in tag.
+
+20 October 2010 Michael Kelsey (kelsey-20101020b)
+-------------------------------------------------
+- G4NucleiModel:  Add scaling factors ("units") for nuclear radii and
+  cross-sections, in order to adapt calculations to the literature.
+
+20 October 2010 Michael Kelsey (kelsey-20101020a)
+-------------------------------------------------
+- G4NucleiModel:  Re-activate changes in worthToPropagate() to use existing
+  nuclear potential values for kaons and hyperons to decide on trapping.
+
+20 October 2010 Michael Kelsey (kelsey-20101020)
+------------------------------------------------
+- G4NucleiModel:  Clean up some bugs in refactoring generateModel().  Back
+  out changes (28 Sep) in worthToPropagate() in order to compare refactored
+  to production code.
+
+20 October 2010 Michael Kelsey (hadr-casc-V09-03-83)
+----------------------------------------------------
+- G4CascadParticle.cc:  Remove debugging output left in by accident.
+
+NOTE:  Tag does not include reorganization of G4NucleiModel (5 Oct 2010)
+
+19 October 2010 Michael Kelsey
+------------------------------
+G4FissionStore, G4FissionConfiguration: Undo migration to integer A and Z.
+
+19 October 2010 Michael Kelsey
+------------------------------
+Clean up CoVerity software analysis reports (three are ignored from obsolete
+code):
+        G4Analyser.cc
+        G4CascadeCheckBalance.cc
+        G4CascadeInterpolator.icc
+        G4CascadeSampler.icc
+        G4CollisionOutput.cc
+        G4ElementaryParticleCollider.cc
+        G4InuclCollider.cc
+        G4NuclWatcher.cc
+        G4NucleiModel.hh,cc
+        G4RegionModel.cc
+
+14 October 2010 Michael Kelsey
+------------------------------
+- G4NucleiModel.cc:  Fix misspelled "G4doulbe" -> G4double.
+
+Pick up several missed bits of the "integer A and Z migration," mostly in
+test and diagnostic code:
+
+        G4Analyser.hh,cc
+        G4FissionConfiguration.hh
+        G4FissionStore.hh,cc
+        G4NuclWatcher.hh,cc
+
+14 October 2010 Michael Kelsey (hadr-casc-V09-03-82)
+----------------------------------------------------
+- G4CascadParticle.cc:  Add protection in getPathToTheNextZone() so that d2
+  is always computed non-negative.  Avoids 10^-5 FPEs.
+
+05 October 2010 Michael Kelsey
+------------------------------
+- G4NucleiModel:  Major reorganization of generateModel() code.
+  1) Rename integration functions from "volNumInt[1]" to "zoneIntegralxxx"
+     with "xxx" = "WoodsSaxon" or "Gaussian".
+  2) Move hardwired constants out to static data members, give clear names.
+  3) Split function into separate pieces to set binding energies, zone
+     radii, volume/density values, and potential heights.
+  4) Within new functions, simplify if-cascades where possible, eliminate
+     duplicated code.
+  5) Move std::vector<> buffers out to data members.
+
+  Also, modify worthToPropagate() to use non-nucleon potential heights to
+  terminate propagation.  Previously, non-nucleons could never fail, as the
+  "getFermiKinetic" always returned zero.
+
+28 September 2010 Michael Kelsey (hadr-casc-V09-03-81)
+------------------------------------------------------
+- G4BigBanger:  Missed some integer A/Z migrations, which caused a bus error
+  in test47.
+
+26 September 2010 Michael Kelsey (hadr-casc-V09-03-80)
+------------------------------------------------------
+- G4VCascadeDeexcitation.hh:  NEW abstract base class for post-cascade
+  processing.  Some common activities will be implemented here from the two
+  actual modules.
+
+- G4CascadeDeexcitation, G4PreCompoundDeexcitation:  Use new base class
+  above, reorganize code in collide() to use new deExcite(G4Fragment*) where
+  appropriate.
+
+- G4InuclCascader.hh:  Reference new base class for de-excitation modules.
+
+25 September 2010 Michael Kelsey
+--------------------------------
+- G4CascadeColliderBase:  Add explosion(G4Fragment*) and explosion(A,Z,Eexc).
+  Move implementation to latter, and write others as call-throughs.  Add
+  Z==0 condition for explosion, regardless of A.
+
+- G4EquilibriumEvaporator.hh:  Remove explosion(G4InuclNuclei*), as
+  unnecessary.  Existing explosion(A,Z,E) now overrides base class.
+
+- G4CollisionOutput:  Add function to process G4ReactionProducts from
+  PreCompound.
+
+- G4PreCompoundDeexcitation.cc:  Move G4ReactionProducts loop to
+  G4CollisionOutput.
+
+24 September 2010 Michael Kelsey
+--------------------------------
+- G4IntraNucleiCascader.cc:  Minor shuffle of post-cascade recoil checking,
+  in preparation for moving fragment construction and processing away.
+
+- G4CollisionOutput.cc:  Replace names "TargetFragment*" and "NucleiFragment*" 
+  with "OutgoingNuclei" and "OutgoingNucleus" for consistency with
+  "OutgoingParticles" names.  This affects the following source files, all
+  of which are modified to match (using |sed|):
+
+        cascade/include/G4Analyser.hh
+        cascade/include/G4CollisionOutput.hh
+        cascade/src/G4Analyser.cc
+        cascade/src/G4CascadeCheckBalance.cc
+        cascade/src/G4CascadeDeexcitation.cc
+        cascade/src/G4CascadeInterface.cc
+        cascade/src/G4CollisionOutput.cc
+        cascade/src/G4EquilibriumEvaporator.cc
+        cascade/src/G4Fissioner.cc
+        cascade/src/G4IntraNucleiCascader.cc
+        cascade/src/G4InuclCollider.cc
+        cascade/src/G4InuclEvaporation.cc
+        cascade/src/G4NonEquilibriumEvaporator.cc
+        cascade/src/G4PreCompoundDeexcitation.cc
+
+- G4CollisionOutput:  Add new G4Fragment data member, with get() and add()
+  accessors.  This fragment is *not* used by setOnShell() or by any of the
+  event-total calculations.  It will be filled by G4IntraNucleiCascader for
+  use by de-excitation modules.  Include G4Fragment in diagnostic output if
+  not empty.
+
+- G4CascadeRecoilMaker.cc:  Remove G4Fragment::SetExcitationEnergy() call.
+
+- G4InuclNuclei:  Add new constructor to create from G4Fragment input, and
+  new makeG4Fragment() to spit out a G4Fragment.  Also provide a casting
+  operator for the latter.
+
+24 September 2010 Michael Kelsey (hadr-casc-V09-03-79)
+------------------------------------------------------
+Fix numerous compilation warnings on Linux (not seen on MacOSX) left over
+from migration to integer A and Z (-73 tag).
+
+        G4BigBanger.cc
+        G4CascadeCheckBalance.hh
+        G4CascadeColliderBase.cc
+        G4CascadeRecoilMaker.cc
+        G4EquilibriumEvaporator.hh,cc
+        G4ExitonConfiguration.hh
+        G4Fissioner.cc
+        G4InuclCollider.cc
+        G4NonEquilibriumEvaporator.cc
+        G4NucleiModel.cc
+
+- G4PreCompoundDeexcitation.cc:  Drop "G4int" casts as no longer required.
+
+23 September 2010 Michael Kelsey (hadr-casc-V09-03-78)
+------------------------------------------------------
+- G4CascadeInterface, G4InuclCollider:  change user selection function name
+  to usePreCompoundDeexcitation(), per Julia Yarba.
+
+- G4PreCompoundCascadeInterface.hh,.cc, G4PreCompountInuclCollider.hh,.cc:
+  REMOVE THESE CLASSES.  These were copied and modified versions of the
+  standard versions, in order to test the new pre-compound de-excitation.
+  That has been integrated with runtime flags, so the parallel versions are
+  not required.
+
+23 September 2010 Michael Kelsey
+--------------------------------
+- G4InuclElementaryParticle.cc:  Drop warning message when converting
+  G4ParticleDefinition to type code.
+
+- G4CascadeInterface:  Add functions for user selection of which
+  post-cascade module to use for de-excitation of the residual nucleus.
+  These call through to the equivalent G4InuclCollider functions.
+
+- G4InuclCollider:  Add functions for user selection of which
+  post-cascade module to use for de-excitation of the residual nucleus.
+  These delete and re-instantiate the appropriate de-excitation "collider".
+
+- G4PreCompoundDeexcitation:  Remove convertFragment() function, remove
+  pass-by-value std::vector<> from getDeExcitedFragments() (use data member
+  instead).  Clean up some obsolete code usage.
+
+- G4CascadeRecoilMaker:  Change fragment-making interface:
+  makeRecoilNuclei() now returns G4InuclNuclei* for standard Bertini colliders.
+  makeRecoilFragment() returns G4Fragment* for external PreCompound models.
+  addExcitonConfiguration() sets local copy of excitons, needed in order to
+    call the individual G4Fragment->Setxxx() functions.
+
+- G4IntraNucleiCascader.cc:  Follow renaming RecoilMaker::makeRecoilNuclei(),
+  use new addExcitonConfiguration().
+
+- G4CascadeRecoilMaker:  Change fragment-making interface:
+  makeRecoilNuclei() now returns G4InuclNuclei*.
+  makeRecoilFragment() returns G4Fragment*.
+  addExcitonConfiguration() sets local copy of excitons, needed in order to
+    call the individual G4Fragment->Setxxx() functions.
+
+- G4IntraNucleiCascader.cc:  Follow renaming RecoilMaker::makeRecoilNuclei(),
+  use new addExcitonConfiguration().
+
+22 September 2010 Michael Kelsey
+--------------------------------
+- G4CascadeMomentum.hh:  This should have been removed from the package long
+  ago.  Apologies for the oversight.
+
+22 September 2010 Julia Yarba 
+-----------------------------
+Added 1st prototype of interface between Bartini cascade and PreCompound model 
+(post-cascade de-excitation).
+
+Changes only to:
+/geant4/source/processes/hadronic/models/cascade/cascade
+
+New developments are those G4PreCompound* classes; the GNUmakefile adapted to 
+include necessary links.
+
+The PreComp wrapper is G4PreCompundDeexcitation class.
+
+The logic, loops, conditions, etc., in the event processing are the same as in 
+the "native" Bertini code. I actually copied it over from G4CascadeInterface,
+and started from there.
+
+However, there're a few technical changes that I'd like to mention.
+
+1. G4PreCompoundInterface inherits directly from G4HadronicInteraction.
+I felt that we'd not need G4VIntra..., because its only advantage is a pointer 
+(and related access methods) to G4VPreCompoundModel.
+This has now moved to the de-excitation wrapper.
+
+2. In the Collider itself, I've made the data member to be
+ G4CascadeColliderBase* theDeexcitation
+This way, it can be transparently set to the native one or the PreComp one 
+(both inherit from G4CascadeColliderBase).
+
+3. There's a setDeExcitation(...) methond in the collider, and also in the interface 
+(which will pass it down to the collider).
+
+19 September 2010 Michael Kelsey (hadr-casc-V09-03-77)
+------------------------------------------------------
+- G4CascadeInterface.cc:  Bug fix to logic in retryInelasticNucleus()
+
+17 September 2010 Michael Kelsey (hadr-casc-V09-03-76)
+------------------------------------------------------
+- G4CascadeInterface:  Add support to pass nucleus projectiles into
+  collider.  Not tested or guaranteed to work, but G4IntraNucleiCascader is
+  supposed to deal with them.
+
+16 September 2010 Michael Kelsey
+--------------------------------
+- G4CascadeInterface:  Extensive reorganization, moving most
+  functionality out of ApplyYourself and into small member functions.
+  Complicated end-of-loop conditions encapsulated.  Data members made into
+  pointers to hide dependences from outside world.
+
+16 September 2010 Michael Kelsey (hadr-casc-V09-03-75)
+------------------------------------------------------
+- G4CascadeInterface.cc:  Add parentheses in end-of-while block to fix
+  compilation warning in GCC 4.3.
+
+- G4IntraNucleiCascader.cc:  In decayTrappedParticle(), check for photon
+  daughters and put them directly onto output list.
+
+16 September 2010 Michael Kelsey
+--------------------------------
+- G4CascadParticle.hh:  Make getGeneration() const.
+
+- G4CascadeDeexcitation.hh:  Fix compiler warning about order of data members.
+
+- G4IntraNucleiCascader:  Intercept hyperons which are "trapped" in nuclear
+  potential, and decay them immediately to produce usable secondaries.  This
+  resolves problem with sub-MeV final-state hyperons (reported by A. Dotti),
+  but is not entire correct.  Properly, the hyperon should be incorporated
+  into a hypernucleus fragment.
+
+- G4InuclElementaryParticle.hh:  Add hyperon() identification function, and
+  constructor to take G4DynamicParticle directly.
+
+- G4InuclParticle.hh, G4InuclNuclei.hh:  Add constructor to take
+  G4DynamicParticle directly.
+
+15 September 2010 Michael Kelsey (hadr-casc-V09-03-74)
+------------------------------------------------------
+- G4CascadeDeexcitation:  NEW collider to handle post-cascade processing of
+  nuclear fragment.  Encapsulates BigBanger and Evaporators.  Eventually
+  will have alternative (configurable or #ifdef) implementation using G4
+  external "pre-compound" models instead of Bertini factories.
+
+- G4InuclCollider:  Remove post-cascade colliders (BigBanger, Evaporators),
+  replace with new G4CascadeDeexcitation.
+
+- G4LorentzConverter:  Move construtor implementations to .cc file, and make
+  sure all constructors have initializers for all data members.
+
+14 September 2010 Michael Kelsey (hadr-casc-V09-03-73)
+------------------------------------------------------
+Migrate to integer A and Z values for nuclear configurations.  Includes both
+explicit arguments for G4InuclNuclei, and function arguments, computational
+parameters involved in nuclear configurations, breakup, etc.
+
+        G4CascadeCheckBalance.cc
+        G4CascadeInterface.cc
+        G4CascadeRecoilMaker.hh
+        G4Fissioner.hh,cc
+        G4InuclEvaporation.cc
+        G4InuclNuclei.hh,cc
+        G4InuclSpecialFunctions.hh
+        G4NonEquilibriumEvaporator.hh,cc
+        G4NucleiModel.hh,cc
+        G4PreCompoundCascadeInterface.cc
+        bindingEnergy.cc
+        bindingEnergyAsymptotic.cc
+        nuclearLevelDensity.cc
+        paraMaker.cc
+
+- G4InuclSpecialFunctions.hh:  Eliminate bindingEnergy functions which are
+  no longer used, along with their .cc files:
+
+        bindingEnergyKummel
+        bindingEnergyExact
+
+Migration validated using 100 events (20 GeV pi+ on lead) with verbose==4.
+Entire log files (9.6M lines, 433 MB) identical before and after migration.
+
+13 September 2010 Michael Kelsey
+--------------------------------
+- G4InuclElementaryParticle:  Move printParticle() implementation to .cc,
+  and add printing of particle name to output.
+
+- G4InuclNuclei.cc:  Move base-class printing to first line of output, to
+  match G4InuclEP.
+
+10 September 2010 Michael Kelsey (hadr-casc-V09-03-72)
+------------------------------------------------------
+- G4CascadeRecoilMaker:  Add new goodNucleus() function with functionality
+  from G4IntraNucleiCascader::goodCase().  Add parameter to set "rounding
+  tolerance" for small/negative mass differences.  Drop getRecoilFragment()
+  in favor of user calling makeRecoilFragment() directly (instead of doing
+  it collide()), which returns a non-const pointer. 
+
+- G4IntraNucleiCascader: Remove previously introduced getResidualMass() and
+  makeResidualFragment() functions, along with G4InuclNuclei object.
+  Replace with new RecoilMaker utility.  Move goodCase() to RecoilMaker, and
+  begin to simplify end-of-cascade processing.
+
+- G4ExitonConfiguration.hh:  Add function to reset values to zero.
+
+- G4InuclNuclei:  Add function to zero out exciton configuration, and use
+  it in fill() functions.
+
+9 September 2010 Michael Kelsey
+-------------------------------
+- G4CascadeCheckBalance:  Add new collide() interface to take lists of both
+  G4InuclElementaryParticles and G4CascadParticles; supports use by
+  G4IntraNucleiCascader and G4CascadeRecoilMaker.
+
+- G4CascadeRecoilMaker:  NEW pseudo-Collider class to compute the nuclear
+  recoil kinematics of a cascade, either in process or completed.  To be
+  used by G4IntraNucleiCascader; uses G4CascadeCheckBalance to do work.
+
+6 September 2010 Michael Kelsey (hadr-casc-V09-03-71)
+-----------------------------------------------------
+- G4NucleiModel:  Add protections in generateInteractionPartners() so that
+  nucleon and quasideuteron interactions are only selected if there are
+  "enough" nucleons of the necessary type still available (xxxNumberCurrent).
+
+- G4InuclNuclei:  Add fill() functions which overwrite the entire particle
+  structure, with same argument lists as constructors.  Will be used mainly
+  by G4IntraNucleiCascader::makeResidualFragment().
+
+- G4IntraNucleiCascader:  Hide non-physical nucleus message behind verbose
+  flag.  Add function and data member to build recoil nucleus (fragment) at
+  end of every iteration; will eventually replace functionality of
+  getResidualMass().
+
+2 September 2010 Michael Kelsey (hadr-casc-V09-03-70)
+-----------------------------------------------------
+- G4NucleiModel.cc:  Restore hadr-casc-V09-03-68 revision 1.71, and remove
+  the two resize(3) actions.
+
+NOTE:  The exercise below has (finally) resolve the cross-section
+discrepancy reported by Sunanda Banerjee back in July.  It turns out that my
+use of "resize(3)" to pre-allocate the qdeutron and acsec buffers on each
+interaction was incorrect.  This doesn't just reserve memory, it fills the
+vectors with three elements, so that subsequent push_back()'s fill [3], [4],
+etc.
+
+2 September 2010 Michael Kelsey
+-------------------------------
+- G4NucleiModel.cc: Reverted to hadr-casc-V09-03-38 revision 1.45, and
+  interface changes incorporated.  See NOTE below from 28 July 2010.
+  Cross-section changes have not been resolved or validated.  Since -38
+  revision has better data/MC matching (ratio near 1), it is restored as the
+  baseline version, and changes will be re-done incrementally.
+
+  Each increment after (0) is tagged "kelsey-20100902x", x=a,b,c,d etc.
+
+        0)  Implement new ctor and generateModel() calls
+            Make cross-section tables statically initialized
+            Include absorptionCrossSection definition
+        1)  Use lookup tables for cross-sections, and new interpolator
+            Change CHC_CHECK to G4CASCADE_DEBUG_CHARGE
+        2)  Use this package's bindingEnergy() function
+        3)  Update diagnostic messages and verbosity levels;
+            Add missing initializers to all constructors;
+            Use data-member buffer for G4CollisionOutput
+        4)  Use generateWithRandomAngles() everywhere appropriate;
+            Collapse if-cascades to use if-return where possible
+        5)  Make generateModel() re-usable with data-member buffers
+        6)  Energy-momentum conservation checking;
+            Implement generateNucleonMomentum;
+            Remove extraneous semicolons and blank lines
+        7)  Use local variables in passFermi();
+            Fix some minor LorentzVector calculations;
+            Use generateNucleonMomentum() for quasi-deuterons
+        8)  Use data-member buffers for quasideuterons, coords, mom.
+        9)  Use enum labels for quasideuteron codes
+        10) Restore Dennis' improved potential (6-zone) calculations
+
+4 August 2010 Michael Kelsey (hadr-casc-V09-03-69)
+--------------------------------------------------
+- G4CascadeData:  Add name data member and optional ctor argument, use
+  when printing tables.
+
+- G4CascadeFunctions.icc:  Pretty-up printing output using name string.
+
+- G4Cascade*Channel.cc:  Add name string to data() ctor, use central part of
+  class name, e.g., "PiMinusP" for G4CascadePiMinusPChannel.
+
+- G4ElementaryParticleCollider.cc:  Add printing of all final-state tables,
+  one time only, protected by G4CASCADE_DEBUG_SAMPLER.
+
+- GNUmakefile:  Add G4CASCADE_DEBUG_SAMPLER preprocessor flag.
+
+3 August 2010 Michael Kelsey (hadr-casc-V09-03-68)
+--------------------------------------------------
+- G4CascadeData.hh: Add print() and printXsec() functions to dump tables.
+- G4CascadeData.icc:  Move implementations out of G4CascadeData.hh file.
+
+- G4CascadeFunctions.hh,icc: Add printTable() to dump tables.
+
+- G4CascadeInterpolator.hh,icc: Add printBins() to dump energy bin edges.
+
+- G4CascadeSampler.hh,icc: Add print() to call through to interpolator.
+
+NOTE:  This should be the last set of functional changes for a while.
+Sunanda Bannerjee reported substantial changes in validation plots for the
+Bertini cascade between 4.9.3-ref-05 (hadr-casc-V09-03-23-01) and
+4.9.3-ref-06/4.9.4-beta-01 (hadr-casc-V09-03-43).
+
+28 July 2010 Michael Kelsey (hadr-casc-V09-03-67)
+-------------------------------------------------
+- G4NucleiModel.cc:  Set deuteron arrays to fixed length (3).  Add A/Z
+  values to generateModel diagnostic message.
+
+- G4EquilibriumEvaporator.cc:  Change some fixed length std::vector<> to
+  simple C arrays.
+
+- G4Fissioner.cc:  Move G4FissionStore to class data member for reuse.
+
+- G4FissionStore:  Add ::clear() function to support reuse, move addConfig()
+  implementation to .cc file.
+
+- G4IntraNucleiCascader:  Move G4NucleiModel to class-level data member
+  (pointer) so it can be reused across collisions.
+
+26 July 2010 Michael Kelsey
+---------------------------
+- G4Analyzer.cc:  Use const-refs for particle lists.
+
+- G4BigBanger, G4ElementaryParticleCollider:  Move std::vector<> buffers to
+  .hh file to allow reuse.
+
+- G4NucleiModel.cc:  Simplify buffers used in generateModel() to reduce
+  G4double allocations.
+
+23 July 2010 Michael Kelsey (hadr-casc-V09-03-66)
+-------------------------------------------------
+- G4CascadeInterface:  Move G4CollisionOutput to .hh file for reuse between
+  events.
+
+- G4IntraNucleiCascader:  Move buffers for cascade output to .hh file, so
+  they can be reused between events.
+
+- G4NucleiModel:  Move G4CollisionOutput to .hh file for reuse; use
+  references for bullet and target args to G4EPCollider.
+
+21 July 2010 Michael Kelsey (hadr-casc-V09-03-65)
+-------------------------------------------------
+- GNUmakefile:  Add new G4CASCADE_CHECK_ECONS to turn off all of the
+  internal conservation checks (reduce std::vector<> memory churn!).
+
+- G4CascadeColliderBase.cc:  Use G4CASCADE_CHECK_ECONS to set the default
+  value of "doConservationChecks": true if defined, false if not defined.
+  This is done with #ifdef-else-endif block, NOT some fancy macro.
+
+- G4CascadeInterface.cc:  Use G4CASCADE_SKIP_ECONS to toggle reporting
+  violations before end of loop (don't need to call balance.okay() twice!).
+
+- G4ElementaryParticleCollider.cc:  Put verbosity check before call to
+  validateOutput() to avoid unnecessary vector manipulations.
+
+- G4IntraNucleiCascader.cc:  Set doConservationChecks = true
+  unconditionally; used to reject bad cascades.
+
+- G4NucleiModel.cc:  Use G4CASCADE_CHECK_ECONS to hide balance.collide().
+
+20 July 2010 Michael Kelsey (hadr-casc-V09-03-64)
+-------------------------------------------------
+- G4CascadeColliderBase: Make G4CascadeCheckBalance a pointer member, move
+  ctor and dtor to .cc file.  Change balance tolerances to 1 per mil and 1
+  MeV.
+
+- G4InuclCollder:  Make all colliders pointer members.
+
+- G4IntraNucleiCascader:  Make EPCollider a pointer member, and call
+  G4CollisionOutput::setOnShell() to balance energy-momentum for null fragments.
+
+- G4InuclNuclei.cc:  Preserve momentum magnitude when setting excitation
+  energy by setting mass and kinetic energy.
+
+19 July 2010 Michael Kelsey (hadr-casc-V09-03-63)
+-------------------------------------------------
+- G4CascadeCheckBalance.cc:  Change zero-tolerance to 10 keV (1e-5) from
+  1 keV (1e-6).
+
+- G4CollisionOutput.hh:  Fix capitalization of "numberOfOutgoingParticles".
+
+- G4EquilibriumEvaporator.cc:  Remove duplicate EEXS setting.
+
+- G4InuclNuclei:  Change implementation of setExitationEnergy() to adjust
+  mass while keeping momentum vector constant.  Currently done with
+  "expensive" LorentzVector manipuation.  Should be done by recomputing
+  kinetic energy given new mass value.
+
+- G4NonEquilibriumEvaporator.cc:  Modify particle-creation block to work as
+  G4EquilibriumEvaporator now does.  Simplify some if-blocks.
+
+16 July 2010 Michael Kelsey (hadr-casc-V09-03-62)
+-------------------------------------------------
+- G4CollisionOutput:  Add interfaces to handle G4CascadParticles.
+
+- G4IntraNucleiCascader:  Eliminate local "inter_case" data member and
+  function.  Use base class G4InteractionCase directly.  Add functionality
+  to test for conservation and recoil before exiting "itry" loop, and redo
+  cascade if unbalanced recoil fragment.
+
+- G4{PreCompound}InuclCollider.cc:  Drop "setInteractionCase()" function.
+
+15 July 2010 Michael Kelsey (hadr-casc-V09-03-61)
+-------------------------------------------------
+- G4CollisionOutput:  New ::add() function to combine two objects.  New
+  accessors to get total charge and baryon number.
+
+- G4CascadeCheckBalance:  Make tempOutput data member, for thread-safety,
+  rather than static variables.  Use new G4CollisionOutput accessors.
+
+- G4IntraNucleiCascader.cc:  Move local output buffers outside while loop,
+  use std::vector<>::insert() for copying.
+
+- G4InuclCollider.cc:  Use new G4CollisionOutput::add() interface.
+
+- G4NucleiModel.cc:  Hide conservation checks behind verbose>2.  Otherwise
+  log files get too big
+
+15 July 2010 Michael Kelsey (hadr-casc-V09-03-60)
+-------------------------------------------------
+- G4CascadeCheckBalance:  Add layer of protection in relativeX() (check
+  initial values) to avoid divide-by-zero errors.  Add another interface to
+  accept list of G4CascadParticles, for use with G4NucleiModel.  Add check
+  of momentum direction balance, not just magnitudes.
+
+- G4IntraNucleiCascader.cc:  Rearrange end-of-loop if-blocks to allow for
+  conservation checking prior to return.  Add reporting of both too-low and
+  too-high residual energy (at verbose=2) for single-nucleon fragments.
+
+- G4InuclCollider.cc:  Remove re-named conservation checks, now all done by
+  colliders themselves.
+
+- G4InuclParticle.hh:  Add setKineticEnergy() function to simplify handling
+  of some kinematics.
+
+- G4NucleiModel:  Add G4InuclNuclei data member (created in generateModel()),
+  for use with conservation checks.  Do check at end of generateParticleFate().
+
+15 July 2010 Michael Kelsey (hadr-casc-V09-03-59)
+-------------------------------------------------
+- G4InuclNuclei:  Use G4DynamicParticle::theDynamicalMass to carry around
+  excitation energy (== theDynamicalMass-thePDGMass).  This allows
+  excitation to be changed on the fly, with all kinematics handled
+  correctly.  It also allows ONLY the ground state nucleus to be created as
+  a G4Ions object (whether global or private), reducing the number of such
+  objects by a factor of thousands.
+
+14 July 2010 Michael Kelsey (hadr-casc-V09-03-58)
+-------------------------------------------------
+- G4CollisionOutput.cc:  For initial four-momentum adjustment, skip over
+  particles which would acquire negative energy after "correction."
+
+- G4ElementaryParticleCollider.cc: Bug fix for two-body final states.
+  Kinematics should not require "rescaling."  That was driven by states
+  which were kinematically forbidden, with m1+m2 > etot_scm (e.g., pi0 p ->
+  pi+ n with just 5 MeV kinetic energy).  Now those states are caught and
+  rejected.
+
+- G4CascadeInterface.cc, G4InuclCollider.cc:  Improve diagnostics reporting
+  number of retries.
+
+- G4EquilibriumEvaporator.cc:  When computing residual excitation
+  (EEXS_new), original thrown energy, S, should be used, not reboosted
+  particle after recoil.  Simplify evaporation kinematics to conserve
+  four-momentum automatically (PEX -= mom).
+
+14 July 2010 Michael Kelsey (hadr-casc-V09-03-57)
+-------------------------------------------------
+- G4CascadeColliderBase: New subclass of G4VCascadeCollider, which takes
+  over all of its concrete functionality, and adds conservation checker
+  (G4CascadeCheckBalance) as data member.  Wrapper functions and control
+  flag allow colliders to user it transparently.
+
+- G4VCascadeCollider:  Move concrete functionality to new base class.
+
+- All colliders:  Change base class to new G4CascadeColliderBase.  Remove
+  local copies of conservation-checker, using new base function instead.
+
+        cascade/include/G4BigBanger.hh
+        cascade/include/G4ElementaryParticleCollider.hh
+        cascade/include/G4EquilibriumEvaporator.hh
+        cascade/include/G4EvaporationInuclCollider.hh
+        cascade/include/G4Fissioner.hh
+        cascade/include/G4IntraNucleiCascader.hh
+        cascade/include/G4InuclCollider.hh
+        cascade/include/G4NonEquilibriumEvaporator.hh
+        cascade/include/G4PreCompoundInuclCollider.hh
+        cascade/src/G4BigBanger.cc
+        cascade/src/G4ElementaryParticleCollider.cc
+        cascade/src/G4EquilibriumEvaporator.cc
+        cascade/src/G4EvaporationInuclCollider.cc
+        cascade/src/G4Fissioner.cc
+        cascade/src/G4IntraNucleiCascader.cc
+        cascade/src/G4InuclCollider.cc
+        cascade/src/G4NonEquilibriumEvaporator.cc
+        cascade/src/G4PreCompoundInuclCollider.cc
+
+- G4EquilibriumEvaporator:  Use new functionality (above) to turn off
+  conservation checks when evaporating fission products.
+
+- G4CascadeInterface.cc:  Report number of retries in verbose messages.
+
+- G4IntraNucleiCascader.cc:  Remove sanity check on afin/zfin vs. model.
+
+13 July 2010 Michael Kelsey (hadr-casc-V09-03-56)
+-------------------------------------------------
+- G4CascadeCheckBalance:  Add zero protection (needed for momentum ratio)
+  using a static tolerance parameter.
+
+- G4EquilibriumEvaporator.cc:  Add dumping of G4CollisionOutput before exit,
+  some fixes for "PEX" adjustments as EEXS changes.
+
+- G4NonEquilibriumEvaporator.cc:  Add conservation checking, dump of output,
+  some fixes for "PEX" adjustments as EEXS changes.
+
+- G4Fissioner.cc:  Add conservation checking, remove addition of excitation
+  energy to input nucleus mass, since already included.
+
+13 July 2010 Michael Kelsey
+---------------------------
+- G4CascadeCheckBalance.cc:  Hide violation reports behind verbose == 1.
+
+- G4ElementaryParticleCollider.cc:  Bump ::collide() message to verbose > 1.
+
+- G4EquilibriumEvaporator.cc:  Add conservation checking, collapse if-cascades
+  by doing failure tests and exits up front, add sanity check on eex_real.
+
+- G4IntraNucleiCascader.cc:  Add diagnostic messages and some FIXME
+  comments, validate afin/zfin against G4NucleiModel result.
+
+- G4InuclCollider.cc:  Move conservation check outside verbosity flag.  This
+  adds to CPU usage, so may need #ifdef protection later.  Adjust a few
+  verbosity levels.
+
+- G4LorentzConvertor.cc:  Bump all diagnostic messages to verbose > 2.
+
+- G4NucleiModel.cc:  Simplify main if-block to put failure up front; add
+  debugging messages when decrementing proton or neutron counts.
+
+11 July 2010 Michael Kelsey
+---------------------------
+- G4CascadeCheckBalance:  Add ctor argument with name of parent class; use
+  this name when reporting results.  Add interface for nuclear fragment
+  list.  Allow parent class name to be changed dynamically.
+
+- G4BigBanger.cc, G4CascadeInterface.cc, G4ElementaryParticleCollider.cc,
+  G4InuclCollider.cc, G4NucleiModel.cc:  Pass name to G4CascadeCheckBalance.
+
+- G4Fissioner.cc:  Use G4CascadeCheckBalance to test energy conservation.
+
+- G4InuclNuclei.hh:  Add model ID as optional ctor argument.
+
+- G4InuclCollider.cc:  Move balance checks on evaporators outside verbosity.
+  Adjust name of G4CascadeCheckBalance for clear reporting.
+
+- G4VCascadeCollider.hh:  Add protected function to change name dynamically.
+
+6 July 2010 Dennis Wright (hadr-casc-V09-03-55)
+-----------------------------------------------
+- G4IntraNucleiCascader.cc: add to final particle list those cascade 
+  particles which were formerly abandoned at cascade termination;
+- G4IntraNucleiCascader.cc: allow meson excitons to pass out of nucleus
+  since they are not proper excitons, and would be lost otherwise 
+  
+4 July 2010 Dennis Wright (hadr-casc-V09-03-54)
+-----------------------------------------------
+- tag changes up to now 
+
+2 July 2010 Michael Kelsey
+--------------------------
+- G4BigBanger.cc:  Moving "generateMomentumModules()" inside loop.  For
+  three-body decays, the angles are not ajustable (so always fails).
+
+- G4CascadeInterface.cc:  Set verbosity on G4CollisionOutput.
+
+- G4EquilibriumEvaporator.cc:  Simplify if-cascades to reduce nesting (test
+  for errors and exit up front, rather than in else-blocks).
+
+- G4InuclCollider.cc:  Set verbosity on G4CollisionOutput.
+
+- G4IntraNucleiCascader.cc:  Throw away casacde which leaves "negative energy
+  recoil" (i.e. rest-frame energy of recoil less than mass of nucleus or
+  nucleon).
+
+1 July 2010 Michael Kelsey
+--------------------------
+- G4InuclCollider.cc:  Used wrong input target for conservation-law checking
+  subsequent to G4EquilibriumEvaporator.
+
+- G4EquilibriumEvaporator.cc, G4NonEquilibriumEvaporator.cc:  Handle
+  excitation energy differently with new G4InuclNuclei interface.  Compute
+  updated values correctly, by taking nuclear mass difference then
+  subtracting energy of evaporated particle.
+
+- G4InuclNuclei:  Reorganization to properly incorporate excited states.
+  Excitation energy included in mass in all cases.  setExcitationEnergy() now
+  replaces G4PartDefn pointer with new excited state.  getNucleiMass() takes
+  (optional) excitation, calls G4NucleiProperties to get mass value.  For
+  weird fragments, multiple isomers allowed, assigned unique non-standard
+  PDG codes.  Include particle name in printout (now in .cc file).
+
+- G4CollisionOutput.cc:  Do not include nuclear excitation energy in output
+  sum; now included with nuclear mass automatically (and properly).
+
+- G4CascadeCheckBalance.cc:  Do not include nucler excitation energy in input
+  sum; now inluded with nuclear mass automatically (and properly).
+
+- G4IntraNucleiCascader.cc:  For recoil nuclear fragments, treat excitation
+  energy properly: compute as difference between recoil "mass" (recoil-frame
+  energy) and ground-state A/Z mass, and pass into G4InuclNuclei ctor.
+  After the cascade, add some recoil energy checks to deal with both
+  floating-point round-off, and quasi-elastic scatters (which change the
+  bullet direction, and leave the nucleus with "missing" energy).
+
+- G4BigBanger.cc:  Simplify boost to "target" nucleus, using "boostVector"
+  instead of G4LorentzConvertor machinery.
+
+30 June 2010 Michael Kelsey
+---------------------------
+- G4BigBanger.cc:  Improve handling of excitation energy as part of nuclear
+  mass, deal with two-body explosion properly (x=0.5, not random).  Some
+  tweaks to diagnostic messages.
+
+- G4CascadeCheckBalance:  Add excitation energy to total energy for bullet
+  and target nuclei.  This is done in G4CollisionOutput for secondaries.  It
+  is NOT relativistically correct (should add to mass, not energy), but is
+  consistent.
+
+- G4CollisionOutput.cc:  Use "getExitationEnergyInGeV()" instead of "0.001*".
+
+- G4IntraNucleiCascader.cc:  More diagnostic messages, improve the afin/zfin
+  sanity check by requiring A>=Z.
+
+28 June 2010 Michael Kelsey
+---------------------------
+- G4CascadeCheckBalance:  Report violations on G4cerr always, regardless
+  of verbosity level.  Add collide(...,<vector>) interface to allow use 
+  from within G4ElementaryParticleCollider.
+
+- G4ElementaryParticleCollider.cc:  Add conservation checks within ::collide().
+
+- G4IntraNucleiCascader.cc:  Add more diagnostic messages.  Test for
+  negative values of afin, zfin, and try again.  Set round-off limit (1e-6),
+  and test nuclear recoil energy for bound (setting to zero) before doing
+  energy-violation test.
+
+- G4InuclNuclei.cc:  Abort job if passed negative Z or A argument.
+
+28 June 2010 Michael Kelsey
+---------------------------
+- G4BigBanger.cc:  Missed one instance of GetBindingEnergy().
+
+26 June 2010 Dennis Wright
+--------------------------
+- G4InuclNuclei:  Fix bug in calculation of PDG code.
+
+25 June 2010 Michael Kelsey (hadr-casc-V09-03-53)
+-------------------------------------------------
+- G4ElementaryParticleCollider.cc:  MAJOR BUG FIX:  Three different places
+  in generating N-body momenta had incorrect "recoil" calculations.  The
+  Nth four-vector was computed as simply -total, which produced negative
+  energy states (obviously).  The correct calculation should have been
+  N = G4LV(0,0,0,etot) - (total) in each instance.
+
+24 June 2010 Michael Kelsey (hadr-casc-V09-03-52)
+-------------------------------------------------
+- G4CascadeCheckBalance.cc:  Bug fixes to limit checks; should be using
+  std::abs(), or negative violations don't get caught!  Also did some
+  cleanup for compiler warnings see on Linux.
+
+23 June 2010 Michael Kelsey
+---------------------------
+- G4CascadeInterface.cc:  Fix infinite-loop bug with energy-violation; must
+  check for (nTries>maxTries) outside parenthetical, with balance.okay()
+  inside.  Do conservation checks all the time, not just with flag set.  If
+  e-violation is "infinite loop", dump last generated event and abort job.
+
+- G4CascadeCheckBalance:  Add checks on baryon number, charge, and kinetic
+  energy.  "Global" sanity check (::okay() function) does not include K.E.
+
+23 June 2010 Michael Kelsey (hadr-casc-V09-03-51)
+-------------------------------------------------
+- bindingEnergy.cc:  Replace entire function with call-through to
+  G4NucleiProperties.  Copy test for valid A/Z from there, and return
+  0. directly on failure; suppresses warning message in 9.4-beta.
+
+- G4EquilibriumEvaporator.cc, G4Fissioner.cc, G4IntraNucleiCascader.cc
+  G4InuclNuclei.cc, G4NonEquilibriumEvaporator.cc, G4NucleiModel.cc
+  G4VCascadeCollider.cc:  Restore calls to local "bindingEnergy()" to
+  get wrapper function above; suppresses warning messages in 9.4-beta.
+
+- G4NucleiModel.cc:  Make class re-usable (had been assumed in -47 tag, but
+  was not true!), by clearing all parameter vectors in generateModel().
+  Values of A,Z are tested first, and previous data kept if same nucleus is
+  being modelled.
+
+23 June 2010 Michael Kelsey
+---------------------------
+- G4CascadeInterface.cc:  Minor fix to rationalize version jumpings.
+- G4NucleiModel.cc:  Fix bug in partners list termination mods below, and
+  remove reflection-recoil from boundaryTransition().
+
+
+21 June 2010 Michael Kelsey (MHK-20100621)
+------------------------------------------
+- G4CascadeInterface.cc:  Bring MHK-20100620 version back to HEAD, make same
+  parenthesis fix for GCC 4.5 compilation.
+
+- G4IntraNucleiCascader.cc:  Bug fix nuclear recoil kinematics introduced in
+  MHK-20100620.
+
+- G4NucleiModel.cc:  Bring MHK-20100620 version back to HEAD, apply ctor
+  argument fix.  Eliminate creating temporaries for terminal "partners" on
+  list; default ctor will take care of it.  Hide negative path-length report
+  behind verbose flag, and don't return an empty list.
+
+  Add block of code in ::boundaryTransition() which recoils the nucleus when
+  the secondary is reflected off the potential, to balance the momentum
+  transfer.  This is implemented by computing the boost for the full nucleus
+  using the computed three-momentum transfer, and applying that boost to the
+  secondary.
+
+23 June 2010 Michael Kelsey (hadr-casc-V09-03-50)
+-------------------------------------------------
+- G4IntraNucleiCascader.cc:  Restore G4NucleiModel as local data member.
+  Not well designed for reuse.
+
+21 June 2010 Michael Kelsey (hadr-casc-V09-03-49)
+-------------------------------------------------
+NOTE:  Changes below were version jumped into V09-03-48.  Tag above ONLY
+       includes these two changes!
+
+- G4CascadeInterface.cc:  Add parentheses in main iteration loop to support
+  mixed && and || (reported by gcc-4.5).
+
+- G4NucleiModel.cc:  Fix mistake in passing ctor args to generateModel().
+
+20 June 2010 Michael Kelsey (MHK-20100620)
+------------------------------------------
+- G4CascadeCheckBalance:  New utility class to do energy and momentum
+  conservation testing.  Implemented as a "collider" to resue interface.
+
+- G4CascadeInterface.cc:  Use new G4CascadeCheckBalance for energy
+  conservation.
+
+- G4CollisionOutput:  Add setVerboseLevel() function, diagnostic reports in
+  ::setOnShell(), and collapse if-cascades to use "if (!xxx) return" instead
+  of multiple nestings.
+
+- G4ElementaryParticleCollider.cc: Collapse if-cascade to use "if (!xxx)
+  return" instead of multiple nestings.
+
+- G4IntraNucleiCascader.cc:  Improve diagnostic messages, change some G4cout
+  to G4cerr.
+
+- G4InuclCollider:  Make output buffers data members to reduce churn;
+  simplify if-cascades to use "if (!xxx) return" instead of nesting; add
+  energy-conservation checks at each stage of cascade.
+
+- G4NucleiModel.cc:  Improve diagnostics, collapse if-cascades to use "if
+  (!xxx) return" instead of multiple nestings.
+
+18 June 2010 Michael Kelsey (hadr-casc-V09-03-48)
+-------------------------------------------------
+- G4CascadeInterface.cc:  Add check on energy conservation between initial
+  and final states; more than 5% imbalance triggers repeat generation of
+  cascade.  This check may be disabled with flag G4CASCADE_SKIP_ECONS.  Also
+  did some minor cleanup of the code.
+
+- cascade/GNUmakefile:  Add support for new flag above.
+
+17 June 2010 Michael Kelsey (hadr-casc-V09-03-47)
+-------------------------------------------------
+- G4CascadeInterface, G4PreCompoundCascadeInterface:  Make G4InuclCollider a
+  local data member, pass verbosity to it in ::ApplyYourself().
+
+- G4InuclCollider, G4ElementaryParticleCollider, G4IntraNucleiCascader,
+  G4EquilibriumEvaporator, G4BigBanger: Make all sub-colliders local data
+  members, pass verbosity to each of them.
+
+- G4NucleiModel:  Add ctors and generateModel() signatures so both can do
+  particle or A/Z initialization.  Add setVerboseLevel() function.
+
+17 June 2010 Michael Kelsey (hadr-casc-V09-03-46)
+-------------------------------------------------
+- cascade/GNUmakefile: Add series of "ifneq (,...)" blocks to map
+  environment variables for debugging onto "-D..." compiler flags:
+
+        G4CASCADE_COULOMB_DEV
+        G4CASCADE_DEBUG_CHARGE
+        G4CASCADE_DEBUG_INTERFACE
+        G4CASCADE_DEBUG_SORT
+        G4CASCADE_WATCHER_HIGHZ
+
+- G4CascadeInterface.cc: Replace preprocessor flags "debug_G4...Interface"
+  with G4CASCADE_DEBUG_INTERFACE, and "BERTDEV" with G4CASCADE_COULOMB_DEV.
+  Fix bug with information access in debugging block.
+
+- G4EquilibriumEvaporator.cc, G4IntraNucleiCascader.cc,
+  G4NonEquilibriumEvaporator.cc:  Eliminate preprocessor flag RUN and all
+  code contained in "#else" blocks.
+
+- G4LorentzConverter.cc:  Add diagnostic output with different levels:  All
+  " >>>" entry messages done for verboseLevel >= 1; all numeric reports for
+  verboseLevel > 2.
+
+- G4NucleiModel: Replace preprocessor flag "CHC_CHECK" with
+  G4CASCADE_DEBUG_CHARGE.
+
+- G4PreCompoundCascadeInterface.cc: Replace preprocessor flag
+  "debug_G4...Interface" with G4CASCADE_DEBUG_INTERFACE.  Fix bug with
+  information access in debugging block.
+
+- G4WatcherGun.cc: Replace preprocessor flag "PB" with G4CASCADE_WATCHER_HIGHZ.
+
+16 June 2010 Michael Kelsey
+---------------------------
+- G4CascadeInterface.hh:  Add setVerboseLevel() function.
+
+- G4IntraNucleiCascader.cc:  Add reporting of momentum content for residual
+  nuclear fragment, list of final-state particles, some cleanup of momentum
+  balancing code.
+
+- G4LorentzConvertor:  Add reporting of bullet and target four-vectors.
+
+- G4NucleiModel.cc:  Add diagnostic output for partner-list generation,
+  replace one random-angle code block with generateWithRandomAngles().
+
+15 June 2010 Michael Kelsey (hadr-casc-V09-03-45)
+-------------------------------------------------
+- G4CascadeInterface.cc, G4PreCompoundCascadeInterface.cc:  When converting
+  output particles from K0/K0bar to K0S/K0L, convert kinetic energy from
+  Bertini internal units (GeV) to GEANT4 units (MeV).
+
+- G4VCascadeCollider:  Make two separate ctors, name and name+verbose, with
+  no default values.  All subclasses *must* provide a literal name string.
+
+12 June 2010 Michael Kelsey (hadr-casc-V09-03-44)
+-------------------------------------------------
+- G4CascadeT33piNChannel.cc:  Swapped 8-body final state tables for pipP
+  vs. pimN.
+
+11 June 2010 Michael Kelsey (hadr-casc-V09-03-43)
+-------------------------------------------------
 - G4CascadeData.hh:  Equally trivial fix for compiler warking on Intel ICC
   for index[] array bounds in ::initialize().  Dimension array as [9], and