G4NeutronElasticXS Class Reference

#include <G4NeutronElasticXS.hh>

Inheritance diagram for G4NeutronElasticXS:

Public Member Functions
	G4NeutronElasticXS ()
	~G4NeutronElasticXS () final
G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *) final
G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *, const G4Material *) final
G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *) final
G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso, const G4Element *elm, const G4Material *mat) final
G4double	ComputeCrossSectionPerElement (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, const G4Element *, const G4Material *) final
G4double	ComputeIsoCrossSection (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, G4int Z, G4int A, const G4Isotope *iso, const G4Element *elm, const G4Material *mat) final
const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE) final
void	BuildPhysicsTable (const G4ParticleDefinition &) final
void	CrossSectionDescription (std::ostream &) const final
G4double	ElementCrossSection (G4double kinEnergy, G4double loge, G4int Z)
G4NeutronElasticXS &	operator= (const G4NeutronElasticXS &right)=delete
	G4NeutronElasticXS (const G4NeutronElasticXS &)=delete
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
G4double	ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual void	SetVerboseLevel (G4int value)
G4double	GetMinKinEnergy () const
void	SetMinKinEnergy (G4double value)
G4double	GetMaxKinEnergy () const
void	SetMaxKinEnergy (G4double value)
bool	ForAllAtomsAndEnergies () const
void	SetForAllAtomsAndEnergies (G4bool val)
const G4String &	GetName () const
void	SetName (const G4String &nam)
G4VCrossSectionDataSet &	operator= (const G4VCrossSectionDataSet &right)=delete
	G4VCrossSectionDataSet (const G4VCrossSectionDataSet &)=delete

Static Public Member Functions
static const char *	Default_Name ()

Additional Inherited Members
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel {0}
G4String	name

Detailed Description

Definition at line 55 of file G4NeutronElasticXS.hh.

Constructor & Destructor Documentation

G4NeutronElasticXS() [1/2]

G4NeutronElasticXS::G4NeutronElasticXS

(

)

Definition at line 66 of file G4NeutronElasticXS.cc.

 : G4VCrossSectionDataSet(Default_Name()),
   neutron(G4Neutron::Neutron())
{
  //  verboseLevel = 0;
  if (verboseLevel > 0){
    G4cout  << "G4NeutronElasticXS::G4NeutronElasticXS Initialise for Z < " 
        << MAXZEL << G4endl;
  }
  ggXsection = 
    G4CrossSectionDataSetRegistry::Instance()->GetComponentCrossSection("Glauber-Gribov");
  if (ggXsection == nullptr)
    ggXsection = new G4ComponentGGHadronNucleusXsc();
  SetForAllAtomsAndEnergies(true);
  FindDirectoryPath();
}

Referenced by G4NeutronElasticXS(), and operator=().

~G4NeutronElasticXS()

G4NeutronElasticXS::~G4NeutronElasticXS ( )

final

Definition at line 83 of file G4NeutronElasticXS.cc.

{
  if (isFirst) {
    for(G4int i=0; i<MAXZEL; ++i) {
      delete data[i];
      data[i] = nullptr;
    }
  }
}

G4NeutronElasticXS() [2/2]

G4NeutronElasticXS::G4NeutronElasticXS ( const G4NeutronElasticXS & )

delete

Member Function Documentation

BuildPhysicsTable()

void G4NeutronElasticXS::BuildPhysicsTable ( const G4ParticleDefinition & p )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 198 of file G4NeutronElasticXS.cc.

{
  if(verboseLevel > 0){
    G4cout << "G4NeutronElasticXS::BuildPhysicsTable for " 
       << p.GetParticleName() << G4endl;
  }
  if(p.GetParticleName() != "neutron") { 
    G4ExceptionDescription ed;
    ed << p.GetParticleName() << " is a wrong particle type -"
       << " only neutron is allowed";
    G4Exception("G4NeutronElasticXS::BuildPhysicsTable(..)","had012",
        FatalException, ed, "");
    return; 
  }
  if (fLock || isFirst) { 
    G4AutoLock l(&nElasticXSMutex);
    if (fLock) { 
      isFirst = true;
      fLock = false;
      FindDirectoryPath();
    }
 
    // Access to elements
    const G4ElementTable* table = G4Element::GetElementTable();
    for ( auto & elm : *table ) {
      G4int Z = std::max( 1, std::min( elm->GetZasInt(), MAXZEL-1) );
      if ( nullptr == data[Z] ) { Initialise(Z); }
    }
    l.unlock();
  }
}

