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 (G4ParticleHPCaptureData *)
G4ParticleHPVector *	GetData (G4ParticleHPElasticData *)
G4ParticleHPVector *	GetData (G4ParticleHPInelasticData *)
G4ParticleHPVector *	MakePhysicsVector (G4Element *theElement, G4ParticleDefinition *projectile, G4ParticleHPFissionData *theSet, char *dataDirVariable)
G4ParticleHPVector *	MakePhysicsVector (G4Element *theElement, G4ParticleDefinition *projectile, G4ParticleHPCaptureData *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 50 of file G4ParticleHPElementData.hh.

Constructor & Destructor Documentation

G4ParticleHPElementData()

G4ParticleHPElementData::G4ParticleHPElementData

(

)

Definition at line 36 of file G4ParticleHPElementData.cc.

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

~G4ParticleHPElementData()

G4ParticleHPElementData::~G4ParticleHPElementData

(

)

Definition at line 47 of file G4ParticleHPElementData.cc.

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

Member Function Documentation

GetData() [1/4]

G4ParticleHPVector * G4ParticleHPElementData::GetData ( G4ParticleHPCaptureData *  )

inline

Definition at line 70 of file G4ParticleHPElementData.hh.

    {return theCaptureData;}

GetData() [2/4]

G4ParticleHPVector * G4ParticleHPElementData::GetData ( G4ParticleHPElasticData *  )

inline

Definition at line 72 of file G4ParticleHPElementData.hh.

    {return theElasticData;}

GetData() [3/4]

G4ParticleHPVector * G4ParticleHPElementData::GetData ( G4ParticleHPFissionData *  )

inline

Definition at line 68 of file G4ParticleHPElementData.hh.

    {return theFissionData;}

Referenced by MakePhysicsVector().

GetData() [4/4]

G4ParticleHPVector * G4ParticleHPElementData::GetData ( G4ParticleHPInelasticData *  )

inline

Definition at line 74 of file G4ParticleHPElementData.hh.

    {return theInelasticData;}

Harmonise()

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

Definition at line 119 of file G4ParticleHPElementData.cc.

  {
    if(theNew == 0) { return; }
    G4int s_tmp = 0, n=0, m_tmp=0;
    G4ParticleHPVector * 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 56 of file G4ParticleHPElementData.cc.

  {
    G4int count = (G4int)theElement->GetNumberOfIsotopes();
      if(count == 0)
        count += theStableOnes.GetNumberOfIsotopes((G4int)theElement->GetZ());
    theIsotopeWiseData = new G4ParticleHPIsoData[count];
    // filename = ein data-set je isotope.
    count = 0;
    G4int nIso = (G4int)theElement->GetNumberOfIsotopes();
    G4int Z = (G4int)theElement->GetZ();
 
    if(nIso!=0)
    {
      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, count++, frac);
        UpdateData(A, Z, M, count++, frac, projectile, dataDirVariable);
      }
    }
    else
    {
      G4int first = theStableOnes.GetFirstIsotope(Z);
      for(G4int i1=0; 
          i1<theStableOnes.GetNumberOfIsotopes((G4int)theElement->GetZ()); ++i1)
      {
        G4int A = theStableOnes.GetIsotopeNucleonCount(first+i1);
        G4double frac = theStableOnes.GetAbundance(first+i1);
        UpdateData(A, Z, count++, frac, projectile, dataDirVariable);
      }
    }
    theElasticData->ThinOut(precision);
    if( projectile == G4Neutron::Neutron() ) theInelasticData->ThinOut(precision);
 
    theCaptureData->ThinOut(precision);
    theFissionData->ThinOut(precision);
  }

Referenced by MakePhysicsVector().

MakePhysicsVector() [1/4]

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

Definition at line 176 of file G4ParticleHPElementData.cc.

  {
    if(projectile != G4Neutron::Neutron()) throw G4HadronicException(__FILE__, __LINE__, "not a neutron");
   Init ( theElement, projectile, dataDirVariable );
   return GetData(theSet);
  }

MakePhysicsVector() [2/4]

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

Definition at line 185 of file G4ParticleHPElementData.cc.

  {
    if(projectile != G4Neutron::Neutron()) throw G4HadronicException(__FILE__, __LINE__, "not a neutron");
   Init ( theElement, projectile, dataDirVariable );
   return GetData(theSet);
  }

MakePhysicsVector() [3/4]

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

Definition at line 167 of file G4ParticleHPElementData.cc.

  {
    if(projectile != G4Neutron::Neutron()) throw G4HadronicException(__FILE__, __LINE__, "not a neutron");
   Init ( theElement, projectile, dataDirVariable );
   return GetData(theSet);
  }

MakePhysicsVector() [4/4]

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

Definition at line 194 of file G4ParticleHPElementData.cc.

  {
    if(projectile != G4Neutron::Neutron()) throw G4HadronicException(__FILE__, __LINE__, "not a neutron");
   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 63 of file G4ParticleHPElementData.hh.

{ G4int M=0; UpdateData( A, Z, M, 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 96 of file G4ParticleHPElementData.cc.

  {
    //Reads in the Data, using G4ParticleHPIsoData[], and its Init
    //theIsotopeWiseData[index].Init(A, Z, abundance);
    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