| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // 26 //
27 /// \file DicomRunAction.cc <<
28 /// \brief Implementation of the DicomRunActio <<
29 // <<
30 <<
31 #include "DicomRunAction.hh" 27 #include "DicomRunAction.hh"
32 <<
33 #include "DicomRun.hh" 28 #include "DicomRun.hh"
34 29
35 //-- In order to obtain detector information. 30 //-- In order to obtain detector information.
36 #include "G4RunManager.hh" <<
37 #include "G4SystemOfUnits.hh" <<
38 #include "G4THitsMap.hh" <<
39 #include "G4UnitsTable.hh" <<
40 <<
41 #include <fstream> 31 #include <fstream>
42 #include <iomanip> 32 #include <iomanip>
>> 33 #include "G4THitsMap.hh"
43 34
44 //....oooOO0OOooo........oooOO0OOooo........oo << 35 #include "G4UnitsTable.hh"
45 /// Constructor << 36 //=======================================================================
46 DicomRunAction::DicomRunAction() : G4UserRunAc << 37 // DicomRunAction
>> 38 //
>> 39 //
>> 40 //
>> 41 //=======================================================================
>> 42 // Constructor
>> 43 DicomRunAction::DicomRunAction():
>> 44 FieldName(15),
>> 45 FieldValue(14)
47 { 46 {
48 // - Prepare data member for DicomRun. 47 // - Prepare data member for DicomRun.
49 // vector represents a list of MultiFuncti 48 // vector represents a list of MultiFunctionalDetector names.
50 fSDName.push_back(G4String("phantomSD")); << 49 theSDName.push_back(G4String("phantomSD"));
51 } 50 }
52 51
53 //....oooOO0OOooo........oooOO0OOooo........oo << 52 // Destructor.
54 /// Destructor. <<
55 DicomRunAction::~DicomRunAction() 53 DicomRunAction::~DicomRunAction()
56 { 54 {
57 fSDName.clear(); << 55 theSDName.clear();
58 } 56 }
59 57
60 //....oooOO0OOooo........oooOO0OOooo........oo << 58 //
>> 59 //==
61 G4Run* DicomRunAction::GenerateRun() 60 G4Run* DicomRunAction::GenerateRun()
62 { 61 {
63 // Generate new RUN object, which is special 62 // Generate new RUN object, which is specially
64 // dedicated for MultiFunctionalDetector sch 63 // dedicated for MultiFunctionalDetector scheme.
65 // Detail description can be found in Dicom 64 // Detail description can be found in DicomRun.hh/cc.
66 // return new DicomRun(fSDName); << 65 return new DicomRun(theSDName);
67 return fDcmrun = new DicomRun(fSDName); <<
68 } 66 }
69 67
70 //....oooOO0OOooo........oooOO0OOooo........oo << 68 //
>> 69 //==
71 void DicomRunAction::BeginOfRunAction(const G4 70 void DicomRunAction::BeginOfRunAction(const G4Run* aRun)
72 { 71 {
73 G4cout << "### Run " << aRun->GetRunID() << 72 G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl;
74 // inform the runManager to save random numb <<
75 G4RunManager::GetRunManager()->SetRandomNumb <<
76 G4RunManager::GetRunManager()->SetRandomNumb <<
77 G4RunManager::GetRunManager()->SetRandomNumb <<
78 <<
79 <<
80 G4int progress = aRun->GetNumberOfEventToBeP <<
81 progress = (progress < 1) ? 1 : progress; <<
82 G4RunManager::GetRunManager()->SetPrintProgr <<
83 } 73 }
84 74
85 //....oooOO0OOooo........oooOO0OOooo........oo << 75 //
>> 76 //==
86 void DicomRunAction::EndOfRunAction(const G4Ru 77 void DicomRunAction::EndOfRunAction(const G4Run* aRun)
87 { 78 {
88 G4int nofEvents = aRun->GetNumberOfEvent(); << 79 G4cout << " ###### EndOfRunAction " <<G4endl;
89 << 80 //- DicomRun object.
90 static G4double local_total_dose = 0; << 81 DicomRun* re02Run = (DicomRun*)aRun;
91 G4double total_dose = 0; <<
92 <<
93 const DicomRun* reRun = static_cast<const Di <<
94 //--- Dump all scored quantities involved in 82 //--- Dump all scored quantities involved in DicomRun.
95 for (G4int i = 0; i < (G4int)fSDName.size(); << 83 for ( G4int i = 0; i < (G4int)theSDName.size(); i++ ){
96 // 84 //
97 //---------------------------------------- 85 //---------------------------------------------
98 // Dump accumulated quantities for this RU 86 // Dump accumulated quantities for this RUN.
99 // (Display only central region of x-y pl 87 // (Display only central region of x-y plane)
100 // 0 ConcreteSD/DoseDeposit 88 // 0 ConcreteSD/DoseDeposit
101 //---------------------------------------- 89 //---------------------------------------------
102 G4THitsMap<G4double>* DoseDeposit = reRun- << 90 G4THitsMap<G4double>* DoseDeposit = re02Run->GetHitsMap(theSDName[i]+"/DoseDeposit");
103 <<
104 if (DoseDeposit && DoseDeposit->GetMap()-> <<
105 std::map<G4int, G4double*>::iterator itr <<
106 for (; itr != DoseDeposit->GetMap()->end <<
107 if (!IsMaster()) { <<
108 local_total_dose += *(itr->second); <<
109 } <<
110 total_dose += *(itr->second); <<
111 } <<
112 } <<
113 } <<
114 91
115 if (IsMaster()) { << 92 G4cout << "=============================================================" <<G4endl;
116 G4cout << "--------------------End of Glob << 93 G4cout << " Number of event processed : "<< aRun->GetNumberOfEvent() << G4endl;
117 G4cout << " The run was " << nofEvents << << 94 G4cout << "=============================================================" <<G4endl;
118 G4cout << "LOCAL TOTAL DOSE : \t" << local << 95
119 G4cout << " TOTAL DOSE : \t" << total << 96 std::ofstream fileout;
120 } << 97 G4String fname = "dicom.out";
121 else { << 98 fileout.open(fname);
122 G4cout << "--------------------End of Loca << 99 G4cout << " opened file " << fname << " for dose output" << G4endl;
123 G4cout << " The run was " << nofEvents << << 100
124 G4cout << "LOCAL TOTAL DOSE : \t" << local << 101
125 G4cout << " TOTAL DOSE : \t" << total << 102 if( DoseDeposit && DoseDeposit->GetMap()->size() != 0 ) {
126 } << 103 std::ostream *myout = &G4cout;
127 << 104 PrintHeader(myout);
128 if (IsMaster()) { << 105 std::map<G4int,G4double*>::iterator itr = DoseDeposit->GetMap()->begin();
129 G4cout << " ###### EndOfRunAction ###### " << 106 for(; itr != DoseDeposit->GetMap()->end(); itr++) {
130 //- DicomRun object. << 107 fileout << itr->first
131 const DicomRun* re02Run = static_cast<cons << 108 << " " << *(itr->second)
132 //--- Dump all scored quantities involved << 109 << G4endl;
133 << 110 G4cout << " " << itr->first
134 for (G4int i = 0; i < (G4int)fSDName.size( << 111 << " " << std::setprecision(6) << *(itr->second) << " Gy"
135 // << 112 << G4endl;
136 //-------------------------------------- <<
137 // Dump accumulated quantities for this <<
138 // (Display only central region of x-y <<
139 // 0 ConcreteSD/DoseDeposit <<
140 //-------------------------------------- <<
141 G4THitsMap<G4double>* DoseDeposit = re02 <<
142 <<
143 G4cout << "============================= <<
144 G4cout << " Number of event processed : <<
145 G4cout << "============================= <<
146 <<
147 std::ofstream fileout; <<
148 G4String fname = "dicom.out"; <<
149 fileout.open(fname); <<
150 G4cout << " opened file " << fname << " <<
151 <<
152 if (DoseDeposit && DoseDeposit->GetMap() <<
153 std::ostream* myout = &G4cout; <<
154 PrintHeader(myout); <<
155 std::map<G4int, G4double*>::iterator i <<
156 for (; itr != DoseDeposit->GetMap()->e <<
157 fileout << itr->first << " " << <<
158 // G4cout << " " << itr->first <<
159 // << " " << std::setprec <<
160 // << *(itr->second)/CLHEP::g <<
161 // << G4endl; <<
162 } <<
163 G4cout << "=========================== <<
164 } 113 }
165 else { << 114 G4cout << "============================================="<<G4endl;
166 G4Exception("DicomRunAction", "000", J <<
167 "DoseDeposit HitsMap is ei <<
168 } <<
169 fileout.close(); <<
170 G4cout << " closed file " << fname << " <<
171 } 115 }
>> 116 fileout.close();
>> 117 G4cout << " closed file " << fname << " for dose output" << G4endl;
>> 118
172 } 119 }
173 <<
174 G4cout << "Finished : End of Run Action " << <<
175 } 120 }
>> 121 //
>> 122 // --
176 123
177 //....oooOO0OOooo........oooOO0OOooo........oo << 124 void DicomRunAction::PrintHeader(std::ostream *out)
178 void DicomRunAction::PrintHeader(std::ostream* <<
179 { 125 {
180 std::vector<G4String> vecScoreName; 126 std::vector<G4String> vecScoreName;
181 vecScoreName.push_back("DoseDeposit"); 127 vecScoreName.push_back("DoseDeposit");
182 128
183 // head line 129 // head line
184 // << 130 //std::string vname = FillString("Volume", ' ', FieldName+1);
>> 131 //*out << vname << '|';
185 std::string vname; 132 std::string vname;
186 *out << std::setw(10) << "Voxel" << 133 *out << std::setw(10) << "Voxel" << " |";
187 << " |"; << 134 for (std::vector<G4String>::iterator it = vecScoreName.begin();
188 for (std::vector<G4String>::iterator it = ve << 135 it != vecScoreName.end(); it++) {
189 // vname = FillString((*it), << 136 //vname = FillString((*it),
190 // ' ', << 137 // ' ',
191 // FieldValue+1, << 138 // FieldValue+1,
192 // false); << 139 // false);
193 // *out << vname << '|'; << 140 // *out << vname << '|';
194 *out << std::setw(fFieldValue) << (*it) << << 141 *out << std::setw(FieldValue) << (*it) << " |";
195 } 142 }
196 *out << G4endl; << 143 *out << G4endl;
197 } 144 }
198 145
199 //....oooOO0OOooo........oooOO0OOooo........oo << 146 std::string DicomRunAction::FillString(const std::string &name,
200 std::string DicomRunAction::FillString(const s << 147 char c, G4int n, G4bool back)
201 { 148 {
202 std::string fname(""); 149 std::string fname("");
203 G4int k = G4int(n - name.size()); << 150 G4int k = n - name.size();
204 if (k > 0) { 151 if (k > 0) {
205 if (back) { 152 if (back) {
206 fname = name; 153 fname = name;
207 fname += std::string(k, c); << 154 fname += std::string(k,c);
208 } 155 }
209 else { 156 else {
210 fname = std::string(k, c); << 157 fname = std::string(k,c);
211 fname += name; 158 fname += name;
212 } 159 }
213 } 160 }
214 else { 161 else {
215 fname = name; 162 fname = name;
216 } 163 }
217 return fname; 164 return fname;
218 } 165 }
219 166