Changeset 1055 for trunk/source/processes/hadronic/models/de_excitation/fission/src/G4FissionProbability.cc

Timestamp:

May 28, 2009, 4:26:57 PM (15 years ago)

Author:

garnier

Message:

maj sur la beta de geant 4.9.3

File:

• trunk/source/processes/hadronic/models/de_excitation/fission/src/G4FissionProbability.cc (modified) (2 diffs)

trunk/source/processes/hadronic/models/de_excitation/fission/src/G4FissionProbability.cc

@@ -25 +25 @@
 //
 //
-// $Id: G4FissionProbability.cc,v 1.8 2007/02/12 09:39:58 ahoward Exp $
-// GEANT4 tag $Name: geant4-09-02 $
+// $Id: G4FissionProbability.cc,v 1.9 2009/02/15 17:03:25 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-03-beta-cand-01 $
 //
 // Hadronic Process: Nuclear De-excitations
 // by V. Lara (Oct 1998)
 //
+//
+// J.M.Quesada (14 february 2009) bug fixed in fission width: missing parenthesis in the denominator 
 
 
@@ -95 +97 @@
     G4double Exp1 = 0.0;
     if (SystemEntropy <= 160.0) Exp1 = std::exp(-SystemEntropy);
-    // @@@@@@@@@@@@@@@@@ hpw changed max to min - cannot notify vicente now since cern mail gave up on me...
+    // @@@@@@@@@@@@@@@@@ hpw changed max to min - cannot notify vicente now
     G4double Exp2 = std::exp( std::min(700.0,Cf-SystemEntropy) ); 
 
+    // JMQ 14/02/09 BUG fixed in fission probability (missing parenthesis at denominator)
     //AH fix from Vincente:    G4double probability = (Exp1 + (1.0-Cf)*Exp2) / 4.0*pi*afission;
-    G4double probability = (Exp1 + (Cf-1.0)*Exp2) / 4.0*pi*afission;
+    //    G4double probability = (Exp1 + (Cf-1.0)*Exp2) / 4.0*pi*afission;
+    G4double probability = (Exp1 + (Cf-1.0)*Exp2) / (4.0*pi*afission);
+
+
     
     return probability;