G4NeutronInelasticCrossSection Class Reference

#include <G4NeutronInelasticCrossSection.hh>

Inheritance diagram for G4NeutronInelasticCrossSection:

Public Member Functions
	G4NeutronInelasticCrossSection ()
	~G4NeutronInelasticCrossSection ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *aPart, G4int Z, const G4Material *)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *)
G4double	GetCrossSection (G4double kineticEnergy, G4int Z, G4int A)
virtual void	CrossSectionDescription (std::ostream &) const
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=nullptr)
virtual G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
G4double	ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=nullptr)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=nullptr, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
virtual const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual void	CrossSectionDescription (std::ostream &) const
virtual G4int	GetVerboseLevel () const
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

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 44 of file G4NeutronInelasticCrossSection.hh.

Constructor & Destructor Documentation

G4NeutronInelasticCrossSection()

G4NeutronInelasticCrossSection::G4NeutronInelasticCrossSection

(

)

Definition at line 40 of file G4NeutronInelasticCrossSection.cc.

  : G4VCrossSectionDataSet("Wellisch-Laidlaw"), 
    minEnergy(19.9*MeV), maxEnergy(19.9*GeV)
{}

~G4NeutronInelasticCrossSection()

G4NeutronInelasticCrossSection::~G4NeutronInelasticCrossSection

(

)

Definition at line 45 of file G4NeutronInelasticCrossSection.cc.

{}

Member Function Documentation

CrossSectionDescription()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 49 of file G4NeutronInelasticCrossSection.cc.

{
  outFile << "G4NeutronInelasticCrossSection calculates the inelastic neutron\n"
          << "scattering cross section for nuclei using the Wellisch-Laidlaw\n"
          << "parameterization between 19.9 MeV and 19.9 GeV.  Above 19.9 GeV\n"
          << "the cross section is assumed to be constant.\n";
}

GetCrossSection()

G4double G4NeutronInelasticCrossSection::GetCrossSection	(	G4double	kineticEnergy,
		G4int	Z,
		G4int	A
	)

Definition at line 74 of file G4NeutronInelasticCrossSection.cc.

{
  if(anEnergy > maxEnergy) { anEnergy = maxEnergy; }
  G4double cross_section = 0.0;
  if(anEnergy < keV) { return cross_section; }
  
  G4Pow* g4pow = G4Pow::GetInstance();
  G4double A13 = g4pow->Z13(A);
  
  G4double elog = std::log10(anEnergy/MeV);
  G4int nOfNeutrons = A - Z;
  G4double atomicNumber = G4double(A);
  static const G4double p1=1.3773;
  G4double p2 = 1. + 10./atomicNumber   - 0.0006*atomicNumber;
  G4double p3 = 0.6+ 13./atomicNumber   - 0.0005*atomicNumber;
  G4double p4 = 7.2449 - 0.018242*atomicNumber;
  G4double p5 = 1.64 - 1.8/atomicNumber - 0.0005*atomicNumber;
  G4double p6 = 1. + 200./atomicNumber + 0.02*atomicNumber;
  G4double p7 = (atomicNumber-70.)*(atomicNumber-200.)/11000.;
 
  G4double logN  = g4pow->logZ(nOfNeutrons);
  G4double part1 = pi*p1*p1*logN;
  G4double part2 = 1.+ A13 - p2*(1.-1./A13);
 
  G4double firstexp = -p4*(elog-p5);
  G4double first    = 1. + G4Exp(firstexp);
  G4double corr     = 1. + p3*(1.-1./first); 
 
  G4double secondexp= -p6*(elog-p7);
  G4double secondv   = 1.+G4Exp(secondexp);
  G4double corr2    = 1./secondv;
 
  G4double xsec = corr*corr2*part1*part2*10.*millibarn;
  if(xsec < 0.0) { xsec = 0.0; }
  return xsec;
}

Referenced by GetElementCrossSection().

GetElementCrossSection()

G4double G4NeutronInelasticCrossSection::GetElementCrossSection ( const G4DynamicParticle *  aPart, G4int  Z, const G4Material *    )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 65 of file G4NeutronInelasticCrossSection.cc.

{
  G4int A = G4int(G4NistManager::Instance()->GetAtomicMassAmu(Z));
  return GetCrossSection(aPart->GetKineticEnergy(), Z, A);
}

IsElementApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 58 of file G4NeutronInelasticCrossSection.cc.

{
  G4double e = part->GetKineticEnergy();
  return (1 < Z && e > minEnergy);    
}

---

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

• G4NeutronInelasticCrossSection.hh
• G4NeutronInelasticCrossSection.cc