86 delete theElementaryParticleCollider;
87 delete theIntraNucleiCascader;
88 delete theDeexcitation;
108 delete theDeexcitation;
114 delete theDeexcitation;
126 const G4int itry_max = 100;
131 G4cout <<
" InuclCollider -> particle on particle collision" <<
G4endl;
133 theElementaryParticleCollider->
collide(bullet, target, globalOutput);
144 G4cerr <<
" InuclCollider -> no collision possible " <<
G4endl;
153 G4cerr <<
" InuclCollider -> ERROR target is not a nucleus " <<
G4endl;
168 G4cerr <<
" InuclCollider -> ERROR bullet is not a hadron " <<
G4endl;
174 G4cerr <<
" InuclCollider -> ERROR can not collide with "
180 btype = pbullet->
type();
185 G4cerr <<
" InuclCollider -> ERROR bullet is not a nucleus " <<
G4endl;
190 ab = nbullet->
getA();
191 zb = nbullet->
getZ();
201 G4cout <<
" InuclCollider -> inelastic interaction is impossible\n"
202 <<
" due to the coulomb barirer " <<
G4endl;
211 G4cout <<
" degenerated? " << convertToTargetRestFrame.
trivial()
226 while (itry < itry_max) {
231 globalOutput.
reset();
234 theIntraNucleiCascader->
collide(zbullet, target, output);
242 G4cout <<
" itry " << itry <<
" finished, moving to lab frame" <<
G4endl;
247 globalOutput.
add(output);
254 G4cout <<
" InuclCollider output after trials " << itry <<
G4endl;
259 G4cerr <<
" InuclCollider setOnShell failed." <<
G4endl;
264 G4cout <<
" InuclCollider -> can not generate acceptable inter. after "
265 << itry_max <<
" attempts " <<
G4endl;
286 globalOutput.
reset();
289 theIntraNucleiCascader->
rescatter(bullet, theSecondaries, theNucleus,
297 globalOutput.
add(output);
300 G4cout <<
" InuclCollider output after trials " << itry <<
G4endl;
308 if (fragment.
GetA() <= 1)
return;
312 const G4int itry_max = 10;
318 theDeexcitation->
deExcite(fragment, DEXoutput);
319 }
while (!
validateOutput(fragment, DEXoutput) && (++itry < itry_max));
322 globalOutput.
add(DEXoutput);
G4bool hadNucleus() const
void trivialise(G4InuclParticle *bullet, G4InuclParticle *target)
void rescatter(G4InuclParticle *bullet, G4KineticTrackVector *theSecondaries, G4V3DNucleus *theNucleus, G4CollisionOutput &globalOutput)
void setVerboseLevel(G4int verbose)
static const G4CascadeChannel * GetTable(G4int initialState)
G4double getTRSMomentum() const
virtual G4bool useEPCollider(G4InuclParticle *bullet, G4InuclParticle *target) const
void removeRecoilFragment(G4int index=-1)
void collide(G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)
virtual G4bool inelasticInteractionPossible(G4InuclParticle *bullet, G4InuclParticle *target, G4double ekin) const
virtual ~G4InuclCollider()
void deexcite(const G4Fragment &fragment, G4CollisionOutput &globalOutput)
const G4String & GetParticleName() const
G4bool acceptable() const
virtual void setVerboseLevel(G4int verbose=0)
virtual void deExcite(const G4Fragment &fragment, G4CollisionOutput &output)=0
void add(const G4CollisionOutput &right)
void useCascadeDeexcitation()
G4GLOB_DLL std::ostream G4cout
void setVerboseLevel(G4int verbose=0)
void rescatter(G4InuclParticle *bullet, G4KineticTrackVector *theSecondaries, G4V3DNucleus *theNucleus, G4CollisionOutput &globalOutput)
void boostToLabFrame(const G4LorentzConvertor &convertor)
void collide(G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &globalOutput)
G4double getKinEnergyInTheTRS() const
void collide(G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &globalOutput)
G4InteractionCase interCase
void set(G4InuclParticle *part1, G4InuclParticle *part2)
virtual G4bool validateOutput(G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)
G4InuclParticle * getBullet() const
virtual void setVerboseLevel(G4int verbose=0)
void setVerboseLevel(G4int verbose=0)
const G4Fragment & getRecoilFragment(G4int index=0) const
void toTheTargetRestFrame()
G4InuclParticle * getTarget() const
G4ParticleDefinition * getDefinition() const
void usePreCompoundDeexcitation()
void setOnShell(G4InuclParticle *bullet, G4InuclParticle *target)
G4GLOB_DLL std::ostream G4cerr