G4NeutronHPElementData Class Reference

#include <G4NeutronHPElementData.hh>

Inheritance diagram for G4NeutronHPElementData:

Public Member Functions
	G4NeutronHPElementData ()
	~G4NeutronHPElementData ()
void	Init (G4Element *theElement)
void	UpdateData (G4int A, G4int Z, G4int index, G4double abundance)
void	UpdateData (G4int A, G4int Z, G4int M, G4int index, G4double abundance)
void	Harmonise (G4NeutronHPVector *&theStore, G4NeutronHPVector *theNew)
G4NeutronHPVector *	GetData (G4NeutronHPFissionData *)
G4NeutronHPVector *	GetData (G4NeutronHPCaptureData *)
G4NeutronHPVector *	GetData (G4NeutronHPElasticData *)
G4NeutronHPVector *	GetData (G4NeutronHPInelasticData *)
G4NeutronHPVector *	MakePhysicsVector (G4Element *theElement, G4ParticleDefinition *theP, G4NeutronHPFissionData *theSet)
G4NeutronHPVector *	MakePhysicsVector (G4Element *theElement, G4ParticleDefinition *theP, G4NeutronHPCaptureData *theSet)
G4NeutronHPVector *	MakePhysicsVector (G4Element *theElement, G4ParticleDefinition *theP, G4NeutronHPElasticData *theSet)
G4NeutronHPVector *	MakePhysicsVector (G4Element *theElement, G4ParticleDefinition *theP, G4NeutronHPInelasticData *theSet)
Public Member Functions inherited from G4HadronCrossSections
	G4HadronCrossSections ()
	~G4HadronCrossSections ()
G4bool	IsApplicable (const G4DynamicParticle *aParticle)
G4double	GetElasticCrossSection (const G4DynamicParticle *, G4int, G4int)
G4double	GetInelasticCrossSection (const G4DynamicParticle *, G4int, G4int)
G4double	GetCaptureCrossSection (const G4DynamicParticle *, G4int)
G4double	GetFissionCrossSection (const G4DynamicParticle *, G4int, G4int)
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel ()

Additional Inherited Members
Static Public Member Functions inherited from G4HadronCrossSections
static G4HadronCrossSections *	Instance ()

Detailed Description

Definition at line 52 of file G4NeutronHPElementData.hh.

Constructor & Destructor Documentation

G4NeutronHPElementData::G4NeutronHPElementData

(

)

Definition at line 35 of file G4NeutronHPElementData.cc.

  {
     precision = 0.02;
     theFissionData = new G4NeutronHPVector;
     theCaptureData = new G4NeutronHPVector;
     theElasticData = new G4NeutronHPVector;
     theInelasticData = new G4NeutronHPVector;
    theIsotopeWiseData = 0;
     theBuffer = NULL;
  }

G4NeutronHPElementData::~G4NeutronHPElementData

(

)

Definition at line 46 of file G4NeutronHPElementData.cc.

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

Member Function Documentation

G4NeutronHPVector* G4NeutronHPElementData::GetData ( G4NeutronHPFissionData *  )

inline

Definition at line 68 of file G4NeutronHPElementData.hh.

Referenced by MakePhysicsVector().

    {return theFissionData;}

G4NeutronHPVector* G4NeutronHPElementData::GetData ( G4NeutronHPCaptureData *  )

inline

Definition at line 70 of file G4NeutronHPElementData.hh.

    {return theCaptureData;}

G4NeutronHPVector* G4NeutronHPElementData::GetData ( G4NeutronHPElasticData *  )

inline

Definition at line 72 of file G4NeutronHPElementData.hh.

    {return theElasticData;}

G4NeutronHPVector* G4NeutronHPElementData::GetData ( G4NeutronHPInelasticData *  )

inline

Definition at line 74 of file G4NeutronHPElementData.hh.

    {return theInelasticData;}

void G4NeutronHPElementData::Harmonise	(	G4NeutronHPVector *&	theStore,
		G4NeutronHPVector *	theNew
	)

Definition at line 142 of file G4NeutronHPElementData.cc.

References test::a, active, G4NeutronHPVector::GetEnergy(), G4NeutronHPVector::GetVectorLength(), G4NeutronHPVector::GetXsec(), G4INCL::Math::max(), n, G4NeutronHPVector::SetData(), and test::x.

