#include <G4NeutronCaptureXS.hh>

Inheritance diagram for G4NeutronCaptureXS:

Public Member Functions
	G4NeutronCaptureXS ()

virtual	~G4NeutronCaptureXS ()

virtual G4bool	IsElementApplicable (const G4DynamicParticle , G4int Z, const G4Material )

virtual G4bool	IsIsoApplicable (const G4DynamicParticle , G4int Z, G4int A, const G4Element , const G4Material *)

virtual G4double	GetElementCrossSection (const G4DynamicParticle , G4int Z, const G4Material mat=0)

virtual G4double	GetIsoCrossSection (const G4DynamicParticle , G4int Z, G4int A, const G4Isotope iso, const G4Element elm, const G4Material mat)

virtual G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy)

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	CrossSectionDescription (std::ostream &) const

Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")

virtual	~G4VCrossSectionDataSet ()

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=0)

G4double	ComputeCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=0)

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

virtual G4int	GetVerboseLevel () const

virtual void	SetVerboseLevel (G4int value)

G4double	GetMinKinEnergy () const

void	SetMinKinEnergy (G4double value)

G4double	GetMaxKinEnergy () const

void	SetMaxKinEnergy (G4double value)

const G4String &	GetName () const

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)

Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Detailed Description

Definition at line 61 of file G4NeutronCaptureXS.hh.

Constructor & Destructor Documentation

G4NeutronCaptureXS::G4NeutronCaptureXS ( )

Definition at line 77 of file G4NeutronCaptureXS.cc.

References G4cout, G4endl, MAXZCAPTURE, and G4VCrossSectionDataSet::verboseLevel.

  : G4VCrossSectionDataSet("G4NeutronCaptureXS"),
    emax(20*MeV),elimit(1.0e-10*eV)
 {
   //  verboseLevel = 0;
   if(verboseLevel > 0){
     G4cout  << "G4NeutronCaptureXS::G4NeutronCaptureXS: Initialise for Z < "
         << MAXZCAPTURE << G4endl;
   }
   data.SetName("NeutronCapture");
   work.resize(13,0);
   temp.resize(13,0.0);
   isInitialized = false;
 }

G4NeutronCaptureXS::~G4NeutronCaptureXS ( )

virtual

Definition at line 92 of file G4NeutronCaptureXS.cc.

93 {}

Member Function Documentation

void G4NeutronCaptureXS::BuildPhysicsTable ( const G4ParticleDefinition & p )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 235 of file G4NeutronCaptureXS.cc.

References FatalException, G4cout, G4endl, G4Exception(), G4Element::GetElementTable(), G4Element::GetNumberOfElements(), G4ParticleDefinition::GetParticleName(), MAXZCAPTURE, and G4VCrossSectionDataSet::verboseLevel.

 {
   if(isInitialized) { return; }
   if(verboseLevel > 0){
     G4cout << "G4NeutronCaptureXS::BuildPhysicsTable for " 
        << p.GetParticleName() << G4endl;
   }
   if(p.GetParticleName() != "neutron") { 
     G4ExceptionDescription ed;
     ed << p.GetParticleName() << " is a wrong particle type -"
        << " only neutron is allowed";
     G4Exception("G4NeutronCaptureXS::BuildPhysicsTable(..)","had012",
         FatalException, ed, "");
     return; 
   }
   isInitialized = true;
 
   // check environment variable 
   // Build the complete string identifying the file with the data set
   char* path = getenv("G4NEUTRONXSDATA");
 
   // Access to elements
   const G4ElementTable* theElmTable = G4Element::GetElementTable();
   size_t numOfElm = G4Element::GetNumberOfElements();
   if(numOfElm > 0) {
     for(size_t i=0; i<numOfElm; ++i) {
       G4int Z = G4int(((*theElmTable)[i])->GetZ());
       if(Z < 1)         { Z = 1; }
       else if(Z >= MAXZCAPTURE) { Z = MAXZCAPTURE-1; }
       //G4cout << "Z= " << Z << G4endl;
       // Initialisation 
       if(!data.GetElementData(Z)) { Initialise(Z, path); }
     }
   }
 }

