#include <G4InuclCollider.hh>

Inheritance diagram for G4InuclCollider:

Public Member Functions
	G4InuclCollider ()

virtual	~G4InuclCollider ()

void	collide (G4InuclParticle bullet, G4InuclParticle target, G4CollisionOutput &globalOutput)

void	rescatter (G4InuclParticle bullet, G4KineticTrackVector theSecondaries, G4V3DNucleus *theNucleus, G4CollisionOutput &globalOutput)

void	setVerboseLevel (G4int verbose=0)

void	useCascadeDeexcitation ()

void	usePreCompoundDeexcitation ()

Public Member Functions inherited from G4CascadeColliderBase
	G4CascadeColliderBase (const char *name, G4int verbose=0)

virtual	~G4CascadeColliderBase ()

Public Member Functions inherited from G4VCascadeCollider
	G4VCascadeCollider (const char *name, G4int verbose=0)

virtual	~G4VCascadeCollider ()

Protected Member Functions
void	deexcite (const G4Fragment &fragment, G4CollisionOutput &globalOutput)

Protected Member Functions inherited from G4CascadeColliderBase
virtual G4bool	useEPCollider (G4InuclParticle bullet, G4InuclParticle target) const

virtual G4bool	inelasticInteractionPossible (G4InuclParticle bullet, G4InuclParticle target, G4double ekin) const

virtual G4bool	validateOutput (G4InuclParticle bullet, G4InuclParticle target, G4CollisionOutput &output)

virtual G4bool	validateOutput (const G4Fragment &fragment, G4CollisionOutput &output)

virtual G4bool	validateOutput (G4InuclParticle bullet, G4InuclParticle target, const std::vector< G4InuclElementaryParticle > &particles)

Protected Member Functions inherited from G4VCascadeCollider
virtual void	setName (const char *name)

Additional Inherited Members
Protected Attributes inherited from G4CascadeColliderBase
G4InteractionCase	interCase

G4CascadeCheckBalance *	balance

Protected Attributes inherited from G4VCascadeCollider
const char *	theName

G4int	verboseLevel

Detailed Description

Definition at line 60 of file G4InuclCollider.hh.

Constructor & Destructor Documentation

G4InuclCollider::G4InuclCollider ( )

Definition at line 79 of file G4InuclCollider.cc.

   : G4CascadeColliderBase("G4InuclCollider"),
     theElementaryParticleCollider(new G4ElementaryParticleCollider),
     theIntraNucleiCascader(new G4IntraNucleiCascader),
     theDeexcitation(new G4CascadeDeexcitation) {}

G4InuclCollider::~G4InuclCollider ( )

virtual

Definition at line 85 of file G4InuclCollider.cc.

                                   {
   delete theElementaryParticleCollider;
   delete theIntraNucleiCascader;
   delete theDeexcitation;
 }

Member Function Documentation

void G4InuclCollider::collide	(	G4InuclParticle *	bullet,
		G4InuclParticle *	target,
		G4CollisionOutput &	globalOutput
	)

virtual

Implements G4VCascadeCollider.

Definition at line 122 of file G4InuclCollider.cc.

References G4CollisionOutput::acceptable(), G4CollisionOutput::add(), G4CollisionOutput::boostToLabFrame(), G4InteractionCase::code(), G4ElementaryParticleCollider::collide(), G4IntraNucleiCascader::collide(), deexcite(), G4cerr, G4cout, G4endl, G4InuclNuclei::getA(), G4InteractionCase::getBullet(), G4InuclParticle::getDefinition(), G4LorentzConvertor::getKinEnergyInTheTRS(), G4ParticleDefinition::GetParticleName(), G4CollisionOutput::getRecoilFragment(), G4CascadeChannelTables::GetTable(), G4InteractionCase::getTarget(), G4LorentzConvertor::getTRSMomentum(), G4InuclNuclei::getZ(), G4InteractionCase::hadNucleus(), G4CascadeColliderBase::inelasticInteractionPossible(), G4CascadeColliderBase::interCase, G4CollisionOutput::removeRecoilFragment(), G4CollisionOutput::reset(), G4InteractionCase::set(), G4CollisionOutput::setOnShell(), CLHEP::HepLorentzVector::setZ(), G4LorentzConvertor::toTheTargetRestFrame(), G4LorentzConvertor::trivial(), G4CollisionOutput::trivialise(), G4InuclElementaryParticle::type(), G4CascadeColliderBase::useEPCollider(), G4InteractionCase::valid(), and G4VCascadeCollider::verboseLevel.

