Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4ChipsKaonMinusInelasticXS Class Reference

#include <G4ChipsKaonMinusInelasticXS.hh>

Inheritance diagram for G4ChipsKaonMinusInelasticXS:
Collaboration diagram for G4ChipsKaonMinusInelasticXS:

Public Member Functions

 G4ChipsKaonMinusInelasticXS ()
 
 ~G4ChipsKaonMinusInelasticXS ()
 
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 G4IsotopeSelectIsotope (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 G4StringGetName () 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 45 of file G4ChipsKaonMinusInelasticXS.hh.

Constructor & Destructor Documentation

G4ChipsKaonMinusInelasticXS::G4ChipsKaonMinusInelasticXS ( )

Definition at line 66 of file G4ChipsKaonMinusInelasticXS.cc.

67 {
68  lastLEN=0; // Pointer to lastArray of LowEn CS
69  lastHEN=0; // Pointer to lastArray of HighEn CS
70  lastN=0; // The last N of calculated nucleus
71  lastZ=0; // The last Z of calculated nucleus
72  lastP=0.; // Last used in CrossSection Momentum
73  lastTH=0.; // Last threshold momentum
74  lastCS=0.; // Last value of the Cross Section
75  lastI=0; // The last position in the DAMDB
76  LEN = new std::vector<G4double*>;
77  HEN = new std::vector<G4double*>;
78 }
G4VCrossSectionDataSet(const G4String &nam="")
Definition: inflate.h:41
G4ChipsKaonMinusInelasticXS::~G4ChipsKaonMinusInelasticXS ( )

Definition at line 80 of file G4ChipsKaonMinusInelasticXS.cc.

81 {
82  G4int lens=LEN->size();
83  for(G4int i=0; i<lens; ++i) delete[] (*LEN)[i];
84  delete LEN;
85 
86  G4int hens=HEN->size();
87  for(G4int i=0; i<hens; ++i) delete[] (*HEN)[i];
88  delete HEN;
89 }
int G4int
Definition: G4Types.hh:78
Definition: inflate.h:41

Member Function Documentation

void G4ChipsKaonMinusInelasticXS::CrossSectionDescription ( std::ostream &  outFile) const
virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 92 of file G4ChipsKaonMinusInelasticXS.cc.

93 {
94  outFile << "G4ChipsKaonMinusInelasticXS provides the inelastic cross\n"
95  << "section for K- nucleus scattering as a function of incident\n"
96  << "momentum. The cross section is calculated using M. Kossov's\n"
97  << "CHIPS parameterization of cross section data.\n";
98 }
static const char* G4ChipsKaonMinusInelasticXS::Default_Name ( )
inlinestatic

Definition at line 54 of file G4ChipsKaonMinusInelasticXS.hh.

54 {return "ChipsKaonMinusInelasticXS";}

Here is the caller graph for this function:

G4double G4ChipsKaonMinusInelasticXS::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 121 of file G4ChipsKaonMinusInelasticXS.cc.

122 {
123  G4bool in=false; // By default the isotope must be found in the AMDB
124  if(tgN!=lastN || tgZ!=lastZ) // The nucleus was not the last used isotope
125  {
126  in = false; // By default the isotope haven't be found in AMDB
127  lastP = 0.; // New momentum history (nothing to compare with)
128  lastN = tgN; // The last N of the calculated nucleus
129  lastZ = tgZ; // The last Z of the calculated nucleus
130  lastI = colN.size(); // Size of the Associative Memory DB in the heap
131  j = 0; // A#0f records found in DB for this projectile
132  if(lastI) for(G4int i=0; i<lastI; i++) // AMDB exists, try to find the (Z,N) isotope
133  {
134  if(colN[i]==tgN && colZ[i]==tgZ) // Try the record "i" in the AMDB
135  {
136  lastI=i; // Remember the index for future fast/last use
137  lastTH =colTH[i]; // The last THreshold (A-dependent)
138  if(pMom<=lastTH)
139  {
140  return 0.; // Energy is below the Threshold value
141  }
142  lastP =colP [i]; // Last Momentum (A-dependent)
143  lastCS =colCS[i]; // Last CrossSect (A-dependent)
144  in = true; // This is the case when the isotop is found in DB
145  // Momentum pMom is in IU ! @@ Units
146  lastCS=CalculateCrossSection(-1,j,-321,lastZ,lastN,pMom); // read & update
147  if(lastCS<=0. && pMom>lastTH) // Correct the threshold (@@ No intermediate Zeros)
148  {
149  lastCS=0.;
150  lastTH=pMom;
151  }
152  break; // Go out of the LOOP
153  }
154  j++; // Increment a#0f records found in DB
155  }
156  if(!in) // This isotope has not been calculated previously
157  {
159  lastCS=CalculateCrossSection(0,j,-321,lastZ,lastN,pMom); //calculate & create
160  //if(lastCS>0.) // It means that the AMBD was initialized
161  //{
162 
163  // lastTH = ThresholdEnergy(tgZ, tgN); // The Threshold Energy which is now the last
164 
165  lastTH = 0; // WP - to be checked!!!
166  colN.push_back(tgN);
167  colZ.push_back(tgZ);
168  colP.push_back(pMom);
169  colTH.push_back(lastTH);
170  colCS.push_back(lastCS);
171  //} // M.K. Presence of H1 with high threshold breaks the syncronization
172  return lastCS*millibarn;
173  } // End of creation of the new set of parameters
174  else
175  {
176  colP[lastI]=pMom;
177  colCS[lastI]=lastCS;
178  }
179  } // End of parameters udate
180  else if(pMom<=lastTH)
181  {
182  return 0.; // Momentum is below the Threshold Value -> CS=0
183  }
184  else // It is the last used -> use the current tables
185  {
186  lastCS=CalculateCrossSection(1,j,-321,lastZ,lastN,pMom); // Only read and UpdateDB
187  lastP=pMom;
188  }
189  return lastCS*millibarn;
190 }
int G4int
Definition: G4Types.hh:78
bool G4bool
Definition: G4Types.hh:79
static constexpr double millibarn
Definition: G4SIunits.hh:106

Here is the caller graph for this function:

G4double G4ChipsKaonMinusInelasticXS::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 110 of file G4ChipsKaonMinusInelasticXS.cc.

114 {
115  G4double pMom=Pt->GetTotalMomentum();
116  G4int tgN = A - tgZ;
117 
118  return GetChipsCrossSection(pMom, tgZ, tgN, -321);
119 }
int G4int
Definition: G4Types.hh:78
G4double GetTotalMomentum() const
double A(double temperature)
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

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

Reimplemented from G4VCrossSectionDataSet.

Definition at line 100 of file G4ChipsKaonMinusInelasticXS.cc.

103 {
104  return true;
105 }

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