#include <G4ParticleHPFissionData.hh>

Inheritance diagram for G4ParticleHPFissionData:

Public Member Functions
	G4ParticleHPFissionData ()

	~G4ParticleHPFissionData ()

G4bool	IsIsoApplicable (const G4DynamicParticle , G4int, G4int, const G4Element , const G4Material *)

G4double	GetIsoCrossSection (const G4DynamicParticle , G4int, G4int, const G4Isotope , const G4Element , const G4Material )

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , G4double aT)

void	BuildPhysicsTable (const G4ParticleDefinition &)

void	DumpPhysicsTable (const G4ParticleDefinition &)

G4int	GetVerboseLevel () const

void	SetVerboseLevel (G4int)

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	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	GetIsoCrossSection (const G4DynamicParticle , G4int Z, G4int A, const G4Isotope iso=nullptr, const G4Element elm=nullptr, 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	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

virtual void	CrossSectionDescription (std::ostream &) 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

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 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 = 0;
   instanceOfWorker = false;
   if ( G4Threading::IsWorkerThread() ) {
      instanceOfWorker = true;
   }
   element_cache = nullptr;
   material_cache = nullptr;
   ke_cache = 0.0; 
   xs_cache = 0.0; 
   //BuildPhysicsTable(*G4Neutron::Neutron());
}

◆ ~G4ParticleHPFissionData()

G4ParticleHPFissionData::~G4ParticleHPFissionData ( )

Definition at line 67 of file G4ParticleHPFissionData.cc.

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

Member Function Documentation

◆ BuildPhysicsTable()

void G4ParticleHPFissionData::BuildPhysicsTable ( const G4ParticleDefinition & aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 105 of file G4ParticleHPFissionData.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()->GetFissionCrossSections();
      return;
   }
 
  std::size_t numberOfElements = G4Element::GetNumberOfElements();
  if ( theCrossSections == nullptr ) 
     theCrossSections = new G4PhysicsTable( numberOfElements );
  else
     theCrossSections->clearAndDestroy();
 
  // make a PhysicsVector for each element
 
  static G4ThreadLocal G4ElementTable *theElementTable  = nullptr ;
  if (!theElementTable) 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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 269 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 & aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 135 of file G4ParticleHPFissionData.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 << "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();
   static G4ThreadLocal G4ElementTable *theElementTable  = nullptr ;
   if (!theElementTable) 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 189 of file G4ParticleHPFissionData.cc.

{
  G4double result = 0;
  if(anE->GetZ()<88) return result;
  G4bool outOfRange;
  G4int 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) 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
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 89 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 259 of file G4ParticleHPFissionData.cc.

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

Referenced by DumpPhysicsTable().

◆ IsIsoApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 76 of file G4ParticleHPFissionData.cc.

{
   G4double eKin = dp->GetKineticEnergy();
   if ( eKin > GetMaxKinEnergy() 
     || eKin < GetMinKinEnergy() 
     || dp->GetDefinition() != G4Neutron::Neutron() ) return false;                                   
 
   return true;
}

◆ SetVerboseLevel()

void G4ParticleHPFissionData::SetVerboseLevel ( G4int newValue )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 264 of file G4ParticleHPFissionData.cc.

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

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4ParticleHPFissionData()

◆ ~G4ParticleHPFissionData()

Member Function Documentation

◆ BuildPhysicsTable()

◆ CrossSectionDescription()

◆ DumpPhysicsTable()

◆ GetCrossSection()

◆ GetIsoCrossSection()

◆ GetVerboseLevel()

◆ IsIsoApplicable()

◆ SetVerboseLevel()