#include <G4ChipsProtonInelasticXS.hh>

Inheritance diagram for G4ChipsProtonInelasticXS:

Public Member Functions
	G4ChipsProtonInelasticXS ()

	~G4ChipsProtonInelasticXS ()

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)

Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")

virtual	~G4VCrossSectionDataSet ()

virtual G4bool	IsElementApplicable (const G4DynamicParticle , G4int Z, const G4Material mat=0)

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=0)

G4double	ComputeCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=0)

virtual G4double	GetElementCrossSection (const G4DynamicParticle , G4int Z, const G4Material mat=0)

virtual G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy)

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

virtual void	CrossSectionDescription (std::ostream &) const

virtual G4int	GetVerboseLevel () const

virtual void	SetVerboseLevel (G4int value)

G4double	GetMinKinEnergy () const

void	SetMinKinEnergy (G4double value)

G4double	GetMaxKinEnergy () const

void	SetMaxKinEnergy (G4double value)

const G4String &	GetName () const

Static Public Member Functions
static const char *	Default_Name ()

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)

Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Detailed Description

Definition at line 44 of file G4ChipsProtonInelasticXS.hh.

Constructor & Destructor Documentation

G4ChipsProtonInelasticXS::G4ChipsProtonInelasticXS ( )

Definition at line 53 of file G4ChipsProtonInelasticXS.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*>;
 }

G4ChipsProtonInelasticXS::~G4ChipsProtonInelasticXS ( )

Definition at line 69 of file G4ChipsProtonInelasticXS.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* G4ChipsProtonInelasticXS::Default_Name ( )

inlinestatic

Definition at line 53 of file G4ChipsProtonInelasticXS.hh.

Referenced by G4ChipsComponentXS::G4ChipsComponentXS().

53 {return "ChipsProtonInelasticXS";}

G4double G4ChipsProtonInelasticXS::GetChipsCrossSection	(	G4double	momentum,
		G4int	Z,
		G4int	N,
		G4int	pdg
	)

virtual

!The slave functions must provide cross-sections in millibarns (mb) !! (not in IU)

Definition at line 104 of file G4ChipsProtonInelasticXS.cc.

References G4ThreadLocal, and python.hepunit::millibarn.

Referenced by GetIsoCrossSection().

 {
   static G4ThreadLocal G4int j;                      // A#0f Z/N-records already tested in AMDB
   static G4ThreadLocal std::vector <G4int>    *colN_G4MT_TLS_ = 0 ; if (!colN_G4MT_TLS_) colN_G4MT_TLS_ = new  std::vector <G4int>     ;  std::vector <G4int>    &colN = *colN_G4MT_TLS_;  // Vector of N for calculated nuclei (isotops)
   static G4ThreadLocal std::vector <G4int>    *colZ_G4MT_TLS_ = 0 ; if (!colZ_G4MT_TLS_) colZ_G4MT_TLS_ = new  std::vector <G4int>     ;  std::vector <G4int>    &colZ = *colZ_G4MT_TLS_;  // Vector of Z for calculated nuclei (isotops)
   static G4ThreadLocal std::vector <G4double> *colP_G4MT_TLS_ = 0 ; if (!colP_G4MT_TLS_) colP_G4MT_TLS_ = new  std::vector <G4double>  ;  std::vector <G4double> &colP = *colP_G4MT_TLS_;  // Vector of last momenta for the reaction
   static G4ThreadLocal std::vector <G4double> *colTH_G4MT_TLS_ = 0 ; if (!colTH_G4MT_TLS_) colTH_G4MT_TLS_ = new  std::vector <G4double>  ;  std::vector <G4double> &colTH = *colTH_G4MT_TLS_; // Vector of energy thresholds for the reaction
   static G4ThreadLocal std::vector <G4double> *colCS_G4MT_TLS_ = 0 ; if (!colCS_G4MT_TLS_) colCS_G4MT_TLS_ = new  std::vector <G4double>  ;  std::vector <G4double> &colCS = *colCS_G4MT_TLS_; // 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 been 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,2212,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
     {
       //!!The slave functions must provide cross-sections in millibarns (mb) !! (not in IU)
       lastCS=CalculateCrossSection(0,j,2212,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,2212,lastZ,lastN,pMom); // Only read and UpdateDB
     lastP=pMom;
   }
   return lastCS*millibarn;
 }

G4double G4ChipsProtonInelasticXS::GetIsoCrossSection	(	const G4DynamicParticle *	Pt,
		G4int	tgZ,
		G4int	A,
		const G4Isotope *	iso = `0`,
		const G4Element *	elm = `0`,
		const G4Material *	mat = `0`
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 93 of file G4ChipsProtonInelasticXS.cc.

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

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

G4bool G4ChipsProtonInelasticXS::IsIsoApplicable	(	const G4DynamicParticle *	Pt,
		G4int	Z,
		G4int	A,
		const G4Element *	elm,
		const G4Material *	mat
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 81 of file G4ChipsProtonInelasticXS.cc.

References G4DynamicParticle::GetDefinition(), and G4Proton::Proton().

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

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

Public Member Functions

Static Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation