Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4Evaporation Class Reference

#include <G4Evaporation.hh>

Inheritance diagram for G4Evaporation:
Collaboration diagram for G4Evaporation:

Public Member Functions

 G4Evaporation (G4VEvaporationChannel *photoEvaporation=nullptr)
 
virtual ~G4Evaporation ()
 
virtual void InitialiseChannels () final
 
virtual void BreakFragment (G4FragmentVector *, G4Fragment *theNucleus) final
 
void SetDefaultChannel ()
 
void SetGEMChannel ()
 
void SetCombinedChannel ()
 
- Public Member Functions inherited from G4VEvaporation
 G4VEvaporation ()
 
virtual ~G4VEvaporation ()
 
virtual void SetPhotonEvaporation (G4VEvaporationChannel *ptr)
 
void SetFermiBreakUp (G4VFermiBreakUp *ptr)
 
G4VFermiBreakUpGetFermiBreakUp () const
 
G4VEvaporationChannelGetPhotonEvaporation ()
 
G4VEvaporationChannelGetFissionChannel ()
 
void SetOPTxs (G4int opt)
 
void UseSICB (G4bool use)
 
size_t GetNumberOfChannels () const
 

Additional Inherited Members

- Protected Member Functions inherited from G4VEvaporation
void CleanChannels ()
 
- Protected Attributes inherited from G4VEvaporation
G4VEvaporationChannelthePhotonEvaporation
 
G4VFermiBreakUptheFBU
 
G4int OPTxs
 
G4bool useSICB
 
std::vector
< G4VEvaporationChannel * > * 
theChannels
 
G4VEvaporationFactorytheChannelFactory
 

Detailed Description

Definition at line 61 of file G4Evaporation.hh.

Constructor & Destructor Documentation

G4Evaporation::G4Evaporation ( G4VEvaporationChannel photoEvaporation = nullptr)
explicit

Definition at line 66 of file G4Evaporation.cc.

67  : G4VEvaporation(),fVerbose(0),nChannels(0),minExcitation(0.1*keV),
68  isInitialised(false)
69 {
70  if(photoEvaporation) { SetPhotonEvaporation(photoEvaporation); }
72 
75  InitialiseChannelFactory();
76 
77  theTableOfIons = G4ParticleTable::GetParticleTable()->GetIonTable();
78  nist = G4NistManager::Instance();
79  unstableBreakUp = new G4UnstableFragmentBreakUp();
80 }
G4VEvaporationChannel * thePhotonEvaporation
static G4NistManager * Instance()
G4IonTable * GetIonTable() const
static G4ParticleTable * GetParticleTable()
virtual void SetPhotonEvaporation(G4VEvaporationChannel *ptr)
G4VEvaporationFactory * theChannelFactory
static constexpr double keV
Definition: G4SIunits.hh:216

Here is the call graph for this function:

G4Evaporation::~G4Evaporation ( )
virtual

Definition at line 82 of file G4Evaporation.cc.

83 {
84  delete unstableBreakUp;
85 }

Member Function Documentation

void G4Evaporation::BreakFragment ( G4FragmentVector theResult,
G4Fragment theNucleus 
)
finalvirtual

Reimplemented from G4VEvaporation.

Definition at line 151 of file G4Evaporation.cc.

