Loading...
Searching...
No Matches
G4ChipsPionMinusInelasticXS Class Reference
#include <G4ChipsPionMinusInelasticXS.hh>
Inheritance diagram for G4ChipsPionMinusInelasticXS:
Static Public Member Functions | |
| static const char * | Default_Name () |
Additional Inherited Members | |
| Protected Attributes inherited from G4VCrossSectionDataSet | |
| G4int | verboseLevel |
| G4String | name |
Detailed Description
Definition at line 47 of file G4ChipsPionMinusInelasticXS.hh.
Constructor & Destructor Documentation
◆ G4ChipsPionMinusInelasticXS()
| G4ChipsPionMinusInelasticXS::G4ChipsPionMinusInelasticXS | ( | ) |
Definition at line 55 of file G4ChipsPionMinusInelasticXS.cc.
56{
57 // Initialization of the
58 lastLEN=0; // Pointer to lastArray of LowEn CS
59 lastHEN=0; // Pointer to lastArray of HighEn CS
60 lastN=0; // The last N of calculated nucleus
61 lastZ=0; // The last Z of calculated nucleus
62 lastP=0.; // Last used cross section Momentum
63 lastTH=0.; // Last threshold momentum
64 lastCS=0.; // Last value of the Cross Section
65 lastI=0; // The last position in the DAMDB
66 LEN = new std::vector<G4double*>;
67 HEN = new std::vector<G4double*>;
68}
◆ ~G4ChipsPionMinusInelasticXS()
| G4ChipsPionMinusInelasticXS::~G4ChipsPionMinusInelasticXS | ( | ) |
Definition at line 70 of file G4ChipsPionMinusInelasticXS.cc.
71{
72 std::size_t lens=LEN->size();
73 for(std::size_t i=0; i<lens; ++i) delete[] (*LEN)[i];
74 delete LEN;
75 std::size_t hens=HEN->size();
76 for(std::size_t i=0; i<hens; ++i) delete[] (*HEN)[i];
77 delete HEN;
78}
Member Function Documentation
◆ CrossSectionDescription()
|
virtual |
Reimplemented from G4VCrossSectionDataSet.
Definition at line 81 of file G4ChipsPionMinusInelasticXS.cc.
82{
83 outFile << "G4ChipsPionMinusInelasticXS provides the inelastic cross\n"
84 << "section for pion- nucleus scattering as a function of incident\n"
85 << "momentum. The cross section is calculated using M. Kossov's\n"
86 << "CHIPS parameterization of cross section data.\n";
87}
◆ Default_Name()
|
inlinestatic |
Definition at line 55 of file G4ChipsPionMinusInelasticXS.hh.
55{return "ChipsPionMinusInelasticXS";}
Referenced by G4ChipsComponentXS::G4ChipsComponentXS().
◆ GetChipsCrossSection()
|
virtual |
!The slave functions must provide cross-sections in millibarns (mb) !! (not in IU)
Definition at line 109 of file G4ChipsPionMinusInelasticXS.cc.
110{
111
113 if(tgN!=lastN || tgZ!=lastZ) // The nucleus was not the last used isotope
114 {
115 in = false; // By default the isotope haven't be found in AMDB
116 lastP = 0.; // New momentum history (nothing to compare with)
117 lastN = tgN; // The last N of the calculated nucleus
118 lastZ = tgZ; // The last Z of the calculated nucleus
120 j = 0; // A#0f records found in DB for this projectile
122 {
123 if(colN[i]==tgN && colZ[i]==tgZ) // Try the record "i" in the AMDB
124 {
125 lastI=i; // Remember the index for future fast/last use
126 lastTH =colTH[i]; // The last THreshold (A-dependent)
127 if(pMom<=lastTH)
128 {
129 return 0.; // Energy is below the Threshold value
130 }
131 lastP =colP [i]; // Last Momentum (A-dependent)
132 lastCS =colCS[i]; // Last CrossSect (A-dependent)
133 in = true; // This is the case when the isotop is found in DB
134 // Momentum pMom is in IU ! @@ Units
135 lastCS=CalculateCrossSection(-1,j,-211,lastZ,lastN,pMom); // read & update
136 if(lastCS<=0. && pMom>lastTH) // Correct the threshold (@@ No intermediate Zeros)
137 {
138 lastCS=0.;
139 lastTH=pMom;
140 }
141 break; // Go out of the LOOP
142 }
143 j++; // Increment a#0f records found in DB
144 }
145 if(!in) // This isotope has not been calculated previously
146 {
147 //!!The slave functions must provide cross-sections in millibarns (mb) !! (not in IU)
148 lastCS=CalculateCrossSection(0,j,-211,lastZ,lastN,pMom); //calculate & create
149 //if(lastCS>0.) // It means that the AMBD was initialized
150 //{
151
152 lastTH = 0; //ThresholdEnergy(tgZ, tgN); // The Threshold Energy which is now the last
153 colN.push_back(tgN);
154 colZ.push_back(tgZ);
155 colP.push_back(pMom);
156 colTH.push_back(lastTH);
157 colCS.push_back(lastCS);
158 //} // M.K. Presence of H1 with high threshold breaks the syncronization
159 return lastCS*millibarn;
160 } // End of creation of the new set of parameters
161 else
162 {
163 colP[lastI]=pMom;
164 colCS[lastI]=lastCS;
165 }
166 } // End of parameters udate
167 else if(pMom<=lastTH)
168 {
169 return 0.; // Momentum is below the Threshold Value -> CS=0
170 }
171 else // It is the last used -> use the current tables
172 {
173 lastCS=CalculateCrossSection(1,j,-211,lastZ,lastN,pMom); // Only read and UpdateDB
174 lastP=pMom;
175 }
176 return lastCS*millibarn;
177}
Referenced by G4ChipsComponentXS::GetInelasticElementCrossSection(), GetIsoCrossSection(), and G4ChipsComponentXS::GetTotalElementCrossSection().
◆ GetIsoCrossSection()
|
virtual |
Reimplemented from G4VCrossSectionDataSet.
Definition at line 98 of file G4ChipsPionMinusInelasticXS.cc.
102{
105
107}
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
Definition: G4ChipsPionMinusInelasticXS.cc:109
G4double GetTotalMomentum() const
◆ IsIsoApplicable()
|
virtual |
Reimplemented from G4VCrossSectionDataSet.
Definition at line 89 of file G4ChipsPionMinusInelasticXS.cc.
92{
93 return true;
94}
The documentation for this class was generated from the following files:
