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 ()
void	useAblaDeexcitation ()
Public Member Functions inherited from G4CascadeColliderBase
	G4CascadeColliderBase (const G4String &name, G4int verbose=0)
virtual	~G4CascadeColliderBase ()
Public Member Functions inherited from G4VCascadeCollider
	G4VCascadeCollider (const G4String &name, G4int verbose=0)
virtual	~G4VCascadeCollider ()

Protected Member Functions
void	deexcite (const G4Fragment &fragment, G4CollisionOutput &globalOutput)
G4bool	photonuclearOkay (G4CollisionOutput &checkOutput) const
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 G4String &name)

Additional Inherited Members
Protected Attributes inherited from G4CascadeColliderBase
G4InteractionCase	interCase
G4CascadeCheckBalance *	balance
Protected Attributes inherited from G4VCascadeCollider
G4String	theName
G4int	verboseLevel

Detailed Description

Definition at line 60 of file G4InuclCollider.hh.

Constructor & Destructor Documentation

G4InuclCollider()

G4InuclCollider::G4InuclCollider

(

)

Definition at line 88 of file G4InuclCollider.cc.

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

~G4InuclCollider()

G4InuclCollider::~G4InuclCollider ( )

virtual

Definition at line 94 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 137 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) { /* Loop checking 08.06.2015 MHK */
    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();
 
    //*** TEMPORARY, USE ENVVAR TO ENABLE/DISABLE THIS TEST ***
    if (std::getenv("G4CASCADE_CHECK_PHOTONUCLEAR"))
      if (!photonuclearOkay(output)) continue;
 
    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;
}

deexcite()

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

protected

Definition at line 323 of file G4InuclCollider.cc.

                                                    {
  if (fragment.GetA_asInt() <= 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 {                  /* Loop checking 08.06.2015 MHK */
    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);
}

Referenced by collide(), and rescatter().

photonuclearOkay()

G4bool G4InuclCollider::photonuclearOkay ( G4CollisionOutput & checkOutput ) const

protected

Definition at line 345 of file G4InuclCollider.cc.

                                                                             {
  if (interCase.twoNuclei()) return true;  // A-A is not photonuclear
 
  G4InuclElementaryParticle* bullet =
    dynamic_cast<G4InuclElementaryParticle*>(interCase.getBullet());
  if (!bullet || !(bullet->isPhoton() || bullet->isElectron())) return true;
 
  if (verboseLevel>1)
    G4cout << " >>> G4InuclCollider::photonuclearOkay" << G4endl;
 
  if (bullet->getKineticEnergy() > 0.050) return true;
 
  if (verboseLevel>2) {
    if (checkOutput.numberOfOutgoingNuclei() > 0) {
      G4cout << " comparing final nucleus with initial target:\n"
             << checkOutput.getOutgoingNuclei()[0] << G4endl
             << *(interCase.getTarget()) << G4endl;
    } else {
      G4cout << " no final nucleus remains when target was "
             << *(interCase.getTarget()) << G4endl;
    }
  }
 
  // Hadron production changes target nucleus
  G4double mfinalNuc = 0.0;
  if (checkOutput.numberOfOutgoingNuclei() > 0)
    mfinalNuc = checkOutput.getOutgoingNuclei()[0].getMass();
  G4double mtargetNuc = interCase.getTarget()->getMass();
  if (mfinalNuc != mtargetNuc) return true; // Mass from G4Ions is fixed
  
  if (verboseLevel>2)
    G4cout << " photonuclear produced only gammas.  Try again." << G4endl;
 
  return false;     // Final state is entirely de-excitation photons
}

Referenced by collide().

rescatter()

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

virtual

Reimplemented from G4CascadeColliderBase.

Definition at line 295 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;
}

setVerboseLevel()

void G4InuclCollider::setVerboseLevel ( G4int verbose = 0 )

virtual

Reimplemented from G4CascadeColliderBase.

Definition at line 102 of file G4InuclCollider.cc.

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

useAblaDeexcitation()

void G4InuclCollider::useAblaDeexcitation

(

)

Definition at line 128 of file G4InuclCollider.cc.

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

useCascadeDeexcitation()

void G4InuclCollider::useCascadeDeexcitation

(

)

Definition at line 116 of file G4InuclCollider.cc.

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

usePreCompoundDeexcitation()

void G4InuclCollider::usePreCompoundDeexcitation

(

)

Definition at line 122 of file G4InuclCollider.cc.

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

---

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

• G4InuclCollider.hh
• G4InuclCollider.cc