Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/DICOM/src/DicomRunAction.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/extended/medical/DICOM/src/DicomRunAction.cc (Version 11.3.0) and /examples/extended/medical/DICOM/src/DicomRunAction.cc (Version 10.1)


  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 84839 2014-10-21 13:44:55Z 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();
 89                                                <<  94   if (nofEvents == 0) return;
 90   static G4double local_total_dose = 0;        <<  95   //print
 91   G4double total_dose = 0;                     <<  96   //
 92                                                <<  97   
                                                   >>  98   static double local_total_dose = 0;
                                                   >>  99   double total_dose = 0;
                                                   >> 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 = 
103                                                << 111       reRun->GetHitsMap(fSDName[i]+"/DoseDeposit");
104     if (DoseDeposit && DoseDeposit->GetMap()-> << 112     
105       std::map<G4int, G4double*>::iterator itr << 113     if( DoseDeposit && DoseDeposit->GetMap()->size() != 0 ) {
106       for (; itr != DoseDeposit->GetMap()->end << 114       std::map<G4int,G4double*>::iterator itr = 
107         if (!IsMaster()) {                     << 115         DoseDeposit->GetMap()->begin();
108           local_total_dose += *(itr->second);  << 116       for(; itr != DoseDeposit->GetMap()->end(); itr++) {
109         }                                      << 117         if(!IsMaster()) { local_total_dose += *(itr->second); }
110         total_dose += *(itr->second);             118         total_dose += *(itr->second);
111       }                                           119       }
112     }                                             120     }
                                                   >> 121     