153 {
154  if(!isInitialised) { InitialiseChannels(); }
155 
156  G4double totprob, prob, oldprob = 0.0;
157  size_t maxchannel, i;
158 
159  G4int Amax = theResidualNucleus->GetA_asInt();
160 
161  // Starts loop over evaporated particles, loop is limited by number
162  // of nucleons
163  for(G4int ia=0; ia<Amax; ++ia) {
164 
165  // g,n,p and light fragments - evaporation is finished
166  G4int Z = theResidualNucleus->GetZ_asInt();
167  G4int A = theResidualNucleus->GetA_asInt();
168  G4double Eex = theResidualNucleus->GetExcitationEnergy();
169 
170  // stop deecitation loop if residual can be deexcited by FBU
171  if(theFBU->IsApplicable(Z, A, Eex)) { break; }
172 
173  // check if it is stable, then finish evaporation
174  G4double abun = nist->GetIsotopeAbundance(Z, A);
175  if(fVerbose > 0) {
176  G4cout << "### G4Evaporation::BreakItUp step " << ia << " Z= " << Z
177  << " A= " << A << " Eex(MeV)= "
178  << theResidualNucleus->GetExcitationEnergy()
179  << " aban= " << abun << G4endl;
180  }
181  // stop deecitation loop in the case of a cold stable fragment
182  if(Eex <= minExcitation && abun > 0.0) { break; }
183 
184  totprob = 0.0;
185  maxchannel = nChannels;
186  if(fVerbose > 1) {
187  G4cout << "### Evaporation loop #" << ia
188  << " Fragment: " << theResidualNucleus << G4endl;
189  }
190  // loop over evaporation channels
191  for(i=0; i<nChannels; ++i) {
192  prob = (*theChannels)[i]->GetEmissionProbability(theResidualNucleus);
193  if(fVerbose > 0 && prob > 0.0) {
194  G4cout << " Channel# " << i << " prob= " << prob << G4endl;
195  }
196  totprob += prob;
197  probabilities[i] = totprob;
198 
199  // if two recent probabilities are near zero stop computations
200  if(i>=8 && prob > 0.0) {
201  if(prob <= totprob*1.e-8 && oldprob <= totprob*1.e-8) {
202  maxchannel = i+1;
203  break;
204  }
205  }
206  oldprob = prob;
207  }
208 
209  // photon evaporation in the case of no other channels available
210  // do evaporation chain and reset total probability
211  if(0.0 < totprob && probabilities[0] == totprob) {
212  if(fVerbose > 1) {
213  G4cout << "Start chain of gamma evaporation" << G4endl;
214  }
215  (*theChannels)[0]->BreakUpChain(theResult, theResidualNucleus);
216  totprob = 0.0;
217  }
218 
219  // stable fragment - evaporation is finished
220  if(0.0 == totprob) {
221 
222  // if fragment is exotic, then force its decay
223  if(0.0 == abun) {
224  if(!unstableBreakUp->BreakUpChain(theResult, theResidualNucleus))
225  { break; }
226  if(fVerbose > 1) { G4cout << "$$$ Decay exotic fragment" << G4endl; }
227  } else {
228  break;
229  }
230  }
231 
232  // select channel
233  totprob *= G4UniformRand();
234  // loop over evaporation channels
235  for(i=0; i<maxchannel; ++i) { if(probabilities[i] >= totprob) { break; } }
236 
237  if(fVerbose > 1) { G4cout << "Channel # " << i << G4endl; }
238  G4Fragment* frag = (*theChannels)[i]->EmittedFragment(theResidualNucleus);
239  if(fVerbose > 1 && frag) { G4cout << " " << *frag << G4endl; }
240 
241  // normaly a fragment should be created
242  if(frag) { theResult->push_back(frag); }
243  else { break; }
244  }
245  // loop is stopped, residual fragment is added to the results
246  theResult->push_back(theResidualNucleus);
247 }
G4VFermiBreakUp * theFBU
virtual G4bool IsApplicable(G4int Z, G4int A, G4double mass) const =0
virtual G4bool BreakUpChain(G4FragmentVector *, G4Fragment *) final
int G4int
Definition: G4Types.hh:78
virtual void InitialiseChannels() final
#define G4UniformRand()
Definition: Randomize.hh:97
G4GLOB_DLL std::ostream G4cout
double A(double temperature)
G4double GetIsotopeAbundance(G4int Z, G4int N) const
#define G4endl
Definition: G4ios.hh:61
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

void G4Evaporation::InitialiseChannels ( )
finalvirtual

Reimplemented from G4VEvaporation.

Definition at line 87 of file G4Evaporation.cc.

88 {
89  if(isInitialised) { return; }
90 
91  G4DeexPrecoParameters* param =
93  minExcitation = param->GetMinExcitation();
94 
95  G4DeexChannelType type = param->GetDeexChannelsType();
96  if(type == fCombined) { SetCombinedChannel(); }
97  else if(type == fGEM) { SetGEMChannel(); }
98 
99  isInitialised = true;
100 }
G4DeexChannelType GetDeexChannelsType() const
void SetGEMChannel()
G4DeexPrecoParameters * GetParameters()
void SetCombinedChannel()
static G4NuclearLevelData * GetInstance()
G4double GetMinExcitation() const

Here is the call graph for this function:

Here is the caller graph for this function:

void G4Evaporation::SetCombinedChannel ( )

Definition at line 140 of file G4Evaporation.cc.

141 {
142  if(fCombined != channelType) {
144  CleanChannels();
145  delete theChannelFactory;
147  InitialiseChannelFactory();
148  }
149 }
G4VEvaporationChannel * thePhotonEvaporation
G4VEvaporationFactory * theChannelFactory

Here is the call graph for this function:

Here is the caller graph for this function:

void G4Evaporation::SetDefaultChannel ( )

Definition at line 118 of file G4Evaporation.cc.

119 {
120  if(fEvaporation != channelType) {
122  CleanChannels();
123  delete theChannelFactory;
125  InitialiseChannelFactory();
126  }
127 }
G4VEvaporationChannel * thePhotonEvaporation
G4VEvaporationFactory * theChannelFactory

Here is the call graph for this function:

Here is the caller graph for this function:

void G4Evaporation::SetGEMChannel ( )

Definition at line 129 of file G4Evaporation.cc.

130 {
131  if(fGEM != channelType) {
133  CleanChannels();
134  delete theChannelFactory;
136  InitialiseChannelFactory();
137  }
138 }
G4VEvaporationChannel * thePhotonEvaporation
G4VEvaporationFactory * theChannelFactory

Here is the call graph for this function:

Here is the caller graph for this function:


The documentation for this class was generated from the following files: