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 // $Id: DicomRunAction.cc 74809 2013-10-22 09:49:26Z gcosmo $ 26 // 27 // 27 /// \file DicomRunAction.cc 28 /// \file DicomRunAction.cc 28 /// \brief Implementation of the DicomRunActio 29 /// \brief Implementation of the DicomRunAction class 29 // 30 // 30 31 31 #include "DicomRunAction.hh" 32 #include "DicomRunAction.hh" 32 << 33 #include "DicomRun.hh" 33 #include "DicomRun.hh" 34 34 35 //-- In order to obtain detector information. 35 //-- In order to obtain detector information. 36 #include "G4RunManager.hh" << 36 #include <fstream> 37 #include "G4SystemOfUnits.hh" << 37 #include <iomanip> 38 #include "G4THitsMap.hh" 38 #include "G4THitsMap.hh" >> 39 39 #include "G4UnitsTable.hh" 40 #include "G4UnitsTable.hh" >> 41 #include "G4SystemOfUnits.hh" 40 42 41 #include <fstream> << 43 #include "G4RunManager.hh" 42 #include <iomanip> << 44 >> 45 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 46 >> 47 DicomRunAction* DicomRunAction::fgInstance = 0; >> 48 >> 49 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 50 >> 51 DicomRunAction* DicomRunAction::Instance() { return fgInstance; } 43 52 44 //....oooOO0OOooo........oooOO0OOooo........oo 53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 45 /// Constructor 54 /// Constructor 46 DicomRunAction::DicomRunAction() : G4UserRunAc << 55 DicomRunAction::DicomRunAction() >> 56 : G4UserRunAction(), dcmrun(0), fFieldValue(14) 47 { 57 { 48 // - Prepare data member for DicomRun. << 58 // - Prepare data member for DicomRun. 49 // vector represents a list of MultiFuncti << 59 // vector represents a list of MultiFunctionalDetector names. 50 fSDName.push_back(G4String("phantomSD")); << 60 fSDName.push_back(G4String("phantomSD")); >> 61 fgInstance = this; 51 } 62 } 52 63 53 //....oooOO0OOooo........oooOO0OOooo........oo 64 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 54 /// Destructor. 65 /// Destructor. 55 DicomRunAction::~DicomRunAction() 66 DicomRunAction::~DicomRunAction() 56 { 67 { 57 fSDName.clear(); << 68 fSDName.clear(); >> 69 fgInstance = 0; 58 } 70 } 59 71 60 //....oooOO0OOooo........oooOO0OOooo........oo 72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 61 G4Run* DicomRunAction::GenerateRun() 73 G4Run* DicomRunAction::GenerateRun() 62 { 74 { 63 // Generate new RUN object, which is special << 75 // Generate new RUN object, which is specially 64 // dedicated for MultiFunctionalDetector sch << 76 // dedicated for MultiFunctionalDetector scheme. 65 // Detail description can be found in Dicom << 77 // Detail description can be found in DicomRun.hh/cc. 66 // return new DicomRun(fSDName); << 78 //return new DicomRun(fSDName); 67 return fDcmrun = new DicomRun(fSDName); << 79 return dcmrun = new DicomRun(fSDName); 68 } 80 } 69 81 70 //....oooOO0OOooo........oooOO0OOooo........oo 82 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 71 void DicomRunAction::BeginOfRunAction(const G4 83 void DicomRunAction::BeginOfRunAction(const G4Run* aRun) 72 { 84 { 73 G4cout << "### Run " << aRun->GetRunID() << << 85 G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl; 74 // inform the runManager to save random numb << 86 //inform the runManager to save random number seed 75 G4RunManager::GetRunManager()->SetRandomNumb << 87 G4RunManager::GetRunManager()->SetRandomNumberStore(false); 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 } 88 } 84 89 85 //....oooOO0OOooo........oooOO0OOooo........oo 90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 86 void DicomRunAction::EndOfRunAction(const G4Ru 91 void DicomRunAction::EndOfRunAction(const G4Run* aRun) 87 { 92 { 88 G4int nofEvents = aRun->GetNumberOfEvent(); << 93 G4int nofEvents = aRun->GetNumberOfEvent(); >> 94 if (nofEvents == 0) return; >> 95 //print >> 96 // 89 97 90 static G4double local_total_dose = 0; << 98 static double local_total_dose = 0; 91 G4double total_dose = 0; << 99 double total_dose = 0; 92 100 93 const DicomRun* reRun = static_cast<const Di << 101 const DicomRun* reRun = static_cast<const DicomRun*>(aRun); 94 //--- Dump all scored quantities involved in << 102 //--- Dump all scored quantities involved in DicomRun. 95 for (G4int i = 0; i < (G4int)fSDName.size(); << 103 for ( G4int i = 0; i < (G4int)fSDName.size(); i++ ){ 96 // << 104 // 97 //---------------------------------------- << 105 //--------------------------------------------- 98 // Dump accumulated quantities for this RU << 106 // Dump accumulated quantities for this RUN. 99 // (Display only central region of x-y pl << 107 // (Display only central region of x-y plane) 100 // 0 ConcreteSD/DoseDeposit << 108 // 0 ConcreteSD/DoseDeposit 101 //---------------------------------------- << 109 //--------------------------------------------- 102 G4THitsMap<G4double>* DoseDeposit = reRun- << 110 G4THitsMap<G4double>* DoseDeposit = reRun->GetHitsMap(fSDName[i]+"/DoseDeposit"); 103 << 111 104 if (DoseDeposit && DoseDeposit->GetMap()-> << 112 if( DoseDeposit && DoseDeposit->GetMap()->size() != 0 ) { 105 std::map<G4int, G4double*>::iterator itr << 113 std::map<G4int,G4double*>::iterator itr = DoseDeposit->GetMap()->begin(); 106 for (; itr != DoseDeposit->GetMap()->end << 114 for(; itr != DoseDeposit->GetMap()->end(); itr++) { 107 if (!IsMaster()) { << 115 if(!IsMaster()) { local_total_dose += *(itr->second); } 108 local_total_dose += *(itr->second); << 116 total_dose += *(itr->second); >> 117 } 109 } 118 } 110 total_dose += *(itr->second); << 119 111 } << 112 } 120 } 113 } << 121 >> 122 if (IsMaster()) >> 123 { >> 124 G4cout >> 125 << "\n--------------------End of Global Run-----------------------" >> 126 << " \n The run was " << nofEvents << " events " << G4endl; >> 127 G4cout << "LOCAL TOTAL DOSE : \t" << local_total_dose/gray << " Gy" << std::endl; >> 128 G4cout << " TOTAL DOSE : \t" << total_dose/gray << " Gy" << std::endl; 114 129 115 if (IsMaster()) { << 130 } 116 G4cout << "--------------------End of Glob << 131 else 117 G4cout << " The run was " << nofEvents << << 132 { 118 G4cout << "LOCAL TOTAL DOSE : \t" << local << 133 G4cout 119 G4cout << " TOTAL DOSE : \t" << total << 134 << "\n--------------------End of Local Run------------------------" 120 } << 135 << " \n The run was " << nofEvents << G4endl; 121 else { << 136 G4cout << "LOCAL TOTAL DOSE : \t" << local_total_dose/gray << " Gy" << std::endl; 122 G4cout << "--------------------End of Loca << 137 G4cout << " TOTAL DOSE : \t" << total_dose/gray << " Gy" << std::endl; 123 G4cout << " The run was " << nofEvents << << 124 G4cout << "LOCAL TOTAL DOSE : \t" << local << 125 G4cout << " TOTAL DOSE : \t" << total << 126 } << 127 138 128 if (IsMaster()) { << 139 } 129 G4cout << " ###### EndOfRunAction ###### " << 130 //- DicomRun object. << 131 const DicomRun* re02Run = static_cast<cons << 132 //--- Dump all scored quantities involved << 133 140 134 for (G4int i = 0; i < (G4int)fSDName.size( << 141 if(IsMaster()) { 135 // << 142 G4cout << " ###### EndOfRunAction ###### " << G4endl; 136 //-------------------------------------- << 143 //- DicomRun object. 137 // Dump accumulated quantities for this << 144 const DicomRun* re02Run = static_cast<const DicomRun*>(aRun); 138 // (Display only central region of x-y << 145 //--- Dump all scored quantities involved in DicomRun. 139 // 0 ConcreteSD/DoseDeposit << 146 140 //-------------------------------------- << 147 for ( G4int i = 0; i < (G4int)fSDName.size(); i++ ){ 141 G4THitsMap<G4double>* DoseDeposit = re02 << 148 // 142 << 149 //--------------------------------------------- 143 G4cout << "============================= << 150 // Dump accumulated quantities for this RUN. 144 G4cout << " Number of event processed : << 151 // (Display only central region of x-y plane) 145 G4cout << "============================= << 152 // 0 ConcreteSD/DoseDeposit 146 << 153 //--------------------------------------------- 147 std::ofstream fileout; << 154 G4THitsMap<G4double>* DoseDeposit = re02Run->GetHitsMap(fSDName[i]+"/DoseDeposit"); 148 G4String fname = "dicom.out"; << 155 149 fileout.open(fname); << 156 G4cout << "=============================================================" <<G4endl; 150 G4cout << " opened file " << fname << " << 157 G4cout << " Number of event processed : "<< aRun->GetNumberOfEvent() << G4endl; 151 << 158 G4cout << "=============================================================" <<G4endl; 152 if (DoseDeposit && DoseDeposit->GetMap() << 159 153 std::ostream* myout = &G4cout; << 160 std::ofstream fileout; 154 PrintHeader(myout); << 161 G4String fname = "dicom.out"; 155 std::map<G4int, G4double*>::iterator i << 162 fileout.open(fname); 156 for (; itr != DoseDeposit->GetMap()->e << 163 G4cout << " opened file " << fname << " for dose output" << G4endl; 157 fileout << itr->first << " " << << 164 158 // G4cout << " " << itr->first << 165 159 // << " " << std::setprec << 166 if( DoseDeposit && DoseDeposit->GetMap()->size() != 0 ) { 160 // << *(itr->second)/CLHEP::g << 167 std::ostream *myout = &G4cout; 161 // << G4endl; << 168 PrintHeader(myout); >> 169 std::map<G4int,G4double*>::iterator itr = DoseDeposit->GetMap()->begin(); >> 170 for(; itr != DoseDeposit->GetMap()->end(); itr++) { >> 171 fileout << itr->first >> 172 << " " << *(itr->second)/CLHEP::gray >> 173 << G4endl; >> 174 G4cout << " " << itr->first >> 175 << " " << std::setprecision(6) << *(itr->second)/CLHEP::gray << " Gy" >> 176 << G4endl; >> 177 } >> 178 G4cout << "============================================="<<G4endl; >> 179 } else { >> 180 G4Exception("DicomRunAction", "000", JustWarning, >> 181 "DoseDeposit HitsMap is either a null pointer of the HitsMap was empty"); >> 182 } >> 183 fileout.close(); >> 184 G4cout << " closed file " << fname << " for dose output" << G4endl; >> 185 162 } 186 } 163 G4cout << "=========================== << 164 } << 165 else { << 166 G4Exception("DicomRunAction", "000", J << 167 "DoseDeposit HitsMap is ei << 168 } << 169 fileout.close(); << 170 G4cout << " closed file " << fname << " << 171 } 187 } 172 } << 173 188 174 G4cout << "Finished : End of Run Action " << << 189 G4cout << "Finished : End of Run Action " << aRun->GetRunID() << G4endl; >> 190 175 } 191 } 176 192 177 //....oooOO0OOooo........oooOO0OOooo........oo 193 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 178 void DicomRunAction::PrintHeader(std::ostream* << 194 void DicomRunAction::PrintHeader(std::ostream *out) 179 { 195 { 180 std::vector<G4String> vecScoreName; 196 std::vector<G4String> vecScoreName; 181 vecScoreName.push_back("DoseDeposit"); 197 vecScoreName.push_back("DoseDeposit"); 182 198 183 // head line 199 // head line 184 // 200 // 185 std::string vname; 201 std::string vname; 186 *out << std::setw(10) << "Voxel" << 202 *out << std::setw(10) << "Voxel" << " |"; 187 << " |"; << 203 for (std::vector<G4String>::iterator it = vecScoreName.begin(); 188 for (std::vector<G4String>::iterator it = ve << 204 it != vecScoreName.end(); it++) { 189 // vname = FillString((*it), << 205 //vname = FillString((*it), 190 // ' ', << 206 // ' ', 191 // FieldValue+1, << 207 // FieldValue+1, 192 // false); << 208 // false); 193 // *out << vname << '|'; << 209 // *out << vname << '|'; 194 *out << std::setw(fFieldValue) << (*it) << << 210 *out << std::setw(fFieldValue) << (*it) << " |"; 195 } 211 } 196 *out << G4endl; 212 *out << G4endl; 197 } 213 } 198 214 199 //....oooOO0OOooo........oooOO0OOooo........oo 215 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 200 std::string DicomRunAction::FillString(const s << 216 std::string DicomRunAction::FillString(const std::string &name, >> 217 char c, G4int n, G4bool back) 201 { 218 { 202 std::string fname(""); 219 std::string fname(""); 203 G4int k = G4int(n - name.size()); << 220 G4int k = n - name.size(); 204 if (k > 0) { 221 if (k > 0) { 205 if (back) { 222 if (back) { 206 fname = name; 223 fname = name; 207 fname += std::string(k, c); << 224 fname += std::string(k,c); 208 } 225 } 209 else { 226 else { 210 fname = std::string(k, c); << 227 fname = std::string(k,c); 211 fname += name; 228 fname += name; 212 } 229 } 213 } 230 } 214 else { 231 else { 215 fname = name; 232 fname = name; 216 } 233 } 217 return fname; 234 return fname; 218 } 235 } 219 236