G4ChipsHyperonInelasticXS Class Reference

#include <G4ChipsHyperonInelasticXS.hh>

Inheritance diagram for G4ChipsHyperonInelasticXS:

Public Member Functions
	G4ChipsHyperonInelasticXS ()
	~G4ChipsHyperonInelasticXS ()
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 G4ChipsHyperonInelasticXS.hh.

---

Constructor & Destructor Documentation

G4ChipsHyperonInelasticXS::G4ChipsHyperonInelasticXS

(

)

Definition at line 53 of file G4ChipsHyperonInelasticXS.cc.

                                                    :G4VCrossSectionDataSet(Default_Name())
{
  // Initialization of the
  lastLEN=0; // Pointer to the lastArray of LowEn CS
  lastHEN=0; // Pointer to the lastArray of HighEn CS
  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
  LEN = new std::vector<G4double*>;
  HEN = new std::vector<G4double*>;
}

G4ChipsHyperonInelasticXS::~G4ChipsHyperonInelasticXS

(

)

Definition at line 68 of file G4ChipsHyperonInelasticXS.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* G4ChipsHyperonInelasticXS::Default_Name

(

)

[inline, static]

Definition at line 58 of file G4ChipsHyperonInelasticXS.hh.

Referenced by G4MiscBuilder::Build(), HadronPhysicsQGSP_FTFP_BERT_95::ConstructProcess(), G4HadronInelasticQBBC::ConstructProcess(), G4ChipsComponentXS::G4ChipsComponentXS(), and G4HyperonFTFPBuilder::G4HyperonFTFPBuilder().

{return "ChipsHyperonInelasticXS";}

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

Definition at line 129 of file G4ChipsHyperonInelasticXS.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,PDG,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,PDG,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,PDG,lastZ,lastN,pMom); // Only read and UpdateDB
    lastP=pMom;
  }
  return lastCS*millibarn;
}

G4double G4ChipsHyperonInelasticXS::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 117 of file G4ChipsHyperonInelasticXS.cc.

References GetChipsCrossSection(), G4DynamicParticle::GetDefinition(), G4ParticleDefinition::GetPDGEncoding(), and G4DynamicParticle::GetTotalMomentum().

{
  G4double pMom=Pt->GetTotalMomentum();
  G4int tgN = A - tgZ;
  G4int pdg = Pt->GetDefinition()->GetPDGEncoding();
  
  return GetChipsCrossSection(pMom, tgZ, tgN, pdg);
}

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

Reimplemented from G4VCrossSectionDataSet.

Definition at line 79 of file G4ChipsHyperonInelasticXS.cc.

References G4DynamicParticle::GetDefinition(), G4Lambda::Lambda(), G4OmegaMinus::OmegaMinus(), G4SigmaMinus::SigmaMinus(), G4SigmaPlus::SigmaPlus(), G4SigmaZero::SigmaZero(), G4XiMinus::XiMinus(), and G4XiZero::XiZero().

{
  G4ParticleDefinition* particle = Pt->GetDefinition();
  if (particle == G4Lambda::Lambda()) 
    {
      return true;
    }
  else if(particle == G4SigmaPlus::SigmaPlus())
    {
    return true;
    }
  else if(particle == G4SigmaMinus::SigmaMinus())
    {
    return true;
    }
  else if(particle == G4SigmaZero::SigmaZero())
    {
      return true;
    }
  else if(particle == G4XiMinus::XiMinus())
    {
      return true;
    }
  else if(particle == G4XiZero::XiZero())
    {
      return true;
    }
  else if(particle == G4OmegaMinus::OmegaMinus())
    {
      return true;
    }
  return false;
}

---

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

• G4ChipsHyperonInelasticXS.hh
• G4ChipsHyperonInelasticXS.cc

---