#include <G4ParaFissionModel.hh>
Inheritance diagram for G4ParaFissionModel:
Public Member Functions | |
G4ParaFissionModel () | |
virtual | ~G4ParaFissionModel () |
virtual G4HadFinalState * | ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &theNucleus) |
Definition at line 44 of file G4ParaFissionModel.hh.
G4ParaFissionModel::G4ParaFissionModel | ( | ) | [inline] |
Definition at line 48 of file G4ParaFissionModel.hh.
References G4HadronicInteraction::SetMaxEnergy(), and G4HadronicInteraction::SetMinEnergy().
00049 { 00050 SetMinEnergy( 0.0 ); 00051 SetMaxEnergy( 60.*MeV ); 00052 }
virtual G4ParaFissionModel::~G4ParaFissionModel | ( | ) | [inline, virtual] |
virtual G4HadFinalState* G4ParaFissionModel::ApplyYourself | ( | const G4HadProjectile & | aTrack, | |
G4Nucleus & | theNucleus | |||
) | [inline, virtual] |
Implements G4HadronicInteraction.
Definition at line 56 of file G4ParaFissionModel.hh.
References G4HadFinalState::AddSecondary(), G4ExcitationHandler::BreakItUp(), G4CompetitiveFission::BreakUp(), G4HadFinalState::Clear(), G4HadProjectile::Get4Momentum(), G4Nucleus::GetA_asInt(), G4ParticleDefinition::GetBaryonNumber(), G4HadProjectile::GetDefinition(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetMomentum(), G4NucleiProperties::GetNuclearMass(), G4Fragment::GetParticleDefinition(), G4ParticleDefinition::GetPDGCharge(), G4Nucleus::GetZ_asInt(), G4HadFinalState::SetEnergyChange(), G4Fragment::SetNumberOfExcitedParticle(), G4Fragment::SetNumberOfHoles(), G4HadFinalState::SetStatusChange(), and stopAndKill.
00058 { 00059 theParticleChange.Clear(); 00060 theParticleChange.SetStatusChange( stopAndKill ); 00061 theParticleChange.SetEnergyChange( 0.0 ); 00062 00063 // prepare the fragment 00064 00065 G4int A = theNucleus.GetA_asInt(); 00066 G4int Z = theNucleus.GetZ_asInt(); 00067 G4double nucMass = G4NucleiProperties::GetNuclearMass(A, Z); 00068 00069 G4int numberOfEx = aTrack.GetDefinition()->GetBaryonNumber(); 00070 G4int numberOfCh = G4int(aTrack.GetDefinition()->GetPDGCharge() + 0.5); 00071 G4int numberOfHoles = 0; 00072 00073 A += numberOfEx; 00074 Z += numberOfCh; 00075 00076 G4LorentzVector v = aTrack.Get4Momentum() + G4LorentzVector(0.0,0.0,0.0,nucMass); 00077 G4Fragment anInitialState(A,Z,v); 00078 anInitialState.SetNumberOfExcitedParticle(numberOfEx,numberOfCh); 00079 anInitialState.SetNumberOfHoles(0,0); 00080 00081 // do the fission 00082 G4FragmentVector * theFissionResult = theFission.BreakUp(anInitialState); 00083 00084 // deexcite the fission fragments and fill result 00085 00086 G4int ll = theFissionResult->size(); 00087 for(G4int i=0; i<ll; i++) 00088 { 00089 G4ReactionProductVector* theExcitationResult = 0; 00090 G4Fragment* aFragment = (*theFissionResult)[i]; 00091 if(aFragment->GetExcitationEnergy() > keV) 00092 { 00093 theExcitationResult = theHandler.BreakItUp(*aFragment); 00094 00095 // add secondaries 00096 for(G4int j = 0; j < G4int(theExcitationResult->size()); j++) 00097 { 00098 G4ReactionProduct* rp0 = (*theExcitationResult)[j]; 00099 G4DynamicParticle* p0 = 00100 new G4DynamicParticle(rp0->GetDefinition(),rp0->GetMomentum()); 00101 theParticleChange.AddSecondary(p0); 00102 delete rp0; 00103 } 00104 delete theExcitationResult; 00105 } 00106 else 00107 { 00108 // add secondary 00109 G4DynamicParticle* p0 = 00110 new G4DynamicParticle(aFragment->GetParticleDefinition(), 00111 aFragment->GetMomentum()); 00112 theParticleChange.AddSecondary(p0); 00113 } 00114 delete aFragment; 00115 } 00116 00117 delete theFissionResult; 00118 00119 return &theParticleChange; 00120 }