90 maxZForFermiBreakUp(9),maxAForFermiBreakUp(17),minEForMultiFrag(4*
GeV),
91 minExcitation(
keV),OPTxs(3),useSICB(false),isEvapLocal(true)
151 }
else if(exEnergy < minExcitation && nist->GetIsotopeAbundance(Z, A) > 0.0) {
167 size_t nsec = theTempResult->size();
175 G4bool deletePrimary =
true;
176 G4FragmentVector::iterator j;
177 for (j = theTempResult->begin(); j != theTempResult->end(); ++j) {
178 if((*j) == theInitialStatePtr) { deletePrimary =
false; }
179 A = (*j)->GetA_asInt();
187 G4double exEnergy1 = (*j)->GetExcitationEnergy();
191 Z = (*j)->GetZ_asInt();
203 if( deletePrimary ) {
delete theInitialStatePtr; }
205 delete theTempResult;
218 std::vector<G4Fragment*>::iterator iList;
237 for(
size_t j=0; j<nsec; ++j) {
238 exEnergy =
results[j]->GetExcitationEnergy();
258 for (
size_t j = 0; j<nsec; ++j) {
267 exEnergy =
results[j]->GetExcitationEnergy();
300 exEnergy = (*iList)->GetExcitationEnergy();
321 theReactionProductVector->reserve(
theResults.size() );
323 G4int theFragmentA, theFragmentZ;
329 theFragmentA = (*iList)->GetA_asInt();
330 theFragmentZ = (*iList)->GetZ_asInt();
331 G4double etot= (*iList)->GetMomentum().e();
334 if (theFragmentA == 0) {
335 theKindOfFragment = (*iList)->GetParticleDefinition();
336 }
else if (theFragmentA == 1 && theFragmentZ == 0) {
338 }
else if (theFragmentA == 1 && theFragmentZ == 1) {
340 }
else if (theFragmentA == 2 && theFragmentZ == 1) {
342 }
else if (theFragmentA == 3 && theFragmentZ == 1) {
344 }
else if (theFragmentA == 3 && theFragmentZ == 2) {
346 }
else if (theFragmentA == 4 && theFragmentZ == 2) {
351 eexc = (*iList)->GetExcitationEnergy();
360 if(theKindOfFragment) {
362 theNew->
SetMomentum((*iList)->GetMomentum().vect());
365 theReactionProductVector->push_back(theNew);
370 if(theKindOfFragment) {
373 if(etot <= ionmass) {
376 G4double ptot = std::sqrt((etot - ionmass)*(etot + ionmass));
377 mom = ((*iList)->GetMomentum().vect().unit())*ptot;
383 theReactionProductVector->push_back(theNew);
397 return theReactionProductVector;
460 outFile <<
"G4ExcitationHandler description\n"
461 <<
"This class samples de-excitation of excited nucleus using\n"
462 <<
"Fermi Break-up model for light fragments (Z < 9, A < 17), "
463 <<
"evaporation, fission, and photo-evaporation models. Evaporated\n"
464 <<
"particle may be proton, neutron, and other light fragment \n"
465 <<
"(Z < 13, A < 29). During photon evaporation produced gamma \n"
466 <<
"or electrons due to internal conversion \n";
static G4Pow * GetInstance()
static G4Triton * TritonDefinition()
static G4He3 * He3Definition()
void ModelDescription(std::ostream &outFile) const
std::vector< G4Fragment * > theResults
CLHEP::Hep3Vector G4ThreeVector
virtual G4bool BreakUpChain(G4FragmentVector *theResult, G4Fragment *theNucleus)
G4VEvaporationChannel * GetPhotonEvaporation()
void SetMomentum(const G4double x, const G4double y, const G4double z)
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
static G4Proton * ProtonDefinition()
virtual void BreakFragment(G4FragmentVector *, G4Fragment *theNucleus)
void SetMinEForMultiFrag(G4double anE)
virtual G4FragmentVector * BreakItUp(const G4Fragment &theNucleus)=0
virtual void Initialise()
static G4NistManager * Instance()
G4ReactionProductVector * BreakItUp(const G4Fragment &theInitialState)
std::vector< G4Fragment * > results
std::vector< G4ReactionProduct * > G4ReactionProductVector
std::vector< G4Fragment * > thePhotoEvapList
G4bool IsApplicable(G4int Z, G4int A, G4double mass) const
G4IonTable * GetIonTable() const
std::vector< G4Fragment * > theEvapList
void SetFermiModel(G4VFermiBreakUp *ptr)
void SetTotalEnergy(const G4double en)
std::vector< G4Fragment * > G4FragmentVector
void SetMultiFragmentation(G4VMultiFragmentation *ptr)
G4VEvaporationChannel * thePhotonEvaporation
G4FermiFragmentsPool * thePool
G4double GetIsotopeAbundance(G4int Z, G4int N) const
G4double GetGroundStateMass() const
G4int maxAForFermiBreakUp
G4VMultiFragmentation * theMultiFragmentation
static const G4double A[nN]
G4VFermiBreakUp * theFermiModel
virtual void BreakFragment(G4FragmentVector *, G4Fragment *theNucleus)=0
G4double minEForMultiFrag
void SetMaxZForFermiBreakUp(G4int aZ)
G4double GetPDGMass() const
static G4ParticleTable * GetParticleTable()
virtual void SetPhotonEvaporation(G4VEvaporationChannel *ptr)
void SetMaxAandZForFermiBreakUp(G4int anA, G4int aZ)
void SetEvaporation(G4VEvaporation *ptr)
G4VEvaporation * theEvaporation
void SetMaxAForFermiBreakUp(G4int anA)
virtual void Initialise()
void SetFormationTime(G4double aTime)
static G4FermiFragmentsPool * Instance()
void SetPhotonEvaporation(G4VEvaporationChannel *ptr)
static G4Deuteron * DeuteronDefinition()
static G4Alpha * AlphaDefinition()
static G4Neutron * NeutronDefinition()
G4int maxZForFermiBreakUp
G4double GetExcitationEnergy() const
G4IonTable * theTableOfIons