G4ParticleHPInelasticData Class Reference

#include <G4ParticleHPInelasticData.hh>

Inheritance diagram for G4ParticleHPInelasticData:

Public Member Functions
	G4ParticleHPInelasticData (G4ParticleDefinition *projectile=G4Neutron::Neutron())
	~G4ParticleHPInelasticData ()
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 &)
G4ParticleDefinition *	GetProjectile ()
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 53 of file G4ParticleHPInelasticData.hh.

Constructor & Destructor Documentation

G4ParticleHPInelasticData()

G4ParticleHPInelasticData::G4ParticleHPInelasticData

(

G4ParticleDefinition *

projectile = G4Neutron::Neutron()

)

Definition at line 48 of file G4ParticleHPInelasticData.cc.

  : G4VCrossSectionDataSet("")
{
  const char* dataDirVariable;
  G4String particleName;
  if( projectile == G4Neutron::Neutron() ) {
    dataDirVariable = "G4NEUTRONHPDATA";
  }else if( projectile == G4Proton::Proton() ) {
    dataDirVariable = "G4PROTONHPDATA";
    particleName = "Proton";
  }else if( projectile == G4Deuteron::Deuteron() ) {
    dataDirVariable = "G4DEUTERONHPDATA";
    particleName = "Deuteron";
  }else if( projectile == G4Triton::Triton() ) {
    dataDirVariable = "G4TRITONHPDATA";
    particleName = "Triton";
  }else if( projectile == G4He3::He3() ) {
    dataDirVariable = "G4HE3HPDATA";
    particleName = "He3";
  }else if( projectile == G4Alpha::Alpha() ) {
    dataDirVariable = "G4ALPHAHPDATA";
    particleName = "Alpha";
  } else {
    G4String message("G4ParticleHPInelasticData may only be called for neutron, proton, deuteron, triton, He3 or alpha, while it is called for " + projectile->GetParticleName());
    throw G4HadronicException(__FILE__, __LINE__,message.c_str());
  }
  G4String dataName = projectile->GetParticleName()+"HPInelasticXS";
  dataName.at(0) = (char)std::toupper(dataName.at(0)) ;
  SetName( dataName );
 
  if ( !G4FindDataDir(dataDirVariable) && !G4FindDataDir( "G4PARTICLEHPDATA" ) ){
    G4String message("Please setenv G4PARTICLEHPDATA (recommended) or, at least setenv " +
                     G4String(dataDirVariable) + " to point to the " + projectile->GetParticleName() + " cross-section files.");
    throw G4HadronicException(__FILE__, __LINE__,message.c_str());
  }
 
  G4String dirName;
  if ( G4FindDataDir(dataDirVariable) ) {
    dirName = G4FindDataDir(dataDirVariable);
  } else {
    G4String baseName = G4FindDataDir( "G4PARTICLEHPDATA" );
    dirName = baseName + "/" + particleName;
  }
  #ifdef G4VERBOSE
  if ( G4HadronicParameters::Instance()->GetVerboseLevel() > 0 ) {
    G4cout << "@@@ G4ParticleHPInelasticData instantiated for particle " << projectile->GetParticleName() << " data directory variable is " << dataDirVariable << " pointing to " << dirName << G4endl;
  }
  #endif
 
  SetMinKinEnergy( 0*CLHEP::MeV );                                   
  SetMaxKinEnergy( 20*CLHEP::MeV );                                   
 
  theCrossSections = 0;
  theProjectile=projectile;
 
  theHPData = nullptr;
  instanceOfWorker = false;
  if ( G4Threading::IsMasterThread() ) {
    theHPData = new G4ParticleHPData( theProjectile ); 
  } else {
    instanceOfWorker = true;
  }
  element_cache = nullptr;
  material_cache = nullptr;
  ke_cache = 0.0; 
  xs_cache = 0.0; 
}

~G4ParticleHPInelasticData()

G4ParticleHPInelasticData::~G4ParticleHPInelasticData

(

)

Definition at line 116 of file G4ParticleHPInelasticData.cc.

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

Member Function Documentation

BuildPhysicsTable()

void G4ParticleHPInelasticData::BuildPhysicsTable ( const G4ParticleDefinition &  projectile )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 159 of file G4ParticleHPInelasticData.cc.

