#include <G4NeutronInelasticXS.hh>

Inheritance diagram for G4NeutronInelasticXS:

Public Member Functions
	G4NeutronInelasticXS ()

virtual	~G4NeutronInelasticXS ()

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 62 of file G4NeutronInelasticXS.hh.

Constructor & Destructor Documentation

G4NeutronInelasticXS::G4NeutronInelasticXS ( )

Definition at line 82 of file G4NeutronInelasticXS.cc.

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

   : G4VCrossSectionDataSet("G4NeutronInelasticXS"),
     proton(G4Proton::Proton())
 {
   //  verboseLevel = 0;
   if(verboseLevel > 0){
     G4cout  << "G4NeutronInelasticXS::G4NeutronInelasticXS Initialise for Z < " 
         << MAXZINEL << G4endl;
   }
   data.SetName("NeutronInelastic");
   work.resize(13,0);
   temp.resize(13,0.0);
   coeff.resize(MAXZINEL, 1.0);
   ggXsection = new G4GlauberGribovCrossSection();
   fNucleon = new G4HadronNucleonXsc();
   isInitialized = false;
 }

G4NeutronInelasticXS::~G4NeutronInelasticXS ( )

virtual

Definition at line 100 of file G4NeutronInelasticXS.cc.

 {
   delete fNucleon;
 }

Member Function Documentation

void G4NeutronInelasticXS::BuildPhysicsTable ( const G4ParticleDefinition & p )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 253 of file G4NeutronInelasticXS.cc.

References FatalException, G4cout, G4endl, G4Exception(), G4lrint(), G4Element::GetElementTable(), G4Element::GetNumberOfElements(), G4ParticleDefinition::GetParticleName(), MAXZINEL, G4Neutron::Neutron(), 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("G4NeutronInelasticXS::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");
 
   G4DynamicParticle* dynParticle = 
     new G4DynamicParticle(G4Neutron::Neutron(),G4ThreeVector(1,0,0),1);
 
   // 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 = G4lrint(((*theElmTable)[i])->GetZ());
       if(Z < 1)         { Z = 1; }
       else if(Z >= MAXZINEL) { Z = MAXZINEL-1; }
       //G4cout << "Z= " << Z << G4endl;
       // Initialisation 
       if(!data.GetElementData(Z)) { Initialise(Z, dynParticle, path); }
     }
   }
   delete dynParticle;
 }

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 105 of file G4NeutronInelasticXS.cc.

 {
   outFile << "G4NeutronInelasticXS calculates the neutron inelastic scattering\n"
           << "cross section on nuclei using data from the high precision\n"
           << "neutron database.  These data are simplified and smoothed over\n"
           << "the resonance region in order to reduce CPU time.\n"
           << "G4NeutronInelasticXS is valid for energies up to 20 MeV, for\n"
           << "nuclei through U.\n";
 }

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 131 of file G4NeutronInelasticXS.cc.

References G4PhysicsVector::Energy(), G4cout, G4endl, G4lrint(), G4HadronNucleonXsc::GetHadronNucleonXscPDG(), G4GlauberGribovCrossSection::GetInelasticGlauberGribovXsc(), G4HadronNucleonXsc::GetInelasticHadronNucleonXsc(), G4GlauberGribovCrossSection::GetIsoCrossSection(), G4DynamicParticle::GetKineticEnergy(), G4PhysicsVector::GetMaxEnergy(), G4NistManager::Instance(), MAXZINEL, G4PhysicsVector::Value(), and G4VCrossSectionDataSet::verboseLevel.

 {
   G4double xs = 0.0;
   G4double ekin = aParticle->GetKineticEnergy();
 
   if(Z < 1 || Z >= MAXZINEL) { return xs; }
   G4int Amean = G4lrint(G4NistManager::Instance()->GetAtomicMassAmu(Z));
 
   G4PhysicsVector* pv = data.GetElementData(Z);
   //  G4cout  << "G4NeutronInelasticXS::GetCrossSection e= " << ekin << " Z= " 
   //    << Z << G4endl;
 
   // element was not initialised
   if(!pv) {
     Initialise(Z);
     pv = data.GetElementData(Z);
     if(!pv) { return xs; }
   }
 
   G4double e1 = pv->Energy(0);
   if(ekin <= e1) { return xs; }
 
   G4double e2 = pv->GetMaxEnergy();
 
   if(ekin <= e2) { 
     xs = pv->Value(ekin); 
   } else if(1 == Z) { 
     fNucleon->GetHadronNucleonXscPDG(aParticle, proton);
     xs = coeff[1]*fNucleon->GetInelasticHadronNucleonXsc();
   } else {          
     ggXsection->GetIsoCrossSection(aParticle, Z, Amean);
     xs = coeff[Z]*ggXsection->GetInelasticGlauberGribovXsc();
   }
 
   if(verboseLevel > 0) {
     G4cout  << "ekin= " << ekin << ",  XSinel= " << xs << G4endl;
   }
   return xs;
 }

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 173 of file G4NeutronInelasticXS.cc.

References G4DynamicParticle::GetKineticEnergy(), and MAXZINEL.

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

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 116 of file G4NeutronInelasticXS.cc.

 {
   return true;
 }

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 123 of file G4NeutronInelasticXS.cc.

 {
   return true;
 }

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 208 of file G4NeutronInelasticXS.cc.

References G4lrint(), G4UniformRand, G4Element::GetIsotopeVector(), G4Element::GetNumberOfIsotopes(), G4Element::GetRelativeAbundanceVector(), G4Element::GetZ(), MAXZINEL, 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 >= MAXZINEL) { Z = MAXZINEL - 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