#include <G4LENDFission.hh>

Inheritance diagram for G4LENDFission:

Public Member Functions
	G4LENDFission (G4ParticleDefinition *pd)

	~G4LENDFission ()

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

Public Member Functions inherited from G4LENDModel
	G4LENDModel (G4String name="LENDModel")

	~G4LENDModel ()

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 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

virtual G4int	GetVerboseLevel () const

virtual 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 46 of file G4LENDFission.hh.

Constructor & Destructor Documentation

G4LENDFission::G4LENDFission ( G4ParticleDefinition * pd )

inline

Definition at line 51 of file G4LENDFission.hh.

References G4LENDModel::create_used_target_map(), and G4LENDModel::proj.

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

G4LENDFission::~G4LENDFission ( )

inline

Definition at line 61 of file G4LENDFission.hh.

61 {;};

Member Function Documentation

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

virtual

Reimplemented from G4LENDModel.

Definition at line 31 of file G4LENDFission.cc.

References G4HadFinalState::AddSecondary(), G4HadFinalState::Clear(), G4Gamma::Gamma(), G4Nucleus::GetA_asInt(), G4GIDI_target::getFissionFinalState(), G4Nucleus::GetIsotope(), G4HadProjectile::GetKineticEnergy(), G4Isotope::Getm(), G4HadProjectile::GetMaterial(), G4LENDManager::GetNucleusEncoding(), G4ParticleTable::GetParticleTable(), G4Material::GetTemperature(), G4Nucleus::GetZ_asInt(), int(), G4LENDModel::lend_manager, python.hepunit::MeV, G4Neutron::Neutron(), python.hepunit::second, G4DynamicParticle::SetDefinition(), G4DynamicParticle::SetMomentum(), G4HadFinalState::SetStatusChange(), stopAndKill, G4HadronicInteraction::theParticleChange, and G4LENDModel::usedTarget_map.

 {
 
    G4double temp = aTrack.GetMaterial()->GetTemperature();
 
    //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();
    //G4int iM = aTarg.GetM_asInt();
    G4int iM = 0;
    if ( aTarg.GetIsotope() != NULL ) {
       iM = aTarg.GetIsotope()->Getm();
    }
 
    G4double ke = aTrack.GetKineticEnergy();
 
    G4HadFinalState* theResult = &theParticleChange;
    theResult->Clear();
 
    G4GIDI_target* aTarget = usedTarget_map.find( lend_manager->GetNucleusEncoding( iZ , iA , iM ) )->second->GetTarget();
    std::vector<G4GIDI_Product>* products = aTarget->getFissionFinalState( ke*MeV, temp, NULL, NULL );
    if ( products != NULL ) 
    {
       for ( G4int j = 0; j < int( products->size() ); j++ ) 
       {
          G4int jZ = (*products)[j].Z; 
          G4int jA = (*products)[j].A; 
 
          //G4cout << "Z = "    << (*products)[j].Z 
          //       << ", A = "  << (*products)[j].A 
          //       << ", EK = " << (*products)[j].kineticEnergy << " [MeV]" 
          //       << ", px = " << (*products)[j].px
          //       << ", py = " << (*products)[j].py
          //       << ", pz = " << (*products)[j].pz
          //       << G4endl;
 
          G4DynamicParticle* theSec = new G4DynamicParticle;
 
          if ( jZ > 0 )
          {
                                                                                         //Ex  j?
             theSec->SetDefinition( G4ParticleTable::GetParticleTable()->FindIon( jZ, jA , 0, 0 ) );
          } 
          else if ( jA == 1 && jZ == 0 )
          {
             theSec->SetDefinition( G4Neutron::Neutron() );
          } 
          else
          {
             theSec->SetDefinition( G4Gamma::Gamma() );
          } 
 
          theSec->SetMomentum( G4ThreeVector( (*products)[j].px*MeV , (*products)[j].py*MeV , (*products)[j].pz*MeV ) );
          //G4cout << theSec->GetDefinition()->GetParticleName() << G4endl;
          theResult->AddSecondary( theSec );
       } 
    }
    delete products;
 
    theResult->SetStatusChange( stopAndKill );
 
    return theResult; 
 
 }

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation