#include <G4ParticleHPIsoProbabilityTable.hh>

Inheritance diagram for G4ParticleHPIsoProbabilityTable:

Public Member Functions
	G4ParticleHPIsoProbabilityTable ()=default

virtual	~G4ParticleHPIsoProbabilityTable ()

virtual void	Init (G4int, G4int, G4int, G4double, const G4String &)

virtual G4double	GetCorrelatedIsoCrossSectionPT (const G4DynamicParticle , G4int, const G4Element , G4double &, G4double &, std::thread::id &)

virtual G4double	GetIsoCrossSectionPT (const G4DynamicParticle , G4int, const G4Element , G4double &, std::map< std::thread::id, G4double > &, std::thread::id &)

Protected Member Functions
G4double	GetDopplerBroadenedElasticXS (const G4DynamicParticle *, G4int, G4int)

G4double	GetDopplerBroadenedCaptureXS (const G4DynamicParticle *, G4int, G4int)

G4double	GetDopplerBroadenedFissionXS (const G4DynamicParticle *, G4int, G4int)

G4double	GetDopplerBroadenedInelasticXS (const G4DynamicParticle *, G4int, G4int)

Protected Attributes
G4int	Z = 0

G4int	A = 0

G4int	m = -1

G4double	T = -1.

G4double	Emin = DBL_MAX

G4double	Emax = 0.

G4int	nEnergies = 0

std::map< std::thread::id, G4double >	energy_cache

std::map< std::thread::id, G4double >	xsela_cache

std::map< std::thread::id, G4double >	xscap_cache

std::map< std::thread::id, G4double >	xsfiss_cache

G4ParticleHPVector *	theEnergies = nullptr

std::vector< std::vector< G4double > * > *	theProbabilities = nullptr

std::vector< std::vector< G4double > * > *	theElasticData = nullptr

std::vector< std::vector< G4double > * > *	theCaptureData = nullptr

std::vector< std::vector< G4double > * > *	theFissionData = nullptr

std::vector< std::vector< G4double > * > *	theInelasticData = nullptr

G4String	filename

Detailed Description

Definition at line 61 of file G4ParticleHPIsoProbabilityTable.hh.

Constructor & Destructor Documentation

◆ G4ParticleHPIsoProbabilityTable()

G4ParticleHPIsoProbabilityTable::G4ParticleHPIsoProbabilityTable ( )

default

◆ ~G4ParticleHPIsoProbabilityTable()

G4ParticleHPIsoProbabilityTable::~G4ParticleHPIsoProbabilityTable ( )

virtual

Definition at line 68 of file G4ParticleHPIsoProbabilityTable.cc.

{
  for ( auto it = theProbabilities->cbegin(); it != theProbabilities->cend(); ++it ) {
    delete* it;
  }
  for ( auto it = theElasticData->cbegin(); it != theElasticData->cend(); ++it ) {
    delete* it;
  }
  for ( auto it = theCaptureData->cbegin(); it != theCaptureData->cend(); ++it ) {
    delete* it;
  }
  for ( auto it = theFissionData->cbegin(); it != theFissionData->cend(); ++it ) {
    delete* it;
  }
  delete theEnergies;
  delete theProbabilities;
  delete theElasticData;
  delete theCaptureData;
  delete theFissionData;
}

Member Function Documentation

◆ GetCorrelatedIsoCrossSectionPT()

G4double G4ParticleHPIsoProbabilityTable::GetCorrelatedIsoCrossSectionPT	(	const G4DynamicParticle *	,
		G4int	,
		const G4Element *	,
		G4double &	,
		G4double &	,
		std::thread::id &	)

virtual

Reimplemented in G4ParticleHPIsoProbabilityTable_CALENDF, and G4ParticleHPIsoProbabilityTable_NJOY.

Definition at line 93 of file G4ParticleHPIsoProbabilityTable.cc.

