G4ParticleHPElementData Class Reference

#include <G4ParticleHPElementData.hh>

Public Member Functions
	G4ParticleHPElementData ()
	~G4ParticleHPElementData ()
void	Init (G4Element *theElement, G4ParticleDefinition *projectile, const char *dataDirVariable)
void	UpdateData (G4int A, G4int Z, G4int index, G4double abundance, G4ParticleDefinition *projectile, const char *dataDirVariable)
void	UpdateData (G4int A, G4int Z, G4int M, G4int index, G4double abundance, G4ParticleDefinition *projectile, const char *dataDirVariable)
void	Harmonise (G4ParticleHPVector *&theStore, G4ParticleHPVector *theNew)
G4ParticleHPVector *	GetData (G4ParticleHPFissionData *)
G4ParticleHPVector *	GetData (G4NeutronHPCaptureData *)
G4ParticleHPVector *	GetData (G4ParticleHPElasticData *)
G4ParticleHPVector *	GetData (G4ParticleHPInelasticData *)
G4ParticleHPVector *	MakePhysicsVector (G4Element *theElement, G4ParticleDefinition *projectile, G4ParticleHPFissionData *theSet, char *dataDirVariable)
G4ParticleHPVector *	MakePhysicsVector (G4Element *theElement, G4ParticleDefinition *projectile, G4NeutronHPCaptureData *theSet, char *dataDirVariable)
G4ParticleHPVector *	MakePhysicsVector (G4Element *theElement, G4ParticleDefinition *projectile, G4ParticleHPElasticData *theSet, char *dataDirVariable)
G4ParticleHPVector *	MakePhysicsVector (G4Element *theElement, G4ParticleDefinition *projectile, G4ParticleHPInelasticData *theSet, char *dataDirVariable)

Detailed Description

Definition at line 51 of file G4ParticleHPElementData.hh.

Constructor & Destructor Documentation

G4ParticleHPElementData()

G4ParticleHPElementData::G4ParticleHPElementData

(

)

Definition at line 37 of file G4ParticleHPElementData.cc.

{
  precision = 0.02;
  theFissionData = new G4ParticleHPVector;
  theCaptureData = new G4ParticleHPVector;
  theElasticData = new G4ParticleHPVector;
  theInelasticData = new G4ParticleHPVector;
  theIsotopeWiseData = nullptr;
  theBuffer = nullptr;
}

~G4ParticleHPElementData()

G4ParticleHPElementData::~G4ParticleHPElementData

(

)

Definition at line 48 of file G4ParticleHPElementData.cc.

{
  delete theFissionData;
  delete theCaptureData;
  delete theElasticData;
  delete theInelasticData;
  delete[] theIsotopeWiseData;
}

Member Function Documentation

GetData() [1/4]

G4ParticleHPVector * G4ParticleHPElementData::GetData ( G4NeutronHPCaptureData * )

inline

Definition at line 75 of file G4ParticleHPElementData.hh.

{ return theCaptureData; }

GetData() [2/4]

G4ParticleHPVector * G4ParticleHPElementData::GetData ( G4ParticleHPElasticData * )

inline

Definition at line 76 of file G4ParticleHPElementData.hh.

{ return theElasticData; }

GetData() [3/4]

G4ParticleHPVector * G4ParticleHPElementData::GetData ( G4ParticleHPFissionData * )

inline

Definition at line 74 of file G4ParticleHPElementData.hh.

{ return theFissionData; }

Referenced by MakePhysicsVector(), MakePhysicsVector(), MakePhysicsVector(), and MakePhysicsVector().

GetData() [4/4]

G4ParticleHPVector * G4ParticleHPElementData::GetData ( G4ParticleHPInelasticData * )

inline

Definition at line 77 of file G4ParticleHPElementData.hh.

{ return theInelasticData; }

Harmonise()

void G4ParticleHPElementData::Harmonise	(	G4ParticleHPVector *&	theStore,
		G4ParticleHPVector *	theNew )

Definition at line 104 of file G4ParticleHPElementData.cc.

{
  if (theNew == nullptr) {
    return;
  }
  G4int s_tmp = 0, n = 0, m_tmp = 0;
  auto theMerge = new G4ParticleHPVector(theStore->GetVectorLength());
  while (theStore->GetEnergy(s_tmp) < theNew->GetEnergy(0)
         && s_tmp < theStore->GetVectorLength())  // Loop checking, 11.05.2015, T. Koi
  {
    theMerge->SetData(m_tmp++, theStore->GetEnergy(s_tmp), theStore->GetXsec(s_tmp));
    ++s_tmp;
  }
  G4ParticleHPVector* active = theStore;
  G4ParticleHPVector* passive = theNew;
  G4ParticleHPVector* tmp;
  G4int a = s_tmp, p = n, t;
  while (a < active->GetVectorLength()
         && p < passive->GetVectorLength())  // Loop checking, 11.05.2015, T. Koi
  {
    if (active->GetEnergy(a) <= passive->GetEnergy(p)) {
      theMerge->SetData(m_tmp, active->GetEnergy(a), active->GetXsec(a));
      G4double x = theMerge->GetEnergy(m_tmp);
      G4double y = std::max(0., passive->GetXsec(x));
      theMerge->SetData(m_tmp, x, theMerge->GetXsec(m_tmp) + y);
      ++m_tmp;
      ++a;
    }
    else {
      tmp = active;
      t = a;
      active = passive;
      a = p;
      passive = tmp;
      p = t;
    }
  }
  while (a != active->GetVectorLength())  // Loop checking, 11.05.2015, T. Koi
  {
    theMerge->SetData(m_tmp++, active->GetEnergy(a), active->GetXsec(a));
    ++a;
  }
  while (p != passive->GetVectorLength())  // Loop checking, 11.05.2015, T. Koi
  {
    G4double x = passive->GetEnergy(p);
    G4double y = std::max(0., active->GetXsec(x));
    theMerge->SetData(m_tmp++, x, passive->GetXsec(p) + y);
    ++p;
  }
  delete theStore;
  theStore = theMerge;
}