void G4NeutronCaptureXS::CrossSectionDescription ( std::ostream & outFile ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 95 of file G4NeutronCaptureXS.cc.

 {
   outFile << "G4NeutronCaptureXS calculates the neutron capture cross sections\n"
           << "on nuclei using data from the high precision neutron database.\n"
           << "These data are simplified and smoothed over the resonance region\n"
           << "in order to reduce CPU time.  G4NeutronCaptureXS is valid up to\n"
           << "20 MeV for all targets through U.\n";
 }

G4double G4NeutronCaptureXS::GetElementCrossSection	(	const G4DynamicParticle *	aParticle,
		G4int	Z,
		const G4Material *	mat = `0`
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 120 of file G4NeutronCaptureXS.cc.

References G4PhysicsVector::Energy(), G4cout, G4endl, G4DynamicParticle::GetKineticEnergy(), MAXZCAPTURE, G4PhysicsVector::Value(), and G4VCrossSectionDataSet::verboseLevel.

 {
   G4double xs = 0.0;
   G4double ekin = aParticle->GetKineticEnergy();
   if(ekin > emax || Z < 1 || Z >= MAXZCAPTURE) { return xs; }
   if(ekin < elimit) { ekin = elimit; }
 
   G4PhysicsVector* pv = data.GetElementData(Z);
 
   // element was not initialised
   if(!pv) {
     Initialise(Z);
     pv = data.GetElementData(Z);
     if(!pv) { return xs; }
   }
 
   G4double e1 = pv->Energy(0);
   if(ekin < e1) { xs = (*pv)[0]*std::sqrt(e1/ekin); }
   else if(ekin <= pv->GetMaxEnergy()) { xs = pv->Value(ekin); }
 
   if(verboseLevel > 0){
     G4cout  << "ekin= " << ekin << ",  xs= " << xs << G4endl;
   }
   return xs;
 }

G4double G4NeutronCaptureXS::GetIsoCrossSection	(	const G4DynamicParticle *	aParticle,
		G4int	Z,
		G4int	A,
		const G4Isotope *	iso,
		const G4Element *	elm,
		const G4Material *	mat
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 148 of file G4NeutronCaptureXS.cc.

References G4DynamicParticle::GetKineticEnergy(), and MAXZCAPTURE.

 {
   G4double xs = 0.0;
   G4double ekin = aParticle->GetKineticEnergy();
   if(ekin <= emax && Z > 0 && Z < MAXZCAPTURE) { 
     xs = IsoCrossSection(ekin, Z, A); 
   }
   return xs;
 }

G4bool G4NeutronCaptureXS::IsElementApplicable	(	const G4DynamicParticle *	,
		G4int	Z,
		const G4Material *
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 105 of file G4NeutronCaptureXS.cc.

 {
   return true;
 }

G4bool G4NeutronCaptureXS::IsIsoApplicable	(	const G4DynamicParticle *	,
		G4int	Z,
		G4int	A,
		const G4Element *	,
		const G4Material *
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 112 of file G4NeutronCaptureXS.cc.

 {
   return true;
 }

G4Isotope * G4NeutronCaptureXS::SelectIsotope	(	const G4Element *	anElement,
		G4double	kinEnergy
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 190 of file G4NeutronCaptureXS.cc.

References G4lrint(), G4UniformRand, G4Element::GetIsotopeVector(), G4Element::GetNumberOfIsotopes(), G4Element::GetRelativeAbundanceVector(), G4Element::GetZ(), MAXZCAPTURE, and nmax.

 {
   G4int nIso = anElement->GetNumberOfIsotopes();
   G4IsotopeVector* isoVector = anElement->GetIsotopeVector();
   G4Isotope* iso = (*isoVector)[0];
 
   // more than 1 isotope
   if(1 < nIso) {
     G4int Z = G4lrint(anElement->GetZ());
     if(Z >= MAXZCAPTURE) { Z = MAXZCAPTURE-1; }
     G4double* abundVector = anElement->GetRelativeAbundanceVector();
     G4double q = G4UniformRand();
     G4double sum = 0.0;
 
     // is there isotope wise cross section?
     if(0 == amin[Z]) {
       for (G4int j = 0; j<nIso; ++j) {
     sum += abundVector[j];
     if(q <= sum) {
       iso = (*isoVector)[j];
       break;
     }
       }
     } else {
       size_t nmax = data.GetNumberOfComponents(Z);
       if(temp.size() < nmax) { temp.resize(nmax,0.0); }
       for (size_t i=0; i<nmax; ++i) {
     G4int A = (*isoVector)[i]->GetN();
         sum += abundVector[i]*IsoCrossSection(kinEnergy, Z, A);
         temp[i] = sum;
       }
       sum *= q;
       for (size_t j = 0; j<nmax; ++j) {
         if(temp[j] >= sum) {
           iso = (*isoVector)[j];
           break;
     }
       }
     }
   }
   return iso;
 }

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