| 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 // MR - 04/04/2012
27 // Based on G4DNACrossSectionDataSet
28 //
29
30
31 #include "G4MicroElecCrossSectionDataSet.hh"
32 #include "G4VDataSetAlgorithm.hh"
33 #include "G4EMDataSet.hh"
34 #include <vector>
35 #include <fstream>
36 #include <sstream>
37
38
39 G4MicroElecCrossSectionDataSet::G4MicroElecCro
40 G4double argUnitEnergies,
41 G4double argUnitData)
42 :
43 algorithm(argAlgorithm), unitEnergies(argUn
44 {;}
45
46 //....oooOO0OOooo........oooOO0OOooo........oo
47 G4MicroElecCrossSectionDataSet::~G4MicroElecCr
48 {
49 CleanUpComponents();
50
51 if (algorithm)
52 delete algorithm;
53 }
54
55 //....oooOO0OOooo........oooOO0OOooo........oo
56 G4bool G4MicroElecCrossSectionDataSet::LoadDat
57 {
58 CleanUpComponents();
59
60 G4String fullFileName(FullFileName(argFileNa
61 std::ifstream in(fullFileName, std::ifstream
62
63 if (!in.is_open())
64 {
65 G4String message("Data file \"");
66 message+=fullFileName;
67 message+="\" not found";
68 G4Exception("G4MicroElecCrossSectionData
69 FatalException,message);
70 return false;
71 }
72
73 std::vector<G4DataVector *> columns;
74 std::vector<G4DataVector *> log_columns;
75
76 std::stringstream *stream(new std::stringstr
77 char c;
78 G4bool comment(false);
79 G4bool space(true);
80 G4bool first(true);
81
82 try
83 {
84 while (!in.eof())
85 {
86 in.get(c);
87
88 switch (c)
89 {
90 case '\r':
91 case '\n':
92 if (!first)
93 {
94 unsigned long i(0);
95 G4double value;
96
97 while (!stream->eof())
98 {
99 (*stream) >> value;
100
101 while (i>=columns.size())
102 {
103 columns.push_back(new G4DataVector);
104 log_columns.push_back
105 }
106
107 columns[i]->push_back(value);
108 // N. A. Karakatsanis
109 // A condition is applied to check i
110 // If yes, then a near-zero value is
111 // If a value is zero, this simplifi
112 // If a value is negative, then it i
113 // logarithmic values should not be
114 //
115 // Therefore, G4LogLogInterpolation
116 // Instead, G4LinInterpolation is sa
117 // SemiLogInterpolation is safe only
118 // G4LinLogInterpolation is safe onl
119
120 if (value <=0.) value =
121 log_columns[i]->push_bac
122
123 ++i;
124 }
125 delete stream;
126 stream=new std::stringstream;
127 }
128 first=true;
129 comment=false;
130 space=true;
131 break;
132
133 case '#':
134 comment=true;
135 break;
136 case '\t':
137 case ' ':
138 space = true;
139 break;
140 default:
141 if (comment) { break; }
142 if (space && (!first)) { (*stream) <<
143 first=false;
144 (*stream) << c;
145 space=false;
146 }
147 }
148 }
149 catch(const std::ios::failure &e)
150 {
151 // some implementations of STL could thr
152 // when read wants read characters after
153 }
154
155 delete stream;
156
157 std::size_t maxI(columns.size());
158
159 if (maxI<2)
160 {
161 G4String message("Data file \"");
162 message+=fullFileName;
163 message+="\" should have at least two co
164 G4Exception("G4MicroElecCrossSectionData
165 FatalException,message);
166 return false;
167 }
168
169 std::size_t i(1);
170 while (i<maxI)
171 {
172 std::size_t maxJ(columns[i]->size());
173
174 if (maxJ!=columns[0]->size())
175 {
176 G4String message("Data file \"");
177 message+=fullFileName;
178 message+="\" has lines with a different nu
179 G4Exception("G4MicroElecCrossSectionDataSe
180 FatalException,message);
181 return false;
182 }
183
184 std::size_t j(0);
185
186 G4DataVector *argEnergies=new G4DataVect
187 G4DataVector *argData=new G4DataVector;
188 G4DataVector *argLogEnergies=new G4DataV
189 G4DataVector *argLogData=new G4DataVecto
190
191 while(j<maxJ)
192 {
193 argEnergies->push_back(columns[0]->operato
194 argData->push_back(columns[i]->operator[]
195 argLogEnergies->push_back(log_columns[0]->
196 argLogData->push_back(log_columns[i]->oper
197 j++;
198 }
199
200 AddComponent(new G4EMDataSet(G4int(i-1),
201 GetAlgorithm()->Clone(), GetUnitEne
202
203 ++i;
204 }
205
206 i=maxI;
207 while (i>0)
208 {
209 --i;
210 delete columns[i];
211 delete log_columns[i];
212 }
213
214 return true;
215 }
216
217 //....oooOO0OOooo........oooOO0OOooo........oo
218
219 G4bool G4MicroElecCrossSectionDataSet::LoadNon
220 {
221 CleanUpComponents();
222
223 G4String fullFileName(FullFileName(argFileNa
224 std::ifstream in(fullFileName, std::ifstream
225
226 if (!in.is_open())
227 {
228 G4String message("Data file \"");
229 message+=fullFileName;
230 message+="\" not found";
231 G4Exception("G4MicroElecCrossSectionData
232 FatalException,message);
233 return false;
234 }
235
236 std::vector<G4DataVector *> columns;
237
238 std::stringstream *stream(new std::stringstr
239 char c;
240 G4bool comment(false);
241 G4bool space(true);
242 G4bool first(true);
243
244 try
245 {
246 while (!in.eof())
247 {
248 in.get(c);
249
250 switch (c)
251 {
252 case '\r':
253 case '\n':
254 if (!first)
255 {
256 unsigned long i(0);
257 G4double value;
258
259 while (!stream->eof())
260 {
261 (*stream) >> value;
262
263 while (i>=columns.size())
264 {
265 columns.push_back(new G4DataVector);
266 }
267
268 columns[i]->push_back(value);
269
270 i++;
271 }
272
273 delete stream;
274 stream=new std::stringstream;
275 }
276
277 first=true;
278 comment=false;
279 space=true;
280 break;
281
282 case '#':
283 comment=true;
284 break;
285
286 case '\t':
287 case ' ':
288 space = true;
289 break;
290
291 default:
292 if (comment) { break; }
293 if (space && (!first)) { (*stream) <<
294
295 first=false;
296 (*stream) << c;
297 space=false;
298 }
299 }
300 }
301 catch(const std::ios::failure &e)
302 {
303 // some implementations of STL could thr
304 // when read wants read characters after
305 }
306
307 delete stream;
308
309 std::size_t maxI(columns.size());
310
311 if (maxI<2)
312 {
313 G4String message("Data file \"");
314 message+=fullFileName;
315 message+="\" should have at least two co
316 G4Exception("G4MicroElecCrossSectionData
317 FatalException,message);
318 return false;
319 }
320
321 std::size_t i(1);
322 while (i<maxI)
323 {
324 std::size_t maxJ(columns[i]->size());
325
326 if (maxJ!=columns[0]->size())
327 {
328 G4String message("Data file \"");
329 message+=fullFileName;
330 message+="\" has lines with a different nu
331 G4Exception("G4MicroElecCrossSection
332 FatalException,message);
333 return false;
334 }
335
336 std::size_t j(0);
337
338 G4DataVector *argEnergies=new G4DataVect
339 G4DataVector *argData=new G4DataVector;
340
341 while(j<maxJ)
342 {
343 argEnergies->push_back(columns[0]->operato
344 argData->push_back(columns[i]->operator[]
345 j++;
346 }
347
348 AddComponent(new G4EMDataSet(G4int(i-1),
349
350 i++;
351 }
352
353 i=maxI;
354 while (i>0)
355 {
356 i--;
357 delete columns[i];
358 }
359
360 return true;
361 }
362
363 //....oooOO0OOooo........oooOO0OOooo........oo
364
365 G4bool G4MicroElecCrossSectionDataSet::SaveDat
366 {
367 const std::size_t n(NumberOfComponents());
368 if (n==0)
369 {
370 G4Exception("G4MicroElecCrossSectionData
371 FatalException,"Expected
372
373 return false;
374 }
375
376 G4String fullFileName(FullFileName(argFileNa
377 std::ofstream out(fullFileName);
378
379 if (!out.is_open())
380 {
381 G4String message("Cannot open \"");
382 message+=fullFileName;
383 message+="\"";
384 G4Exception("G4MicroElecCrossSectionData
385 FatalException,message);
386 return false;
387 }
388
389 G4DataVector::const_iterator iEnergies(GetCo
390 G4DataVector::const_iterator iEnergiesEnd(Ge
391 G4DataVector::const_iterator * iData(new G4D
392
393 std::size_t k(n);
394
395 while (k>0)
396 {
397 k--;
398 iData[k]=GetComponent((G4int)k)->GetData
399 }
400
401 while (iEnergies!=iEnergiesEnd)
402 {
403 out.precision(10);
404 out.width(15);
405 out.setf(std::ofstream::left);
406 out << ((*iEnergies)/GetUnitEnergies());
407
408 k=0;
409
410 while (k<n)
411 {
412 out << ' ';
413 out.precision(10);
414 out.width(15);
415 out.setf(std::ofstream::left);
416 out << ((*(iData[k]))/GetUnitData());
417
418 iData[k]++;
419 k++;
420 }
421 out << std::endl;
422 iEnergies++;
423 }
424 delete[] iData;
425
426 return true;
427 }
428
429 //....oooOO0OOooo........oooOO0OOooo........oo
430
431 G4String G4MicroElecCrossSectionDataSet::FullF
432 {
433 const char* path = G4FindDataDir("G4LEDATA")
434 if (!path)
435 {
436 G4Exception("G4MicroElecCrossSectionData
437 FatalException,"G4LEDATA
438
439 return "";
440 }
441
442 std::ostringstream fullFileName;
443
444 fullFileName << path << "/" << argFileName <
445
446 return G4String(fullFileName.str().c_str());
447 }
448
449 //....oooOO0OOooo........oooOO0OOooo........oo
450
451 G4double G4MicroElecCrossSectionDataSet::FindV
452 {
453 // Returns the sum over the shells correspon
454 G4double value = 0.;
455
456 std::vector<G4VEMDataSet *>::const_iterator
457 std::vector<G4VEMDataSet *>::const_iterator
458
459 while (i!=end)
460 {
461 value+=(*i)->FindValue(argEnergy);
462 i++;
463 }
464
465 return value;
466 }
467
468 //....oooOO0OOooo........oooOO0OOooo........oo
469
470 void G4MicroElecCrossSectionDataSet::PrintData
471 {
472 const std::size_t n(NumberOfComponents());
473
474 G4cout << "The data set has " << n << " comp
475 G4cout << G4endl;
476
477 G4int i(0);
478
479 while (i<(G4int)n)
480 {
481 G4cout << "--- Component " << i << " ---
482 GetComponent(i)->PrintData();
483 i++;
484 }
485 }
486
487 //....oooOO0OOooo........oooOO0OOooo........oo
488
489 void G4MicroElecCrossSectionDataSet::SetEnergi
490 G4DataVector* argData,
491 G4int argComponentId)
492 {
493 G4VEMDataSet * component(components[argCompo
494
495 if (component)
496 {
497 component->SetEnergiesData(argEnergies,
498 return;
499 }
500
501 std::ostringstream message;
502 message << "Component " << argComponentId <<
503
504 G4Exception("G4MicroElecCrossSectionDataSet:
505 FatalException,message.str().
506
507 }
508
509 //....oooOO0OOooo........oooOO0OOooo........oo
510
511 void G4MicroElecCrossSectionDataSet::SetLogEne
512 G4DataVector* argData,
513 G4Dat
514 G4Dat
515 G4int argComponentId)
516 {
517 G4VEMDataSet * component(components[argCompo
518
519 if (component)
520 {
521 component->SetLogEnergiesData(argEnergie
522 return;
523 }
524
525 std::ostringstream message;
526 message << "Component " << argComponentId <<
527
528 G4Exception("G4MicroElecCrossSectionDataSet:
529 FatalException,message.str().
530
531 }
532
533 //....oooOO0OOooo........oooOO0OOooo........oo
534
535 void G4MicroElecCrossSectionDataSet::CleanUpCo
536 {
537 while (!components.empty())
538 {
539 if (components.back()) delete components
540 components.pop_back();
541 }
542 }
543
544
545