G4ParticleHPFissionData Class Reference

#include <G4ParticleHPFissionData.hh>

Inheritance diagram for G4ParticleHPFissionData:

Public Member Functions
	G4ParticleHPFissionData ()
	~G4ParticleHPFissionData () override
G4bool	IsIsoApplicable (const G4DynamicParticle *, G4int, G4int, const G4Element *, const G4Material *) override
G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int, G4int, 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
G4int	GetVerboseLevel () const
void	SetVerboseLevel (G4int) 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)
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 48 of file G4ParticleHPFissionData.hh.

Constructor & Destructor Documentation

G4ParticleHPFissionData()

G4ParticleHPFissionData::G4ParticleHPFissionData

(

)

Definition at line 49 of file G4ParticleHPFissionData.cc.

                                                 : G4VCrossSectionDataSet("NeutronHPFissionXS")
{
  SetMinKinEnergy(0 * MeV);
  SetMaxKinEnergy(20 * MeV);
 
  theCrossSections = nullptr;
  instanceOfWorker = false;
  if (G4Threading::IsWorkerThread()) {
    instanceOfWorker = true;
  }
  element_cache = nullptr;
  material_cache = nullptr;
  ke_cache = 0.0;
  xs_cache = 0.0;
}

~G4ParticleHPFissionData()

G4ParticleHPFissionData::~G4ParticleHPFissionData ( )

override

Definition at line 65 of file G4ParticleHPFissionData.cc.

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

Member Function Documentation

BuildPhysicsTable()

void G4ParticleHPFissionData::BuildPhysicsTable ( const G4ParticleDefinition & )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 99 of file G4ParticleHPFissionData.cc.

{
  if (G4Threading::IsWorkerThread()) {
    theCrossSections = G4ParticleHPManager::GetInstance()->GetFissionCrossSections();
    return;
  }
 
  std::size_t numberOfElements = G4Element::GetNumberOfElements();
  if (theCrossSections == nullptr)
    theCrossSections = new G4PhysicsTable(numberOfElements);
  else
    theCrossSections->clearAndDestroy();
 
  // make a PhysicsVector for each element
 
  auto theElementTable = G4Element::GetElementTable();
  for (std::size_t i = 0; i < numberOfElements; ++i) {
    G4PhysicsVector* physVec = G4ParticleHPData::Instance(G4Neutron::Neutron())
                                 ->MakePhysicsVector((*theElementTable)[i], this);
    theCrossSections->push_back(physVec);
  }
 
  G4ParticleHPManager::GetInstance()->RegisterFissionCrossSections(theCrossSections);
}

CrossSectionDescription()

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

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 253 of file G4ParticleHPFissionData.cc.

{
  outFile << "High Precision cross data based on Evaluated Nuclear Data Files (ENDF)\n"
          << "for induced fission reaction of neutrons below 20MeV\n";
}

DumpPhysicsTable()

void G4ParticleHPFissionData::DumpPhysicsTable ( const G4ParticleDefinition & )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 124 of file G4ParticleHPFissionData.cc.

{
#ifdef G4VERBOSE
  if (G4HadronicParameters::Instance()->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 << "Fission 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;
 
    if ((*((*theCrossSections)(i))).GetVectorLength() == 0) {
      G4cout << "The cross-section data of the fission of this element is not available." << G4endl;
      G4cout << G4endl;
      continue;
    }
 
    for (G4int ie = 0; ie < 130; ++ie) {
      G4double eKinetic = 1.0e-5 * G4Pow::GetInstance()->powA(10.0, ie / 10.0) * eV;
      G4bool outOfRange = false;
 
      if (eKinetic < 20 * MeV) {
        G4cout << eKinetic / eV << " "
               << (*((*theCrossSections)(i))).GetValue(eKinetic, outOfRange) / barn << G4endl;
      }
    }
 
    G4cout << G4endl;
  }
#endif
}

GetCrossSection()

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

Definition at line 171 of file G4ParticleHPFissionData.cc.

{
  G4double result = 0;
  if (anE->GetZ() < 88) return result;
  G4bool outOfRange;
  auto index = (G4int)anE->GetIndex();
 
  if (((*theCrossSections)(index))->GetVectorLength() == 0) return result;
 
  // prepare neutron
  G4double eKinetic = aP->GetKineticEnergy();
  G4ReactionProduct theNeutronRP(aP->GetDefinition());
  theNeutronRP.SetMomentum(aP->GetMomentum());
  theNeutronRP.SetKineticEnergy(eKinetic);
 
  if (G4ParticleHPManager::GetInstance()->GetNeglectDoppler()) {
    // NEGLECT_DOPPLER
    G4double factor = 1.0;
    if (eKinetic < aT * k_Boltzmann) {
      // below 0.1 eV neutrons
      // Have to do some, but now just igonre.
      // Will take care after performance check.
      // factor = factor * targetV;
    }
    return ((*((*theCrossSections)(index))).GetValue(eKinetic, outOfRange)) * factor;
  }
 
  // prepare thermal nucleus
  G4Nucleus aNuc;
  G4double eps = 0.0001;
  G4double theA = anE->GetN();
  G4double theZ = anE->GetZ();
  G4double eleMass;
  eleMass = (G4NucleiProperties::GetNuclearMass(static_cast<G4int>(theA + eps),
                                                static_cast<G4int>(theZ + eps)))
            / G4Neutron::Neutron()->GetPDGMass();
 
  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() * theNeutronRP.GetMomentum();
  G4double neutronVMag = neutronVelocity.mag();
 
  while (counter == 0
         || std::abs(buffer - result / std::max(1, counter))
              > 0.01 * 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(theNeutronRP, aThermalNuc);
      G4double theEkin = boosted.GetKineticEnergy();
      aXsection = (*((*theCrossSections)(index))).GetValue(theEkin, outOfRange);
      // velocity correction.
      G4ThreeVector targetVelocity = 1. / aThermalNuc.GetMass() * aThermalNuc.GetMomentum();
      aXsection *= (targetVelocity - neutronVelocity).mag() / neutronVMag;
      result += aXsection;
    }
    size += size;
  }
  result /= counter;
  return result;
}

Referenced by GetIsoCrossSection().

GetIsoCrossSection()

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

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 83 of file G4ParticleHPFissionData.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;
}

GetVerboseLevel()

G4int G4ParticleHPFissionData::GetVerboseLevel

(

)

const

Definition at line 243 of file G4ParticleHPFissionData.cc.

{
  return G4ParticleHPManager::GetInstance()->GetVerboseLevel();
}

Referenced by DumpPhysicsTable().

IsIsoApplicable()

G4bool G4ParticleHPFissionData::IsIsoApplicable ( const G4DynamicParticle * dp, G4int , G4int , const G4Element * , const G4Material *  )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 74 of file G4ParticleHPFissionData.cc.

{
  G4double eKin = dp->GetKineticEnergy();
  return eKin <= GetMaxKinEnergy() && eKin >= GetMinKinEnergy()
         && dp->GetDefinition() == G4Neutron::Neutron();
}

SetVerboseLevel()

void G4ParticleHPFissionData::SetVerboseLevel ( G4int newValue )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 248 of file G4ParticleHPFissionData.cc.

{
  G4ParticleHPManager::GetInstance()->SetVerboseLevel(newValue);
}

---

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

• G4ParticleHPFissionData.hh
• G4ParticleHPFissionData.cc