117 G4cout <<
"### G4Evaporation::InitialiseChannelFactory for "
121 (*theChannels)[i]->SetOPTxs(
OPTxs);
122 (*theChannels)[i]->Initialise();
184 G4double totprob, prob, oldprob = 0.0;
185 size_t maxchannel, i;
189 G4cout <<
"### G4Evaporation::BreakItUp loop" <<
G4endl;
194 for(
G4int ia=0; ia<Amax; ++ia) {
199 if(A <= 1) {
break; }
209 (abun > 0.0 || (A == 3 && (Z == 1 || Z == 2)))) {
break; }
214 G4cout <<
"Evaporation# " << ia <<
" Z= " << Z <<
" A= " << A
216 <<
" aban= " << abun <<
G4endl;
220 prob = (*theChannels)[i]->GetEmissionProbability(theResidualNucleus);
222 G4cout <<
" Channel# " << i <<
" prob= " << prob <<
G4endl;
228 if(i>=8 && prob > 0.0) {
229 if(prob <= totprob*1.
e-8 && oldprob <= totprob*1.
e-8) {
241 G4cout <<
"$$$ Start chain of gamma evaporation" <<
G4endl;
243 (*theChannels)[0]->BreakUpChain(theResult, theResidualNucleus);
268 for(i=0; i<maxchannel; ++i) {
if(
probabilities[i] >= totprob) {
break; } }
271 G4Fragment* frag = (*theChannels)[i]->EmittedFragment(theResidualNucleus);
275 if(frag) { theResult->push_back(frag); }