{
  G4Exception( "G4ParticleHPIsoProbabilityTable::GetCorrelatedIsoCrossSectionPT", "hadhp01",
           FatalException, "The base method G4ParticleHPIsoProbabilityTable::GetCorrelatedIsoCrossSectionPT"
           "is trying to be accessed, which is not allowed,"
           "the derived class for NJOY or CALENDF was not properly declared." );
  return 0.0;
}

◆ GetDopplerBroadenedCaptureXS()

G4double G4ParticleHPIsoProbabilityTable::GetDopplerBroadenedCaptureXS	(	const G4DynamicParticle *	dP,
		G4int	indexEl,
		G4int	isotopeJ )

protected

Definition at line 150 of file G4ParticleHPIsoProbabilityTable.cc.

                                                                                                                                   {
  G4double result = 0.0;
  G4ReactionProduct theNeutron( dP->GetDefinition() );
  theNeutron.SetMomentum( dP->GetMomentum() );
  theNeutron.SetKineticEnergy( dP->GetKineticEnergy() );
  // prepare thermal nucleus
  G4Nucleus aNuc;
  G4double eleMass = G4NucleiProperties::GetNuclearMass(A, Z) / G4Neutron::Neutron()->GetPDGMass();
  G4ReactionProduct boosted;
  G4double aXsection;
  G4int counter = 0;
  G4double buffer = 0.0;
  G4int size = G4int( std::max( 10.0, T / 60.0 * kelvin));
  G4ThreeVector neutronVelocity = 1.0 / G4Neutron::Neutron()->GetPDGMass() * theNeutron.GetMomentum();
  G4double neutronVMag = neutronVelocity.mag();
  while ( counter == 0  ||  std::abs( buffer - result / std::max( 1, counter ) ) > 0.03 * buffer ) {
    if ( counter ) buffer = result / counter;
    while ( counter < size ) {
      ++counter;
      G4ReactionProduct aThermalNuc = aNuc.GetThermalNucleus( eleMass, T );
      boosted.Lorentz( theNeutron, aThermalNuc );
      G4double theEkin = boosted.GetKineticEnergy();
      aXsection = (*G4ParticleHPManager::GetInstance()->GetCaptureFinalStates())[indexEl]->GetWeightedXsec( theEkin, isotopeJ );
      // velocity correction
      G4ThreeVector targetVelocity = 1.0 / aThermalNuc.GetMass() * aThermalNuc.GetMomentum();
      aXsection *= (targetVelocity - neutronVelocity).mag() / neutronVMag;
      result += aXsection;
    }
    size += size;
  }
  result /= counter;
  return result;
}

Referenced by G4ParticleHPIsoProbabilityTable_CALENDF::GetCorrelatedIsoCrossSectionPT(), G4ParticleHPIsoProbabilityTable_NJOY::GetCorrelatedIsoCrossSectionPT(), G4ParticleHPIsoProbabilityTable_CALENDF::GetIsoCrossSectionPT(), and G4ParticleHPIsoProbabilityTable_NJOY::GetIsoCrossSectionPT().

◆ GetDopplerBroadenedElasticXS()

G4double G4ParticleHPIsoProbabilityTable::GetDopplerBroadenedElasticXS	(	const G4DynamicParticle *	dP,
		G4int	indexEl,
		G4int	isotopeJ )

protected

