Changeset 962 for trunk/source/processes/hadronic/models/de_excitation/multifragmentation/src/G4StatMFChannel.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/G4StatMFChannel.cc (modified) (3 diffs)

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

@@ -25 +25 @@
 //
 //
-// $Id: G4StatMFChannel.cc,v 1.6.2.1 2008/04/23 16:19:50 gcosmo Exp $
-// GEANT4 tag $Name: geant4-09-01-patch-02 $
+// $Id: G4StatMFChannel.cc,v 1.10 2008/11/19 14:33:31 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // Hadronic Process: Nuclear De-excitations
 // by V. Lara
+//
+// Modified:
+// 25.07.08 I.Pshenichnov (in collaboration with Alexander Botvina and Igor 
+//          Mishustin (FIAS, Frankfurt, INR, Moscow and Kurchatov Institute, 
+//          Moscow, pshenich@fias.uni-frankfurt.de) fixed semi-infinite loop 
+
 
 #include "G4StatMFChannel.hh"
@@ -92 +98 @@
         G4int A = static_cast<G4int>((*i)->GetA());
         G4int Z = static_cast<G4int>((*i)->GetZ());
-        if ( (A > 1 && (Z >= A || Z <= 0)) || (A==1 && Z > A) || A <= 0 ) return false;
+        if ( (A > 1 && (Z > A || Z <= 0)) || (A==1 && Z > A) || A <= 0 ) return false;
     }
     
@@ -315 +321 @@
         G4double CosTheta1;
         G4double Sign;
-        do 
-          {
-            do 
-              {
-                CosTheta1 = 1.0 - 2.0*G4UniformRand();
-              } 
-            while (CosTheta1*CosTheta1 < CTM12);
-          }
-        while (CTM12 >= 0.0 && CosTheta1 < 0.0);
+
+        if (CTM12 > 0.9999) {CosTheta1 = 1.;} 
+        else {
+         do 
+           {
+             do 
+               {
+                 CosTheta1 = 1.0 - 2.0*G4UniformRand();
+               } 
+             while (CosTheta1*CosTheta1 < CTM12);
+           }
+          while (CTM12 >= 0.0 && CosTheta1 < 0.0);
+        }
 
         if (CTM12 < 0.0) Sign = 1.0;