G4ParticleHPCaptureData Class Reference

#include <G4ParticleHPCaptureData.hh>

Inheritance diagram for G4ParticleHPCaptureData:

Public Member Functions
	G4ParticleHPCaptureData ()
	~G4ParticleHPCaptureData ()
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 &)
void	IgnoreOnFlightDopplerBroadening ()
void	EnableOnFlightDopplerBroadening ()
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	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 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 G4int	GetVerboseLevel () 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

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Detailed Description

Definition at line 51 of file G4ParticleHPCaptureData.hh.

Constructor & Destructor Documentation

G4ParticleHPCaptureData()

G4ParticleHPCaptureData::G4ParticleHPCaptureData

(

)

Definition at line 51 of file G4ParticleHPCaptureData.cc.

:G4VCrossSectionDataSet("NeutronHPCaptureXS")
{
   SetMinKinEnergy( 0*MeV );                                   
   SetMaxKinEnergy( 20*MeV );                                   
 
   theCrossSections = 0;
   onFlightDB = true;
 
   instanceOfWorker = false;
   if ( G4Threading::IsWorkerThread() ) {
      instanceOfWorker = true;
   }
 
   element_cache = NULL;
   material_cache = NULL;
   ke_cache = 0.0; 
   xs_cache = 0.0; 
    
   //BuildPhysicsTable(*G4Neutron::Neutron());
}

~G4ParticleHPCaptureData()

G4ParticleHPCaptureData::~G4ParticleHPCaptureData

(

)

Definition at line 73 of file G4ParticleHPCaptureData.cc.

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

Member Function Documentation

BuildPhysicsTable()

void G4ParticleHPCaptureData::BuildPhysicsTable ( const G4ParticleDefinition &  aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 121 of file G4ParticleHPCaptureData.cc.

{
  if(&aP!=G4Neutron::Neutron()) 
     throw G4HadronicException(__FILE__, __LINE__, "Attempt to use NeutronHP data for particles other than neutrons!!!");  
 
//080428
   if ( G4ParticleHPManager::GetInstance()->GetNeglectDoppler() )
   {
      onFlightDB = false;
      #ifdef G4VERBOSE
      if ( G4HadronicParameters::Instance()->GetVerboseLevel() > 0 ) {
    G4cout << "Find a flag of \"G4NEUTRONHP_NEGLECT_DOPPLER\"." << G4endl;
    G4cout << "On the fly Doppler broadening will be neglect in the cross section calculation of capture reaction of neutrons (<20MeV)." << G4endl;
      }
      #endif
   }
 
   if ( G4Threading::IsWorkerThread() ) {
      theCrossSections = G4ParticleHPManager::GetInstance()->GetCaptureCrossSections();
      return;
   }
  
  size_t numberOfElements = G4Element::GetNumberOfElements();
  // G4cout << "CALLED G4ParticleHPCaptureData::BuildPhysicsTable "<<numberOfElements<<G4endl;
   // TKDB
   //if ( theCrossSections == 0 ) theCrossSections = new G4PhysicsTable( numberOfElements );
   if ( theCrossSections == NULL ) 
      theCrossSections = new G4PhysicsTable( numberOfElements );
   else
      theCrossSections->clearAndDestroy();
 
  // make a PhysicsVector for each element
 
  static G4ThreadLocal G4ElementTable *theElementTable  = 0 ; if (!theElementTable) theElementTable= G4Element::GetElementTable();
  for( size_t i=0; i<numberOfElements; ++i )
  {
     #ifdef G4VERBOSE
     if(std::getenv("CaptureDataIndexDebug"))
     {
       G4int index_debug = ((*theElementTable)[i])->GetIndex();
       if ( G4HadronicParameters::Instance()->GetVerboseLevel() > 0 ) G4cout << "IndexDebug "<< i <<" "<<index_debug<<G4endl;
     }
     #endif
     G4PhysicsVector* physVec = G4ParticleHPData::
      Instance(G4Neutron::Neutron())->MakePhysicsVector((*theElementTable)[i], this);
     theCrossSections->push_back(physVec);
  }
 
  G4ParticleHPManager::GetInstance()->RegisterCaptureCrossSections( theCrossSections );
}

CrossSectionDescription()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 308 of file G4ParticleHPCaptureData.cc.

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

DumpPhysicsTable()

void G4ParticleHPCaptureData::DumpPhysicsTable ( const G4ParticleDefinition &  aP )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 172 of file G4ParticleHPCaptureData.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 << "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;
 
   size_t numberOfElements = G4Element::GetNumberOfElements();
   static G4ThreadLocal G4ElementTable *theElementTable  = 0 ; if (!theElementTable) theElementTable= G4Element::GetElementTable();
 
   for ( 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;
   }
 
   //G4cout << "G4ParticleHPCaptureData::DumpPhysicsTable still to be implemented"<<G4endl;
   #endif
}

EnableOnFlightDopplerBroadening()

void G4ParticleHPCaptureData::EnableOnFlightDopplerBroadening ( )

inline

Definition at line 83 of file G4ParticleHPCaptureData.hh.

{ onFlightDB = true; };

GetCrossSection()

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

Definition at line 227 of file G4ParticleHPCaptureData.cc.

{
  G4double result = 0;
  G4bool outOfRange;
  G4int index = anE->GetIndex();
 
  // prepare neutron
  G4double eKinetic = aP->GetKineticEnergy();
 
  if ( !onFlightDB )
  {
     //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; 
  }
 
  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( 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()*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) 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, 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))).GetValue(eKinetic, outOfRange) << " " 
         << (*((*theCrossSections)(index))).GetValue(eKinetic, outOfRange) /result << G4endl;
*/
  return result;
}

Referenced by GetIsoCrossSection().

GetIsoCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 95 of file G4ParticleHPCaptureData.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 G4ParticleHPCaptureData::GetVerboseLevel ( ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 300 of file G4ParticleHPCaptureData.cc.

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

Referenced by BuildPhysicsTable(), and DumpPhysicsTable().

IgnoreOnFlightDopplerBroadening()

void G4ParticleHPCaptureData::IgnoreOnFlightDopplerBroadening ( )

inline

Definition at line 82 of file G4ParticleHPCaptureData.hh.

{ onFlightDB = false; };

IsIsoApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 82 of file G4ParticleHPCaptureData.cc.

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

SetVerboseLevel()

void G4ParticleHPCaptureData::SetVerboseLevel ( G4int  newValue )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 304 of file G4ParticleHPCaptureData.cc.

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

---

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

• G4ParticleHPCaptureData.hh
• G4ParticleHPCaptureData.cc