Definition at line 115 of file G4ParticleHPIsoProbabilityTable.cc.

                                                                                                                                   {
  G4double result = 0.0;
  G4ReactionProduct theNeutron( dP->GetDefinition() );
  theNeutron.SetMomentum( dP->GetMomentum() );
  theNeutron.SetKineticEnergy( dP->GetKineticEnergy() );
  // prepare thermal nucleus
  G4Nucleus aNuc;
  G4double eleMass = G4NucleiProperties::GetNuclearMass(A, Z) / G4Neutron::Neutron()->GetPDGMass();
  G4ReactionProduct boosted;
  G4double aXsection;
  G4int counter = 0;
  G4double buffer = 0.0;
  G4int size = G4int( std::max( 10.0, T / 60.0 * kelvin ) );
  G4ThreeVector neutronVelocity = 1.0 / G4Neutron::Neutron()->GetPDGMass() * theNeutron.GetMomentum();
  G4double neutronVMag = neutronVelocity.mag();
  while ( counter == 0  ||  std::abs( buffer - result / std::max(1, counter) ) > 0.03 * buffer ) {
    if ( counter ) buffer = result / counter;
    while ( counter < size ) {
      ++counter;
      G4ReactionProduct aThermalNuc = aNuc.GetThermalNucleus( eleMass, T );
      boosted.Lorentz( theNeutron, aThermalNuc );
      G4double theEkin = boosted.GetKineticEnergy();
      aXsection = (*G4ParticleHPManager::GetInstance()->GetElasticFinalStates())[indexEl]->GetWeightedXsec( theEkin, isotopeJ );
      // velocity correction
      G4ThreeVector targetVelocity = 1.0 / aThermalNuc.GetMass() * aThermalNuc.GetMomentum();
      aXsection *= (targetVelocity - neutronVelocity).mag() / neutronVMag;
      result += aXsection;
    }
    size += size;
  }
  result /= counter;
  return result;
}

Referenced by G4ParticleHPIsoProbabilityTable_CALENDF::GetCorrelatedIsoCrossSectionPT(), G4ParticleHPIsoProbabilityTable_NJOY::GetCorrelatedIsoCrossSectionPT(), G4ParticleHPIsoProbabilityTable_CALENDF::GetIsoCrossSectionPT(), and G4ParticleHPIsoProbabilityTable_NJOY::GetIsoCrossSectionPT().

◆ GetDopplerBroadenedFissionXS()

G4double G4ParticleHPIsoProbabilityTable::GetDopplerBroadenedFissionXS	(	const G4DynamicParticle *	dP,
		G4int	indexEl,
		G4int	isotopeJ )

protected

Definition at line 185 of file G4ParticleHPIsoProbabilityTable.cc.

                                                                                                                                   {
  G4double result = 0.0;
  G4ReactionProduct theNeutron( dP->GetDefinition() );
  theNeutron.SetMomentum( dP->GetMomentum() );
  theNeutron.SetKineticEnergy( dP->GetKineticEnergy() );
  // prepare thermal nucleus
  G4Nucleus aNuc;
  G4double eleMass;
  eleMass = G4NucleiProperties::GetNuclearMass(A, Z) / G4Neutron::Neutron()->GetPDGMass();
  G4ReactionProduct boosted;
  G4double aXsection;
  G4int counter = 0;
  G4double buffer = 0.0;
  G4int size = G4int( std::max( 10.0, T / 60.0 * kelvin ) );
  G4ThreeVector neutronVelocity = 1.0 / G4Neutron::Neutron()->GetPDGMass() * theNeutron.GetMomentum();
  G4double neutronVMag = neutronVelocity.mag();
  while ( counter == 0  ||  std::abs( buffer - result / std::max( 1, counter ) ) > 0.01 * buffer ) {
    if ( counter ) buffer = result / counter;
    while ( counter < size ) {
      ++counter;
      G4ReactionProduct aThermalNuc = aNuc.GetThermalNucleus( eleMass, T );
      boosted.Lorentz( theNeutron, aThermalNuc );
      G4double theEkin = boosted.GetKineticEnergy();
      aXsection = (*G4ParticleHPManager::GetInstance()->GetFissionFinalStates())[indexEl]->GetWeightedXsec( theEkin, isotopeJ );
      // velocity correction
      G4ThreeVector targetVelocity = 1.0 / aThermalNuc.GetMass() * aThermalNuc.GetMomentum();
      aXsection *= (targetVelocity - neutronVelocity).mag() / neutronVMag;
      result += aXsection;
    }
    size += size;
  }
  result /= counter;
  return result;
}

Referenced by G4ParticleHPIsoProbabilityTable_CALENDF::GetCorrelatedIsoCrossSectionPT(), G4ParticleHPIsoProbabilityTable_NJOY::GetCorrelatedIsoCrossSectionPT(), G4ParticleHPIsoProbabilityTable_CALENDF::GetIsoCrossSectionPT(), and G4ParticleHPIsoProbabilityTable_NJOY::GetIsoCrossSectionPT().