ComputeCrossSectionPerElement()

G4double G4NeutronElasticXS::ComputeCrossSectionPerElement ( G4double kinEnergy, G4double loge, const G4ParticleDefinition * , const G4Element * elm, const G4Material *  )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 125 of file G4NeutronElasticXS.cc.

{
  return ElementCrossSection(ekin, loge, elm->GetZasInt());
}

ComputeIsoCrossSection()

G4double G4NeutronElasticXS::ComputeIsoCrossSection ( G4double kinEnergy, G4double loge, const G4ParticleDefinition * , G4int Z, G4int A, const G4Isotope * iso, const G4Element * elm, const G4Material * mat )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 153 of file G4NeutronElasticXS.cc.

{
  return ElementCrossSection(ekin, loge, Z)*A/aeff[Z];
}

CrossSectionDescription()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 93 of file G4NeutronElasticXS.cc.

{
  outFile << "G4NeutronElasticXS calculates the neutron elastic 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"
          << "For high energies Glauber-Gribiv cross section is used.\n";
}

Default_Name()

static const char * G4NeutronElasticXS::Default_Name ( )

inlinestatic

Definition at line 63 of file G4NeutronElasticXS.hh.

{return "G4NeutronElasticXS";}

Referenced by G4NeutronElasticXS().

ElementCrossSection()

G4double G4NeutronElasticXS::ElementCrossSection	(	G4double	kinEnergy,
		G4double	loge,
		G4int	Z )

Definition at line 133 of file G4NeutronElasticXS.cc.

{
  G4int Z = (ZZ >= MAXZEL) ? MAXZEL - 1 : ZZ;
  auto pv = GetPhysicsVector(Z);
 
  G4double xs = (ekin <= pv->GetMaxEnergy()) ? pv->LogVectorValue(ekin, loge) 
    : coeff[Z]*ggXsection->GetElasticElementCrossSection(neutron, ekin,
                                                         Z, aeff[Z]);
 
#ifdef G4VERBOSE
  if(verboseLevel > 1) {
    G4cout  << "Z= " << Z << " Ekin(MeV)= " << ekin/CLHEP::MeV 
        << ",  nElmXSel(b)= " << xs/CLHEP::barn 
        << G4endl;
  }
#endif
  return xs;
}

Referenced by ComputeCrossSectionPerElement(), ComputeIsoCrossSection(), GetElementCrossSection(), and GetIsoCrossSection().

GetElementCrossSection()

G4double G4NeutronElasticXS::GetElementCrossSection ( const G4DynamicParticle * aParticle, G4int Z, const G4Material *  )

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 117 of file G4NeutronElasticXS.cc.

{
  return ElementCrossSection(aParticle->GetKineticEnergy(),
                 aParticle->GetLogKineticEnergy(), Z);
}

GetIsoCrossSection()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 163 of file G4NeutronElasticXS.cc.

{
  return ElementCrossSection(aParticle->GetKineticEnergy(),
                 aParticle->GetLogKineticEnergy(), Z)*A/aeff[Z];
 
}

IsElementApplicable()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 103 of file G4NeutronElasticXS.cc.

{
  return true;
}

IsIsoApplicable()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 109 of file G4NeutronElasticXS.cc.

{
  return false;
}

operator=()

G4NeutronElasticXS & G4NeutronElasticXS::operator= ( const G4NeutronElasticXS & right )

delete

SelectIsotope()

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

finalvirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 173 of file G4NeutronElasticXS.cc.

{
  G4int nIso = (G4int)anElement->GetNumberOfIsotopes();
  const G4Isotope* iso = anElement->GetIsotope(0);
 
  //G4cout << "SelectIsotope NIso= " << nIso << G4endl;
  if(1 == nIso) { return iso; }
 
  const G4double* abundVector = anElement->GetRelativeAbundanceVector();
  G4double q = G4UniformRand();
  G4double sum = 0.0;
 
  // isotope wise cross section not used
  for (G4int j=0; j<nIso; ++j) {
    sum += abundVector[j];
    if(q <= sum) {
      iso = anElement->GetIsotope(j);
      break;
    }
  }
  return iso;
}

---

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

• G4NeutronElasticXS.hh
• G4NeutronElasticXS.cc