G4LEDeuteronInelastic Class Reference

#include <G4LEDeuteronInelastic.hh>

Inheritance diagram for G4LEDeuteronInelastic:

Public Member Functions
	G4LEDeuteronInelastic ()
	~G4LEDeuteronInelastic ()
G4HadFinalState *	ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
virtual void	ModelDescription (std::ostream &outFile) const

---

Detailed Description

Definition at line 45 of file G4LEDeuteronInelastic.hh.

---

Constructor & Destructor Documentation

G4LEDeuteronInelastic::G4LEDeuteronInelastic

(

)

[inline]

Definition at line 49 of file G4LEDeuteronInelastic.hh.

References G4cout, G4endl, G4HadronicInteraction::SetMaxEnergy(), and G4HadronicInteraction::SetMinEnergy().

                            : G4InelasticInteraction("G4LEDeuteronInelastic")
    {
      SetMinEnergy( 0.0 );
      // SetMaxEnergy( 100.*CLHEP::MeV );  // NUCREC only worked for energies < 100MeV
      // Work around to avoid exception in G4EnergyRangeManager
      SetMaxEnergy( 10.*CLHEP::TeV );  // NUCREC only worked for energies < 100MeV
      G4cout << "WARNING: model G4LEDeuteronInelastic is being deprecated and will\n"
             << "disappear in Geant4 version 10.0"  << G4endl;
    }

G4LEDeuteronInelastic::~G4LEDeuteronInelastic

(

)

[inline]

Definition at line 59 of file G4LEDeuteronInelastic.hh.

{}

---

Member Function Documentation

G4HadFinalState * G4LEDeuteronInelastic::ApplyYourself	(	const G4HadProjectile &	aTrack,
		G4Nucleus &	targetNucleus
	)			[virtual]

Implements G4HadronicInteraction.

Definition at line 51 of file G4LEDeuteronInelastic.cc.

References G4HadFinalState::AddSecondary(), G4Nucleus::AtomicMass(), G4HadFinalState::Clear(), G4InelasticInteraction::DoIsotopeCounting(), G4cout, G4endl, G4HadProjectile::Get4Momentum(), G4Nucleus::GetA_asInt(), G4HadProjectile::GetKineticEnergy(), G4HadProjectile::GetMaterial(), G4Material::GetName(), G4Nucleus::GetZ_asInt(), G4FastVector< Type, N >::Initialize(), isAlive, G4InelasticInteraction::isotopeProduction, G4ReactionDynamics::NuclearReaction(), G4DynamicParticle::SetDefinition(), G4HadFinalState::SetEnergyChange(), G4DynamicParticle::SetMomentum(), G4HadFinalState::SetMomentumChange(), G4HadFinalState::SetStatusChange(), G4HadronicInteraction::theParticleChange, G4InelasticInteraction::theReactionDynamics, and G4HadronicInteraction::verboseLevel.

{ 
  theParticleChange.Clear();
  const G4HadProjectile* originalIncident = &aTrack;
    
  if (verboseLevel > 1) {
    const G4Material *targetMaterial = aTrack.GetMaterial();
    G4cout << "G4LEDeuteronInelastic::ApplyYourself called" << G4endl;
    G4cout << "kinetic energy = " << originalIncident->GetKineticEnergy()/MeV << "MeV, ";
    G4cout << "target material = " << targetMaterial->GetName() << ", ";
  }
    
  // Work-around for lack of model above 100 MeV
  if (originalIncident->GetKineticEnergy()/MeV > 100. ||
      originalIncident->GetKineticEnergy() <= 0.1*MeV) {
    theParticleChange.SetStatusChange(isAlive);
    theParticleChange.SetEnergyChange(aTrack.GetKineticEnergy());
    theParticleChange.SetMomentumChange(aTrack.Get4Momentum().vect().unit()); 
    return &theParticleChange;      
  }

  G4double A = targetNucleus.GetA_asInt();
  G4double Z = targetNucleus.GetZ_asInt();
  G4double theAtomicMass = targetNucleus.AtomicMass(A, Z);
  G4double massVec[9];
  massVec[0] = targetNucleus.AtomicMass( A+2.0, Z+1.0 );
  massVec[1] = targetNucleus.AtomicMass( A+1.0, Z+1.0 );
  massVec[2] = targetNucleus.AtomicMass( A+1.0, Z     );
  massVec[3] = theAtomicMass;
  massVec[4] = 0.;
  if (A > 1.0 && A-1.0 > Z) 
    massVec[4] = targetNucleus.AtomicMass(A-1.0, Z);
  massVec[5] = 0.;
  if (A > 2.0 && Z > 1.0 && A-2.0 > Z-1.0) 
    massVec[5] = targetNucleus.AtomicMass(A-2.0, Z-1.0);
  massVec[6] = 0.;
  if (A > Z+1.0) 
    massVec[6] = targetNucleus.AtomicMass(A, Z+1.0);
  massVec[7] = massVec[3];
  massVec[8] = 0.;
  if (Z > 1.0) massVec[8] = targetNucleus.AtomicMass(A,Z-1.0);
    
  G4FastVector<G4ReactionProduct,4> vec;  // vec will contain the secondary particles
  G4int vecLen = 0;
  vec.Initialize( 0 );
    
  theReactionDynamics.NuclearReaction(vec, vecLen, originalIncident,
                                      targetNucleus, theAtomicMass, massVec);

  G4double p = vec[0]->GetMomentum().mag();
  theParticleChange.SetMomentumChange( vec[0]->GetMomentum() * (1.0/p)  );
  theParticleChange.SetEnergyChange( vec[0]->GetKineticEnergy() );
  delete vec[0];

  if (vecLen <= 1) 
    {
      theParticleChange.SetStatusChange(isAlive);
      theParticleChange.SetEnergyChange(aTrack.GetKineticEnergy());
      theParticleChange.SetMomentumChange(aTrack.Get4Momentum().vect().unit());
      if (isotopeProduction) DoIsotopeCounting(originalIncident, targetNucleus);
      return &theParticleChange;      
    }

  G4DynamicParticle* pd;
  for (G4int i=1; i<vecLen; ++i) {
    pd = new G4DynamicParticle();
    pd->SetDefinition( vec[i]->GetDefinition() );
    pd->SetMomentum( vec[i]->GetMomentum() );
    theParticleChange.AddSecondary( pd );
    delete vec[i];
  }

  if (isotopeProduction) DoIsotopeCounting(originalIncident, targetNucleus); 
  return &theParticleChange;
}

void G4LEDeuteronInelastic::ModelDescription

(

std::ostream &

outFile

)

const [virtual]

Reimplemented from G4HadronicInteraction.

Definition at line 36 of file G4LEDeuteronInelastic.cc.

{
  outFile << "G4LEDeuteronInelastic is one of the Low Energy Parameterized\n"
          << "(LEP) models used to implement inelastic deuteron scattering\n"
          << "from nuclei.  It is a re-engineered version of the GHEISHA\n"
          << "code of H. Fesefeldt.  It divides the initial collision\n"
          << "products into backward- and forward-going clusters which are\n"
          << "then decayed into final state hadrons.  The model does not\n"
          << "conserve energy on an event-by-event basis.  It may be\n"
          << "applied to deuterons with initial energies between 0 and 10\n"
          << "TeV.\n";
}

---

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

• G4LEDeuteronInelastic.hh
• G4LEDeuteronInelastic.cc

---