Referenced by G4CascadeInterface::ApplyYourself().

                                                    {
   if (verboseLevel) G4cout << " >>> G4InuclCollider::collide" << G4endl;
 
   const G4int itry_max = 100;
 
   // Particle-on-particle collision; no nucleus involved
   if (useEPCollider(bullet,target)) {
     if (verboseLevel > 2)
       G4cout << " InuclCollider -> particle on particle collision" << G4endl;
  
     theElementaryParticleCollider->collide(bullet, target, globalOutput);
     return;
   }
   
   interCase.set(bullet,target);     // Classify collision type
   if (verboseLevel > 2) {
     G4cout << " InuclCollider -> inter case " << interCase.code() << G4endl;
   }
 
   if (!interCase.valid()) {
     if (verboseLevel > 1)
       G4cerr << " InuclCollider -> no collision possible " << G4endl;
 
     globalOutput.trivialise(bullet, target);
     return;
   }
 
   // Target must be a nucleus
   G4InuclNuclei* ntarget = dynamic_cast<G4InuclNuclei*>(interCase.getTarget());
   if (!ntarget) {
     G4cerr << " InuclCollider -> ERROR target is not a nucleus " << G4endl;
 
     globalOutput.trivialise(bullet, target);
     return;
   }
 
   G4int btype = 0;
   G4int ab = 0;
   G4int zb = 0;
   
   if (interCase.hadNucleus()) {     // particle with nuclei
     G4InuclElementaryParticle* pbullet = 
       dynamic_cast<G4InuclElementaryParticle*>(interCase.getBullet());
 
     if (!pbullet) {
       G4cerr << " InuclCollider -> ERROR bullet is not a hadron " << G4endl;
       globalOutput.trivialise(bullet, target);
       return;
     }
 
     if (!G4CascadeChannelTables::GetTable(pbullet->type())) {
       G4cerr << " InuclCollider -> ERROR can not collide with "
          << pbullet->getDefinition()->GetParticleName() << G4endl;
       globalOutput.trivialise(bullet, target);
       return;
     }
 
     btype = pbullet->type();
   } else {              // nuclei with nuclei
     G4InuclNuclei* nbullet = 
       dynamic_cast<G4InuclNuclei*>(interCase.getBullet());
     if (!nbullet) {
       G4cerr << " InuclCollider -> ERROR bullet is not a nucleus " << G4endl;
       globalOutput.trivialise(bullet, target);
       return;
     }
     
     ab = nbullet->getA();
     zb = nbullet->getZ();
   }
 
   G4LorentzConvertor convertToTargetRestFrame(bullet, ntarget);
   G4double ekin = convertToTargetRestFrame.getKinEnergyInTheTRS();
   
   if (verboseLevel > 3) G4cout << " ekin in trs " << ekin << G4endl;
 
   if (!inelasticInteractionPossible(bullet, target, ekin)) {
     if (verboseLevel > 3)
       G4cout << " InuclCollider -> inelastic interaction is impossible\n"
          << " due to the coulomb barirer " << G4endl;
 
     globalOutput.trivialise(bullet, target);
     return;
   }
 
   // Generate interaction secondaries in rest frame of target nucleus
   convertToTargetRestFrame.toTheTargetRestFrame();
   if (verboseLevel > 3) {
     G4cout << " degenerated? " << convertToTargetRestFrame.trivial()
        << G4endl;
   }
   
   G4LorentzVector bmom;         // Bullet is along local Z
   bmom.setZ(convertToTargetRestFrame.getTRSMomentum());
 
   // Need to make copy of bullet with momentum realigned
   G4InuclParticle* zbullet = 0;
   if (interCase.hadNucleus())
     zbullet = new G4InuclElementaryParticle(bmom, btype);
   else
     zbullet = new G4InuclNuclei(bmom, ab, zb);
 
   G4int itry = 0;
   while (itry < itry_max) {
     itry++;
     if (verboseLevel > 2)
       G4cout << " InuclCollider itry " << itry << G4endl;
 
     globalOutput.reset();       // Clear buffers for this attempt
     output.reset();
 
     theIntraNucleiCascader->collide(zbullet, target, output);
     
     if (verboseLevel > 1) G4cout << " After Cascade " << G4endl;
 
     deexcite(output.getRecoilFragment(), output);
     output.removeRecoilFragment();
 
     if (verboseLevel > 2)
       G4cout << " itry " << itry << " finished, moving to lab frame" << G4endl;
 
     // convert to the LAB frame and add to final result
     output.boostToLabFrame(convertToTargetRestFrame);
 
     globalOutput.add(output);
 
     // Adjust final state particles to balance momentum and energy
     // FIXME:  This should no longer be necessary!
     globalOutput.setOnShell(bullet, target);
     if (globalOutput.acceptable()) {
       if (verboseLevel) 
     G4cout << " InuclCollider output after trials " << itry << G4endl;
       delete zbullet;
       return;
     } else {
       if (verboseLevel>2)
     G4cerr << " InuclCollider setOnShell failed." << G4endl;
     }
   } // while (itry < itry_max)
   
   if (verboseLevel) {
     G4cout << " InuclCollider -> can not generate acceptable inter. after " 
        << itry_max << " attempts " << G4endl;
   }
   
   globalOutput.trivialise(bullet, target);
 
   delete zbullet;
   return;
 }

