G4ChipsNeutronInelasticXS Class Reference

#include <G4ChipsNeutronInelasticXS.hh>

Inheritance diagram for G4ChipsNeutronInelasticXS:

Public Member Functions
	G4ChipsNeutronInelasticXS ()
	~G4ChipsNeutronInelasticXS ()
virtual G4bool	IsIsoApplicable (const G4DynamicParticle *Pt, G4int Z, G4int A, const G4Element *elm, const G4Material *mat)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int tgZ, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
virtual G4double	GetChipsCrossSection (G4double momentum, G4int Z, G4int N, G4int pdg)
Static Public Member Functions
static const char *	Default_Name ()

---

Detailed Description

Definition at line 49 of file G4ChipsNeutronInelasticXS.hh.

---

Constructor & Destructor Documentation

G4ChipsNeutronInelasticXS::G4ChipsNeutronInelasticXS

(

)

Definition at line 54 of file G4ChipsNeutronInelasticXS.cc.

                                                    :G4VCrossSectionDataSet(Default_Name())
{
  lastLEN=0; // Pointer to the lastArray of LowEn CS
  lastHEN=0; // Pointer to the lastArray of HighEnCS
  lastN=0;   // The last N of calculated nucleus
  lastZ=0;   // The last Z of calculated nucleus
  lastP=0.;  // Last used in cross section Momentum
  lastTH=0.; // Last threshold momentum
  lastCS=0.; // Last value of the Cross Section
  lastI=0;   // The last position in the DAMDB
  HEthresh=0.;// HE threshold for the CS calculation
  LEN = new std::vector<G4double*>;
  HEN = new std::vector<G4double*>;
}

G4ChipsNeutronInelasticXS::~G4ChipsNeutronInelasticXS

(

)

Definition at line 69 of file G4ChipsNeutronInelasticXS.cc.

{
  G4int lens=LEN->size();
  for(G4int i=0; i<lens; ++i) delete[] (*LEN)[i];
  delete LEN;
  G4int hens=HEN->size();
  for(G4int i=0; i<hens; ++i) delete[] (*HEN)[i];
  delete HEN;
}

---

Member Function Documentation

static const char* G4ChipsNeutronInelasticXS::Default_Name

(

)

[inline, static]

Definition at line 59 of file G4ChipsNeutronInelasticXS.hh.

Referenced by G4ChipsComponentXS::G4ChipsComponentXS().

{return "ChipsNeutronInelasticXS";}

G4double G4ChipsNeutronInelasticXS::GetChipsCrossSection	(	G4double	momentum,
		G4int	Z,
		G4int	N,
		G4int	pdg
	)			[virtual]

Definition at line 102 of file G4ChipsNeutronInelasticXS.cc.

Referenced by GetIsoCrossSection().

{
  static G4int j;                      // A#0f Z/N-records already tested in AMDB
  static std::vector <G4int>    colN;  // Vector of N for calculated nuclei (isotops)
  static std::vector <G4int>    colZ;  // Vector of Z for calculated nuclei (isotops)
  static std::vector <G4double> colP;  // Vector of last momenta for the reaction
  static std::vector <G4double> colTH; // Vector of energy thresholds for the reaction
  static std::vector <G4double> colCS; // Vector of last cross sections for the reaction
  // ***---*** End of the mandatory Static Definitions of the Associative Memory ***---***


  G4bool in=false;                     // By default the isotope must be found in the AMDB
  if(tgN!=lastN || tgZ!=lastZ)         // The nucleus was not the last used isotope
  {
    in = false;                        // By default the isotope haven't be found in AMDB  
    lastP   = 0.;                      // New momentum history (nothing to compare with)
    lastN   = tgN;                     // The last N of the calculated nucleus
    lastZ   = tgZ;                     // The last Z of the calculated nucleus
    lastI   = colN.size();             // Size of the Associative Memory DB in the heap
    j  = 0;                            // A#0f records found in DB for this projectile
    if(lastI) for(G4int i=0; i<lastI; i++) // AMDB exists, try to find the (Z,N) isotope
    {
      if(colN[i]==tgN && colZ[i]==tgZ) // Try the record "i" in the AMDB
      {
        lastI=i;                       // Remember the index for future fast/last use
        lastTH =colTH[i];              // The last THreshold (A-dependent)
        if(pMom<=lastTH)
        {
          return 0.;                   // Energy is below the Threshold value
        }
        lastP  =colP [i];              // Last Momentum  (A-dependent)
        lastCS =colCS[i];              // Last CrossSect (A-dependent)
        in = true;                     // This is the case when the isotop is found in DB
        // Momentum pMom is in IU ! @@ Units
        lastCS=CalculateCrossSection(-1,j,2112,lastZ,lastN,pMom); // read & update
        if(lastCS<=0. && pMom>lastTH)  // Correct the threshold (@@ No intermediate Zeros)
        {
          lastCS=0.;
          lastTH=pMom;
        }
        break;                         // Go out of the LOOP
      }
      j++;                             // Increment a#0f records found in DB
    }
    if(!in)                            // This isotope has not been calculated previously
    {
      lastCS=CalculateCrossSection(0,j,2112,lastZ,lastN,pMom); //calculate & create
      //if(lastCS>0.)                   // It means that the AMBD was initialized
      //{

      lastTH = 0; //ThresholdEnergy(tgZ, tgN); // The Threshold Energy which is now the last
        colN.push_back(tgN);
        colZ.push_back(tgZ);
        colP.push_back(pMom);
        colTH.push_back(lastTH);
        colCS.push_back(lastCS);
        //} // M.K. Presence of H1 with high threshold breaks the syncronization
      return lastCS*millibarn;
    } // End of creation of the new set of parameters
    else
    {
      colP[lastI]=pMom;
      colCS[lastI]=lastCS;
    }
  } // End of parameters udate
  else if(pMom<=lastTH)
  {
    return 0.;                         // Momentum is below the Threshold Value -> CS=0
  }
  else                                 // It is the last used -> use the current tables
  {
    lastCS=CalculateCrossSection(1,j,2112,lastZ,lastN,pMom); // Only read and UpdateDB
    lastP=pMom;
  }
  return lastCS*millibarn;
}

G4double G4ChipsNeutronInelasticXS::GetIsoCrossSection	(	const G4DynamicParticle *	,
		G4int	tgZ,
		G4int	A,
		const G4Isotope *	iso = 0,
		const G4Element *	elm = 0,
		const G4Material *	mat = 0
	)			[virtual]

Reimplemented from G4VCrossSectionDataSet.

Definition at line 89 of file G4ChipsNeutronInelasticXS.cc.

References GetChipsCrossSection(), and G4DynamicParticle::GetTotalMomentum().

{
  G4double pMom=Pt->GetTotalMomentum();
  G4int tgN = A - tgZ;

  return GetChipsCrossSection(pMom, tgZ, tgN, 2112);
}

G4bool G4ChipsNeutronInelasticXS::IsIsoApplicable	(	const G4DynamicParticle *	Pt,
		G4int	Z,
		G4int	A,
		const G4Element *	elm,
		const G4Material *	mat
	)			[virtual]

Reimplemented from G4VCrossSectionDataSet.

Definition at line 79 of file G4ChipsNeutronInelasticXS.cc.

References G4DynamicParticle::GetDefinition(), and G4Neutron::Neutron().

{
  G4ParticleDefinition* particle = Pt->GetDefinition();
  if (particle == G4Neutron::Neutron()      ) return true;
  return false;
}

---

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

• G4ChipsNeutronInelasticXS.hh
• G4ChipsNeutronInelasticXS.cc

---