Changeset 962 for trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/src/G4VPreCompoundFragment.cc

Timestamp:

Apr 6, 2009, 12:30:29 PM (15 years ago)

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/src/G4VPreCompoundFragment.cc (modified) (4 diffs)

trunk/source/processes/hadronic/models/pre_equilibrium/exciton_model/src/G4VPreCompoundFragment.cc

@@ -24 +24 @@
 // ********************************************************************
 //
+// $Id: G4VPreCompoundFragment.cc,v 1.12 2009/02/10 16:01:37 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
-// $Id: G4VPreCompoundFragment.cc,v 1.4 2007/07/23 09:56:40 ahoward Exp $
-// GEANT4 tag $Name:  $
+// J. M. Quesada (August 2008).  
+// Based  on previous work by V. Lara
 //
-// by V. Lara
  
 #include "G4VPreCompoundFragment.hh"
@@ -142 +143 @@
   if ((theRestNucleusA < theRestNucleusZ) ||
       (theRestNucleusA < theA) ||
-      (theRestNucleusZ < theZ) ||
-      (aFragment.GetNumberOfCharged() < theZ)) // AH last argument from JMQ
+      (theRestNucleusZ < theZ)) 
     {
       // In order to be sure that emission probability will be 0.
@@ -155 +155 @@
     GetCoulombBarrier(static_cast<G4int>(theRestNucleusA),static_cast<G4int>(theRestNucleusZ),
                       aFragment.GetExcitationEnergy());
-  
+
+
   // Compute Binding Energies for fragments 
   // (needed to separate a fragment from the nucleus)
@@ -164 +165 @@
   
   // Compute Maximal Kinetic Energy which can be carried by fragments after separation
+  // This is the true (assimptotic) maximal kinetic energy
   G4double m = aFragment.GetMomentum().m();
   G4double rm = GetRestNuclearMass();
   G4double em = GetNuclearMass();
   theMaximalKineticEnergy = ((m - rm)*(m + rm) + em*em)/(2.0*m) - em;
-  // - theCoulombBarrier;
+ 
   
   return;