G4NeutronIsotopeProduction Class Reference

#include <G4NeutronIsotopeProduction.hh>

Inheritance diagram for G4NeutronIsotopeProduction:

Public Member Functions
	G4NeutronIsotopeProduction ()
	G4NeutronIsotopeProduction (const G4NeutronIsotopeProduction &nip)
virtual	~G4NeutronIsotopeProduction ()
G4IsoResult *	GetIsotope (const G4HadProjectile *aTrack, const G4Nucleus &aNucleus)
Public Member Functions inherited from G4VIsotopeProduction
	G4VIsotopeProduction ()
virtual	~G4VIsotopeProduction ()
virtual G4IsoResult *	GetIsotope (const G4HadProjectile *aTrack, const G4Nucleus &aNucleus)=0
G4bool	operator== (const G4VIsotopeProduction &aProd)

Detailed Description

Definition at line 44 of file G4NeutronIsotopeProduction.hh.

Constructor & Destructor Documentation

G4NeutronIsotopeProduction() [1/2]

G4NeutronIsotopeProduction::G4NeutronIsotopeProduction

(

)

Definition at line 29 of file G4NeutronIsotopeProduction.cc.

{
  numberOfElements = G4Element::GetNumberOfElements();
  theData = new G4ElementIsoCrossSections<G4NeutronIsoIsoCrossSections>* [numberOfElements];
  for (G4int i = 0; i < numberOfElements; i++) {
    theData[i] = new G4ElementIsoCrossSections<G4NeutronIsoIsoCrossSections>;
    if((*(G4Element::GetElementTable()))[i]->GetZ() > 9 &&
       (*(G4Element::GetElementTable()))[i]->GetZ() < 84) {
      // Workaround to be fixed in G4NeutronHPNames.
      theData[i]->Init((*(G4Element::GetElementTable()))[i]);
    }
  }
  G4cout << "WARNING: G4NeutronIsotopeProduction is deprecated and will be removed with Geant4 version 10"
         << G4endl;
}

G4NeutronIsotopeProduction() [2/2]

G4NeutronIsotopeProduction::G4NeutronIsotopeProduction

(

const G4NeutronIsotopeProduction &

nip

)

~G4NeutronIsotopeProduction()

G4NeutronIsotopeProduction::~G4NeutronIsotopeProduction ( )

virtual

Definition at line 46 of file G4NeutronIsotopeProduction.cc.

{
  for (G4int i = 0; i < numberOfElements; i++) delete theData[i];
  if (theData) delete [] theData;
}

Member Function Documentation

GetIsotope()

G4IsoResult * G4NeutronIsotopeProduction::GetIsotope ( const G4HadProjectile *  aTrack, const G4Nucleus &  aNucleus  )

virtual

Implements G4VIsotopeProduction.

Definition at line 54 of file G4NeutronIsotopeProduction.cc.

{
  G4IsoResult* result = 0;
  // is applicable?
  if (aTrack->GetDefinition() != G4Neutron::Neutron()) return result;
  G4double incidentKE = aTrack->GetKineticEnergy();
  if (incidentKE > 100*MeV) return result;
 
  // get the isotope
  const G4Material* theMaterial = aTrack->GetMaterial();
  G4int nEleInMat = theMaterial->GetNumberOfElements();
  for (G4int check = 0; check < nEleInMat; check++) {
    // Workaround to be fixed in G4NeutronHPNames
    if (theMaterial->GetElement(check)->GetZ() < 10) return result;
    // Workaround to be fixed in G4NeutronHPNames.
    if (theMaterial->GetElement(check)->GetZ() > 83) return result;
  }
  G4int index(0);
  G4double* xSec = new G4double[nEleInMat];
  G4double sum = 0;
 
  for (G4int i = 0; i < nEleInMat; i++) {
    index = theMaterial->GetElement(i)->GetIndex();
    xSec[i] = theData[index]->GetCrossSection(incidentKE);
    sum += xSec[i];
  }
 
  G4double random = G4UniformRand();
  G4double running = 0;
 
  for (G4int i = 0; i < nEleInMat; i++) {
    running += xSec[i];
    index = theMaterial->GetElement(i)->GetIndex();
    if (random <= running/sum) break;
  }
 
  delete [] xSec;
  result = theData[index]->GetProductIsotope(incidentKE);
  return result;
}

---

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

• G4NeutronIsotopeProduction.hh
• G4NeutronIsotopeProduction.cc