G4NeutronHPCaptureData Class Reference

#include <G4NeutronHPCaptureData.hh>

Inheritance diagram for G4NeutronHPCaptureData:

Public Member Functions
	G4NeutronHPCaptureData ()
	~G4NeutronHPCaptureData () override
G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int Z, G4int A, const G4Element *, const G4Material *) override
G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *, const G4Element *, const G4Material *) override
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, G4double aT)
void	BuildPhysicsTable (const G4ParticleDefinition &) override
void	DumpPhysicsTable (const G4ParticleDefinition &) override
void	CrossSectionDescription (std::ostream &) const override
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, 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	ComputeCrossSectionPerElement (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, const G4Element *, const G4Material *mat=nullptr)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=nullptr)
virtual G4double	ComputeIsoCrossSection (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, 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	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

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

Detailed Description

Definition at line 55 of file G4NeutronHPCaptureData.hh.

Constructor & Destructor Documentation

G4NeutronHPCaptureData()

G4NeutronHPCaptureData::G4NeutronHPCaptureData

(

)

Definition at line 65 of file G4NeutronHPCaptureData.cc.

  : G4VCrossSectionDataSet("NeutronHPCaptureXS")
{
  emax = 20.*CLHEP::MeV;
  fManager = G4ParticleHPManager::GetInstance();
}

~G4NeutronHPCaptureData()

G4NeutronHPCaptureData::~G4NeutronHPCaptureData ( )

override

Definition at line 72 of file G4NeutronHPCaptureData.cc.

{
  if (isFirst) {
    if (nullptr != theCrossSections)
      theCrossSections->clearAndDestroy();
    delete theCrossSections;
    theCrossSections = nullptr;
  }
}

Member Function Documentation

BuildPhysicsTable()

void G4NeutronHPCaptureData::BuildPhysicsTable ( const G4ParticleDefinition & p )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 108 of file G4NeutronHPCaptureData.cc.

{
  // the choice of the first instance
  if (fLock) {
    G4AutoLock l(&theHPCaptureData);
    if (fLock) {
      isFirst = true;
      fLock = false;
    }
    l.unlock();
  }
  if (!isFirst) { return; }
  if (p.GetParticleName() != "neutron") { 
    G4ExceptionDescription ed;
    ed << p.GetParticleName() << " is a wrong particle type -"
       << " only neutron is allowed";
    G4Exception("G4NeutronHPCaptureData::BuildPhysicsTable(..)","had012",
                FatalException, ed, "");
    return; 
  }
 
  // initialisation for the first instance, others are locked
  G4AutoLock l(&theHPCaptureData);
  if (theCrossSections != nullptr) {
    theCrossSections->clearAndDestroy();
    delete theCrossSections;
  }
  std::size_t numberOfElements = G4Element::GetNumberOfElements();
  theCrossSections = new G4PhysicsTable(numberOfElements);
  
  // make a PhysicsVector for each element
  auto theElementTable = G4Element::GetElementTable();
  for (std::size_t i = 0; i < numberOfElements; ++i) {
    auto elm = (*theElementTable)[i];
#ifdef G4VERBOSE
    if (fManager->GetDEBUG()) {
      G4cout << "ElementIndex " << elm->GetIndex() << "  " 
             << elm->GetName() << G4endl;
    }
#endif
    G4PhysicsVector* physVec =
      G4ParticleHPData::Instance(G4Neutron::Neutron())->MakePhysicsVector(elm, this);
    theCrossSections->push_back(physVec);
  }
  fManager->RegisterCaptureCrossSections(theCrossSections);
  l.unlock();
}

CrossSectionDescription()

void G4NeutronHPCaptureData::CrossSectionDescription ( std::ostream & outF ) const

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 254 of file G4NeutronHPCaptureData.cc.

{
  outF << "High Precision cross data based on Evaluated Nuclear Data Files"
       << " (ENDF) for radiative capture reaction of neutrons below 20 MeV";
}

DumpPhysicsTable()

void G4NeutronHPCaptureData::DumpPhysicsTable ( const G4ParticleDefinition & )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 156 of file G4NeutronHPCaptureData.cc.

{
#ifdef G4VERBOSE
  if (fManager->GetVerboseLevel() == 0) return;
 
  //
  // Dump element based cross section
  // range 10e-5 eV to 20 MeV
  // 10 point per decade
  // in barn
  //
 
  G4cout << G4endl;
  G4cout << G4endl;
  G4cout << "Capture Cross Section of Neutron HP" << G4endl;
  G4cout << "(Pointwise cross-section at 0 Kelvin.)" << G4endl;
  G4cout << G4endl;
  G4cout << "Name of Element" << G4endl;
  G4cout << "Energy[eV]  XS[barn]" << G4endl;
  G4cout << G4endl;
 
  std::size_t numberOfElements = G4Element::GetNumberOfElements();
  auto theElementTable = G4Element::GetElementTable();
 
  for (std::size_t i = 0; i < numberOfElements; ++i) {
    G4cout << (*theElementTable)[i]->GetName() << G4endl;
    G4cout << *((*theCrossSections)(i)) << G4endl;
  }
#endif
}

GetCrossSection()

G4double G4NeutronHPCaptureData::GetCrossSection	(	const G4DynamicParticle *	aP,
		const G4Element *	anE,
		G4double	aT )

Definition at line 187 of file G4NeutronHPCaptureData.cc.

{
  G4double result = 0;
  auto idx = (G4int)anE->GetIndex();
 
  // prepare neutron
  G4double eKinetic = aP->GetKineticEnergy();
  if (eKinetic >= emax) { return 0.0; }
 
  // NEGLECT_DOPPLER
  if (fManager->GetNeglectDoppler()) {
    return (*((*theCrossSections)(idx))).Value(eKinetic);
  }
 
  G4ReactionProduct theNeutron(aP->GetDefinition());
  theNeutron.SetMomentum(aP->GetMomentum());
  theNeutron.SetKineticEnergy(eKinetic);
 
  // prepare thermal nucleus
  G4Nucleus aNuc;
  G4int theA = anE->GetN();
  G4int theZ = anE->GetZasInt();
  G4double eleMass = G4NucleiProperties::GetNuclearMass(theA, theZ)
    / CLHEP::neutron_mass_c2;
 
  G4ReactionProduct boosted;
  G4double aXsection;
 
  // MC integration loop
  G4int counter = 0;
  G4double buffer = 0;
  G4int size = G4int(std::max(10., aT / 60 * kelvin));
  G4ThreeVector neutronVelocity =
    1. / G4Neutron::Neutron()->GetPDGMass() * theNeutron.GetMomentum();
  G4double neutronVMag = neutronVelocity.mag();
 
  while (counter == 0 ||
         std::abs(buffer - result / std::max(1, counter)) > 0.03 * buffer)
  // Loop checking, 11.05.2015, T. Koi
  {
    if (counter != 0) buffer = result / counter;
    while (counter < size)  // Loop checking, 11.05.2015, T. Koi
    {
      ++counter;
      G4ReactionProduct aThermalNuc = aNuc.GetThermalNucleus(eleMass, aT);
      boosted.Lorentz(theNeutron, aThermalNuc);
      G4double theEkin = boosted.GetKineticEnergy();
      aXsection = (*((*theCrossSections)(idx))).Value(theEkin);
      // velocity correction, or luminosity factor...
      G4ThreeVector targetVelocity = 1. / aThermalNuc.GetMass() * aThermalNuc.GetMomentum();
      aXsection *= (targetVelocity - neutronVelocity).mag() / neutronVMag;
      result += aXsection;
    }
    size += size;
  }
  result /= counter;
  /*
    // Checking impact of  G4NEUTRONHP_NEGLECT_DOPPLER
    G4cout << " result " << result << " "
           << (*((*theCrossSections)(index))).Value(eKinetic) << " "
           << (*((*theCrossSections)(index))).Value(eKinetic) /result << G4endl;
  */
  return result;
}

Referenced by GetIsoCrossSection().

GetIsoCrossSection()

G4double G4NeutronHPCaptureData::GetIsoCrossSection ( const G4DynamicParticle * dp, G4int Z, G4int A, const G4Isotope * , const G4Element * element, const G4Material * material )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 90 of file G4NeutronHPCaptureData.cc.

{
  if (dp->GetKineticEnergy() == ke_cache && element == element_cache &&
      material == material_cache)
    return xs_cache;
 
  ke_cache = dp->GetKineticEnergy();
  element_cache = element;
  material_cache = material;
  G4double xs = GetCrossSection(dp, element, material->GetTemperature());
  xs_cache = xs;
  return xs;
}

IsIsoApplicable()

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

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 82 of file G4NeutronHPCaptureData.cc.

{
  return true;
}

---

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

• G4NeutronHPCaptureData.hh
• G4NeutronHPCaptureData.cc