G4InuclCollider Class Reference

#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 ()
virtual void	rescatter (G4InuclParticle *, G4KineticTrackVector *, G4V3DNucleus *, G4CollisionOutput &)
virtual void	setVerboseLevel (G4int verbose=0)
virtual void	setConservationChecks (G4bool doBalance=true)
Public Member Functions inherited from G4VCascadeCollider
	G4VCascadeCollider (const char *name, G4int verbose=0)
virtual	~G4VCascadeCollider ()
virtual void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)=0
virtual void	setVerboseLevel (G4int verbose=0)

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	explosion (G4InuclNuclei *target) const
virtual G4bool	explosion (G4Fragment *target) const
virtual G4bool	explosion (G4int A, G4int Z, G4double excitation) const
virtual G4bool	inelasticInteractionPossible (G4InuclParticle *bullet, G4InuclParticle *target, G4double ekin) const
virtual G4bool	validateOutput (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)
virtual G4bool	validateOutput (G4InuclParticle *bullet, G4InuclParticle *target, const std::vector< G4InuclElementaryParticle > &particles)
virtual G4bool	validateOutput (G4InuclParticle *bullet, G4InuclParticle *target, const std::vector< G4InuclNuclei > &fragments)
Protected Member Functions inherited from G4VCascadeCollider
virtual void	setName (const char *name)

Additional Inherited Members
Protected Attributes inherited from G4CascadeColliderBase
G4InteractionCase	interCase
G4bool	doConservationChecks
G4CascadeCheckBalance *	balance
Protected Attributes inherited from G4VCascadeCollider
const char *	theName
G4int	verboseLevel

Detailed Description

Definition at line 59 of file G4InuclCollider.hh.

Constructor & Destructor Documentation

G4InuclCollider()

G4InuclCollider::G4InuclCollider

(

)

Definition at line 78 of file G4InuclCollider.cc.

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

~G4InuclCollider()

G4InuclCollider::~G4InuclCollider ( )

virtual

Definition at line 84 of file G4InuclCollider.cc.

                                  {
  delete theElementaryParticleCollider;
  delete theIntraNucleiCascader;
  delete theDeexcitation;
}

Member Function Documentation

collide()

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

virtual

Implements G4VCascadeCollider.

Definition at line 121 of file G4InuclCollider.cc.

                                                   {
  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;
}

Referenced by G4CascadeInterface::ApplyYourself().

deexcite()

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

protected

Definition at line 305 of file G4InuclCollider.cc.

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

Referenced by collide(), and rescatter().

rescatter()

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

virtual

Reimplemented from G4CascadeColliderBase.

Definition at line 276 of file G4InuclCollider.cc.

                                                 {
  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;
}

Referenced by G4CascadeInterface::Propagate().

setVerboseLevel()

void G4InuclCollider::setVerboseLevel ( G4int  verbose = 0 )

virtual

Reimplemented from G4CascadeColliderBase.

Definition at line 92 of file G4InuclCollider.cc.

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

Referenced by G4CascadeInterface::SetVerboseLevel().

useCascadeDeexcitation()

void G4InuclCollider::useCascadeDeexcitation

(

)

Definition at line 106 of file G4InuclCollider.cc.

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

Referenced by G4CascadeInterface::useCascadeDeexcitation().

usePreCompoundDeexcitation()

void G4InuclCollider::usePreCompoundDeexcitation

(

)

Definition at line 112 of file G4InuclCollider.cc.

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

Referenced by G4CascadeInterface::usePreCompoundDeexcitation().

---

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

• G4InuclCollider.hh
• G4InuclCollider.cc