◆ GetDopplerBroadenedInelasticXS()

G4double G4ParticleHPIsoProbabilityTable::GetDopplerBroadenedInelasticXS	(	const G4DynamicParticle *	dP,
		G4int	indexEl,
		G4int	isotopeJ )

protected

Definition at line 221 of file G4ParticleHPIsoProbabilityTable.cc.

                                                                                                                                     {
  G4double result = 0.0;
  G4ReactionProduct theNeutron( dP->GetDefinition() );
  theNeutron.SetMomentum( dP->GetMomentum() );
  theNeutron.SetKineticEnergy( dP->GetKineticEnergy() );
  // prepare thermal nucleus
  G4Nucleus aNuc;
  G4double eleMass = G4NucleiProperties::GetNuclearMass(A, Z) / G4Neutron::Neutron()->GetPDGMass();
  G4ReactionProduct boosted;
  G4double aXsection;
  G4int counter = 0;
  G4int failCount = 0;
  G4double buffer = 0.0;
  G4int size = G4int( std::max( 10.0, T / 60.0 * kelvin ) );
  G4ThreeVector neutronVelocity = 1.0 / G4Neutron::Neutron()->GetPDGMass() * theNeutron.GetMomentum();
  G4double neutronVMag = neutronVelocity.mag();
  #ifndef G4PHP_DOPPLER_LOOP_ONCE
  while ( counter == 0  ||  std::abs( buffer - result / std::max( 1, counter ) ) > 0.01 * buffer ) {
    if ( counter ) buffer = result / counter;
    while ( counter < size ) {
      ++counter;
  #endif
      G4ReactionProduct aThermalNuc = aNuc.GetThermalNucleus( eleMass / G4Neutron::Neutron()->GetPDGMass(), T );
      boosted.Lorentz( theNeutron, aThermalNuc );
      G4double theEkin = boosted.GetKineticEnergy();
      aXsection = ((*G4ParticleHPManager::GetInstance()->GetInelasticFinalStates( dP->GetDefinition() ))[indexEl]
                  ->GetWeightedXsec( theEkin, isotopeJ ) ) * barn;
      if ( aXsection < 0.0 ) {
    if ( failCount < 1000 ) {
      ++failCount;
          #ifndef G4PHP_DOPPLER_LOOP_ONCE
      --counter;
      continue;
          #endif
    } else {
      aXsection = 0.0;
    }
      }
      // velocity correction
      G4ThreeVector targetVelocity = 1.0 / aThermalNuc.GetMass() * aThermalNuc.GetMomentum();
      aXsection *= (targetVelocity - neutronVelocity).mag() / neutronVMag;
      result += aXsection;
    }
#ifndef G4PHP_DOPPLER_LOOP_ONCE
    size += size;
  }
  result /= counter;
#endif
  return result;
}

Referenced by G4ParticleHPIsoProbabilityTable_CALENDF::GetCorrelatedIsoCrossSectionPT(), and G4ParticleHPIsoProbabilityTable_CALENDF::GetIsoCrossSectionPT().

◆ GetIsoCrossSectionPT()

G4double G4ParticleHPIsoProbabilityTable::GetIsoCrossSectionPT	(	const G4DynamicParticle *	,
		G4int	,
		const G4Element *	,
		G4double &	,
		std::map< std::thread::id, G4double > &	,
		std::thread::id &	)

virtual

Reimplemented in G4ParticleHPIsoProbabilityTable_CALENDF, and G4ParticleHPIsoProbabilityTable_NJOY.

Definition at line 104 of file G4ParticleHPIsoProbabilityTable.cc.

{
  G4Exception( "G4ParticleHPIsoProbabilityTable::GetIsoCrossSectionPT", "hadhp01",
           FatalException, "The base method G4ParticleHPIsoProbabilityTable::GetIsoCrossSectionPT"
           "is trying to be accessed, which is not allowed,"
           "the derived class for NJOY or CALENDF was not properly declared." );
  return 0.0;
}