113   }                                               122   }
114                                                << 123   
115   if (IsMaster()) {                            << 124   if (IsMaster())
116     G4cout << "--------------------End of Glob << 125     {
117     G4cout << " The run was " << nofEvents <<  << 126       G4cout
118     G4cout << "LOCAL TOTAL DOSE : \t" << local << 127         << "\n--------------------End of Global Run-----------------------"
119     G4cout << "      TOTAL DOSE : \t" << total << 128         << " \n The run was " << nofEvents << " events " << G4endl;
120   }                                            << 129       G4cout << "LOCAL TOTAL DOSE : \t" << local_total_dose/gray 
121   else {                                       << 130              << " Gy" << std::endl;
122     G4cout << "--------------------End of Loca << 131       G4cout << "      TOTAL DOSE : \t" << total_dose/gray 
123     G4cout << " The run was " << nofEvents <<  << 132              << " Gy" << std::endl;
124     G4cout << "LOCAL TOTAL DOSE : \t" << local << 133       
125     G4cout << "      TOTAL DOSE : \t" << total << 134     }
126   }                                            << 135   else
127                                                << 136     {
128   if (IsMaster()) {                            << 137       G4cout
                                                   >> 138         << "\n--------------------End of Local Run------------------------"
                                                   >> 139         << " \n The run was " << nofEvents << G4endl;
                                                   >> 140       G4cout << "LOCAL TOTAL DOSE : \t" << local_total_dose/gray 
                                                   >> 141              << " Gy" << std::endl;
                                                   >> 142       G4cout << "      TOTAL DOSE : \t" << total_dose/gray 
                                                   >> 143              << " Gy" << std::endl;
                                                   >> 144       
                                                   >> 145     }
                                                   >> 146   
                                                   >> 147   if(IsMaster()) {
129     G4cout << " ###### EndOfRunAction ###### "    148     G4cout << " ###### EndOfRunAction ###### " << G4endl;
130     //- DicomRun object.                          149     //- DicomRun object.
131     const DicomRun* re02Run = static_cast<cons    150     const DicomRun* re02Run = static_cast<const DicomRun*>(aRun);
132     //--- Dump all scored quantities involved     151     //--- Dump all scored quantities involved in DicomRun.
133                                                << 152     
134     for (G4int i = 0; i < (G4int)fSDName.size( << 153     for ( G4int i = 0; i < (G4int)fSDName.size(); i++ ){
135       //                                          154       //
136       //--------------------------------------    155       //---------------------------------------------
137       // Dump accumulated quantities for this     156       // Dump accumulated quantities for this RUN.
138       //  (Display only central region of x-y     157       //  (Display only central region of x-y plane)
139       //      0       ConcreteSD/DoseDeposit      158       //      0       ConcreteSD/DoseDeposit
140       //--------------------------------------    159       //---------------------------------------------
141       G4THitsMap<G4double>* DoseDeposit = re02 << 160       G4THitsMap<G4double>* DoseDeposit = 
142                                                << 161         re02Run->GetHitsMap(fSDName[i]+"/DoseDeposit");
143       G4cout << "============================= << 162       
144       G4cout << " Number of event processed :  << 163       G4cout << "============================================================="
145       G4cout << "============================= << 164              <<G4endl;
146                                                << 165       G4cout << " Number of event processed : "
                                                   >> 166              << aRun->GetNumberOfEvent() << G4endl;
                                                   >> 167       G4cout << "============================================================="
                                                   >> 168              <<G4endl;
                                                   >> 169       
147       std::ofstream fileout;                      170       std::ofstream fileout;
148       G4String fname = "dicom.out";               171       G4String fname = "dicom.out";
149       fileout.open(fname);                        172       fileout.open(fname);
150       G4cout << " opened file " << fname << "     173       G4cout << " opened file " << fname << " for dose output" << G4endl;
151                                                << 174       
152       if (DoseDeposit && DoseDeposit->GetMap() << 175       
153         std::ostream* myout = &G4cout;         << 176       if( DoseDeposit && DoseDeposit->GetMap()->size() != 0 ) {
                                                   >> 177         std::ostream *myout = &G4cout;
154         PrintHeader(myout);                       178         PrintHeader(myout);
155         std::map<G4int, G4double*>::iterator i << 179         std::map<G4int,G4double*>::iterator itr =DoseDeposit->GetMap()->begin();
156         for (; itr != DoseDeposit->GetMap()->e << 180         for(; itr != DoseDeposit->GetMap()->end(); itr++) {
157           fileout << itr->first << "     " <<  << 181           fileout <<  itr->first
158           // G4cout << "    " << itr->first    << 182                   << "     "  << *(itr->second)/CLHEP::gray
159           //        << "     " << std::setprec << 183                   << G4endl;
160           //        << *(itr->second)/CLHEP::g << 184           G4cout << "    " << itr->first
161           //        << G4endl;                 << 185                  << "     " << std::setprecision(6) 
                                                   >> 186                  << *(itr->second)/CLHEP::gray << " Gy"
                                                   >> 187                  << G4endl;
162         }                                         188         }
163         G4cout << "=========================== << 189         G4cout << "============================================="<<G4endl;
164       }                                        << 190             } else {
165       else {                                   << 191         G4Exception("DicomRunAction", "000", JustWarning, 
166         G4Exception("DicomRunAction", "000", J << 192       "DoseDeposit HitsMap is either a null pointer of the HitsMap was empty");
167                     "DoseDeposit HitsMap is ei << 
168       }                                           193       }
169       fileout.close();                            194       fileout.close();
170       G4cout << " closed file " << fname << "     195       G4cout << " closed file " << fname << " for dose output" << G4endl;
                                                   >> 196       
171     }                                             197     }
172   }                                               198   }
173                                                << 199   
174   G4cout << "Finished : End of Run Action " <<    200   G4cout << "Finished : End of Run Action " << aRun->GetRunID() << G4endl;
                                                   >> 201   
175 }                                                 202 }
176                                                   203 
177 //....oooOO0OOooo........oooOO0OOooo........oo    204 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
178 void DicomRunAction::PrintHeader(std::ostream* << 205 void DicomRunAction::PrintHeader(std::ostream *out)
179 {                                                 206 {
180   std::vector<G4String> vecScoreName;             207   std::vector<G4String> vecScoreName;
181   vecScoreName.push_back("DoseDeposit");          208   vecScoreName.push_back("DoseDeposit");
182                                                << 209   
183   // head line                                    210   // head line
184   //                                              211   //
185   std::string vname;                              212   std::string vname;
186   *out << std::setw(10) << "Voxel"             << 213   *out << std::setw(10) << "Voxel" << " |";
187        << " |";                                << 214   for (std::vector<G4String>::iterator it = vecScoreName.begin();
188   for (std::vector<G4String>::iterator it = ve << 215        it != vecScoreName.end(); it++) {
189     // vname = FillString((*it),               << 216     //vname = FillString((*it),
190     //                        ' ',             << 217     //                       ' ',
191     //                        FieldValue+1,    << 218     //                       FieldValue+1,
192     //                        false);          << 219     //                       false);
193     //     *out << vname << '|';               << 220     //    *out << vname << '|';
194     *out << std::setw(fFieldValue) << (*it) <<    221     *out << std::setw(fFieldValue) << (*it) << "  |";
195   }                                               222   }
196   *out << G4endl;                                 223   *out << G4endl;
197 }                                                 224 }
198                                                   225 
199 //....oooOO0OOooo........oooOO0OOooo........oo    226 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
200 std::string DicomRunAction::FillString(const s << 227 std::string DicomRunAction::FillString(const std::string &name,
                                                   >> 228                                        char c, G4int n, G4bool back)
201 {                                                 229 {
202   std::string fname("");                          230   std::string fname("");
203   G4int k = G4int(n - name.size());            << 231   G4int k = n - name.size();
204   if (k > 0) {                                    232   if (k > 0) {
205     if (back) {                                   233     if (back) {
206       fname = name;                               234       fname = name;
207       fname += std::string(k, c);              << 235       fname += std::string(k,c);
208     }                                             236     }
209     else {                                        237     else {
210       fname = std::string(k, c);               << 238       fname = std::string(k,c);
211       fname += name;                              239       fname += name;
212     }                                             240     }
213   }                                               241   }
214   else {                                          242   else {
215     fname = name;                                 243     fname = name;
216   }                                               244   }
217   return fname;                                   245   return fname;
218 }                                                 246 }
219                                                   247