Referenced by UpdateData().

Init()

void G4ParticleHPElementData::Init	(	G4Element *	theElement,
		G4ParticleDefinition *	projectile,
		const char *	dataDirVariable )

Definition at line 57 of file G4ParticleHPElementData.cc.

{
  auto nIso = (G4int)theElement->GetNumberOfIsotopes();
  auto Z = theElement->GetZasInt();
  const std::size_t dsize = nIso > 0 ? nIso : 1;
  theIsotopeWiseData = new G4ParticleHPIsoData[dsize];
 
  for (G4int i1 = 0; i1 < nIso; ++i1) {
    G4int A = theElement->GetIsotope(i1)->GetN();
    G4int M = theElement->GetIsotope(i1)->Getm();
    G4double frac = theElement->GetRelativeAbundanceVector()[i1] / CLHEP::perCent;
    UpdateData(A, Z, M, i1, frac, projectile, dataDirVariable);
  }
  theElasticData->ThinOut(precision);
  if (projectile == G4Neutron::Neutron()) theInelasticData->ThinOut(precision);
 
  theCaptureData->ThinOut(precision);
  theFissionData->ThinOut(precision);
}

Referenced by MakePhysicsVector(), MakePhysicsVector(), MakePhysicsVector(), and MakePhysicsVector().

MakePhysicsVector() [1/4]

G4ParticleHPVector * G4ParticleHPElementData::MakePhysicsVector	(	G4Element *	theElement,
		G4ParticleDefinition *	projectile,
		G4NeutronHPCaptureData *	theSet,
		char *	dataDirVariable )

Definition at line 168 of file G4ParticleHPElementData.cc.

{
  Init(theElement, projectile, dataDirVariable);
  return GetData(theSet);
}

MakePhysicsVector() [2/4]

G4ParticleHPVector * G4ParticleHPElementData::MakePhysicsVector	(	G4Element *	theElement,
		G4ParticleDefinition *	projectile,
		G4ParticleHPElasticData *	theSet,
		char *	dataDirVariable )

Definition at line 178 of file G4ParticleHPElementData.cc.

{
  Init(theElement, projectile, dataDirVariable);
  return GetData(theSet);
}

MakePhysicsVector() [3/4]

G4ParticleHPVector * G4ParticleHPElementData::MakePhysicsVector	(	G4Element *	theElement,
		G4ParticleDefinition *	projectile,
		G4ParticleHPFissionData *	theSet,
		char *	dataDirVariable )

Definition at line 158 of file G4ParticleHPElementData.cc.

{
  Init(theElement, projectile, dataDirVariable);
  return GetData(theSet);
}

MakePhysicsVector() [4/4]

G4ParticleHPVector * G4ParticleHPElementData::MakePhysicsVector	(	G4Element *	theElement,
		G4ParticleDefinition *	projectile,
		G4ParticleHPInelasticData *	theSet,
		char *	dataDirVariable )

Definition at line 187 of file G4ParticleHPElementData.cc.

{
  Init(theElement, projectile, dataDirVariable);
  return GetData(theSet);
}

UpdateData() [1/2]

void G4ParticleHPElementData::UpdateData ( G4int A, G4int Z, G4int index, G4double abundance, G4ParticleDefinition * projectile, const char * dataDirVariable )

inline

Definition at line 61 of file G4ParticleHPElementData.hh.

    {
      UpdateData(A, Z, 0, index, abundance, projectile, dataDirVariable);
    };

Referenced by Init(), and UpdateData().

UpdateData() [2/2]

void G4ParticleHPElementData::UpdateData	(	G4int	A,
		G4int	Z,
		G4int	M,
		G4int	index,
		G4double	abundance,
		G4ParticleDefinition *	projectile,
		const char *	dataDirVariable )

Definition at line 79 of file G4ParticleHPElementData.cc.

{
  // Reads in the Data, using G4ParticleHPIsoData[], and its Init
  theIsotopeWiseData[index].Init(A, Z, M, abundance, projectile, dataDirVariable);
 
  theBuffer = theIsotopeWiseData[index].MakeElasticData();
  Harmonise(theElasticData, theBuffer);
  delete theBuffer;
 
  theBuffer = theIsotopeWiseData[index].MakeInelasticData();
  Harmonise(theInelasticData, theBuffer);
  delete theBuffer;
 
  theBuffer = theIsotopeWiseData[index].MakeCaptureData();
  Harmonise(theCaptureData, theBuffer);
  delete theBuffer;
 
  theBuffer = theIsotopeWiseData[index].MakeFissionData();
  Harmonise(theFissionData, theBuffer);
  delete theBuffer;
}

---

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

• G4ParticleHPElementData.hh
• G4ParticleHPElementData.cc