Changeset 962 for trunk/source/processes/hadronic/models/de_excitation/multifragmentation/src/G4StatMFMacroNucleon.cc

Timestamp:

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

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/hadronic/models/de_excitation/multifragmentation/src/G4StatMFMacroNucleon.cc (modified) (3 diffs)

trunk/source/processes/hadronic/models/de_excitation/multifragmentation/src/G4StatMFMacroNucleon.cc

@@ -25 +25 @@
 //
 //
-// $Id: G4StatMFMacroNucleon.cc,v 1.5 2006/06/29 20:25:19 gunter Exp $
-// GEANT4 tag $Name:  $
+// $Id: G4StatMFMacroNucleon.cc,v 1.6 2008/07/25 11:20:47 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // Hadronic Process: Nuclear De-excitations
@@ -72 +72 @@
     G4double exponent_neutron = mu/T;
 
-    if (exponent_neutron > 700.0) exponent_proton = 700.0;
+    if (exponent_neutron > 700.0) exponent_neutron = 700.0;
     if (exponent_proton > 700.0) exponent_proton = 700.0;
 
-    _NeutronMeanMultiplicity = (degeneracy*FreeVol/lambda3)*std::exp(mu/T);
+    _NeutronMeanMultiplicity = (degeneracy*FreeVol/lambda3)*std::exp(exponent_neutron);
         
-    _ProtonMeanMultiplicity = (degeneracy*FreeVol/lambda3)*
-        std::exp((mu+nu-Coulomb)/T);
+    _ProtonMeanMultiplicity = (degeneracy*FreeVol/lambda3)*std::exp(exponent_proton);
 
         
@@ -92 +91 @@
         (1.0 - 1.0/std::pow(1.0+G4StatMFParameters::GetKappaCoulomb(),1./3.));
                                                                         
-    return _Energy = Coulomb * theZARatio + (3./2.) * T;
+    return _Energy = Coulomb * theZARatio * theZARatio + (3./2.) * T;
                                                         
 }