#include <G4NeutronInelasticXS.hh>

Inheritance diagram for G4NeutronInelasticXS:

Collaboration 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

Static Public Member Functions
static const char *	Default_Name ()

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 93 of file G4NeutronInelasticXS.cc.

   : G4VCrossSectionDataSet(Default_Name()),
     proton(G4Proton::Proton())
 {
   //  verboseLevel = 0;
   if(verboseLevel > 0){
     G4cout << "G4NeutronInelasticXS::G4NeutronInelasticXS Initialise for Z < " 
         << MAXZINEL << G4endl;
   }
   ggXsection = new G4ComponentGGHadronNucleusXsc();
   fNucleon = new G4HadronNucleonXsc();
   isMaster = false;
 }

G4NeutronInelasticXS::~G4NeutronInelasticXS ( )

virtual

Definition at line 107 of file G4NeutronInelasticXS.cc.

 {
   //G4cout << "G4NeutronInelasticXS::~G4NeutronInelasticXS() " 
   // << " isMaster= " << isMaster << "  data: " << data << G4endl;
   delete fNucleon;
   delete ggXsection;
   if(isMaster) { delete data; data = 0; }
 }

Member Function Documentation

void G4NeutronInelasticXS::BuildPhysicsTable ( const G4ParticleDefinition & p )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 279 of file G4NeutronInelasticXS.cc.

 {
   if(verboseLevel > 0){
     G4cout << "G4NeutronInelasticXS::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; 
   }
 
   if(!data) { 
     isMaster = true;
     data = new G4ElementData(); 
     data->SetName("NeutronInelastic");
     temp.resize(13,0.0);
   }
 
   // it is possible re-initialisation for the new run
   if(isMaster) {
 
     // 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;
   }
 }

Here is the call graph for this function:

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 116 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";
 }

static const char* G4NeutronInelasticXS::Default_Name ( )

inlinestatic

Definition at line 70 of file G4NeutronInelasticXS.hh.

70 {return "G4NeutronInelasticXS";}

Here is the caller graph for this function:

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 141 of file G4NeutronInelasticXS.cc.

 {
   G4double xs = 0.0;
   G4double ekin = aParticle->GetKineticEnergy();
 
   if(Z >= MAXZINEL) { Z = MAXZINEL - 1; }
   else if(Z < 1)    { Z = 1; }
 
   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;
 }

Here is the call graph for this function:

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 185 of file G4NeutronInelasticXS.cc.

 {
   if(Z >= MAXZINEL) { Z = MAXZINEL - 1; }
   else if(Z < 1)    { Z = 1; }
 
   return IsoCrossSection(aParticle->GetKineticEnergy(), Z, A);
 }

Here is the call graph for this function:

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 127 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 134 of file G4NeutronInelasticXS.cc.

 {
   return true;
 }

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 224 of file G4NeutronInelasticXS.cc.

 {
   size_t nIso = anElement->GetNumberOfIsotopes();
   G4IsotopeVector* isoVector = anElement->GetIsotopeVector();
   G4Isotope* iso = (*isoVector)[0];
 
   //G4cout << "SelectIsotope NIso= " << nIso << G4endl;
 
   // more than 1 isotope
   if(1 < nIso) {
     G4int Z = G4lrint(anElement->GetZ());
     //G4cout << "SelectIsotope Z= " << Z << G4endl;
 
     G4double* abundVector = anElement->GetRelativeAbundanceVector();
     G4double q = G4UniformRand();
     G4double sum = 0.0;
 
     // is there isotope wise cross section?
     size_t j;
     if(0 == amin[Z] || Z >= MAXZINEL) {
       for (j = 0; j<nIso; ++j) {
     sum += abundVector[j];
     if(q <= sum) {
       iso = (*isoVector)[j];
       break;
     }
       }
     } else {
 
       // element may be not initialised in unit test
       if(!data->GetElementData(Z)) { Initialise(Z); }
       size_t nn = temp.size();
       if(nn < nIso) { temp.resize(nIso, 0.); }
 
       for (j=0; j<nIso; ++j) {
     //G4cout << j << "-th isotope " << (*isoVector)[j]->GetN() 
     //       <<  " abund= " << abundVector[j] << G4endl;
         sum += abundVector[j]*IsoCrossSection(kinEnergy, Z, 
                           (*isoVector)[j]->GetN());
         temp[j] = sum;
       }
       sum *= q;
       for (j = 0; j<nIso; ++j) {
         if(temp[j] >= sum) {
           iso = (*isoVector)[j];
           break;
     }
       }
     }
   }
   return iso;
 }

Here is the call graph for this function:

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

geant4.10.03.p01/source/processes/hadronic/cross_sections/include/G4NeutronInelasticXS.hh
geant4.10.03.p01/source/processes/hadronic/cross_sections/src/G4NeutronInelasticXS.cc

Public Member Functions

Static Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation