Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/DICOM/src/DicomRun.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/DicomRun.cc (Version 11.3.0) and /examples/extended/medical/DICOM/src/DicomRun.cc (Version 10.6.p1)


  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 medical/DICOM/src/DicomRun.cc            27 /// \file medical/DICOM/src/DicomRun.cc
 28 /// \brief Implementation of the DicomRun clas     28 /// \brief Implementation of the DicomRun class
 29                                                    29 
 30 //============================================     30 //=====================================================================
 31 ///                                                31 ///
 32 ///  (Description)                                 32 ///  (Description)
 33 ///    DicomRun Class is for accumulating scor <<  33 ///    DicomRun Class is for accumulating scored quantities which is 
 34 ///  scored using G4MutiFunctionalDetector and     34 ///  scored using G4MutiFunctionalDetector and G4VPrimitiveScorer.
 35 ///  Accumulation is done using G4THitsMap obj     35 ///  Accumulation is done using G4THitsMap object.
 36 ///                                                36 ///
 37 ///    The constructor DicomRun(const std::vec     37 ///    The constructor DicomRun(const std::vector<G4String> mfdName)
 38 ///  needs a vector filled with MultiFunctiona     38 ///  needs a vector filled with MultiFunctionalDetector names which
 39 ///  was assigned at instantiation of MultiFun     39 ///  was assigned at instantiation of MultiFunctionalDetector(MFD).
 40 ///  Then DicomRun constructor automatically s     40 ///  Then DicomRun constructor automatically scans primitive scorers
 41 ///  in the MFD, and obtains collectionIDs of      41 ///  in the MFD, and obtains collectionIDs of all collections associated
 42 ///  to those primitive scorers. Futhermore, t <<  42 ///  to those primitive scorers. Futhermore, the G4THitsMap objects 
 43 ///  for accumulating during a RUN are automat     43 ///  for accumulating during a RUN are automatically created too.
 44 ///  (*) Collection Name is same as primitive      44 ///  (*) Collection Name is same as primitive scorer name.
 45 ///                                            <<  45 /// 
 46 ///    The resultant information is kept insid     46 ///    The resultant information is kept inside DicomRun objects as
 47 ///  data members.                                 47 ///  data members.
 48 ///  std::vector<G4String> fCollName;              48 ///  std::vector<G4String> fCollName;            // Collection Name,
 49 ///  std::vector<G4int> fCollID;                   49 ///  std::vector<G4int> fCollID;                 // Collection ID,
 50 ///  std::vector<G4THitsMap<G4double>*> fRunMa     50 ///  std::vector<G4THitsMap<G4double>*> fRunMap; // HitsMap for RUN.
 51 ///                                                51 ///
 52 ///  The resualtant HitsMap objects are obtain     52 ///  The resualtant HitsMap objects are obtain using access method,
 53 ///  GetHitsMap(..).                               53 ///  GetHitsMap(..).
 54 ///                                                54 ///
 55 //============================================     55 //=====================================================================
 56                                                    56 
 57 #include "DicomRun.hh"                             57 #include "DicomRun.hh"
                                                   >>  58 #include "G4SDManager.hh"
 58                                                    59 
 59 #include "G4MultiFunctionalDetector.hh"            60 #include "G4MultiFunctionalDetector.hh"
 60 #include "G4SDManager.hh"                      << 
 61 #include "G4VPrimitiveScorer.hh"                   61 #include "G4VPrimitiveScorer.hh"
 62                                                    62 
 63 //....oooOO0OOooo........oooOO0OOooo........oo     63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 64 //                                                 64 //
 65 //  Constructor.                               <<  65 //  Constructor. 
 66 DicomRun::DicomRun() : G4Run() {}              <<  66 DicomRun::DicomRun()
                                                   >>  67 : G4Run()
                                                   >>  68 { }
 67                                                    69 
 68 //....oooOO0OOooo........oooOO0OOooo........oo     70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 69 //                                                 71 //
 70 //  Constructor.                                   72 //  Constructor.
 71 //   (The vector of MultiFunctionalDetector na     73 //   (The vector of MultiFunctionalDetector name has to given.)
 72 DicomRun::DicomRun(const std::vector<G4String> <<  74 DicomRun::DicomRun(const std::vector<G4String> mfdName): G4Run()
 73 {                                                  75 {
 74   ConstructMFD(mfdName);                       <<  76     ConstructMFD(mfdName);
 75 }                                                  77 }
 76                                                    78 
 77 //....oooOO0OOooo........oooOO0OOooo........oo     79 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 78 //                                                 80 //
 79 // Destructor                                      81 // Destructor
 80 //    clear all data members.                      82 //    clear all data members.
 81 DicomRun::~DicomRun()                              83 DicomRun::~DicomRun()
 82 {                                                  84 {
 83   //--- Clear HitsMap for RUN                      85   //--- Clear HitsMap for RUN
 84   size_t Nmap = fRunMap.size();                    86   size_t Nmap = fRunMap.size();
 85   for (size_t i = 0; i < Nmap; ++i) {          <<  87   for ( size_t i = 0; i < Nmap; ++i)
 86     if (fRunMap[i]) fRunMap[i]->clear();       <<  88   {
                                                   >>  89     if(fRunMap[i] ) fRunMap[i]->clear();
 87   }                                                90   }
 88   fCollName.clear();                               91   fCollName.clear();
 89   fCollID.clear();                                 92   fCollID.clear();
 90   fRunMap.clear();                                 93   fRunMap.clear();
 91 }                                                  94 }
 92                                                    95 
 93 //....oooOO0OOooo........oooOO0OOooo........oo     96 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 94 //                                                 97 //
 95 // Destructor                                      98 // Destructor
 96 //    clear all data members.                      99 //    clear all data members.
 97 void DicomRun::ConstructMFD(const std::vector<    100 void DicomRun::ConstructMFD(const std::vector<G4String>& mfdName)
 98 {                                                 101 {
                                                   >> 102 
 99   G4SDManager* SDman = G4SDManager::GetSDMpoin    103   G4SDManager* SDman = G4SDManager::GetSDMpointer();
100   //==========================================    104   //=================================================
101   //  Initalize RunMaps for accumulation.         105   //  Initalize RunMaps for accumulation.
102   //  Get CollectionIDs for HitCollections.       106   //  Get CollectionIDs for HitCollections.
103   //==========================================    107   //=================================================
104   size_t Nmfd = mfdName.size();                   108   size_t Nmfd = mfdName.size();
105   for (size_t idet = 0; idet < Nmfd; ++idet)   << 109   for ( size_t idet = 0; idet < Nmfd ; ++idet)  // Loop for all MFD.
106   {                                               110   {
107     G4String detName = mfdName[idet];             111     G4String detName = mfdName[idet];
108     //--- Seek and Obtain MFD objects from SDm    112     //--- Seek and Obtain MFD objects from SDmanager.
109     G4MultiFunctionalDetector* mfd =              113     G4MultiFunctionalDetector* mfd =
110       (G4MultiFunctionalDetector*)(SDman->Find    114       (G4MultiFunctionalDetector*)(SDman->FindSensitiveDetector(detName));
111     //                                            115     //
112     if (mfd) {                                 << 116     if ( mfd )
                                                   >> 117     {
113       //--- Loop over the registered primitive    118       //--- Loop over the registered primitive scorers.
114       for (G4int icol = 0; icol < mfd->GetNumb << 119       for (G4int icol = 0; icol < mfd->GetNumberOfPrimitives(); ++icol)
                                                   >> 120       {
115         // Get Primitive Scorer object.           121         // Get Primitive Scorer object.
116         G4VPrimitiveScorer* scorer = mfd->GetP    122         G4VPrimitiveScorer* scorer = mfd->GetPrimitive(icol);
117         // collection name and collectionID fo    123         // collection name and collectionID for HitsCollection,
118         // where type of HitsCollection is G4T << 124         // where type of HitsCollection is G4THitsMap in case 
119         // of primitive scorer.                   125         // of primitive scorer.
120         // The collection name is given by <MF << 126         // The collection name is given by <MFD name>/<Primitive 
121         // Scorer name>.                          127         // Scorer name>.
122         G4String collectionName = scorer->GetN    128         G4String collectionName = scorer->GetName();
123         G4String fullCollectionName = detName  << 129         G4String fullCollectionName = detName+"/"+collectionName;
124         G4int collectionID = SDman->GetCollect << 130         G4int    collectionID = SDman->GetCollectionID(fullCollectionName);
125         //                                        131         //
126         if (collectionID >= 0) {               << 132         if ( collectionID >= 0 )
127           G4cout << "++ " << fullCollectionNam << 133         {
                                                   >> 134           G4cout << "++ "<<fullCollectionName<< " id " << collectionID 
                                                   >> 135                  << G4endl;
128           // Store obtained HitsCollection inf    136           // Store obtained HitsCollection information into data members.
129           // And, creates new G4THitsMap for a    137           // And, creates new G4THitsMap for accumulating quantities during RUN.
130           fCollName.push_back(fullCollectionNa    138           fCollName.push_back(fullCollectionName);
131           fCollID.push_back(collectionID);        139           fCollID.push_back(collectionID);
132           fRunMap.push_back(new G4THitsMap<G4d << 140           fRunMap.push_back(new G4THitsMap<G4double>(detName,collectionName));
133         }                                         141         }
134         else {                                 << 142         else
135           G4cout << "** collection " << fullCo << 143         {
                                                   >> 144           G4cout << "** collection " << fullCollectionName << " not found. "
                                                   >> 145                  <<G4endl;
136         }                                         146         }
137       }                                           147       }
138     }                                             148     }
139   }                                               149   }
140 }                                                 150 }
141                                                   151 
142 //....oooOO0OOooo........oooOO0OOooo........oo    152 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
143 //                                                153 //
144 //  RecordEvent is called at end of event.        154 //  RecordEvent is called at end of event.
145 //  For scoring purpose, the resultant quantit    155 //  For scoring purpose, the resultant quantity in a event,
146 //  is accumulated during a Run.                  156 //  is accumulated during a Run.
147 void DicomRun::RecordEvent(const G4Event* aEve    157 void DicomRun::RecordEvent(const G4Event* aEvent)
148 {                                                 158 {
149   // disable below because this is done in G4R    159   // disable below because this is done in G4Run::RecordEvent(aEvent);
150   // numberOfEvent++;  // This is an original  << 160   //numberOfEvent++;  // This is an original line.
151                                                << 161   
152   // G4cout << "Dicom Run :: Recording event " << 162   //G4cout << "Dicom Run :: Recording event " << aEvent->GetEventID() 
153   //<< "..." << G4endl;                           163   //<< "..." << G4endl;
154   //=============================                 164   //=============================
155   //  HitsCollection of This Event             << 165   // HitsCollection of This Event
156   //============================                  166   //============================
157   G4HCofThisEvent* HCE = aEvent->GetHCofThisEv    167   G4HCofThisEvent* HCE = aEvent->GetHCofThisEvent();
158   if (!HCE) return;                               168   if (!HCE) return;
159                                                << 169   
160   //==========================================    170   //=======================================================
161   // Sum up HitsMap of this Event  into HitsMa    171   // Sum up HitsMap of this Event  into HitsMap of this RUN
162   //==========================================    172   //=======================================================
163   size_t Ncol = fCollID.size();                   173   size_t Ncol = fCollID.size();
164   for (size_t i = 0; i < Ncol; ++i)  // Loop o << 174   for ( size_t i = 0; i < Ncol ; ++i )  // Loop over HitsCollection
165   {                                               175   {
166     G4THitsMap<G4double>* EvtMap = nullptr;       176     G4THitsMap<G4double>* EvtMap = nullptr;
167     if (fCollID[i] >= 0)  // Collection is att << 177     if ( fCollID[i] >= 0 )           // Collection is attached to HCE
168     {                                             178     {
169       EvtMap = static_cast<G4THitsMap<G4double    179       EvtMap = static_cast<G4THitsMap<G4double>*>(HCE->GetHC(fCollID[i]));
170     }                                             180     }
171     else {                                     << 181     else
172       G4cout << " Error EvtMap Not Found " <<  << 182     {
                                                   >> 183       G4cout <<" Error EvtMap Not Found "<< i << G4endl;
173     }                                             184     }
174     if (EvtMap) {                              << 185     if ( EvtMap )
                                                   >> 186     {
175       //=== Sum up HitsMap of this event to Hi    187       //=== Sum up HitsMap of this event to HitsMap of RUN.===
176       *fRunMap[i] += *EvtMap;                     188       *fRunMap[i] += *EvtMap;
177       // G4cout << "Summing EvtMap into RunMap << 189       //G4cout << "Summing EvtMap into RunMap at " << i << "..." << G4endl;
178       //======================================    190       //======================================================
179     }  // else {                               << 191     } //else {
180     // G4cout << "Null pointer to EvtMap at "  << 192     //G4cout << "Null pointer to EvtMap at " << i << "..." << G4endl;
181     // }                                       << 193     //}
182   }                                               194   }
183                                                << 195   
184   G4Run::RecordEvent(aEvent);                     196   G4Run::RecordEvent(aEvent);
                                                   >> 197   
185 }                                                 198 }
186                                                   199 
187 //....oooOO0OOooo........oooOO0OOooo........oo    200 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
188 // Merge hits map from threads                    201 // Merge hits map from threads
189 void DicomRun::Merge(const G4Run* aRun)           202 void DicomRun::Merge(const G4Run* aRun)
190 {                                                 203 {
191   const DicomRun* localRun = static_cast<const    204   const DicomRun* localRun = static_cast<const DicomRun*>(aRun);
192   Copy(fCollName, localRun->fCollName);           205   Copy(fCollName, localRun->fCollName);
193   Copy(fCollID, localRun->fCollID);               206   Copy(fCollID, localRun->fCollID);
194   size_t ncopies = Copy(fRunMap, localRun->fRu    207   size_t ncopies = Copy(fRunMap, localRun->fRunMap);
195   // copy function returns the fRunMap size if    208   // copy function returns the fRunMap size if all data is copied
196   // so this loop isn't executed the first tim    209   // so this loop isn't executed the first time around
197   for (size_t i = ncopies; i < fRunMap.size(); << 210   for(size_t i = ncopies; i < fRunMap.size(); ++i)
                                                   >> 211   {
198     *fRunMap[i] += *localRun->fRunMap[i];         212     *fRunMap[i] += *localRun->fRunMap[i];
199   }                                               213   }
200   G4Run::Merge(aRun);                             214   G4Run::Merge(aRun);
201 }                                                 215 }
202                                                   216 
203 //....oooOO0OOooo........oooOO0OOooo........oo    217 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
204 //============================================    218 //=================================================================
205 //  Access method for HitsMap of the RUN          219 //  Access method for HitsMap of the RUN
206 //                                                220 //
207 //-----                                           221 //-----
208 // Access HitsMap.                                222 // Access HitsMap.
209 //  By  MultiFunctionalDetector name and Colle    223 //  By  MultiFunctionalDetector name and Collection Name.
210 G4THitsMap<G4double>* DicomRun::GetHitsMap(con << 224 G4THitsMap<G4double>* DicomRun::GetHitsMap(const G4String& detName,
                                                   >> 225                                            const G4String& colName) const
211 {                                                 226 {
212   G4String fullName = detName + "/" + colName; << 227   G4String fullName = detName+"/"+colName;
213   return GetHitsMap(fullName);                    228   return GetHitsMap(fullName);
214 }                                                 229 }
215                                                   230 
216 //....oooOO0OOooo........oooOO0OOooo........oo    231 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
217 // Access HitsMap.                                232 // Access HitsMap.
218 //  By full description of collection name, th    233 //  By full description of collection name, that is
219 //    <MultiFunctional Detector Name>/<Primiti    234 //    <MultiFunctional Detector Name>/<Primitive Scorer Name>
220 G4THitsMap<G4double>* DicomRun::GetHitsMap(con    235 G4THitsMap<G4double>* DicomRun::GetHitsMap(const G4String& fullName) const
221 {                                                 236 {
222   // G4THitsMap<G4double>* hitsmap = 0;        << 237   
223                                                << 238   //G4THitsMap<G4double>* hitsmap = 0;
                                                   >> 239   
224   size_t Ncol = fCollName.size();                 240   size_t Ncol = fCollName.size();
225   for (size_t i = 0; i < Ncol; ++i) {          << 241   for ( size_t i = 0; i < Ncol; ++i)
226     if (fCollName[i] == fullName) {            << 242   {
                                                   >> 243     if ( fCollName[i] == fullName )
                                                   >> 244     {
227       return fRunMap[i];                          245       return fRunMap[i];
228       // if(hitsmap) { *hitsmap += *fRunMap[i] << 246             //if(hitsmap) { *hitsmap += *fRunMap[i]; }
229       // if(!hitsmap) { hitsmap = fRunMap[i];  << 247             //if(!hitsmap) { hitsmap = fRunMap[i]; }
230     }                                             248     }
231   }                                               249   }
232                                                << 250   
233   G4Exception("DicomRun", fullName.c_str(), Ju    251   G4Exception("DicomRun", fullName.c_str(), JustWarning,
234               "GetHitsMap failed to locate the    252               "GetHitsMap failed to locate the requested HitsMap");
235   return nullptr;                                 253   return nullptr;
236 }                                                 254 }
237                                                   255