Geant4.10
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Static Public Member Functions
G4ChipsPionPlusInelasticXS Class Reference

#include <G4ChipsPionPlusInelasticXS.hh>

Inheritance diagram for G4ChipsPionPlusInelasticXS:
G4VCrossSectionDataSet

Public Member Functions

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

Constructor & Destructor Documentation

G4ChipsPionPlusInelasticXS::G4ChipsPionPlusInelasticXS ( )

Definition at line 51 of file G4ChipsPionPlusInelasticXS.cc.

51  :G4VCrossSectionDataSet(Default_Name())
52 {
53  // Initialization of the
54  lastLEN=0; // Pointer to lastArray of LowEn CS
55  lastHEN=0; // Pointer to lastArray of HighEn CS
56  lastN=0; // The last N of calculated nucleus
57  lastZ=0; // The last Z of calculated nucleus
58  lastP=0.; // Last used in cross section Momentum
59  lastTH=0.; // Last threshold momentum
60  lastCS=0.; // Last value of the Cross Section
61  lastI=0; // The last position in the DAMDB
62  LEN = new std::vector<G4double*>;
63  HEN = new std::vector<G4double*>;
64 }
G4VCrossSectionDataSet::G4VCrossSectionDataSet
G4VCrossSectionDataSet(const G4String &nam="")
Definition: G4VCrossSectionDataSet.cc:54
LEN
Definition: inflate.h:41
G4ChipsPionPlusInelasticXS::~G4ChipsPionPlusInelasticXS ( )

Definition at line 67 of file G4ChipsPionPlusInelasticXS.cc.

68 {
69  G4int lens=LEN->size();
70  for(G4int i=0; i<lens; ++i) delete[] (*LEN)[i];
71  delete LEN;
72  G4int hens=HEN->size();
73  for(G4int i=0; i<hens; ++i) delete[] (*HEN)[i];
74  delete HEN;
75 }
G4int
int G4int
Definition: G4Types.hh:78
LEN
Definition: inflate.h:41

Member Function Documentation

static const char* G4ChipsPionPlusInelasticXS::Default_Name ( )
inlinestatic

Definition at line 57 of file G4ChipsPionPlusInelasticXS.hh.

Referenced by G4ChipsComponentXS::G4ChipsComponentXS().

57 {return "ChipsPionPlusInelasticXS";}
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 100 of file G4ChipsPionPlusInelasticXS.cc.

References G4ThreadLocal, and python.hepunit::millibarn.

Referenced by GetIsoCrossSection().

101 {
102  static G4ThreadLocal G4int j; // A#0f Z/N-records already tested in AMDB
103  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)
104  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)
105  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
106  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
107  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
108  // ***---*** End of the mandatory Static Definitions of the Associative Memory ***---***
109 
110  G4bool in=false; // By default the isotope must be found in the AMDB
111  if(tgN!=lastN || tgZ!=lastZ) // The nucleus was not the last used isotope
112  {
113  in = false; // By default the isotope haven't be found in AMDB
114  lastP = 0.; // New momentum history (nothing to compare with)
115  lastN = tgN; // The last N of the calculated nucleus
116  lastZ = tgZ; // The last Z of the calculated nucleus
117  lastI = colN.size(); // Size of the Associative Memory DB in the heap
118  j = 0; // A#0f records found in DB for this projectile
119  if(lastI) for(G4int i=0; i<lastI; i++) // AMDB exists, try to find the (Z,N) isotope
120  {
121  if(colN[i]==tgN && colZ[i]==tgZ) // Try the record "i" in the AMDB
122  {
123  lastI=i; // Remember the index for future fast/last use
124  lastTH =colTH[i]; // The last THreshold (A-dependent)
125  if(pMom<=lastTH)
126  {
127  return 0.; // Energy is below the Threshold value
128  }
129  lastP =colP [i]; // Last Momentum (A-dependent)
130  lastCS =colCS[i]; // Last CrossSect (A-dependent)
131  in = true; // This is the case when the isotop is found in DB
132  // Momentum pMom is in IU ! @@ Units
133  lastCS=CalculateCrossSection(-1,j,211,lastZ,lastN,pMom); // read & update
134  if(lastCS<=0. && pMom>lastTH) // Correct the threshold (@@ No intermediate Zeros)
135  {
136  lastCS=0.;
137  lastTH=pMom;
138  }
139  break; // Go out of the LOOP
140  }
141  j++; // Increment a#0f records found in DB
142  }
143  if(!in) // This isotope has not been calculated previously
144  {
145  //!!The slave functions must provide cross-sections in millibarns (mb) !! (not in IU)
146  lastCS=CalculateCrossSection(0,j,211,lastZ,lastN,pMom); //calculate & create
147  //if(lastCS>0.) // It means that the AMBD was initialized
148  //{
149 
150  lastTH = 0; //ThresholdEnergy(tgZ, tgN); // The Threshold Energy which is now the last
151  colN.push_back(tgN);
152  colZ.push_back(tgZ);
153  colP.push_back(pMom);
154  colTH.push_back(lastTH);
155  colCS.push_back(lastCS);
156  //} // M.K. Presence of H1 with high threshold breaks the syncronization
157  return lastCS*millibarn;
158  } // End of creation of the new set of parameters
159  else
160  {
161  colP[lastI]=pMom;
162  colCS[lastI]=lastCS;
163  }
164  } // End of parameters udate
165  else if(pMom<=lastTH)
166  {
167  return 0.; // Momentum is below the Threshold Value -> CS=0
168  }
169  else // It is the last used -> use the current tables
170  {
171  lastCS=CalculateCrossSection(1,j,211,lastZ,lastN,pMom); // Only read and UpdateDB
172  lastP=pMom;
173  }
174  return lastCS*millibarn;
175 }
G4ThreadLocal
#define G4ThreadLocal
Definition: tls.hh:52
G4int
int G4int
Definition: G4Types.hh:78
python.hepunit.millibarn
int millibarn
Definition: hepunit.py:40
G4bool
bool G4bool
Definition: G4Types.hh:79
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 88 of file G4ChipsPionPlusInelasticXS.cc.

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

92 {
93  G4double pMom=Pt->GetTotalMomentum();
94  G4int tgN = A - tgZ;
95 
96  return GetChipsCrossSection(pMom, tgZ, tgN, 211);
97 }
G4int
int G4int
Definition: G4Types.hh:78
G4ChipsPionPlusInelasticXS::GetChipsCrossSection
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
Definition: G4ChipsPionPlusInelasticXS.cc:100
G4DynamicParticle::GetTotalMomentum
G4double GetTotalMomentum() const
G4double
double G4double
Definition: G4Types.hh:76
G4bool G4ChipsPionPlusInelasticXS::IsIsoApplicable ( const G4DynamicParticle *  Pt,
G4int  Z,
G4int  A,
const G4Element *  elm,
const G4Material *  mat 
)
virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 77 of file G4ChipsPionPlusInelasticXS.cc.

References G4DynamicParticle::GetDefinition(), and G4PionPlus::PionPlus().

80 {
81  G4ParticleDefinition* particle = Pt->GetDefinition();
82  if (particle == G4PionPlus::PionPlus() ) return true;
83  return false;
84 }
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4PionPlus::PionPlus
static G4PionPlus * PionPlus()
Definition: G4PionPlus.cc:98
The documentation for this class was generated from the following files:
Generated on Wed Apr 30 2014 15:56:16 for Geant4.10 by   doxygen 1.8.7