| 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 // G4Region << 26 //
>> 27 //
>> 28 // class G4Region
27 // 29 //
28 // Class description: 30 // Class description:
29 // 31 //
30 // Defines a region or a group of regions in t 32 // Defines a region or a group of regions in the detector geometry
31 // setup, sharing properties associated to mat 33 // setup, sharing properties associated to materials or production
32 // cuts which may affect or bias specific phys 34 // cuts which may affect or bias specific physics processes.
33 35
34 // 18.09.02, G.Cosmo - Initial version << 36 // History:
>> 37 // 18.09.02 G.Cosmo Initial version
35 // ------------------------------------------- 38 // --------------------------------------------------------------------
36 #ifndef G4REGION_HH 39 #ifndef G4REGION_HH
37 #define G4REGION_HH 1 << 40 #define G4REGION_HH
38 41
39 #include <vector> 42 #include <vector>
40 #include <map> 43 #include <map>
41 #include <algorithm> 44 #include <algorithm>
42 45
43 #include "G4Types.hh" 46 #include "G4Types.hh"
44 #include "G4String.hh" 47 #include "G4String.hh"
45 #include "G4GeomSplitter.hh" 48 #include "G4GeomSplitter.hh"
46 49
47 class G4ProductionCuts; 50 class G4ProductionCuts;
48 class G4LogicalVolume; 51 class G4LogicalVolume;
49 class G4Material; 52 class G4Material;
50 class G4VUserRegionInformation; 53 class G4VUserRegionInformation;
51 class G4MaterialCutsCouple; 54 class G4MaterialCutsCouple;
52 class G4UserLimits; 55 class G4UserLimits;
53 class G4FieldManager; 56 class G4FieldManager;
54 class G4FastSimulationManager; 57 class G4FastSimulationManager;
55 class G4VPhysicalVolume; 58 class G4VPhysicalVolume;
56 class G4UserSteppingAction; 59 class G4UserSteppingAction;
57 60
58 class G4RegionData 61 class G4RegionData
59 { 62 {
60 // Encapsulates the fields associated to the 63 // Encapsulates the fields associated to the class
61 // G4Region that may not be read-only. 64 // G4Region that may not be read-only.
62 65
63 public: 66 public:
64 67
65 void initialize() 68 void initialize()
66 { 69 {
67 fFastSimulationManager = nullptr; << 70 fFastSimulationManager = 0;
68 fRegionalSteppingAction = nullptr; << 71 fRegionalSteppingAction = 0;
69 } 72 }
70 73
71 G4FastSimulationManager* fFastSimulationMa 74 G4FastSimulationManager* fFastSimulationManager;
72 G4UserSteppingAction* fRegionalSteppingAct 75 G4UserSteppingAction* fRegionalSteppingAction;
73 }; 76 };
74 77
75 // The type G4RegionManager is introduced to e 78 // The type G4RegionManager is introduced to encapsulate the methods used by
76 // both the master thread and worker threads t 79 // both the master thread and worker threads to allocate memory space for
77 // the fields encapsulated by the class G4Regi 80 // the fields encapsulated by the class G4RegionData. When each thread
78 // initializes the value for these fields, it 81 // initializes the value for these fields, it refers to them using a macro
79 // definition defined below. For every G4Regio 82 // definition defined below. For every G4Region instance, there is a
80 // corresponding G4RegionData instance. All G4 83 // corresponding G4RegionData instance. All G4RegionData instances are
81 // organized by the class G4RegionManager as a 84 // organized by the class G4RegionManager as an array.
82 // The field "int instanceID" is added to the 85 // The field "int instanceID" is added to the class G4Region.
83 // The value of this field in each G4Region in 86 // The value of this field in each G4Region instance is the subscript
84 // of the corresponding G4RegionData instance. 87 // of the corresponding G4RegionData instance.
85 // In order to use the class G4RegionManager, 88 // In order to use the class G4RegionManager, we add a static member in
86 // the class G4Region as follows: "static G4Re 89 // the class G4Region as follows: "static G4RegionManager subInstanceManager".
87 // For the master thread, the array for G4Regi 90 // For the master thread, the array for G4RegionData instances grows
88 // dynamically along with G4Region instances a 91 // dynamically along with G4Region instances are created. For each worker
89 // thread, it copies the array of G4RegionData 92 // thread, it copies the array of G4RegionData instances from the master thread.
90 // In addition, it invokes a method similiar t 93 // In addition, it invokes a method similiar to the constructor explicitly
91 // to achieve the partial effect for each inst 94 // to achieve the partial effect for each instance in the array.
92 // 95 //
93 using G4RegionManager = G4GeomSplitter<G4Regio << 96 typedef G4GeomSplitter<G4RegionData> G4RegionManager;
94 97
95 class G4Region 98 class G4Region
96 { 99 {
97 public: << 100 typedef std::vector<G4LogicalVolume*> G4RootLVList;
>> 101 typedef std::vector<G4Material*> G4MaterialList;
>> 102 typedef std::pair<G4Material*,G4MaterialCutsCouple*> G4MaterialCouplePair;
>> 103 typedef std::map<G4Material*,G4MaterialCutsCouple*> G4MaterialCoupleMap;
>> 104
>> 105 public: // with description
98 106
99 G4Region(const G4String& name); 107 G4Region(const G4String& name);
100 virtual ~G4Region(); 108 virtual ~G4Region();
101 109
102 G4Region(const G4Region&) = delete; <<
103 G4Region& operator=(const G4Region&) = del <<
104 // Copy constructor and assignment opera <<
105 <<
106 inline G4bool operator==(const G4Region& r 110 inline G4bool operator==(const G4Region& rg) const;
107 // Equality defined by address only. 111 // Equality defined by address only.
108 112
109 void AddRootLogicalVolume(G4LogicalVolume* 113 void AddRootLogicalVolume(G4LogicalVolume* lv, G4bool search=true);
110 void RemoveRootLogicalVolume(G4LogicalVolu 114 void RemoveRootLogicalVolume(G4LogicalVolume* lv, G4bool scan=true);
111 // Add/remove root logical volumes and s 115 // Add/remove root logical volumes and set/reset their
112 // daughters flags as regions. They also 116 // daughters flags as regions. They also recompute the
113 // materials list for the region. Flag f 117 // materials list for the region. Flag for scanning the subtree
114 // always enabled by default. Search in 118 // always enabled by default. Search in the tree can be turned off
115 // when adding, assuming the user guaran 119 // when adding, assuming the user guarantees the logical volume is
116 // NOT already inserted, in which case s 120 // NOT already inserted, in which case significant speedup can be
117 // achieved in very complex flat geometr 121 // achieved in very complex flat geometry setups.
118 122
119 void SetName(const G4String& name); << 123 inline void SetName(const G4String& name);
120 inline const G4String& GetName() const; 124 inline const G4String& GetName() const;
121 // Set/get region's name. 125 // Set/get region's name.
122 126
123 inline void RegionModified(G4bool flag); 127 inline void RegionModified(G4bool flag);
124 inline G4bool IsModified() const; 128 inline G4bool IsModified() const;
125 // Accessors to flag identifying if a re 129 // Accessors to flag identifying if a region has been modified
126 // (and still cuts needs to be computed) 130 // (and still cuts needs to be computed) or not.
127 131
128 inline void SetProductionCuts(G4Production 132 inline void SetProductionCuts(G4ProductionCuts* cut);
129 inline G4ProductionCuts* GetProductionCuts 133 inline G4ProductionCuts* GetProductionCuts() const;
130 134
131 inline std::vector<G4LogicalVolume*>::iter 135 inline std::vector<G4LogicalVolume*>::iterator
132 GetRootLogicalVolumeIterator(); 136 GetRootLogicalVolumeIterator();
133 inline std::vector<G4Material*>::const_ite 137 inline std::vector<G4Material*>::const_iterator
134 GetMaterialIterator() const; 138 GetMaterialIterator() const;
135 // Return iterators to lists of root log 139 // Return iterators to lists of root logical volumes and materials.
136 140
137 inline std::size_t GetNumberOfMaterials() << 141 inline size_t GetNumberOfMaterials() const;
138 inline std::size_t GetNumberOfRootVolumes( << 142 inline size_t GetNumberOfRootVolumes() const;
139 // Return the number of elements in the 143 // Return the number of elements in the lists of materials and
140 // root logical volumes. 144 // root logical volumes.
141 145
142 void UpdateMaterialList(); 146 void UpdateMaterialList();
143 // Clears material list and recomputes i 147 // Clears material list and recomputes it looping through
144 // each root logical volume in the regio 148 // each root logical volume in the region.
145 149
146 void ClearMaterialList(); 150 void ClearMaterialList();
147 // Clears the material list. 151 // Clears the material list.
148 152
149 void ScanVolumeTree(G4LogicalVolume* lv, G 153 void ScanVolumeTree(G4LogicalVolume* lv, G4bool region);
150 // Scans recursively the 'lv' logical vo 154 // Scans recursively the 'lv' logical volume tree, retrieves
151 // and places all materials in the list 155 // and places all materials in the list if becoming a region.
152 156
153 inline void SetUserInformation(G4VUserRegi 157 inline void SetUserInformation(G4VUserRegionInformation* ui);
154 inline G4VUserRegionInformation* GetUserIn 158 inline G4VUserRegionInformation* GetUserInformation() const;
155 // Set and Get methods for user informat 159 // Set and Get methods for user information.
156 160
157 inline void SetUserLimits(G4UserLimits* ul 161 inline void SetUserLimits(G4UserLimits* ul);
158 inline G4UserLimits* GetUserLimits() const 162 inline G4UserLimits* GetUserLimits() const;
159 // Set and Get methods for userL-limits 163 // Set and Get methods for userL-limits associated to a region.
160 // Once user-limits are set, it will pro 164 // Once user-limits are set, it will propagate to daughter volumes.
161 165
162 inline void ClearMap(); 166 inline void ClearMap();
163 // Reset G4MaterialCoupleMap 167 // Reset G4MaterialCoupleMap
164 168
165 inline void RegisterMaterialCouplePair(G4M 169 inline void RegisterMaterialCouplePair(G4Material* mat,
166 G4M 170 G4MaterialCutsCouple* couple);
167 // Method invoked by G4ProductionCutsTab 171 // Method invoked by G4ProductionCutsTable to register the pair.
168 172
169 inline G4MaterialCutsCouple* FindCouple(G4 173 inline G4MaterialCutsCouple* FindCouple(G4Material* mat);
170 // Find a G4MaterialCutsCouple which cor 174 // Find a G4MaterialCutsCouple which corresponds to the material
171 // in this region. 175 // in this region.
172 176
173 void SetFastSimulationManager(G4FastSimula 177 void SetFastSimulationManager(G4FastSimulationManager* fsm);
174 G4FastSimulationManager* GetFastSimulation 178 G4FastSimulationManager* GetFastSimulationManager() const;
175 // Set and Get methods for G4FastSimulat 179 // Set and Get methods for G4FastSimulationManager.
176 // The root logical volume that has the 180 // The root logical volume that has the region with G4FastSimulationManager
177 // becomes an envelope of fast simulatio 181 // becomes an envelope of fast simulation.
178 182
179 void ClearFastSimulationManager(); 183 void ClearFastSimulationManager();
180 // Set G4FastSimulationManager pointer t 184 // Set G4FastSimulationManager pointer to the one for the parent region
181 // if it exists. Otherwise set to null. 185 // if it exists. Otherwise set to null.
182 186
183 inline void SetFieldManager(G4FieldManager 187 inline void SetFieldManager(G4FieldManager* fm);
184 inline G4FieldManager* GetFieldManager() c 188 inline G4FieldManager* GetFieldManager() const;
185 // Set and Get methods for G4FieldManage 189 // Set and Get methods for G4FieldManager.
186 // The region with assigned field-manage 190 // The region with assigned field-manager sets the field to the
187 // geometrical area associated with it; 191 // geometrical area associated with it; priority is anyhow given
188 // to local fields eventually set to log 192 // to local fields eventually set to logical volumes.
189 193
190 inline G4VPhysicalVolume* GetWorldPhysical 194 inline G4VPhysicalVolume* GetWorldPhysical() const;
191 // Get method for the world physical vol 195 // Get method for the world physical volume which this region
192 // belongs to. A valid pointer will be a 196 // belongs to. A valid pointer will be assigned by G4RunManagerKernel
193 // through G4RegionStore when the geomet 197 // through G4RegionStore when the geometry is to be closed. Thus, this
194 // pointer may be incorrect at PreInit a 198 // pointer may be incorrect at PreInit and Idle state. If the pointer
195 // is null at the proper state, this par 199 // is null at the proper state, this particular region does not belong
196 // to any world (maybe not assigned to a 200 // to any world (maybe not assigned to any volume, etc.).
197 201
198 void SetWorld(G4VPhysicalVolume* wp); 202 void SetWorld(G4VPhysicalVolume* wp);
199 // Set the world physical volume if this 203 // Set the world physical volume if this region belongs to this world.
200 // If wp is null, reset the pointer. 204 // If wp is null, reset the pointer.
201 205
202 G4bool BelongsTo(G4VPhysicalVolume* thePhy 206 G4bool BelongsTo(G4VPhysicalVolume* thePhys) const;
203 // Returns whether this region belongs t 207 // Returns whether this region belongs to the given physical volume
204 // (recursively scanned to the bottom of 208 // (recursively scanned to the bottom of the hierarchy).
205 209
206 G4Region* GetParentRegion(G4bool& unique) 210 G4Region* GetParentRegion(G4bool& unique) const;
207 // Returns a region that contains this r 211 // Returns a region that contains this region. Otherwise null returned.
208 // Flag 'unique' is true if there is onl 212 // Flag 'unique' is true if there is only one parent region containing
209 // the current region. 213 // the current region.
210 214
211 void SetRegionalSteppingAction(G4UserStepp 215 void SetRegionalSteppingAction(G4UserSteppingAction* rusa);
212 G4UserSteppingAction* GetRegionalSteppingA 216 G4UserSteppingAction* GetRegionalSteppingAction() const;
213 // Set/Get method of the regional user s 217 // Set/Get method of the regional user stepping action
214 218
215 public: << 219 public: // without description
216 220
217 G4Region(__void__&); 221 G4Region(__void__&);
218 // Fake default constructor for usage re 222 // Fake default constructor for usage restricted to direct object
219 // persistency for clients requiring pre 223 // persistency for clients requiring preallocation of memory for
220 // persistifiable objects. 224 // persistifiable objects.
221 225
222 inline G4int GetInstanceID() const; 226 inline G4int GetInstanceID() const;
223 // Returns the instance ID. 227 // Returns the instance ID.
224 228
225 static const G4RegionManager& GetSubInstan 229 static const G4RegionManager& GetSubInstanceManager();
226 // Returns the private data instance man 230 // Returns the private data instance manager.
227 231
228 static void Clean(); 232 static void Clean();
229 // Clear memory allocated by sub-instanc 233 // Clear memory allocated by sub-instance manager.
230 234
231 inline void UsedInMassGeometry(G4bool val << 235 inline void UsedInMassGeometry(G4bool val=true);
232 inline void UsedInParallelGeometry(G4bool << 236 inline void UsedInParallelGeometry(G4bool val=true);
233 inline G4bool IsInMassGeometry() const; 237 inline G4bool IsInMassGeometry() const;
234 inline G4bool IsInParallelGeometry() const 238 inline G4bool IsInParallelGeometry() const;
235 // Utility methods to identify if region 239 // Utility methods to identify if region is part of the main mass
236 // geometry for tracking or a parallel g 240 // geometry for tracking or a parallel geometry.
237 241
238 private: 242 private:
239 243
>> 244 G4Region(const G4Region&);
>> 245 G4Region& operator=(const G4Region&);
>> 246 // Private copy constructor and assignment operator.
>> 247
240 inline void AddMaterial (G4Material* aMate 248 inline void AddMaterial (G4Material* aMaterial);
241 // Searchs the specified material in the 249 // Searchs the specified material in the material table and
242 // if not present adds it. 250 // if not present adds it.
243 251
244 private: 252 private:
245 253
246 using G4RootLVList = std::vector<G4Logical <<
247 using G4MaterialList = std::vector<G4Mater <<
248 using G4MaterialCouplePair = std::pair<G4M <<
249 using G4MaterialCoupleMap = std::map<G4Mat <<
250 <<
251 G4String fName; 254 G4String fName;
252 255
253 G4RootLVList fRootVolumes; 256 G4RootLVList fRootVolumes;
254 G4MaterialList fMaterials; 257 G4MaterialList fMaterials;
255 G4MaterialCoupleMap fMaterialCoupleMap; 258 G4MaterialCoupleMap fMaterialCoupleMap;
256 259
257 G4bool fRegionMod = true; << 260 G4bool fRegionMod;
258 G4ProductionCuts* fCut = nullptr; << 261 G4ProductionCuts* fCut;
259 262
260 G4VUserRegionInformation* fUserInfo = null << 263 G4VUserRegionInformation* fUserInfo;
261 G4UserLimits* fUserLimits = nullptr; << 264 G4UserLimits* fUserLimits;
262 G4FieldManager* fFieldManager = nullptr; << 265 G4FieldManager* fFieldManager;
263 266
264 G4VPhysicalVolume* fWorldPhys = nullptr; << 267 G4VPhysicalVolume* fWorldPhys;
265 268
266 G4bool fInMassGeometry = false; << 269 G4bool fInMassGeometry;
267 G4bool fInParallelGeometry = false; << 270 G4bool fInParallelGeometry;
268 271
269 G4int instanceID; 272 G4int instanceID;
270 // This field is used as instance ID. 273 // This field is used as instance ID.
271 G4GEOM_DLL static G4RegionManager subInsta 274 G4GEOM_DLL static G4RegionManager subInstanceManager;
272 // This field helps to use the class G4R 275 // This field helps to use the class G4RegionManager introduced above.
273 }; 276 };
274 277
275 #include "G4Region.icc" 278 #include "G4Region.icc"
276 279
277 #endif 280 #endif
278 281