| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 // -------------------------------------------
28 //
29 // GEANT4 Class file
30 //
31 // File name: G4EmDataHandler
32 //
33 // Author: V. Ivanchenko
34 //
35 // Creation date: 16 August 2017
36 //
37 // Modifications:
38 //
39 // -------------------------------------------
40 //
41 //
42
43 #include "G4EmDataHandler.hh"
44 #include "G4EmDataRegistry.hh"
45 #include "G4ParticleDefinition.hh"
46 #include "G4EmParameters.hh"
47 #include "G4PhysicsTableHelper.hh"
48 #include "G4VEmProcess.hh"
49 #include "G4VEnergyLossProcess.hh"
50
51 //....oooOO0OOooo........oooOO0OOooo........oo
52
53 G4EmDataHandler::G4EmDataHandler(std::size_t n
54 : tLength(n), fName(nam)
55 {
56 data.resize(n, nullptr);
57 fMaxXS = new std::vector<G4double>;
58 fXSpeaks = new std::vector<G4TwoPeaksXS*>;
59 G4EmDataRegistry::Instance()->Register(this)
60 }
61
62 //....oooOO0OOooo........oooOO0OOooo........oo
63
64 G4EmDataHandler::~G4EmDataHandler()
65 {
66 if (!fUseBaseParticleTable) {
67 for (std::size_t i=0; i<tLength; ++i) {
68 CleanTable(i);
69 }
70 delete fMaxXS;
71 delete fXSpeaks;
72 }
73 if (!fElemSelectors.empty()) {
74 for (auto const & ptr : fElemSelectors) {
75 if (nullptr != ptr) {
76 for (auto const & p : *ptr) { delete p; }
77 }
78 }
79 }
80
81 }
82
83 //....oooOO0OOooo........oooOO0OOooo........oo
84
85 std::size_t G4EmDataHandler::SetTable(G4Physic
86 {
87 for (std::size_t i=0; i<tLength; ++i) {
88 if (ptr == data[i]) { return i; }
89 }
90 data.push_back(ptr);
91 ++tLength;
92 return tLength-1;
93 }
94
95 //....oooOO0OOooo........oooOO0OOooo........oo
96
97 void G4EmDataHandler::UpdateTable(G4PhysicsTab
98 {
99 // update table pointer but not delete previ
100 if (idx < tLength) {
101 if (ptr != data[idx]) { data[idx] = ptr; }
102 data[idx] = G4PhysicsTableHelper::PrepareP
103 } else {
104 G4cout << "### G4EmDataHandler::UpdateTabl
105 << " length=" << tLength << G4endl;
106 }
107 }
108
109 //....oooOO0OOooo........oooOO0OOooo........oo
110
111 void G4EmDataHandler::SaveTable(G4PhysicsTable
112 {
113 if (idx < tLength) {
114 data[idx] = ptr;
115 }
116 }
117
118 //....oooOO0OOooo........oooOO0OOooo........oo
119
120 G4PhysicsTable* G4EmDataHandler::MakeTable(std
121 {
122 std::size_t idx = i;
123 if (idx >= tLength) {
124 data.push_back(nullptr);
125 idx = tLength;
126 ++tLength;
127 }
128 data[idx] = G4PhysicsTableHelper::PreparePhy
129 return data[idx];
130 }
131
132 //....oooOO0OOooo........oooOO0OOooo........oo
133
134 G4PhysicsTable* G4EmDataHandler::MakeTable(G4P
135 {
136 std::size_t idx = i;
137 // create new table only if index correspond
138 // position in the vector
139 if (idx < tLength) {
140 if (ptr != data[idx]) {
141 CleanTable(idx);
142 data[idx] = ptr;
143 }
144 } else {
145 data.push_back(ptr);
146 idx = tLength;
147 ++tLength;
148 }
149 data[idx] = G4PhysicsTableHelper::PreparePhy
150 return data[idx];
151 }
152
153 //....oooOO0OOooo........oooOO0OOooo........oo
154
155 void G4EmDataHandler::CleanTable(std::size_t i
156 {
157 if (i < tLength && nullptr != data[i]) {
158 auto ptr = data[i];
159 ptr->clearAndDestroy();
160 delete ptr;
161 for (std::size_t j=0; j<tLength; ++j) {
162 if (ptr == data[j]) { data[j] = nullptr;
163 }
164 }
165 }
166
167 //....oooOO0OOooo........oooOO0OOooo........oo
168
169 G4bool G4EmDataHandler::StorePhysicsTable(std:
170 const G4ParticleDef
171 const G4String& fname,
172 G4bool ascii)
173 {
174 G4bool yes = true;
175 if(nullptr != data[idx]) {
176 yes = data[idx]->StorePhysicsTable(fname,
177
178 if ( yes ) {
179 G4cout << "### Physics table is stored f
180 << part->GetParticleName()
181 << " <" << fname << "> " << G4end
182 } else {
183 G4cout << "### Fail to store Physics Tab
184 << part->GetParticleName()
185 << " <" << fname << "> " << G4end
186 }
187 }
188 return yes;
189 }
190
191 //....oooOO0OOooo........oooOO0OOooo........oo
192
193 G4bool G4EmDataHandler::RetrievePhysicsTable(s
194 const G4ParticleDef
195 const G4String& fname,
196 G4bool ascii, G4bool spline)
197 {
198 G4PhysicsTable* table = Table(idx);
199 G4bool yes = G4PhysicsTableHelper::RetrieveP
200 G4EmParameters* param = G4EmParameters::Inst
201 if ( yes ) {
202 if (0 < param->Verbose()) {
203 G4cout << "### Physics table " << idx <<
204 << part->GetParticleName()
205 << " is retrieved from <" << fname << "
206 << G4endl;
207 }
208 } else if (1 < param->Verbose()) {
209 G4cout << "### Fail to retrieve physics ta
210 << part->GetParticleName() << " from <"
211 << fname << ">" << G4endl;
212 }
213 return yes;
214 }
215
216 //....oooOO0OOooo........oooOO0OOooo........oo
217
218 void G4EmDataHandler::SetMasterProcess(const G
219 {
220 masterProcess.push_back(ptr);
221 }
222
223 //....oooOO0OOooo........oooOO0OOooo........oo
224
225 const G4VEmProcess* G4EmDataHandler::GetMaster
226 {
227 return (idx < masterProcess.size()) ? master
228 }
229
230 //....oooOO0OOooo........oooOO0OOooo........oo
231
232 void G4EmDataHandler::SetElementSelectors(std:
233 std::size_t i)
234 {
235 if (i < eLength) {
236 if (fElemSelectors[i] != p) {
237 delete fElemSelectors[i];
238 }
239 fElemSelectors[i] = p;
240 } else {
241 fElemSelectors.push_back(p);
242 ++eLength;
243 }
244 }
245
246 //....oooOO0OOooo........oooOO0OOooo........oo
247