G4LENDInelastic Class Reference

#include <G4LENDInelastic.hh>

Inheritance diagram for G4LENDInelastic:

Public Member Functions
	G4LENDInelastic (G4ParticleDefinition *pd)
	~G4LENDInelastic ()
G4HadFinalState *	ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
Public Member Functions inherited from G4LENDModel
	G4LENDModel (G4String name="LENDModel")
	~G4LENDModel ()
virtual G4HadFinalState *	ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
void	ChangeDefaultEvaluation (G4String name)
void	AllowNaturalAbundanceTarget ()
void	AllowAnyCandidateTarget ()
void	BuildPhysicsTable (const G4ParticleDefinition &)
Public Member Functions inherited from G4HadronicInteraction
	G4HadronicInteraction (const G4String &modelName="HadronicModel")
virtual	~G4HadronicInteraction ()
virtual G4HadFinalState *	ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)=0
virtual G4double	SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
virtual G4bool	IsApplicable (const G4HadProjectile &, G4Nucleus &)
G4double	GetMinEnergy () const
G4double	GetMinEnergy (const G4Material *aMaterial, const G4Element *anElement) const
void	SetMinEnergy (G4double anEnergy)
void	SetMinEnergy (G4double anEnergy, const G4Element *anElement)
void	SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)
G4double	GetMaxEnergy () const
G4double	GetMaxEnergy (const G4Material *aMaterial, const G4Element *anElement) const
void	SetMaxEnergy (const G4double anEnergy)
void	SetMaxEnergy (G4double anEnergy, const G4Element *anElement)
void	SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)
const G4HadronicInteraction *	GetMyPointer () const
G4int	GetVerboseLevel () const
void	SetVerboseLevel (G4int value)
const G4String &	GetModelName () const
void	DeActivateFor (const G4Material *aMaterial)
void	ActivateFor (const G4Material *aMaterial)
void	DeActivateFor (const G4Element *anElement)
void	ActivateFor (const G4Element *anElement)
G4bool	IsBlocked (const G4Material *aMaterial) const
G4bool	IsBlocked (const G4Element *anElement) const
void	SetRecoilEnergyThreshold (G4double val)
G4double	GetRecoilEnergyThreshold () const
G4bool	operator== (const G4HadronicInteraction &right) const
G4bool	operator!= (const G4HadronicInteraction &right) const
virtual const std::pair< G4double, G4double >	GetFatalEnergyCheckLevels () const
virtual std::pair< G4double, G4double >	GetEnergyMomentumCheckLevels () const
void	SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
virtual void	ModelDescription (std::ostream &outFile) const

Additional Inherited Members
Protected Member Functions inherited from G4LENDModel
void	create_used_target_map ()
void	recreate_used_target_map ()
Protected Member Functions inherited from G4HadronicInteraction
void	SetModelName (const G4String &nam)
G4bool	IsBlocked () const
void	Block ()
Protected Attributes inherited from G4LENDModel
G4ParticleDefinition *	proj
G4LENDManager *	lend_manager
std::map< G4int, G4LENDUsedTarget * >	usedTarget_map
Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState	theParticleChange
G4int	verboseLevel
G4double	theMinEnergy
G4double	theMaxEnergy
G4bool	isBlocked

Detailed Description

Definition at line 44 of file G4LENDInelastic.hh.

Constructor & Destructor Documentation

G4LENDInelastic()

G4LENDInelastic::G4LENDInelastic ( G4ParticleDefinition *  pd )

inline

Definition at line 49 of file G4LENDInelastic.hh.

     :G4LENDModel( "LENDInelastic" ) 
     { 
        proj = pd; 
       
//        theModelName = "LENDInelastic for "; 
//        theModelName += proj->GetParticleName(); 
        create_used_target_map();
     };

~G4LENDInelastic()

G4LENDInelastic::~G4LENDInelastic ( )

inline

Definition at line 59 of file G4LENDInelastic.hh.

{;};

Member Function Documentation

ApplyYourself()

G4HadFinalState * G4LENDInelastic::ApplyYourself ( const G4HadProjectile &  aTrack, G4Nucleus &  aTargetNucleus  )

virtual

Reimplemented from G4LENDModel.

Definition at line 31 of file G4LENDInelastic.cc.

{
 
   G4ThreeVector proj_p = aTrack.Get4Momentum().vect();
 
   G4double temp = aTrack.GetMaterial()->GetTemperature();
 
   //G4int iZ = int ( aTarg.GetZ() );
   //G4int iA = int ( aTarg.GetN() );
   //migrate to integer A and Z (GetN_asInt returns number of neutrons in the nucleus since this) 
   G4int iZ = aTarg.GetZ_asInt();
   G4int iA = aTarg.GetA_asInt();
   //G4cout << "target: Z = " << iZ << " N = " << iA << G4endl;
 
   G4double ke = aTrack.GetKineticEnergy();
   //G4cout << "projectile: KE = " << ke/MeV << " [MeV]" << G4endl;
 
   G4HadFinalState* theResult = &theParticleChange;
   theResult->Clear();
 
   G4GIDI_target* aTarget = usedTarget_map.find( lend_manager->GetNucleusEncoding( iZ , iA ) )->second->GetTarget();
   std::vector<G4GIDI_Product>* products = aTarget->getOthersFinalState( ke*MeV, temp, NULL, NULL );
   if ( products != NULL ) 
   {
 
      G4ThreeVector psum(0);
 
      int totN = 0;
      for ( G4int j = 0; j < int( products->size() ); j++ ) 
      {
 
         G4int jZ = (*products)[j].Z; 
         G4int jA = (*products)[j].A; 
 
         //G4cout << "ZA = " << 1000 * (*products)[j].Z + (*products)[j].A << "  EK = "
         //     << (*products)[j].kineticEnergy
         //     << " px  " <<  (*products)[j].px
         //     << " py  " <<  (*products)[j].py
         //     << " pz  " <<  (*products)[j].pz
         //     << G4endl;
 
         G4DynamicParticle* theSec = new G4DynamicParticle;
 
         if ( jA == 1 && jZ == 1 )
         {
            theSec->SetDefinition( G4Proton::Proton() );
            totN += 1;
         }
         else if ( jA == 1 && jZ == 0 )
         {
            theSec->SetDefinition( G4Neutron::Neutron() );
            totN += 1;
         } 
         else if ( jZ > 0 )
         {
            if ( jA != 0 )
            {
               theSec->SetDefinition( G4ParticleTable::GetParticleTable()->FindIon( jZ , jA , 0 , 0 ) );
               totN += jA;
            }
            else 
            {
               theSec->SetDefinition( G4ParticleTable::GetParticleTable()->FindIon( jZ , iA+1-totN , 0 , 0 ) );
            }
         } 
         else
         {
            theSec->SetDefinition( G4Gamma::Gamma() );
         } 
 
         G4ThreeVector p( (*products)[j].px*MeV , (*products)[j].py*MeV , (*products)[j].pz*MeV ); 
         psum += p; 
         if ( p.mag() == 0 ) p = proj_p - psum;
 
         theSec->SetMomentum( p );
 
         theResult->AddSecondary( theSec );
      } 
   }
   delete products;
 
   theResult->SetStatusChange( stopAndKill );
 
   return theResult; 
 
}

---

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

• G4LENDInelastic.hh
• G4LENDInelastic.cc