| 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 // Derived from
27 // https://twiki.cern.ch/twiki/bin/view/Geant
28 // Courtesy of A. Dotti
29 //
30 // This example is provided by the Geant4-DNA
31 // Any report or published results obtained us
32 // shall cite the following Geant4-DNA collabo
33 // Med. Phys. 45 (2018) e722-e739
34 // Phys. Med. 31 (2015) 861-874
35 // Med. Phys. 37 (2010) 4692-4708
36 // Int. J. Model. Simul. Sci. Comput. 1 (2010)
37 //
38 // The Geant4-DNA web site is available at htt
39 //
40 /// \file medical/dna/svalue/src/MyFileReader.
41 /// \brief Implementation of the MyFileReader
42
43 #include "MyFileReader.hh"
44
45 //....oooOO0OOooo........oooOO0OOooo........oo
46
47 MyFileReader::MyFileReader()
48 {
49 //******************************************
50 //*** Specify file name containing list of i
51 //******************************************
52 G4String fileName = "spectrum.txt";
53 //
54
55 fInputFile.open(fileName.data());
56 }
57
58 //....oooOO0OOooo........oooOO0OOooo........oo
59
60 MyFileReader::~MyFileReader()
61 {
62 fInputFile.close();
63 }
64
65 //....oooOO0OOooo........oooOO0OOooo........oo
66
67 G4double MyFileReader::GetAnEvent()
68 {
69 //******************************************
70 //*** Specify number of lines to read at onc
71 //******************************************
72 G4int numberOfLinesToRead = 100;
73 //
74
75 if (fEvList.size() == 0) {
76 for (G4int i = 0; i < numberOfLinesToRead;
77 G4double nrj;
78 fInputFile >> nrj;
79 fEvList.push_back(nrj);
80 }
81 }
82
83 // Get first element from list
84 G4double ev = fEvList.front();
85
86 // Remove this element from list
87 fEvList.pop_front();
88
89 return ev;
90 }
91