71 char* env = getenv(
"G4UsePhotonEvaporationOLD");
119 (*theChannels)[i]->SetOPTxs(
OPTxs);
120 (*theChannels)[i]->UseSICB(
useSICB);
121 (*theChannels)[i]->Initialise();
157 G4double totprob, prob, oldprob = 0.0;
158 size_t maxchannel, i;
164 for(
G4int ia=0; ia<Amax; ++ia) {
175 theResult->push_back(theResidualNucleus);
188 if(Eex <= minExcitation && abun > 0.0) {
189 theResult->push_back(theResidualNucleus);
201 prob = (*theChannels)[i]->GetEmissionProbability(theResidualNucleus);
209 if(prob <= totprob*1.e-8 && oldprob <= totprob*1.e-8) {
221 (*theChannels)[0]->BreakUpChain(theResult, theResidualNucleus);
229 if(0.0 == abun && Z < 20) {
233 theResult->push_back(theResidualNucleus);
241 for(i=0; i<maxchannel; ++i) {
if(
probabilities[i] >= totprob) {
break; } }
244 G4Fragment* frag = (*theChannels)[i]->EmittedFragment(theResidualNucleus);
246 if(frag) { theResult->push_back(frag); }
249 theResult->push_back(theResidualNucleus);
255 theResult->push_back(theResidualNucleus);
G4FragmentVector * BreakItUp(const G4Fragment &theNucleus)
void BreakFragment(G4FragmentVector *, G4Fragment *theNucleus)
virtual std::vector< G4VEvaporationChannel * > * GetChannel()=0
std::vector< G4double > probabilities
G4VEvaporationChannel * thePhotonEvaporation
static G4NistManager * Instance()
std::vector< G4VEvaporationChannel * > * theChannels
virtual void InitialiseChannels()
G4bool IsApplicable(G4int Z, G4int A, G4double mass) const
G4IonTable * GetIonTable() const
double A(double temperature)
std::vector< G4Fragment * > G4FragmentVector
G4IonTable * theTableOfIons
G4double GetIsotopeAbundance(G4int Z, G4int N) const
G4double GetGroundStateMass() const
void SetCombinedChannel()
G4FermiFragmentsPool * thePool
static G4ParticleTable * GetParticleTable()
void InitialiseEvaporation()
virtual void SetPhotonEvaporation(G4VEvaporationChannel *ptr)
G4UnstableFragmentBreakUp unstableBreakUp
virtual G4bool BreakUpChain(G4FragmentVector *theResult, G4Fragment *theNucleus)
static G4FermiFragmentsPool * Instance()
G4VEvaporationFactory * theChannelFactory
G4double GetExcitationEnergy() const