G4ParticleHPElasticData Class Reference

#include <G4ParticleHPElasticData.hh>

Inheritance diagram for G4ParticleHPElasticData:

Public Member Functions
	G4ParticleHPElasticData ()
	~G4ParticleHPElasticData () 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
void	SetVerboseLevel (G4int) override
G4int	GetVerboseLevel () const
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
G4String	name

Detailed Description

Definition at line 51 of file G4ParticleHPElasticData.hh.

Constructor & Destructor Documentation

G4ParticleHPElasticData()

G4ParticleHPElasticData::G4ParticleHPElasticData

(

)

Definition at line 53 of file G4ParticleHPElasticData.cc.

                                                 : G4VCrossSectionDataSet("NeutronHPElasticXS")
{
  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;
  // BuildPhysicsTable( *G4Neutron::Neutron() );
}

~G4ParticleHPElasticData()

G4ParticleHPElasticData::~G4ParticleHPElasticData ( )

override

Definition at line 70 of file G4ParticleHPElasticData.cc.

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

Member Function Documentation

BuildPhysicsTable()

void G4ParticleHPElasticData::BuildPhysicsTable ( const G4ParticleDefinition & aP )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 114 of file G4ParticleHPElasticData.cc.

{
  if (&aP != G4Neutron::Neutron())
    throw G4HadronicException(__FILE__, __LINE__,
                              "Attempt to use NeutronHP data for particles other than neutrons!!!");
 
  if (G4Threading::IsWorkerThread()) {
    theCrossSections = G4ParticleHPManager::GetInstance()->GetElasticCrossSections();
    return;
  }
 
  std::size_t numberOfElements = G4Element::GetNumberOfElements();
  // TKDB
  // if ( theCrossSections == 0 ) theCrossSections = new G4PhysicsTable( numberOfElements );
  if (theCrossSections == nullptr)
    theCrossSections = new G4PhysicsTable(numberOfElements);
  else
    theCrossSections->clearAndDestroy();
 
  // make a PhysicsVector for each element
 
  static G4ThreadLocal G4ElementTable* theElementTable = nullptr;
  if (theElementTable == nullptr) 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()->RegisterElasticCrossSections(theCrossSections);
}

CrossSectionDescription()

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

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 279 of file G4ParticleHPElasticData.cc.

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

DumpPhysicsTable()

void G4ParticleHPElasticData::DumpPhysicsTable ( const G4ParticleDefinition & aP )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 146 of file G4ParticleHPElasticData.cc.

{
  if (&aP != G4Neutron::Neutron())
    throw G4HadronicException(__FILE__, __LINE__,
                              "Attempt to use NeutronHP data for particles other than neutrons!!!");
 
#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 << "Elastic 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();
  static G4ThreadLocal G4ElementTable* theElementTable = nullptr;
  if (theElementTable == nullptr) theElementTable = G4Element::GetElementTable();
 
  for (std::size_t i = 0; i < numberOfElements; ++i) {
    G4cout << (*theElementTable)[i]->GetName() << G4endl;
    G4int ie = 0;
 
    for (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 G4ParticleHPElasticData::GetCrossSection	(	const G4DynamicParticle *	aP,
		const G4Element *	anE,
		G4double	aT )

Definition at line 193 of file G4ParticleHPElasticData.cc.

{
  G4double result = 0;
  G4bool outOfRange;
  auto index = (G4int)anE->GetIndex();
 
  // prepare neutron
  G4double eKinetic = aP->GetKineticEnergy();
 
  if (G4ParticleHPManager::GetInstance()->GetNeglectDoppler()) {
    // NEGLECT_DOPPLER_B.
    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;
  }
 
  G4ReactionProduct theNeutron(aP->GetDefinition());
  theNeutron.SetMomentum(aP->GetMomentum());
  theNeutron.SetKineticEnergy(eKinetic);
 
  // prepare thermal nucleus
  G4Nucleus aNuc;
  G4double eps = 0.0001;
  G4double theA = anE->GetN();
  G4double theZ = anE->GetZ();
  G4double eleMass;
 
  eleMass = (G4NucleiProperties::GetNuclearMass(G4int(theA + eps), 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() * 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)(index))).GetValue(theEkin, outOfRange);
      // velocity correction.
      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))).GetValue(eKinetic, outOfRange) << " "
           << (*((*theCrossSections)(index))).GetValue(eKinetic, outOfRange) /result << G4endl;
  */
  return result;
}

Referenced by GetIsoCrossSection().

GetIsoCrossSection()

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

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 88 of file G4ParticleHPElasticData.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 G4ParticleHPElasticData::GetVerboseLevel

(

)

const

Definition at line 269 of file G4ParticleHPElasticData.cc.

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

Referenced by DumpPhysicsTable().

IsIsoApplicable()

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

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 79 of file G4ParticleHPElasticData.cc.

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

SetVerboseLevel()

void G4ParticleHPElasticData::SetVerboseLevel ( G4int newValue )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 274 of file G4ParticleHPElasticData.cc.

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

---

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

• G4ParticleHPElasticData.hh
• G4ParticleHPElasticData.cc