G4ParticleHPJENDLHEData Class Reference

#include <G4ParticleHPJENDLHEData.hh>

Inheritance diagram for G4ParticleHPJENDLHEData:

Public Member Functions
	G4ParticleHPJENDLHEData ()
	G4ParticleHPJENDLHEData (G4String, G4ParticleDefinition *)
	~G4ParticleHPJENDLHEData () override
G4bool	IsApplicable (const G4DynamicParticle *, const G4Element *)
G4bool	IsZAApplicable (const G4DynamicParticle *, G4double, G4double)
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, G4double aT)
void	BuildPhysicsTable (const G4ParticleDefinition &) override
void	DumpPhysicsTable (const G4ParticleDefinition &) override
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	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 49 of file G4ParticleHPJENDLHEData.hh.

Constructor & Destructor Documentation

G4ParticleHPJENDLHEData() [1/2]

G4ParticleHPJENDLHEData::G4ParticleHPJENDLHEData

(

)

Definition at line 57 of file G4ParticleHPJENDLHEData.cc.

{
  for (auto& itZ : mIsotope) {
    std::map<G4int, G4PhysicsVector*>* pointer_map = itZ.second;
    if (pointer_map != nullptr) {
      for (auto& itA : *pointer_map) {
        G4PhysicsVector* pointerPhysicsVector = itA.second;
        if (pointerPhysicsVector != nullptr) {
          delete pointerPhysicsVector;
          itA.second = NULL;
        }
      }
      delete pointer_map;
      itZ.second = NULL;
    }
  }
  mIsotope.clear();
}

G4ParticleHPJENDLHEData() [2/2]

G4ParticleHPJENDLHEData::G4ParticleHPJENDLHEData	(	G4String	reaction,
		G4ParticleDefinition *	pd )

Definition at line 76 of file G4ParticleHPJENDLHEData.cc.

  : G4VCrossSectionDataSet("JENDLHE" + reaction + "CrossSection")
{
  reactionName = reaction;
  BuildPhysicsTable(*pd);
}

~G4ParticleHPJENDLHEData()

G4ParticleHPJENDLHEData::~G4ParticleHPJENDLHEData ( )

overridedefault

Member Function Documentation

BuildPhysicsTable()

void G4ParticleHPJENDLHEData::BuildPhysicsTable ( const G4ParticleDefinition & aP )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 85 of file G4ParticleHPJENDLHEData.cc.

{
  particleName = aP.GetParticleName();
 
  G4String baseName = G4FindDataDir("G4NEUTRONHPDATA");
  G4String dirName = baseName + "/JENDL_HE/" + particleName + "/" + reactionName;
  G4String aFSType = "/CrossSection/";
  G4ParticleHPNames theNames;
 
  G4String filename;
 
  // Create JENDL_HE data
  // Create map element or isotope
 
  std::size_t numberOfElements = G4Element::GetNumberOfElements();
 
  // make a PhysicsVector for each element
 
  auto theElementTable = G4Element::GetElementTable();
  vElement.clear();
  vElement.resize(numberOfElements);
  for (std::size_t i = 0; i < numberOfElements; ++i) {
    G4Element* theElement = (*theElementTable)[i];
    vElement[i] = false;
 
    // isotope
    auto nIso = (G4int)(*theElementTable)[i]->GetNumberOfIsotopes();
    auto Z = (G4int)(*theElementTable)[i]->GetZ();
    for (G4int i1 = 0; i1 < nIso; ++i1) {
      G4int A = theElement->GetIsotope(i1)->GetN();
 
      if (isThisNewIsotope(Z, A)) {
    std::stringstream ss;
    ss << dirName << aFSType << Z << "_" << A << "_" << theNames.GetName(Z - 1);
    filename = ss.str();
    std::fstream file;
    file.open(filename, std::fstream::in);
    G4int dummy;
    file >> dummy;
    if (file.good()) {
      vElement[i] = true;
 
      // read the file
      G4PhysicsVector* aPhysVec = readAFile(&file);
      registAPhysicsVector(Z, A, aPhysVec);
    }
    file.close();
      }
    }
  }
}

Referenced by G4ParticleHPJENDLHEData().

DumpPhysicsTable()

void G4ParticleHPJENDLHEData::DumpPhysicsTable ( const G4ParticleDefinition & )

overridevirtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 137 of file G4ParticleHPJENDLHEData.cc.

{}

GetCrossSection()

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

Definition at line 140 of file G4ParticleHPJENDLHEData.cc.

{
  // Primary energy >20MeV
  // Thus not taking into account of Doppler broadening
  // also not taking into account of Target thermal motions
 
  G4double result = 0;
 
  G4double ek = aP->GetKineticEnergy();
 
  auto nIso = (G4int)anE->GetNumberOfIsotopes();
  auto Z = (G4int)anE->GetZ();
  for (G4int i1 = 0; i1 < nIso; ++i1) {
    G4int A = anE->GetIsotope(i1)->GetN();
    G4double frac = anE->GetRelativeAbundanceVector()[i1];
    // This case does NOT request "*perCent".
    result += frac * getXSfromThisIsotope(Z, A, ek);
  }
  return result;
}

IsApplicable()

G4bool G4ParticleHPJENDLHEData::IsApplicable	(	const G4DynamicParticle *	aP,
		const G4Element *	anE )

Definition at line 42 of file G4ParticleHPJENDLHEData.cc.

{
  G4bool result = true;
  G4double eKin = aP->GetKineticEnergy();
  // if(eKin>20*MeV||aP->GetDefinition()!=G4Neutron::Neutron()) result = false;
  if (eKin < 20 * MeV || 3 * GeV < eKin || aP->GetDefinition() != G4Neutron::Neutron()) {
    result = false;
  }
  // Element Check
  else if (!(vElement[anE->GetIndex()]))
    result = false;
 
  return result;
}

IsZAApplicable()

G4bool G4ParticleHPJENDLHEData::IsZAApplicable ( const G4DynamicParticle * , G4double , G4double  )

inline

Definition at line 59 of file G4ParticleHPJENDLHEData.hh.

    {
      return false;
    }

---

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

• G4ParticleHPJENDLHEData.hh
• G4ParticleHPJENDLHEData.cc