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

Timestamp:

Dec 22, 2010, 3:52:27 PM (13 years ago)

Author:

garnier

Message:

geant4 tag 9.4

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.9 2009/02/15 17:03:25 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-03-ref-09 $
+// $Id: G4FissionProbability.cc,v 1.10 2010/11/17 20:22:46 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-04-ref-00 $
 //
 // Hadronic Process: Nuclear De-excitations
@@ -38 +38 @@
 #include "G4PairingCorrection.hh"
 
+G4FissionProbability::G4FissionProbability()
+{}
+
+G4FissionProbability::~G4FissionProbability()
+{}
 
 
-G4FissionProbability::G4FissionProbability(const G4FissionProbability &) : G4VEmissionProbability()
+G4double 
+G4FissionProbability::EmissionProbability(const G4Fragment & fragment, 
+                                          G4double MaximalKineticEnergy)
+  // Compute integrated probability of fission channel
 {
-    throw G4HadronicException(__FILE__, __LINE__, "G4FissionProbability::copy_constructor meant to not be accessable");
+  if (MaximalKineticEnergy <= 0.0) return 0.0;
+  G4int A = fragment.GetA_asInt();
+  G4int Z = fragment.GetZ_asInt();
+  G4double U = fragment.GetExcitationEnergy();
+  
+  G4double Ucompound = U - 
+    G4PairingCorrection::GetInstance()->GetPairingCorrection(A,Z);
+
+  G4double Ufission = U - 
+    G4PairingCorrection::GetInstance()->GetFissionPairingCorrection(A,Z);
+  
+  G4double SystemEntropy = 
+    2.0*std::sqrt(theEvapLDP.LevelDensityParameter(A,Z,Ucompound)*Ucompound);
+        
+  G4double afission = theFissLDP.LevelDensityParameter(A,Z,Ufission);
+
+  G4double Cf = 2.0*std::sqrt(afission*MaximalKineticEnergy);
+
+  //    G4double Q1 = 1.0 + (Cf - 1.0)*std::exp(Cf);
+  //    G4double Q2 = 4.0*pi*afission*std::exp(SystemEntropy);
+  
+  //    G4double probability = Q1/Q2;
+   
+  G4double Exp1 = 0.0;
+  if (SystemEntropy <= 160.0) { Exp1 = std::exp(-SystemEntropy); }
+  // @@@@@@@@@@@@@@@@@ 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);
+
+  return probability;
 }
 
-
-
-
-const G4FissionProbability & G4FissionProbability::operator=(const G4FissionProbability &)
-{
-    throw G4HadronicException(__FILE__, __LINE__, "G4FissionProbability::operator= meant to not be accessable");
-    return *this;
-}
-
-
-G4bool G4FissionProbability::operator==(const G4FissionProbability &) const
-{
-    return false;
-}
-
-G4bool G4FissionProbability::operator!=(const G4FissionProbability &) const
-{
-    return true;
-}
-
-
-G4double G4FissionProbability::EmissionProbability(const G4Fragment & fragment, const G4double MaximalKineticEnergy)
-    // Compute integrated probability of fission channel
-{
-    if (MaximalKineticEnergy <= 0.0) return 0.0;
-    G4double A = fragment.GetA();
-    G4double Z = fragment.GetZ();
-    G4double U = fragment.GetExcitationEnergy();
-  
-    G4double Ucompound = U - G4PairingCorrection::GetInstance()->GetPairingCorrection(static_cast<G4int>(A),
-                                                                                      static_cast<G4int>(Z));
-    G4double Ufission = U - G4PairingCorrection::GetInstance()->GetFissionPairingCorrection(static_cast<G4int>(A),
-                                                                                            static_cast<G4int>(Z));
-  
-    G4double SystemEntropy = 2.0*std::sqrt(theEvapLDP.LevelDensityParameter(static_cast<G4int>(A),
-                                                                       static_cast<G4int>(Z),
-                                                                       Ucompound)*Ucompound);
-        
-    G4double afission = theFissLDP.LevelDensityParameter(static_cast<G4int>(A),
-                                                         static_cast<G4int>(Z),
-                                                         Ufission);
-
-    G4double Cf = 2.0*std::sqrt(afission*MaximalKineticEnergy);
-
-    //    G4double Q1 = 1.0 + (Cf - 1.0)*std::exp(Cf);
-    //    G4double Q2 = 4.0*pi*afission*std::exp(SystemEntropy);
-        
-    //    G4double probability = Q1/Q2;
- 
-    
-    G4double Exp1 = 0.0;
-    if (SystemEntropy <= 160.0) Exp1 = std::exp(-SystemEntropy);
-    // @@@@@@@@@@@@@@@@@ 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);
-
-
-    
-    return probability;
-}
-