Referenced by UpdateData().

  {
    if(theNew == 0) { return; }
    G4int s_tmp = 0, n=0, m_tmp=0;
    G4NeutronHPVector * theMerge = new G4NeutronHPVector(theStore->GetVectorLength());
//    G4cout << "Harmonise 1: "<<theStore->GetEnergy(s)<<" "<<theNew->GetEnergy(0)<<G4endl;
    while ( theStore->GetEnergy(s_tmp)<theNew->GetEnergy(0)&&s_tmp<theStore->GetVectorLength() )
    {
      theMerge->SetData(m_tmp++, theStore->GetEnergy(s_tmp), theStore->GetXsec(s_tmp));
      s_tmp++;
    }
    G4NeutronHPVector *active = theStore;
    G4NeutronHPVector * passive = theNew;
    G4NeutronHPVector * tmp;
    G4int a = s_tmp, p = n, t;
//    G4cout << "Harmonise 2: "<<active->GetVectorLength()<<" "<<passive->GetVectorLength()<<G4endl;
    while (a<active->GetVectorLength()&&p<passive->GetVectorLength())
    {
      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 {
//        G4cout << "swapping in Harmonise"<<G4endl;
        tmp = active; t=a;
        active = passive; a=p;
        passive = tmp; p=t;
      }
    }
//    G4cout << "Harmonise 3: "<< a <<" "<<active->GetVectorLength()<<" "<<m<<G4endl;
    while (a!=active->GetVectorLength())
    {
      theMerge->SetData(m_tmp++, active->GetEnergy(a), active->GetXsec(a));
      a++;
    }
//    G4cout << "Harmonise 4: "<< p <<" "<<passive->GetVectorLength()<<" "<<m<<G4endl;
    while (p!=passive->GetVectorLength())
    {
      // Modified by T. KOI
      //theMerge->SetData(m++, passive->GetEnergy(p), passive->GetXsec(p));
      G4double x = passive->GetEnergy(p);
      G4double y = std::max(0., active->GetXsec(x));
      theMerge->SetData(m_tmp++, x, passive->GetXsec(p)+y);
      p++;
    }
//    G4cout <<"Harmonise 5: "<< theMerge->GetVectorLength() << " " << m << G4endl;
    delete theStore;
    theStore = theMerge;
//    G4cout <<"Harmonise 6: "<< theStore->GetVectorLength() << " " << m << G4endl;
  }

void G4NeutronHPElementData::Init

(

G4Element *

theElement

)

Definition at line 55 of file G4NeutronHPElementData.cc.

References G4StableIsotopes::GetAbundance(), G4StableIsotopes::GetFirstIsotope(), G4Element::GetIsotope(), G4StableIsotopes::GetIsotopeNucleonCount(), G4Isotope::Getm(), G4Isotope::GetN(), G4StableIsotopes::GetNumberOfIsotopes(), G4Element::GetNumberOfIsotopes(), G4Element::GetRelativeAbundanceVector(), G4Element::GetZ(), python.hepunit::perCent, G4NeutronHPVector::ThinOut(), and UpdateData().

Referenced by MakePhysicsVector().

  {
    G4int count = theElement->GetNumberOfIsotopes();
      if(count == 0) count +=
         theStableOnes.GetNumberOfIsotopes(static_cast<G4int>(theElement->GetZ()));
    theIsotopeWiseData = new G4NeutronHPIsoData[count];
    // filename = ein data-set je isotope.
    count = 0;
    G4int nIso = theElement->GetNumberOfIsotopes();
    G4int Z = static_cast<G4int> (theElement->GetZ());
    //G4int i1;
    if(nIso!=0)
    {
      for (G4int i1=0; i1<nIso; i1++)
      {
//        G4cout <<" Init: normal case"<<G4endl;
        G4int A = theElement->GetIsotope(i1)->GetN();
        G4int M = theElement->GetIsotope(i1)->Getm();
        G4double frac = theElement->GetRelativeAbundanceVector()[i1]/perCent;
        //UpdateData(A, Z, count++, frac);
        UpdateData(A, Z, M, count++, frac);
      }
    }else{
//      G4cout <<" Init: theStableOnes case: Z="<<Z<<G4endl;
      G4int first = theStableOnes.GetFirstIsotope(Z);
//      G4cout <<"first="<<first<<" "<<theStableOnes.GetNumberOfIsotopes(theElement->GetZ())<<G4endl;
      for(G4int i1=0; 
        i1<theStableOnes.GetNumberOfIsotopes(static_cast<G4int>(theElement->GetZ()) );
        i1++)
      {
//        G4cout <<" Init: theStableOnes in the loop"<<G4endl;
        G4int A = theStableOnes.GetIsotopeNucleonCount(first+i1);
        G4double frac = theStableOnes.GetAbundance(first+i1);
//        G4cout <<" Init: theStableOnes in the loop: "<<A<<G4endl;
        UpdateData(A, Z, count++, frac);
      }
    }
    theElasticData->ThinOut(precision);
    theInelasticData->ThinOut(precision);
    theCaptureData->ThinOut(precision);
    theFissionData->ThinOut(precision);
  }

