Changeset 1347 for trunk/source/processes/hadronic/models/de_excitation/evaporation/src/G4Evaporation.cc

Timestamp:

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

Author:

garnier

Message:

geant4 tag 9.4

File:

• trunk/source/processes/hadronic/models/de_excitation/evaporation/src/G4Evaporation.cc (modified) (6 diffs)

trunk/source/processes/hadronic/models/de_excitation/evaporation/src/G4Evaporation.cc

@@ -25 +25 @@
 //
 //
-// $Id: G4Evaporation.cc,v 1.25 2010/06/09 11:56:47 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-03-ref-09 $
+// $Id: G4Evaporation.cc,v 1.26 2010/11/23 18:10:10 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-04-ref-00 $
 //
 // Hadronic Process: Nuclear De-excitations
@@ -142 +142 @@
     G4int Z = theResidualNucleus->GetZ_asInt();
     G4double abun = nist->GetIsotopeAbundance(Z, A); 
-    /*
-    G4cout << "### G4Evaporation::BreakItUp step " << ia << " Z= " << Z
-           << " A= " << A << " Eex(MeV)= " 
-           << theResidualNucleus->GetExcitationEnergy()
-           << " aban= " << abun << G4endl;
-    */
+    
+    // G4cout << "### G4Evaporation::BreakItUp step " << ia << " Z= " << Z
+    //     << " A= " << A << " Eex(MeV)= " 
+    //     << theResidualNucleus->GetExcitationEnergy()
+    //     << " aban= " << abun << G4endl;
+   
     if(theResidualNucleus->GetExcitationEnergy() <= minExcitation && 
        (abun > 0.0)) 
@@ -184 +184 @@
     // do evaporation chain and reset total probability
     if(0.0 < totprob && probabilities[0] == totprob) {
+      //G4cout << "Start gamma evaporation" << G4endl;
       theTempResult = (*theChannels)[0]->BreakUpFragment(theResidualNucleus);
       if(theTempResult) {
@@ -228 +229 @@
     // single photon evaporation, primary pointer is kept
     if(0 == i) {
+      //G4cout << "Single gamma" << G4endl;
       G4Fragment* gamma = (*theChannels)[0]->EmittedFragment(theResidualNucleus);
       if(gamma) { theResult->push_back(gamma); }
@@ -233 +235 @@
       // fission, return results to the main loop if fission is succesful
     } else if(1 == i) {
+      //G4cout << "Fission" << G4endl;
       theTempResult = (*theChannels)[1]->BreakUp(*theResidualNucleus);
       if(theTempResult) {
@@ -248 +251 @@
       // other channels
     } else {
+      //G4cout << "Channel # " << i << G4endl;
       theTempResult = (*theChannels)[i]->BreakUp(*theResidualNucleus);
       if(theTempResult) {