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 // G4RegionStore implementation << 27 // 26 // 28 // 18.09.02, G.Cosmo - Initial version << 27 // $Id: G4RegionStore.cc,v 1.14 2008/07/10 09:46:01 gcosmo Exp $ >> 28 // GEANT4 tag $Name: geant4-09-03-patch-02 $ >> 29 // >> 30 // G4RegionStore >> 31 // >> 32 // Implementation for singleton container >> 33 // >> 34 // History: >> 35 // 18.09.02 G.Cosmo Initial version 29 // ------------------------------------------- 36 // -------------------------------------------------------------------- 30 37 31 #include "G4Region.hh" 38 #include "G4Region.hh" 32 #include "G4RegionStore.hh" 39 #include "G4RegionStore.hh" 33 #include "G4GeometryManager.hh" 40 #include "G4GeometryManager.hh" 34 #include "G4VPhysicalVolume.hh" 41 #include "G4VPhysicalVolume.hh" 35 #include "G4PhysicalVolumeStore.hh" 42 #include "G4PhysicalVolumeStore.hh" 36 43 37 #include "G4ios.hh" 44 #include "G4ios.hh" 38 #include "G4AutoLock.hh" << 39 << 40 namespace << 41 { << 42 G4Mutex mapMutex = G4MUTEX_INITIALIZER; << 43 } << 44 45 45 // ******************************************* 46 // *************************************************************************** 46 // Static class variables 47 // Static class variables 47 // ******************************************* 48 // *************************************************************************** 48 // 49 // 49 G4RegionStore* G4RegionStore::fgInstance = nul << 50 G4RegionStore* G4RegionStore::fgInstance = 0; 50 G4ThreadLocal G4VStoreNotifier* G4RegionStore: << 51 G4VStoreNotifier* G4RegionStore::fgNotifier = 0; 51 G4ThreadLocal G4bool G4RegionStore::locked = f << 52 G4bool G4RegionStore::locked = false; 52 53 53 // ******************************************* 54 // *************************************************************************** 54 // Protected constructor: Construct underlying 55 // Protected constructor: Construct underlying container with 55 // initial size of 20 entries 56 // initial size of 20 entries 56 // ******************************************* 57 // *************************************************************************** 57 // 58 // 58 G4RegionStore::G4RegionStore() 59 G4RegionStore::G4RegionStore() 59 << 60 : std::vector<G4Region*>() 60 { 61 { 61 reserve(20); 62 reserve(20); 62 } 63 } 63 64 64 // ******************************************* 65 // *************************************************************************** 65 // Destructor 66 // Destructor 66 // ******************************************* 67 // *************************************************************************** 67 // 68 // 68 G4RegionStore::~G4RegionStore() 69 G4RegionStore::~G4RegionStore() 69 { 70 { 70 Clean(); // Delete all regions in the store << 71 Clean(); 71 G4Region::Clean(); // Delete allocated sub- << 72 } 72 } 73 73 74 // ******************************************* 74 // *************************************************************************** 75 // Delete all regions from the store except fo 75 // Delete all regions from the store except for the world region 76 // ******************************************* 76 // *************************************************************************** 77 // 77 // 78 void G4RegionStore::Clean() 78 void G4RegionStore::Clean() 79 { 79 { 80 // Do nothing if geometry is closed 80 // Do nothing if geometry is closed 81 // 81 // 82 if (G4GeometryManager::GetInstance()->IsGeom 82 if (G4GeometryManager::GetInstance()->IsGeometryClosed()) 83 { 83 { 84 G4cout << "WARNING - Attempt to delete the 84 G4cout << "WARNING - Attempt to delete the region store" 85 << " while geometry closed !" << G4 85 << " while geometry closed !" << G4endl; 86 return; 86 return; 87 } 87 } 88 88 89 // Locks store for deletion of regions. De-r 89 // Locks store for deletion of regions. De-registration will be 90 // performed at this stage. G4Regions will n 90 // performed at this stage. G4Regions will not de-register themselves. 91 // 91 // 92 locked = true; 92 locked = true; 93 93 >> 94 size_t i=0; 94 G4RegionStore* store = GetInstance(); 95 G4RegionStore* store = GetInstance(); 95 96 96 for(auto pos=store->cbegin(); pos!=store->ce << 97 #ifdef G4GEOMETRY_VOXELDEBUG >> 98 G4cout << "Deleting Regions ... "; >> 99 #endif >> 100 >> 101 for(iterator pos=store->begin(); pos!=store->end(); ++pos) 97 { 102 { 98 if (fgNotifier != nullptr) { fgNotifier->N << 103 if (fgNotifier) { fgNotifier->NotifyDeRegistration(); } 99 delete *pos; << 104 if (*pos) { delete *pos; } >> 105 i++; 100 } 106 } 101 107 102 store->bmap.clear(); store->mvalid = false; << 108 #ifdef G4GEOMETRY_VOXELDEBUG >> 109 if (store->size() < i-1) >> 110 { G4cout << "No regions deleted. Already deleted by user ?" << G4endl; } >> 111 else >> 112 { G4cout << i-1 << " regions deleted !" << G4endl; } >> 113 #endif >> 114 103 locked = false; 115 locked = false; 104 store->clear(); 116 store->clear(); 105 } 117 } 106 118 107 // ******************************************* 119 // *************************************************************************** 108 // Associate user notifier to the store 120 // Associate user notifier to the store 109 // ******************************************* 121 // *************************************************************************** 110 // 122 // 111 void G4RegionStore::SetNotifier(G4VStoreNotifi 123 void G4RegionStore::SetNotifier(G4VStoreNotifier* pNotifier) 112 { 124 { 113 GetInstance(); 125 GetInstance(); 114 fgNotifier = pNotifier; 126 fgNotifier = pNotifier; 115 } 127 } 116 128 117 // ******************************************* 129 // *************************************************************************** 118 // Bring contents of internal map up to date a << 119 // ******************************************* << 120 // << 121 void G4RegionStore::UpdateMap() << 122 { << 123 G4AutoLock l(&mapMutex); // to avoid thread << 124 if (mvalid) return; << 125 bmap.clear(); << 126 for(auto pos=GetInstance()->cbegin(); pos!=G << 127 { << 128 const G4String& reg_name = (*pos)->GetName << 129 auto it = bmap.find(reg_name); << 130 if (it != bmap.cend()) << 131 { << 132 it->second.push_back(*pos); << 133 } << 134 else << 135 { << 136 std::vector<G4Region*> reg_vec { *pos }; << 137 bmap.insert(std::make_pair(reg_name, reg << 138 } << 139 } << 140 mvalid = true; << 141 l.unlock(); << 142 } << 143 << 144 // ******************************************* << 145 // Add Region to container 130 // Add Region to container 146 // ******************************************* 131 // *************************************************************************** 147 // 132 // 148 void G4RegionStore::Register(G4Region* pRegion 133 void G4RegionStore::Register(G4Region* pRegion) 149 { 134 { 150 G4RegionStore* store = GetInstance(); << 135 GetInstance()->push_back(pRegion); 151 store->push_back(pRegion); << 136 if (fgNotifier) { fgNotifier->NotifyRegistration(); } 152 const G4String& reg_name = pRegion->GetName( << 153 auto it = store->bmap.find(reg_name); << 154 if (it != store->bmap.cend()) << 155 { << 156 it->second.push_back(pRegion); << 157 } << 158 else << 159 { << 160 std::vector<G4Region*> reg_vec { pRegion } << 161 store->bmap.insert(std::make_pair(reg_name << 162 } << 163 if (fgNotifier != nullptr) { fgNotifier->Not << 164 store->mvalid = true; << 165 } 137 } 166 138 167 // ******************************************* 139 // *************************************************************************** 168 // Remove Region from container 140 // Remove Region from container 169 // ******************************************* 141 // *************************************************************************** 170 // 142 // 171 void G4RegionStore::DeRegister(G4Region* pRegi 143 void G4RegionStore::DeRegister(G4Region* pRegion) 172 { 144 { 173 G4RegionStore* store = GetInstance(); << 174 if (!locked) // Do not de-register if loc 145 if (!locked) // Do not de-register if locked ! 175 { 146 { 176 if (fgNotifier != nullptr) { fgNotifier-> << 147 if (fgNotifier) { fgNotifier->NotifyDeRegistration(); } 177 for (auto i=store->cbegin(); i!=store->cen << 148 for (iterator i=GetInstance()->begin(); i!=GetInstance()->end(); i++) 178 { 149 { 179 if (**i==*pRegion) 150 if (**i==*pRegion) 180 { 151 { 181 store->erase(i); << 152 GetInstance()->erase(i); 182 break; 153 break; 183 } 154 } 184 } 155 } 185 const G4String& reg_name = pRegion->GetNam << 186 auto it = store->bmap.find(reg_name); << 187 if (it != store->bmap.cend()) << 188 { << 189 if (it->second.size() > 1) << 190 { << 191 for (auto i=it->second.cbegin(); i!=it << 192 { << 193 if (**i==*pRegion) << 194 { << 195 it->second.erase(i); << 196 break; << 197 } << 198 } << 199 } << 200 else << 201 { << 202 store->bmap.erase(it); << 203 } << 204 } << 205 } 156 } 206 } 157 } 207 158 208 // ******************************************* 159 // *************************************************************************** 209 // Return ptr to Store, setting if necessary 160 // Return ptr to Store, setting if necessary 210 // ******************************************* 161 // *************************************************************************** 211 // 162 // 212 G4RegionStore* G4RegionStore::GetInstance() 163 G4RegionStore* G4RegionStore::GetInstance() 213 { 164 { 214 static G4RegionStore worldStore; 165 static G4RegionStore worldStore; 215 if (fgInstance == nullptr) << 166 if (!fgInstance) 216 { 167 { 217 fgInstance = &worldStore; 168 fgInstance = &worldStore; 218 } 169 } 219 return fgInstance; 170 return fgInstance; 220 } 171 } 221 172 222 // ******************************************* 173 // *************************************************************************** 223 // Loops through all regions to verify if a re 174 // Loops through all regions to verify if a region has been modified. 224 // It returns TRUE if just one region is modif 175 // It returns TRUE if just one region is modified. 225 // ******************************************* 176 // *************************************************************************** 226 // 177 // 227 G4bool G4RegionStore::IsModified() const 178 G4bool G4RegionStore::IsModified() const 228 { 179 { 229 for (auto i=GetInstance()->cbegin(); i!=GetI << 180 for (iterator i=GetInstance()->begin(); i!=GetInstance()->end(); i++) 230 { 181 { 231 if ((*i)->IsModified()) { return true; } 182 if ((*i)->IsModified()) { return true; } 232 } 183 } 233 return false; 184 return false; 234 } 185 } 235 186 236 // ******************************************* 187 // *************************************************************************** 237 // Loops through all regions to reset flag for 188 // Loops through all regions to reset flag for modification to FALSE. 238 // Used by the run manager to notify that the 189 // Used by the run manager to notify that the physics table has been updated. 239 // ******************************************* 190 // *************************************************************************** 240 // 191 // 241 void G4RegionStore::ResetRegionModified() 192 void G4RegionStore::ResetRegionModified() 242 { 193 { 243 for (auto i=GetInstance()->cbegin(); i!=GetI << 194 for (iterator i=GetInstance()->begin(); i!=GetInstance()->end(); i++) 244 { 195 { 245 (*i)->RegionModified(false); 196 (*i)->RegionModified(false); 246 } 197 } 247 } 198 } 248 199 249 // ******************************************* 200 // *************************************************************************** 250 // Forces recomputation of material lists in a 201 // Forces recomputation of material lists in all regions in the store. 251 // ******************************************* 202 // *************************************************************************** 252 // 203 // 253 void G4RegionStore::UpdateMaterialList(G4VPhys 204 void G4RegionStore::UpdateMaterialList(G4VPhysicalVolume* currentWorld) 254 { 205 { 255 for (auto i=GetInstance()->cbegin(); i!=GetI << 206 for (iterator i=GetInstance()->begin(); i!=GetInstance()->end(); i++) 256 { 207 { 257 if((*i)->IsInMassGeometry() || (*i)->IsInP << 208 if((*i)->GetWorldPhysical()==currentWorld) { (*i)->UpdateMaterialList(); } 258 || (currentWor << 259 { (*i)->UpdateMaterialList(); } << 260 } 209 } 261 } 210 } 262 211 263 // ******************************************* 212 // *************************************************************************** 264 // Returns a region through its name specifica 213 // Returns a region through its name specification. 265 // ******************************************* 214 // *************************************************************************** 266 // 215 // 267 G4Region* G4RegionStore::GetRegion(const G4Str 216 G4Region* G4RegionStore::GetRegion(const G4String& name, G4bool verbose) const 268 { 217 { 269 G4RegionStore* store = GetInstance(); << 218 for (iterator i=GetInstance()->begin(); i!=GetInstance()->end(); i++) 270 if (!store->mvalid) { store->UpdateMap(); } << 271 auto pos = store->bmap.find(name); << 272 if(pos != store->bmap.cend()) << 273 { 219 { 274 if ((verbose) && (pos->second.size()>1)) << 220 if ((*i)->GetName() == name) { return *i; } 275 { << 276 std::ostringstream message; << 277 message << "There exists more than ONE r << 278 << name << "!" << G4endl << 279 << "Returning the first found."; << 280 G4Exception("G4RegionStore::GetSolid()", << 281 "GeomMgt1001", JustWarning, << 282 } << 283 return pos->second[0]; << 284 } 221 } 285 if (verbose) 222 if (verbose) 286 { 223 { 287 std::ostringstream message; << 224 G4cerr << "ERROR - G4RegionStore::GetRegion()" << G4endl 288 message << "Region NOT found in store !" < << 225 << " Region " << name << " NOT found in store !" << G4endl 289 << " Region " << name << " << 226 << " Returning NULL pointer." << G4endl; 290 << " Returning NULL pointer << 227 G4Exception("G4RegionStore::GetRegion()", "InvalidQuery", 291 G4Exception("G4RegionStore::GetRegion()", << 228 JustWarning, "Region NOT found in store !"); 292 "GeomMgt1001", JustWarning, me << 293 } 229 } 294 return nullptr; << 230 return 0; 295 } 231 } 296 232 297 // ******************************************* 233 // *************************************************************************** 298 // Returns a region through its name specifica 234 // Returns a region through its name specification, if it exists. 299 // If it does not exist it will allocate a new 235 // If it does not exist it will allocate a new region with the given 300 // name, delegating the ownership to the calle 236 // name, delegating the ownership to the caller client. 301 // ******************************************* 237 // *************************************************************************** 302 // 238 // 303 G4Region* G4RegionStore::FindOrCreateRegion(co 239 G4Region* G4RegionStore::FindOrCreateRegion(const G4String& name) 304 { 240 { 305 G4Region* target = GetRegion(name,false); 241 G4Region* target = GetRegion(name,false); 306 if (target == nullptr) << 242 if (!target) 307 { 243 { 308 target = new G4Region(name); 244 target = new G4Region(name); 309 } 245 } 310 return target; 246 return target; 311 } 247 } 312 248 313 // ******************************************* 249 // ************************************************************************** 314 // Set a world physical volume pointer to a re 250 // Set a world physical volume pointer to a region that belongs to it. 315 // Scan over all world volumes. 251 // Scan over all world volumes. 316 // ******************************************* 252 // ************************************************************************** 317 // 253 // 318 void G4RegionStore::SetWorldVolume() 254 void G4RegionStore::SetWorldVolume() 319 { 255 { 320 // Reset all pointers first 256 // Reset all pointers first 321 // 257 // 322 for (auto i=GetInstance()->cbegin(); i!=GetI << 258 for (iterator i=GetInstance()->begin(); i!=GetInstance()->end(); i++) 323 { (*i)->SetWorld(nullptr); } << 259 { (*i)->SetWorld(0); } 324 260 325 // Find world volumes 261 // Find world volumes 326 // 262 // 327 G4PhysicalVolumeStore* fPhysicalVolumeStore 263 G4PhysicalVolumeStore* fPhysicalVolumeStore 328 = G4PhysicalVolumeStore::GetInstance(); 264 = G4PhysicalVolumeStore::GetInstance(); 329 std::size_t nPhys = fPhysicalVolumeStore->si << 265 size_t nPhys = fPhysicalVolumeStore->size(); 330 for(std::size_t iPhys=0; iPhys<nPhys; ++iPhy << 266 for(size_t iPhys=0; iPhys<nPhys; iPhys++) 331 { 267 { 332 G4VPhysicalVolume* fPhys = (*fPhysicalVolu 268 G4VPhysicalVolume* fPhys = (*fPhysicalVolumeStore)[iPhys]; 333 if(fPhys->GetMotherLogical() != nullptr) { << 269 if(fPhys->GetMotherLogical()) { continue; } // not a world volume 334 270 335 // Now 'fPhys' is a world volume, set it t 271 // Now 'fPhys' is a world volume, set it to regions that belong to it. 336 // 272 // 337 for (auto i=GetInstance()->cbegin(); i!=Ge << 273 for (iterator i=GetInstance()->begin(); i!=GetInstance()->end(); i++) 338 { (*i)->SetWorld(fPhys); } 274 { (*i)->SetWorld(fPhys); } 339 } 275 } 340 } 276 } 341 277 342 278