G4NeutronHPVector * G4NeutronHPElementData::MakePhysicsVector	(	G4Element *	theElement,
		G4ParticleDefinition *	theP,
		G4NeutronHPFissionData *	theSet
	)

Definition at line 197 of file G4NeutronHPElementData.cc.

References GetData(), Init(), and G4Neutron::Neutron().

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

G4NeutronHPVector * G4NeutronHPElementData::MakePhysicsVector	(	G4Element *	theElement,
		G4ParticleDefinition *	theP,
		G4NeutronHPCaptureData *	theSet
	)

Definition at line 205 of file G4NeutronHPElementData.cc.

References GetData(), Init(), and G4Neutron::Neutron().

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

G4NeutronHPVector * G4NeutronHPElementData::MakePhysicsVector	(	G4Element *	theElement,
		G4ParticleDefinition *	theP,
		G4NeutronHPElasticData *	theSet
	)

Definition at line 213 of file G4NeutronHPElementData.cc.

References GetData(), Init(), and G4Neutron::Neutron().

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

G4NeutronHPVector * G4NeutronHPElementData::MakePhysicsVector	(	G4Element *	theElement,
		G4ParticleDefinition *	theP,
		G4NeutronHPInelasticData *	theSet
	)

Definition at line 221 of file G4NeutronHPElementData.cc.

References GetData(), Init(), and G4Neutron::Neutron().

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

void G4NeutronHPElementData::UpdateData ( G4int  A, G4int  Z, G4int  index, G4double  abundance  )

inline

Definition at line 63 of file G4NeutronHPElementData.hh.

References UpdateData().

Referenced by Init(), and UpdateData().

{ G4int M=0; UpdateData( A, Z, M, index, abundance); };

void G4NeutronHPElementData::UpdateData	(	G4int	A,
		G4int	Z,
		G4int	M,
		G4int	index,
		G4double	abundance
	)

Definition at line 99 of file G4NeutronHPElementData.cc.

References Harmonise(), G4NeutronHPIsoData::Init(), G4NeutronHPIsoData::MakeCaptureData(), G4NeutronHPIsoData::MakeElasticData(), G4NeutronHPIsoData::MakeFissionData(), and G4NeutronHPIsoData::MakeInelasticData().

  {
    //Reads in the Data, using G4NeutronHPIsoData[], and its Init
//    G4cout << "entered: ElementWiseData::UpdateData"<<G4endl;
    //theIsotopeWiseData[index].Init(A, Z, abundance);
    theIsotopeWiseData[index].Init(A, Z, M, abundance);
//    G4cout << "ElementWiseData::UpdateData Init finished"<<G4endl;

    theBuffer = theIsotopeWiseData[index].MakeElasticData();
//    G4cout << "ElementWiseData::UpdateData MakeElasticData finished: "
//         <<theBuffer->GetVectorLength()<<G4endl;
    Harmonise(theElasticData, theBuffer);
//    G4cout << "ElementWiseData::UpdateData Harmonise finished: "
//         <<theElasticData->GetVectorLength()<<G4endl;
    delete theBuffer;
    
    theBuffer = theIsotopeWiseData[index].MakeInelasticData();
//    G4cout << "ElementWiseData::UpdateData MakeInelasticData finished: "
//         <<theBuffer->GetVectorLength()<<G4endl;
    Harmonise(theInelasticData, theBuffer);
//    G4cout << "ElementWiseData::UpdateData Harmonise finished: "
//         <<theInelasticData->GetVectorLength()<<G4endl;
    delete theBuffer;
    
    theBuffer = theIsotopeWiseData[index].MakeCaptureData();
//    G4cout << "ElementWiseData::UpdateData MakeCaptureData finished: "
//         <<theBuffer->GetVectorLength()<<G4endl;
    Harmonise(theCaptureData, theBuffer);
//    G4cout << "ElementWiseData::UpdateData Harmonise finished: "
//         <<theCaptureData->GetVectorLength()<<G4endl;
    delete theBuffer;
    
    theBuffer = theIsotopeWiseData[index].MakeFissionData();
//    G4cout << "ElementWiseData::UpdateData MakeFissionData finished: "
//         <<theBuffer->GetVectorLength()<<G4endl;
    Harmonise(theFissionData, theBuffer);
//    G4cout << "ElementWiseData::UpdateData Harmonise finished: "
//         <<theFissionData->GetVectorLength()<<G4endl;
    delete theBuffer;
    
//    G4cout << "ElementWiseData::UpdateData finished"<endl;
  }

---

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

• G4NeutronHPElementData.hh
• G4NeutronHPElementData.cc