Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/geometry/management/src/G4PhysicalVolumeStore.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 /geometry/management/src/G4PhysicalVolumeStore.cc (Version 11.3.0) and /geometry/management/src/G4PhysicalVolumeStore.cc (Version 9.2.p4)


  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 // G4PhysicalVolumeStore implementation        << 
 27 //                                                 26 //
 28 // 25.07.95, P.Kent, G.Cosmo                   <<  27 // $Id: G4PhysicalVolumeStore.cc,v 1.20 2008/07/10 09:41:20 gcosmo Exp $
                                                   >>  28 // GEANT4 tag $Name: geant4-09-02-patch-04 $
                                                   >>  29 //
                                                   >>  30 // G4PhysicalVolumeStore
                                                   >>  31 //
                                                   >>  32 // Implementation for singleton container
                                                   >>  33 //
                                                   >>  34 // History:
                                                   >>  35 // 25.07.95 P.Kent Initial version
 29 // -------------------------------------------     36 // --------------------------------------------------------------------
 30                                                    37 
 31 #include "G4Types.hh"                              38 #include "G4Types.hh"
 32 #include "G4PhysicalVolumeStore.hh"                39 #include "G4PhysicalVolumeStore.hh"
 33 #include "G4GeometryManager.hh"                    40 #include "G4GeometryManager.hh"
 34 #include "G4LogicalVolume.hh"                      41 #include "G4LogicalVolume.hh"
 35                                                    42 
 36 #include "G4AutoLock.hh"                       << 
 37                                                << 
 38 namespace                                      << 
 39 {                                              << 
 40   G4Mutex mapMutex = G4MUTEX_INITIALIZER;      << 
 41 }                                              << 
 42                                                << 
 43 // *******************************************     43 // ***************************************************************************
 44 // Static class variables                          44 // Static class variables
 45 // *******************************************     45 // ***************************************************************************
 46 //                                                 46 //
 47 G4PhysicalVolumeStore* G4PhysicalVolumeStore:: <<  47 G4PhysicalVolumeStore* G4PhysicalVolumeStore::fgInstance = 0;
 48 G4ThreadLocal G4VStoreNotifier* G4PhysicalVolu <<  48 G4VStoreNotifier* G4PhysicalVolumeStore::fgNotifier = 0;
 49 G4ThreadLocal G4bool G4PhysicalVolumeStore::lo <<  49 G4bool G4PhysicalVolumeStore::locked = false;
 50                                                    50 
 51 // *******************************************     51 // ***************************************************************************
 52 // Protected constructor: Construct underlying     52 // Protected constructor: Construct underlying container with
 53 // initial size of 100 entries                     53 // initial size of 100 entries
 54 // *******************************************     54 // ***************************************************************************
 55 //                                                 55 //
 56 G4PhysicalVolumeStore::G4PhysicalVolumeStore()     56 G4PhysicalVolumeStore::G4PhysicalVolumeStore()
 57                                                <<  57   : std::vector<G4VPhysicalVolume*>()
 58 {                                                  58 {
 59   reserve(100);                                    59   reserve(100);
 60 }                                                  60 }
 61                                                    61 
 62 // *******************************************     62 // ***************************************************************************
 63 // Destructor                                      63 // Destructor
 64 // *******************************************     64 // ***************************************************************************
 65 //                                                 65 //
 66 G4PhysicalVolumeStore::~G4PhysicalVolumeStore(     66 G4PhysicalVolumeStore::~G4PhysicalVolumeStore()
 67 {                                                  67 {
 68   Clean();  // Delete all volumes in the store <<  68   Clean();
 69   G4VPhysicalVolume::Clean();  // Delete alloc << 
 70 }                                                  69 }
 71                                                    70 
 72 // *******************************************     71 // ***************************************************************************
 73 // Delete all elements from the store              72 // Delete all elements from the store
 74 // *******************************************     73 // ***************************************************************************
 75 //                                                 74 //
 76 void G4PhysicalVolumeStore::Clean()                75 void G4PhysicalVolumeStore::Clean()
 77 {                                                  76 {
 78   // Do nothing if geometry is closed              77   // Do nothing if geometry is closed
 79   //                                               78   //
 80   if (G4GeometryManager::GetInstance()->IsGeom     79   if (G4GeometryManager::GetInstance()->IsGeometryClosed())
 81   {                                                80   {
 82     G4cout << "WARNING - Attempt to delete the     81     G4cout << "WARNING - Attempt to delete the physical volume store"
 83            << " while geometry closed !" << G4     82            << " while geometry closed !" << G4endl;
 84     return;                                        83     return;
 85   }                                                84   }
 86                                                    85 
 87   // Locks store for deletion of volumes. De-r     86   // Locks store for deletion of volumes. De-registration will be
 88   // performed at this stage. G4VPhysicalVolum     87   // performed at this stage. G4VPhysicalVolumes will not de-register
 89   // themselves.                                   88   // themselves.
 90   //                                               89   //
 91   locked = true;                                   90   locked = true;
 92                                                    91 
                                                   >>  92   size_t i=0;
 93   G4PhysicalVolumeStore* store = GetInstance()     93   G4PhysicalVolumeStore* store = GetInstance();
                                                   >>  94   std::vector<G4VPhysicalVolume*>::iterator pos;
                                                   >>  95 
                                                   >>  96 #ifdef G4GEOMETRY_VOXELDEBUG
                                                   >>  97   G4cout << "Deleting Physical Volumes ... ";
                                                   >>  98 #endif
 94                                                    99 
 95   for(auto pos=store->cbegin(); pos!=store->ce << 100   for(pos=store->begin(); pos!=store->end(); pos++)
 96   {                                               101   {
 97     if (fgNotifier != nullptr) { fgNotifier->N << 102     if (fgNotifier) { fgNotifier->NotifyDeRegistration(); }
 98     delete *pos;                               << 103     if (*pos) { delete *pos; }
                                                   >> 104     i++;
 99   }                                               105   }
100                                                   106 
101   store->bmap.clear(); store->mvalid = false;  << 107 #ifdef G4GEOMETRY_VOXELDEBUG
                                                   >> 108   if (store->size() < i-1)
                                                   >> 109     { G4cout << "No volumes deleted. Already deleted by user ?" << G4endl; }
                                                   >> 110   else
                                                   >> 111     { G4cout << i-1 << " volumes deleted !" << G4endl; }
                                                   >> 112 #endif
                                                   >> 113 
102   locked = false;                                 114   locked = false;
103   store->clear();                                 115   store->clear();
104 }                                                 116 }
105                                                   117 
106 // *******************************************    118 // ***************************************************************************
107 // Associate user notifier to the store           119 // Associate user notifier to the store
108 // *******************************************    120 // ***************************************************************************
109 //                                                121 //
110 void G4PhysicalVolumeStore::SetNotifier(G4VSto    122 void G4PhysicalVolumeStore::SetNotifier(G4VStoreNotifier* pNotifier)
111 {                                                 123 {
112   GetInstance();                                  124   GetInstance();
113   fgNotifier = pNotifier;                         125   fgNotifier = pNotifier;
114 }                                                 126 }
115                                                   127 
116 // *******************************************    128 // ***************************************************************************
117 // Bring contents of internal map up to date a << 
118 // ******************************************* << 
119 //                                             << 
120 void G4PhysicalVolumeStore::UpdateMap()        << 
121 {                                              << 
122   G4AutoLock l(&mapMutex);  // to avoid thread << 
123   if (mvalid) return;                          << 
124   bmap.clear();                                << 
125   for(auto pos=GetInstance()->cbegin(); pos!=G << 
126   {                                            << 
127     const G4String& vol_name = (*pos)->GetName << 
128     auto it = bmap.find(vol_name);             << 
129     if (it != bmap.cend())                     << 
130     {                                          << 
131       it->second.push_back(*pos);              << 
132     }                                          << 
133     else                                       << 
134     {                                          << 
135       std::vector<G4VPhysicalVolume*> vol_vec  << 
136       bmap.insert(std::make_pair(vol_name, vol << 
137     }                                          << 
138   }                                            << 
139   mvalid = true;                               << 
140   l.unlock();                                  << 
141 }                                              << 
142                                                << 
143 // ******************************************* << 
144 // Add Volume to container                        129 // Add Volume to container
145 // *******************************************    130 // ***************************************************************************
146 //                                                131 //
147 void G4PhysicalVolumeStore::Register(G4VPhysic    132 void G4PhysicalVolumeStore::Register(G4VPhysicalVolume* pVolume)
148 {                                                 133 {
149   G4PhysicalVolumeStore* store = GetInstance() << 134   GetInstance()->push_back(pVolume);
150   store->push_back(pVolume);                   << 135   if (fgNotifier) { fgNotifier->NotifyRegistration(); }
151   const G4String& vol_name = pVolume->GetName( << 
152   auto it = store->bmap.find(vol_name);        << 
153   if (it != store->bmap.cend())                << 
154   {                                            << 
155     it->second.push_back(pVolume);             << 
156   }                                            << 
157   else                                         << 
158   {                                            << 
159     std::vector<G4VPhysicalVolume*> vol_vec {  << 
160     store->bmap.insert(std::make_pair(vol_name << 
161   }                                            << 
162   if (fgNotifier != nullptr) { fgNotifier->Not << 
163   store->mvalid = true;                        << 
164 }                                                 136 }
165                                                   137 
166 // *******************************************    138 // ***************************************************************************
167 // Remove Volume from container and update the    139 // Remove Volume from container and update the list of daughters
168 // of the mother's logical volume                 140 // of the mother's logical volume
169 // *******************************************    141 // ***************************************************************************
170 //                                                142 //
171 void G4PhysicalVolumeStore::DeRegister(G4VPhys    143 void G4PhysicalVolumeStore::DeRegister(G4VPhysicalVolume* pVolume)
172 {                                                 144 {
173   G4PhysicalVolumeStore* 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->N << 147     if (fgNotifier) { fgNotifier->NotifyDeRegistration(); }
177     G4LogicalVolume* motherLogical = pVolume->    148     G4LogicalVolume* motherLogical = pVolume->GetMotherLogical();
178     if (motherLogical != nullptr) { motherLogi << 149     if (motherLogical) { motherLogical->RemoveDaughter(pVolume); }
179     for (auto i=store->cbegin(); i!=store->cen << 150     for (iterator i=GetInstance()->begin(); i!=GetInstance()->end(); i++)
180     {                                             151     {
181       if (**i==*pVolume)                          152       if (**i==*pVolume)
182       {                                           153       {
183         store->erase(i);                       << 154         GetInstance()->erase(i);
184         break;                                    155         break;
185       }                                           156       }
186     }                                             157     }
187     const G4String& vol_name = pVolume->GetNam << 
188     auto it = store->bmap.find(vol_name);      << 
189     if (it != store->bmap.cend())              << 
190     {                                          << 
191       if (it->second.size() > 1)               << 
192       {                                        << 
193         for (auto i=it->second.cbegin(); i!=it << 
194         {                                      << 
195           if (**i==*pVolume)                   << 
196           {                                    << 
197             it->second.erase(i);               << 
198             break;                             << 
199           }                                    << 
200         }                                      << 
201       }                                        << 
202       else                                     << 
203       {                                        << 
204         store->bmap.erase(it);                 << 
205       }                                        << 
206     }                                          << 
207   }                                               158   }
208 }                                                 159 }
209                                                   160 
210 // *******************************************    161 // ***************************************************************************
211 // Retrieve the first or last volume pointer i << 162 // Retrieve the first volume pointer in the container having that name
212 // *******************************************    163 // ***************************************************************************
213 //                                                164 //
214 G4VPhysicalVolume*                                165 G4VPhysicalVolume*
215 G4PhysicalVolumeStore::GetVolume(const G4Strin << 166 G4PhysicalVolumeStore::GetVolume(const G4String& name, G4bool verbose) const
216                                  G4bool revers << 
217 {                                                 167 {
218   G4PhysicalVolumeStore* store = GetInstance() << 168   for (iterator i=GetInstance()->begin(); i!=GetInstance()->end(); i++)
219   if (!store->mvalid)  { store->UpdateMap(); } << 
220   auto pos = store->bmap.find(name);           << 
221   if(pos != store->bmap.cend())                << 
222   {                                               169   {
223     if ((verbose) && (pos->second.size()>1))   << 170     if ((*i)->GetName() == name) { return *i; }
224     {                                          << 
225       std::ostringstream message;              << 
226       message << "There exists more than ONE p << 
227               << name << "!" << G4endl         << 
228               << "Returning the first found."; << 
229       G4Exception("G4PhysicalVolumeStore::GetV << 
230                   "GeomMgt1001", JustWarning,  << 
231     }                                          << 
232     if(reverseSearch)                          << 
233     {                                          << 
234       return pos->second[pos->second.size()-1] << 
235     }                                          << 
236     else                                       << 
237     {                                          << 
238       return pos->second[0];                   << 
239     }                                          << 
240   }                                               171   }
241   if (verbose)                                    172   if (verbose)
242   {                                               173   {
243      std::ostringstream message;               << 174      G4cerr << "ERROR - G4PhysicalVolumeStore::GetVolume()" << G4endl
244      message << "Volume NOT found in store !"  << 
245             << "        Volume " << name << "     175             << "        Volume " << name << " NOT found in store !" << G4endl
246             << "        Returning NULL pointer << 176             << "        Returning NULL pointer." << G4endl;
247      G4Exception("G4PhysicalVolumeStore::GetVo << 177      G4Exception("G4PhysicalVolumeStore::GetVolume()", "InvalidQuery",
248                  "GeomMgt1001", JustWarning, m << 178                  JustWarning, "Volume NOT found in store !");
249   }                                               179   }
250   return nullptr;                              << 180   return 0;
251 }                                                 181 }
252                                                   182 
253 // *******************************************    183 // ***************************************************************************
254 // Return ptr to Store, setting if necessary      184 // Return ptr to Store, setting if necessary
255 // *******************************************    185 // ***************************************************************************
256 //                                                186 //
257 G4PhysicalVolumeStore* G4PhysicalVolumeStore::    187 G4PhysicalVolumeStore* G4PhysicalVolumeStore::GetInstance()
258 {                                                 188 {
259   static G4PhysicalVolumeStore worldStore;        189   static G4PhysicalVolumeStore worldStore;
260   if (fgInstance == nullptr)                   << 190   if (!fgInstance)
261   {                                               191   {
262     fgInstance = &worldStore;                     192     fgInstance = &worldStore;
263   }                                               193   }
264   return fgInstance;                              194   return fgInstance;
265 }                                                 195 }
266                                                   196