{
   if ( G4Threading::IsWorkerThread() ) {
      theCrossSections = G4ParticleHPManager::GetInstance()->GetInelasticCrossSections( &projectile );
      return;
   } else {
      if ( theHPData == nullptr ) theHPData = G4ParticleHPData::Instance( const_cast<G4ParticleDefinition*> ( &projectile ) ); 
   }
 
  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 = theHPData->MakePhysicsVector((*theElementTable)[i], this);
    theCrossSections->push_back(physVec);
  }
  G4ParticleHPManager::GetInstance()->RegisterInelasticCrossSections( &projectile , theCrossSections );
}

CrossSectionDescription()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 332 of file G4ParticleHPInelasticData.cc.

{
   outFile << "Extension of High Precision cross section for inelastic reaction of proton, deuteron, triton, He3 and alpha below 20MeV\n";
}

DumpPhysicsTable()

void G4ParticleHPInelasticData::DumpPhysicsTable ( const G4ParticleDefinition &  projectile )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 186 of file G4ParticleHPInelasticData.cc.

{
  if(&projectile!=theProjectile) 
     throw G4HadronicException(__FILE__, __LINE__, "Attempt to use ParticleHP data for a wrong projectile!!!");  
 
#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 << "Inelastic 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;
 
     G4int ie = 0;
 
     for ( ie = 0 ; ie < 130 ; ie++ )
     {
       G4double eKinetic = 1.0e-5 * G4Pow::GetInstance()->powA ( 10.0 , ie/10.0 ) *CLHEP::eV;
       G4bool outOfRange = false;
 
       if ( eKinetic < 20*CLHEP::MeV )
       {
         G4cout << eKinetic/CLHEP::eV << " "
                << (*((*theCrossSections)(i))).GetValue(eKinetic, outOfRange)/CLHEP::barn
                << G4endl;
       }
     }
     G4cout << G4endl;
  }
#endif
}

GetCrossSection()

G4double G4ParticleHPInelasticData::GetCrossSection	(	const G4DynamicParticle *	projectile,
		const G4Element *	anE,
		G4double	aT
	)

Definition at line 235 of file G4ParticleHPInelasticData.cc.

{
  G4double result = 0;
  G4bool outOfRange;
  std::size_t index = anE->GetIndex();
 
  // prepare neutron
  G4double eKinetic = projectile->GetKineticEnergy();
 
  if ( G4ParticleHPManager::GetInstance()->GetNeglectDoppler() )
  {
     //NEGLECT_DOPPLER
     G4double factor = 1.0;
     if ( eKinetic < aT * CLHEP::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( projectile->GetDefinition() );
  theNeutron.SetMomentum( projectile->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) );
  
  G4ReactionProduct boosted;
  G4double aXsection;
  
  // MC integration loop
  G4int counter = 0;
  G4int failCount = 0;
  G4double buffer = 0; G4int size = G4int(std::max(10., aT/60*CLHEP::kelvin));
  G4ThreeVector neutronVelocity = 1./theProjectile->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) // Loop checking, 11.05.2015, T. Koi
  {
    if(counter) buffer = result/counter;
    while (counter<size) // Loop checking, 11.05.2015, T. Koi
    {
      ++counter;
#endif
      G4ReactionProduct aThermalNuc = aNuc.GetThermalNucleus( eleMass/G4Neutron::Neutron()->GetPDGMass(), aT );
      boosted.Lorentz(theNeutron, aThermalNuc);
      G4double theEkin = boosted.GetKineticEnergy();
      aXsection = (*((*theCrossSections)(index))).GetValue(theEkin, outOfRange);
     if(aXsection <0) 
      {
        if(failCount<1000)
    {
      ++failCount;
#ifndef G4PHP_DOPPLER_LOOP_ONCE
      --counter;
      continue;
#endif
    }
    else
    {
      aXsection = 0;
    }
      }
      // velocity correction.
      G4ThreeVector targetVelocity = 1./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 GetIsoCrossSection().

GetIsoCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

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

GetProjectile()

G4ParticleDefinition * G4ParticleHPInelasticData::GetProjectile ( )

inline

Definition at line 85 of file G4ParticleHPInelasticData.hh.

{return theProjectile;}

GetVerboseLevel()

G4int G4ParticleHPInelasticData::GetVerboseLevel

(

)

const

Definition at line 322 of file G4ParticleHPInelasticData.cc.

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

Referenced by DumpPhysicsTable(), and G4ParticleHPInelasticData().

IsIsoApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 129 of file G4ParticleHPInelasticData.cc.

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

SetVerboseLevel()

void G4ParticleHPInelasticData::SetVerboseLevel ( G4int  newValue )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 327 of file G4ParticleHPInelasticData.cc.

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

---

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

• G4ParticleHPInelasticData.hh
• G4ParticleHPInelasticData.cc