Changeset 962 for trunk/source/processes/hadronic/models/de_excitation/fission/src/G4FissionProbability.cc
- Timestamp:
- Apr 6, 2009, 12:30:29 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/de_excitation/fission/src/G4FissionProbability.cc
r819 r962 25 25 // 26 26 // 27 // $Id: G4FissionProbability.cc,v 1. 8 2007/02/12 09:39:58 ahowardExp $28 // GEANT4 tag $Name: $27 // $Id: G4FissionProbability.cc,v 1.9 2009/02/15 17:03:25 vnivanch Exp $ 28 // GEANT4 tag $Name: geant4-09-02-ref-02 $ 29 29 // 30 30 // Hadronic Process: Nuclear De-excitations 31 31 // by V. Lara (Oct 1998) 32 32 // 33 // 34 // J.M.Quesada (14 february 2009) bug fixed in fission width: missing parenthesis in the denominator 33 35 34 36 … … 95 97 G4double Exp1 = 0.0; 96 98 if (SystemEntropy <= 160.0) Exp1 = std::exp(-SystemEntropy); 97 // @@@@@@@@@@@@@@@@@ hpw changed max to min - cannot notify vicente now since cern mail gave up on me...99 // @@@@@@@@@@@@@@@@@ hpw changed max to min - cannot notify vicente now 98 100 G4double Exp2 = std::exp( std::min(700.0,Cf-SystemEntropy) ); 99 101 102 // JMQ 14/02/09 BUG fixed in fission probability (missing parenthesis at denominator) 100 103 //AH fix from Vincente: G4double probability = (Exp1 + (1.0-Cf)*Exp2) / 4.0*pi*afission; 101 G4double probability = (Exp1 + (Cf-1.0)*Exp2) / 4.0*pi*afission; 104 // G4double probability = (Exp1 + (Cf-1.0)*Exp2) / 4.0*pi*afission; 105 G4double probability = (Exp1 + (Cf-1.0)*Exp2) / (4.0*pi*afission); 106 107 102 108 103 109 return probability;
Note: See TracChangeset
for help on using the changeset viewer.