void G4InuclCollider::deexcite	(	const G4Fragment &	fragment,
		G4CollisionOutput &	globalOutput
	)

protected

Definition at line 306 of file G4InuclCollider.cc.

References G4CollisionOutput::add(), G4VCascadeDeexcitation::deExcite(), G4cout, G4endl, G4Fragment::GetA(), G4CollisionOutput::reset(), G4CascadeColliderBase::validateOutput(), and G4VCascadeCollider::verboseLevel.

Referenced by collide(), and rescatter().

                                                     {
   if (fragment.GetA() <= 1) return; // Nothing real to be de-excited
 
   if (verboseLevel) G4cout << " >>> G4InuclCollider::deexcite" << G4endl;
 
   const G4int itry_max = 10;        // Maximum number of attempts
   G4int itry = 0;
   do {
     if (verboseLevel > 2) G4cout << " deexcite itry " << itry << G4endl;
 
     DEXoutput.reset();
     theDeexcitation->deExcite(fragment, DEXoutput);
   } while (!validateOutput(fragment, DEXoutput) && (++itry < itry_max));
 
   // Add de-excitation products to output buffer
   globalOutput.add(DEXoutput);
 }

void G4InuclCollider::rescatter	(	G4InuclParticle *	bullet,
		G4KineticTrackVector *	theSecondaries,
		G4V3DNucleus *	theNucleus,
		G4CollisionOutput &	globalOutput
	)

virtual

Reimplemented from G4CascadeColliderBase.

Definition at line 277 of file G4InuclCollider.cc.

References G4CollisionOutput::add(), deexcite(), G4cout, G4endl, G4CollisionOutput::getRecoilFragment(), G4CollisionOutput::removeRecoilFragment(), G4IntraNucleiCascader::rescatter(), G4CollisionOutput::reset(), and G4VCascadeCollider::verboseLevel.

Referenced by G4CascadeInterface::Propagate().

                                                  {
   if (verboseLevel) G4cout << " >>> G4InuclCollider::rescatter" << G4endl;
 
   G4int itry=1;     // For diagnostic post-processing only
   if (verboseLevel > 2) G4cout << " InuclCollider itry " << itry << G4endl;
 
   globalOutput.reset();     // Clear buffers for this attempt
   output.reset();
 
   theIntraNucleiCascader->rescatter(bullet, theSecondaries, theNucleus, 
                     output);
 
   if (verboseLevel > 1) G4cout << " After Rescatter" << G4endl;
 
   deexcite(output.getRecoilFragment(), output);
   output.removeRecoilFragment();
 
   globalOutput.add(output); // Add local results to global output
 
   if (verboseLevel) 
     G4cout << " InuclCollider output after trials " << itry << G4endl;
 }

void G4InuclCollider::setVerboseLevel ( G4int verbose = 0 )

virtual

Reimplemented from G4CascadeColliderBase.

Definition at line 93 of file G4InuclCollider.cc.

References G4VCascadeCollider::setVerboseLevel(), G4CollisionOutput::setVerboseLevel(), G4CascadeColliderBase::setVerboseLevel(), G4IntraNucleiCascader::setVerboseLevel(), and G4VCascadeCollider::verboseLevel.

Referenced by G4CascadeInterface::SetVerboseLevel().

                                                    {
   G4CascadeColliderBase::setVerboseLevel(verbose);
 
   theElementaryParticleCollider->setVerboseLevel(verboseLevel);
   theIntraNucleiCascader->setVerboseLevel(verboseLevel);
   theDeexcitation->setVerboseLevel(verboseLevel);
 
   output.setVerboseLevel(verboseLevel);
   DEXoutput.setVerboseLevel(verboseLevel);
 }

void G4InuclCollider::useCascadeDeexcitation ( )

Definition at line 107 of file G4InuclCollider.cc.

References G4VCascadeCollider::setVerboseLevel(), and G4VCascadeCollider::verboseLevel.

Referenced by G4CascadeInterface::useCascadeDeexcitation().

                                              {
   delete theDeexcitation;
   theDeexcitation = new G4CascadeDeexcitation;
   theDeexcitation->setVerboseLevel(verboseLevel);
 }

void G4InuclCollider::usePreCompoundDeexcitation ( )

Definition at line 113 of file G4InuclCollider.cc.

References G4VCascadeCollider::setVerboseLevel(), and G4VCascadeCollider::verboseLevel.

Referenced by G4CascadeInterface::usePreCompoundDeexcitation().

                                                  {
   delete theDeexcitation;
   theDeexcitation = new G4PreCompoundDeexcitation;
   theDeexcitation->setVerboseLevel(verboseLevel);
 }

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

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation