Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/geometry/management/src/G4SolidStore.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/G4SolidStore.cc (Version 11.3.0) and /geometry/management/src/G4SolidStore.cc (Version 7.0.p1)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                    <<   3 // * DISCLAIMER                                                       *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th <<   5 // * The following disclaimer summarizes all the specific disclaimers *
  6 // * the Geant4 Collaboration.  It is provided <<   6 // * of contributors to this software. The specific disclaimers,which *
  7 // * conditions of the Geant4 Software License <<   7 // * govern, are listed with their locations in:                      *
  8 // * LICENSE and available at  http://cern.ch/ <<   8 // *   http://cern.ch/geant4/license                                  *
  9 // * include a list of copyright holders.      << 
 10 // *                                                9 // *                                                                  *
 11 // * Neither the authors of this software syst     10 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     11 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     12 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     13 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  <<  14 // * use.                                                             *
 16 // * for the full disclaimer and the limitatio << 
 17 // *                                               15 // *                                                                  *
 18 // * This  code  implementation is the result  <<  16 // * This  code  implementation is the  intellectual property  of the *
 19 // * technical work of the GEANT4 collaboratio <<  17 // * GEANT4 collaboration.                                            *
 20 // * By using,  copying,  modifying or  distri <<  18 // * By copying,  distributing  or modifying the Program (or any work *
 21 // * any work based  on the software)  you  ag <<  19 // * based  on  the Program)  you indicate  your  acceptance of  this *
 22 // * use  in  resulting  scientific  publicati <<  20 // * statement, and all its terms.                                    *
 23 // * acceptance of all terms of the Geant4 Sof << 
 24 // *******************************************     21 // ********************************************************************
 25 //                                                 22 //
 26 // G4SolidStore implementation                 << 
 27 //                                                 23 //
 28 // 10.07.95, P.Kent, G.Cosmo                   <<  24 // $Id: G4SolidStore.cc,v 1.13 2004/09/03 08:17:58 gcosmo Exp $
                                                   >>  25 // GEANT4 tag $Name: geant4-07-00-patch-01 $
                                                   >>  26 //
                                                   >>  27 // G4SolidStore
                                                   >>  28 //
                                                   >>  29 // Implementation for singleton container
                                                   >>  30 //
                                                   >>  31 // History:
                                                   >>  32 // 10.07.95 P.Kent Initial version
 29 // -------------------------------------------     33 // --------------------------------------------------------------------
 30                                                    34 
 31 #include "globals.hh"                              35 #include "globals.hh"
 32 #include "G4SolidStore.hh"                         36 #include "G4SolidStore.hh"
 33 #include "G4GeometryManager.hh"                    37 #include "G4GeometryManager.hh"
 34                                                <<  38 #include "G4VStoreNotifier.hh"
 35 #include "G4AutoLock.hh"                       << 
 36                                                << 
 37 namespace                                      << 
 38 {                                              << 
 39   G4Mutex mapMutex = G4MUTEX_INITIALIZER;      << 
 40 }                                              << 
 41                                                    39 
 42 // *******************************************     40 // ***************************************************************************
 43 // Static class variables                          41 // Static class variables
 44 // *******************************************     42 // ***************************************************************************
 45 //                                                 43 //
 46 G4SolidStore* G4SolidStore::fgInstance = nullp <<  44 G4SolidStore* G4SolidStore::fgInstance = 0;
 47 G4ThreadLocal G4VStoreNotifier* G4SolidStore:: <<  45 G4VStoreNotifier* G4SolidStore::fgNotifier = 0;
 48 G4ThreadLocal G4bool G4SolidStore::locked = fa <<  46 G4bool G4SolidStore::locked = false;
 49                                                    47 
 50 // *******************************************     48 // ***************************************************************************
 51 // Protected constructor: Construct underlying     49 // Protected constructor: Construct underlying container with
 52 // initial size of 100 entries                     50 // initial size of 100 entries
 53 // *******************************************     51 // ***************************************************************************
 54 //                                                 52 //
 55 G4SolidStore::G4SolidStore()                       53 G4SolidStore::G4SolidStore()
 56                                                <<  54   : std::vector<G4VSolid*>()
 57 {                                                  55 {
 58   reserve(100);                                    56   reserve(100);
 59 }                                                  57 }
 60                                                    58 
 61 // *******************************************     59 // ***************************************************************************
 62 // Destructor                                      60 // Destructor
 63 // *******************************************     61 // ***************************************************************************
 64 //                                                 62 //
 65 G4SolidStore::~G4SolidStore()                      63 G4SolidStore::~G4SolidStore() 
 66 {                                                  64 {
                                                   >>  65   // NOTE: destruction of solids is client responsibility !
                                                   >>  66   // clear();
 67   Clean();                                         67   Clean();
 68 }                                                  68 }
 69                                                    69 
 70 // *******************************************     70 // ***************************************************************************
 71 // Delete all elements from the store              71 // Delete all elements from the store
 72 // *******************************************     72 // ***************************************************************************
 73 //                                                 73 //
 74 void G4SolidStore::Clean()                         74 void G4SolidStore::Clean()
 75 {                                                  75 {
 76   // Do nothing if geometry is closed              76   // Do nothing if geometry is closed
 77   //                                               77   //
 78   if (G4GeometryManager::GetInstance()->IsGeom     78   if (G4GeometryManager::GetInstance()->IsGeometryClosed())
 79   {                                                79   {
 80     G4cout << "WARNING - Attempt to delete the     80     G4cout << "WARNING - Attempt to delete the solid store"
 81            << " while geometry closed !" << G4     81            << " while geometry closed !" << G4endl;
 82     return;                                        82     return;
 83   }                                                83   }
 84                                                    84 
 85   // Locks store for deletion of solids. De-re     85   // Locks store for deletion of solids. De-registration will be
 86   // performed at this stage. G4VSolids will n     86   // performed at this stage. G4VSolids will not de-register themselves.
 87   //                                               87   //
 88   locked = true;                                   88   locked = true;  
 89                                                    89 
                                                   >>  90   size_t i=0;
 90   G4SolidStore* store = GetInstance();             91   G4SolidStore* store = GetInstance();
                                                   >>  92   std::vector<G4VSolid*>::iterator pos;
 91                                                    93 
 92   for(auto pos=store->cbegin(); pos!=store->ce <<  94 #ifdef G4GEOMETRY_VOXELDEBUG
                                                   >>  95   G4cout << "Deleting Solids ... ";
                                                   >>  96 #endif
                                                   >>  97 
                                                   >>  98   for(pos=store->begin(); pos!=store->end(); pos++)
 93   {                                                99   {
 94     if (fgNotifier != nullptr) { fgNotifier->N << 100     if (fgNotifier) fgNotifier->NotifyDeRegistration();
 95     delete *pos;                               << 101     if (*pos) delete *pos; i++;
 96   }                                               102   }
 97                                                   103 
 98   store->bmap.clear(); store->mvalid = false;  << 104 #ifdef G4GEOMETRY_VOXELDEBUG
                                                   >> 105   if (store->size() < i-1)
                                                   >> 106     { G4cout << "No solids deleted. Already deleted by user ?" << G4endl; }
                                                   >> 107   else
                                                   >> 108     { G4cout << i-1 << " solids deleted !" << G4endl; }
                                                   >> 109 #endif
                                                   >> 110 
 99   locked = false;                                 111   locked = false;
100   store->clear();                                 112   store->clear();
101 }                                                 113 }
102                                                   114 
103 // *******************************************    115 // ***************************************************************************
104 // Associate user notifier to the store           116 // Associate user notifier to the store
105 // *******************************************    117 // ***************************************************************************
106 //                                                118 //
107 void G4SolidStore::SetNotifier(G4VStoreNotifie    119 void G4SolidStore::SetNotifier(G4VStoreNotifier* pNotifier)
108 {                                                 120 {
109   GetInstance();                                  121   GetInstance();
110   fgNotifier = pNotifier;                         122   fgNotifier = pNotifier;
111 }                                                 123 }
112                                                   124 
113 // *******************************************    125 // ***************************************************************************
114 // Bring contents of internal map up to date a << 
115 // ******************************************* << 
116 //                                             << 
117 void G4SolidStore::UpdateMap()                 << 
118 {                                              << 
119   G4AutoLock l(&mapMutex);  // to avoid thread << 
120   if (mvalid) return;                          << 
121   bmap.clear();                                << 
122   for(auto pos=GetInstance()->cbegin(); pos!=G << 
123   {                                            << 
124     const G4String& sol_name = (*pos)->GetName << 
125     auto it = bmap.find(sol_name);             << 
126     if (it != bmap.cend())                     << 
127     {                                          << 
128       it->second.push_back(*pos);              << 
129     }                                          << 
130     else                                       << 
131     {                                          << 
132       std::vector<G4VSolid*> sol_vec { *pos }; << 
133       bmap.insert(std::make_pair(sol_name, sol << 
134     }                                          << 
135   }                                            << 
136   mvalid = true;                               << 
137   l.unlock();                                  << 
138 }                                              << 
139                                                << 
140 // ******************************************* << 
141 // Add Solid to container                         126 // Add Solid to container
142 // *******************************************    127 // ***************************************************************************
143 //                                                128 //
144 void G4SolidStore::Register(G4VSolid* pSolid)     129 void G4SolidStore::Register(G4VSolid* pSolid)
145 {                                                 130 {
146   G4SolidStore* store = GetInstance();         << 131   GetInstance()->push_back(pSolid);
147   store->push_back(pSolid);                    << 132   if (fgNotifier) fgNotifier->NotifyRegistration();
148   const G4String& sol_name = pSolid->GetName() << 
149   auto it = store->bmap.find(sol_name);        << 
150   if (it != store->bmap.cend())                << 
151   {                                            << 
152     it->second.push_back(pSolid);              << 
153   }                                            << 
154   else                                         << 
155   {                                            << 
156     std::vector<G4VSolid*> sol_vec { pSolid }; << 
157     store->bmap.insert(std::make_pair(sol_name << 
158   }                                            << 
159   if (fgNotifier != nullptr) { fgNotifier->Not << 
160   store->mvalid = true;                        << 
161 }                                                 133 }
162                                                   134 
163 // *******************************************    135 // ***************************************************************************
164 // Remove Solid from container                    136 // Remove Solid from container
165 // *******************************************    137 // ***************************************************************************
166 //                                                138 //
167 void G4SolidStore::DeRegister(G4VSolid* pSolid    139 void G4SolidStore::DeRegister(G4VSolid* pSolid)
168 {                                                 140 {
169   G4SolidStore* store = GetInstance();         << 
170   if (!locked)    // Do not de-register if loc    141   if (!locked)    // Do not de-register if locked !
171   {                                               142   {
172     if (fgNotifier != nullptr) { fgNotifier->N << 143     if (fgNotifier) fgNotifier->NotifyDeRegistration();
173     for (auto i=store->crbegin(); i!=store->cr << 144     for (iterator i=GetInstance()->begin(); i!=GetInstance()->end(); i++)
174     {                                             145     {
175       if (**i==*pSolid)                           146       if (**i==*pSolid)
176       {                                           147       {
177         store->erase(std::next(i).base());     << 148         GetInstance()->erase(i);
178         store->mvalid = false;                 << 
179         break;                                    149         break;
180       }                                           150       }
181     }                                             151     }
182     const G4String& sol_name = pSolid->GetName << 
183     auto it = store->bmap.find(sol_name);      << 
184     if (it != store->bmap.cend())              << 
185     {                                          << 
186       if (it->second.size() > 1)               << 
187       {                                        << 
188         for (auto i=it->second.cbegin(); i!=it << 
189         {                                      << 
190           if (**i==*pSolid)                    << 
191           {                                    << 
192             it->second.erase(i);               << 
193             break;                             << 
194           }                                    << 
195         }                                      << 
196       }                                        << 
197       else                                     << 
198       {                                        << 
199         store->bmap.erase(it);                 << 
200       }                                        << 
201     }                                          << 
202   }                                            << 
203 }                                              << 
204                                                << 
205 // ******************************************* << 
206 // Retrieve the first or last solid pointer in << 
207 // ******************************************* << 
208 //                                             << 
209 G4VSolid* G4SolidStore::GetSolid(const G4Strin << 
210                                  G4bool revers << 
211 {                                              << 
212   G4SolidStore* store = GetInstance();         << 
213   if (!store->mvalid)  { store->UpdateMap(); } << 
214   auto pos = store->bmap.find(name);           << 
215   if(pos != store->bmap.cend())                << 
216   {                                            << 
217     if ((verbose) && (pos->second.size()>1))   << 
218     {                                          << 
219       std::ostringstream message;              << 
220       message << "There exists more than ONE s << 
221               << name << "!" << G4endl         << 
222               << "Returning the first found."; << 
223       G4Exception("G4SolidStore::GetSolid()",  << 
224                   "GeomMgt1001", JustWarning,  << 
225     }                                          << 
226     if(reverseSearch)                          << 
227     {                                          << 
228       return pos->second[pos->second.size()-1] << 
229     }                                          << 
230     else                                       << 
231     {                                          << 
232       return pos->second[0];                   << 
233     }                                          << 
234   }                                            << 
235   if (verbose)                                 << 
236   {                                            << 
237     std::ostringstream message;                << 
238     message << "Solid " << name << " not found << 
239             << "Returning NULL pointer.";      << 
240     G4Exception("G4SolidStore::GetSolid()",    << 
241                 "GeomMgt1001", JustWarning, me << 
242   }                                               152   }
243   return nullptr;                              << 
244 }                                                 153 }
245                                                   154 
246 // *******************************************    155 // ***************************************************************************
247 // Return ptr to Store, setting if necessary      156 // Return ptr to Store, setting if necessary
248 // *******************************************    157 // ***************************************************************************
249 //                                                158 //
250 G4SolidStore* G4SolidStore::GetInstance()         159 G4SolidStore* G4SolidStore::GetInstance()
251 {                                                 160 {
252   static G4SolidStore worldStore;                 161   static G4SolidStore worldStore;
253   if (fgInstance == nullptr)                   << 162   if (!fgInstance)
254   {                                               163   {
255     fgInstance = &worldStore;                     164     fgInstance = &worldStore;
256   }                                               165   }
257   return fgInstance;                              166   return fgInstance;
258 }                                                 167 }
259                                                   168