72 useGNASHTransition(false), OPTxs(3), useSICB(false),
73 useNGB(false), useSCO(false), useCEMtr(true), maxZ(3), maxA(5)
114 ed <<
"G4PreCompoundModel is used for ";
123 if(primary ==
proton) { Zp = 1; }
148 for(G4ReactionProductVector::iterator i= result->begin();
149 i != result->end(); ++i)
153 (*i)->GetTotalEnergy(),
154 (*i)->GetMomentum());
177 if ((Z <
maxZ && A <
maxA) || Eex < MeV /*|| Eex > 3.*
MeV*A*/) {
184 const G4int countmax = 1000;
188 G4int EquilibriumExcitonNumber =
199 G4bool ThereIsTransition =
false;
215 if (test <= EquilibriumExcitonNumber) { go_ahead=
true; }
226 G4double TotalTransitionProbability =
237 if(!go_ahead || P1 <= P2+P3 ||
263 G4double TotalProbability = TotalEmissionProbability
264 + TotalTransitionProbability;
267 if (TotalProbability*
G4UniformRand() > TotalEmissionProbability)
271 ThereIsTransition =
true;
279 ThereIsTransition =
false;
284 }
while (ThereIsTransition);
287 if(count >= countmax) {
289 ed <<
"G4PreCompoundModel loop over " << countmax <<
" iterations; "
290 <<
"current G4Fragment: \n" << aFragment;
365 outFile <<
"The GEANT4 precompound model is considered as an extension of the\n"
366 <<
"hadron kinetic model. It gives a possibility to extend the low energy range\n"
367 <<
"of the hadron kinetic model for nucleon-nucleus inelastic collision and it \n"
368 <<
"provides a ”smooth” transition from kinetic stage of reaction described by the\n"
369 <<
"hadron kinetic model to the equilibrium stage of reaction described by the\n"
370 <<
"equilibrium deexcitation models.\n"
371 <<
"The initial information for calculation of pre-compound nuclear stage\n"
372 <<
"consists of the atomic mass number A, charge Z of residual nucleus, its\n"
373 <<
"four momentum P0 , excitation energy U and number of excitons n, which equals\n"
374 <<
"the sum of the number of particles p (from them p_Z are charged) and the number of\n"
376 <<
"At the preequilibrium stage of reaction, we follow the exciton model approach in ref. [1],\n"
377 <<
"taking into account the competition among all possible nuclear transitions\n"
378 <<
"with ∆n = +2, −2, 0 (which are defined by their associated transition probabilities) and\n"
379 <<
"the emission of neutrons, protons, deutrons, thritium and helium nuclei (also defined by\n"
380 <<
"their associated emission probabilities according to exciton model)\n"
382 <<
"[1] K.K. Gudima, S.G. Mashnik, V.D. Toneev, Nucl. Phys. A401 329 (1983)\n"
388 outFile <<
"description of precompound model as used with DeExcite()"
virtual void PerformTransition(G4Fragment &aFragment)=0
void UseDefaultEmission()
void PerformEquilibriumEmission(const G4Fragment &aFragment, G4ReactionProductVector *theResult) const
static G4double GetNuclearMass(const G4double A, const G4double Z)
virtual ~G4PreCompoundModel()
static const G4double * P1[nN]
std::ostringstream G4ExceptionDescription
G4PreCompoundEmission * theEmission
G4PreCompoundModel(G4ExcitationHandler *ptr=0)
const G4ParticleDefinition * neutron
void SetNumberOfHoles(G4int valueTot, G4int valueP=0)
G4bool useGNASHTransition
G4double GetTransitionProb2() const
void UseGNASHTransition()
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &thePrimary, G4Nucleus &theNucleus)
const G4String & GetParticleName() const
G4HadFinalState theResult
void SetStatusChange(G4HadFinalStateStatus aS)
void SetExcitationHandler(G4ExcitationHandler *ptr)
std::vector< G4ReactionProduct * > G4ReactionProductVector
void SetNumberOfExcitedParticle(G4int valueTot, G4int valueP)
G4ExcitationHandler * GetExcitationHandler() const
double A(double temperature)
void UseDefaultTransition()
const G4ParticleDefinition * GetDefinition() const
G4VPreCompoundTransitions * theTransition
void UseCEMtr(G4bool use)
G4double GetGlobalTime() const
static G4Proton * Proton()
static G4Neutron * Neutron()
G4double GetTransitionProb1() const
const G4LorentzVector & Get4Momentum() const
virtual G4ReactionProductVector * DeExcite(G4Fragment &aFragment)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
virtual G4double CalculateProbability(const G4Fragment &aFragment)=0
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4ReactionProduct * PerformEmission(G4Fragment &aFragment)
void SetCreationTime(G4double time)
G4int GetNumberOfExcitons() const
const G4double x[NPOINTSGL]
virtual void ModelDescription(std::ostream &outFile) const
static const G4double * P2[nN]
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
G4double GetTransitionProb3() const
const G4ParticleDefinition * proton
virtual void DeExciteModelDescription(std::ostream &outFile) const
G4double GetLevelDensity()
G4double GetExcitationEnergy() const
G4double GetTotalProbability(const G4Fragment &aFragment)
CLHEP::HepLorentzVector G4LorentzVector