Changeset 962 for trunk/source/processes/hadronic/models/de_excitation/evaporation/src/G4Evaporation.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/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.7 2007/02/14 13:37:49 ahoward Exp $
-// GEANT4 tag $Name:  $
+// $Id: G4Evaporation.cc,v 1.12 2008/12/09 17:57:36 ahoward Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // Hadronic Process: Nuclear De-excitations
@@ -33 +33 @@
 // Alex Howard - added protection for negative probabilities in the sum, 14/2/07
 //
+// Modif (03 September 2008) by J. M. Quesada for external choice of inverse 
+// cross section option
+// JMQ (06 September 2008) Also external choices have been added for 
+// superimposed Coulomb barrier (if useSICBis set true, by default is false) 
+
 #include "G4Evaporation.hh"
 #include "G4EvaporationFactory.hh"
@@ -89 +94 @@
 }
 
-
 G4FragmentVector * G4Evaporation::BreakItUp(const G4Fragment &theNucleus)
 {
@@ -105 +109 @@
     // Number of channels
     G4int TotNumberOfChannels = theChannels->size();  
-        
 
     // Starts loop over evaporated particles
     for (;;) 
-      {
+ 
+   {    
         // loop over evaporation channels
         std::vector<G4VEvaporationChannel*>::iterator i;
         for (i=theChannels->begin(); i != theChannels->end(); i++) 
           {
+  // for inverse cross section choice
+            (*i)->SetOPTxs(OPTxs);
+  // for superimposed Coulomb Barrier for inverse cross sections
+            (*i)->UseSICB(useSICB);
+
             (*i)->Initialize(theResidualNucleus);
           }
@@ -172 +181 @@
             if( j >= TotNumberOfChannels ) 
               {
+                G4cerr << " Residual A: " << theResidualNucleus.GetA() << " Residual Z: " << theResidualNucleus.GetZ() << " Excitation Energy: " << theResidualNucleus.GetExcitationEnergy() << G4endl;
+                G4cerr << " j has not chosen a channel, j = " << j << " TotNumberOfChannels " << TotNumberOfChannels << " Total Probability: " << TotalProbability << G4endl;
+                for (j=0; j < TotNumberOfChannels; j++) 
+                  {
+                    G4cerr << " j: " << j << " EmissionProbChannel: " << EmissionProbChannel[j] << " and shoot: " << shoot << " (<ProbChannel?) " << G4endl;
+                  }             
                 throw G4HadronicException(__FILE__, __LINE__,  "G4Evaporation::BreakItUp: Can't define emission probability of the channels!" );
               } 
@@ -216 +231 @@
                   theResidualNucleus.SetCreatorModel(G4String("ResidualNucleus"));
 #endif
+
                 }
               }