#include <G4ChipsPionPlusInelasticXS.hh>

Inheritance diagram for G4ChipsPionPlusInelasticXS:

Collaboration diagram for G4ChipsPionPlusInelasticXS:

Public Member Functions
	G4ChipsPionPlusInelasticXS ()

	~G4ChipsPionPlusInelasticXS ()

virtual void	CrossSectionDescription (std::ostream &) const

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 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 48 of file G4ChipsPionPlusInelasticXS.hh.

Constructor & Destructor Documentation

G4ChipsPionPlusInelasticXS::G4ChipsPionPlusInelasticXS ( )

Definition at line 55 of file G4ChipsPionPlusInelasticXS.cc.

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

G4ChipsPionPlusInelasticXS::~G4ChipsPionPlusInelasticXS ( )

Definition at line 71 of file G4ChipsPionPlusInelasticXS.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

void G4ChipsPionPlusInelasticXS::CrossSectionDescription ( std::ostream & outFile ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 82 of file G4ChipsPionPlusInelasticXS.cc.

 {
     outFile << "G4ChipsPionPlusInelasticXS provides the inelastic cross\n"
             << "section for pion+ nucleus scattering as a function of incident\n"
             << "momentum. The cross section is calculated using M. Kossov's\n"
             << "CHIPS parameterization of cross section data.\n";
 }

static const char* G4ChipsPionPlusInelasticXS::Default_Name ( )

inlinestatic

Definition at line 57 of file G4ChipsPionPlusInelasticXS.hh.

57 {return "ChipsPionPlusInelasticXS";}

Here is the caller graph for this function:

G4double G4ChipsPionPlusInelasticXS::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 111 of file G4ChipsPionPlusInelasticXS.cc.

 {
 
   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,211,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,211,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,211,lastZ,lastN,pMom); // Only read and UpdateDB
     lastP=pMom;
   }
   return lastCS*millibarn;
 }

Here is the caller graph for this function:

G4double G4ChipsPionPlusInelasticXS::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 99 of file G4ChipsPionPlusInelasticXS.cc.

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

Here is the call graph for this function:

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 90 of file G4ChipsPionPlusInelasticXS.cc.

 {
   return true;
 }

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

source/geant4.10.03.p02/source/processes/hadronic/cross_sections/include/G4ChipsPionPlusInelasticXS.hh
source/geant4.10.03.p02/source/processes/hadronic/cross_sections/src/G4ChipsPionPlusInelasticXS.cc

Public Member Functions

Static Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation