| 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 //
>> 26 //
>> 27 // $Id: G4LogicalVolume.icc,v 1.25 2008-07-10 09:40:08 gcosmo Exp $
>> 28 // GEANT4 tag $Name: geant4-09-04-patch-01 $
>> 29 //
25 // 30 //
26 // class G4LogicalVolume inline implementation << 31 // class G4LogicalVolume Inline Implementation file
27 // 32 //
28 // 15.01.13 - G.Cosmo, A.Dotti: Modified for t << 33 // 10.20.97 - P. MoraDeFreitas : Added SetFastSimulation method.
29 // 10.20.97 - P. MoraDeFreitas : Added SetFast <<
30 // 05.11.98 - M. Verderi: Add Get/Set methods 34 // 05.11.98 - M. Verderi: Add Get/Set methods for fBiasWeight
31 // 09.11.98 - J.Apostolakis: Changed Magnetic << 35 // 09.11.98 - J. Apostolakis: Changed MagneticField to FieldManager
32 // 12.02.99 - S.Giani: Added set/get methods f 36 // 12.02.99 - S.Giani: Added set/get methods for voxelization quality
33 // 18.04.01 - G.Cosmo: Migrated to STL vector 37 // 18.04.01 - G.Cosmo: Migrated to STL vector
34 // 17.05.02 - G.Cosmo: Added IsToOptimise() me 38 // 17.05.02 - G.Cosmo: Added IsToOptimise() method
35 // ------------------------------------------- 39 // --------------------------------------------------------------------
36 40
37 // ******************************************* 41 // ********************************************************************
38 // GetName 42 // GetName
39 // ******************************************* 43 // ********************************************************************
40 // 44 //
41 inline 45 inline
42 const G4String& G4LogicalVolume::GetName() con << 46 G4String G4LogicalVolume::GetName() const
43 { 47 {
44 return fName; 48 return fName;
45 } 49 }
46 50
47 // ******************************************* 51 // ********************************************************************
48 // GetInstanceID << 52 // SetName
49 // ******************************************* 53 // ********************************************************************
50 // 54 //
51 inline 55 inline
52 G4int G4LogicalVolume::GetInstanceID() const << 56 void G4LogicalVolume::SetName(const G4String& pName)
53 { 57 {
54 return instanceID; << 58 fName = pName;
55 } 59 }
56 60
57 // ******************************************* 61 // ********************************************************************
58 // GetMasterFieldManager << 62 // GetFieldManager
59 // ******************************************* 63 // ********************************************************************
60 // 64 //
61 inline 65 inline
62 G4FieldManager* G4LogicalVolume::GetMasterFiel << 66 G4FieldManager* G4LogicalVolume::GetFieldManager() const
63 { 67 {
64 return fFieldManager; 68 return fFieldManager;
65 } 69 }
66 70
67 // ******************************************* 71 // ********************************************************************
68 // GetNoDaughters 72 // GetNoDaughters
69 // ******************************************* 73 // ********************************************************************
70 // 74 //
71 inline 75 inline
72 std::size_t G4LogicalVolume::GetNoDaughters() << 76 G4int G4LogicalVolume::GetNoDaughters() const
73 { 77 {
74 return fDaughters.size(); 78 return fDaughters.size();
75 } 79 }
76 80
77 // ******************************************* 81 // ********************************************************************
78 // GetDaughter 82 // GetDaughter
79 // ******************************************* 83 // ********************************************************************
80 // 84 //
81 inline 85 inline
82 G4VPhysicalVolume* G4LogicalVolume::GetDaughte << 86 G4VPhysicalVolume* G4LogicalVolume::GetDaughter(const G4int i) const
83 { 87 {
84 return fDaughters[i]; 88 return fDaughters[i];
85 } 89 }
86 90
87 // ******************************************* 91 // ********************************************************************
88 // GetFastSimulationManager 92 // GetFastSimulationManager
89 // ******************************************* 93 // ********************************************************************
90 // 94 //
91 inline 95 inline
92 G4FastSimulationManager* G4LogicalVolume::GetF 96 G4FastSimulationManager* G4LogicalVolume::GetFastSimulationManager () const
93 { 97 {
94 G4FastSimulationManager* fFSM = nullptr; << 98 G4FastSimulationManager* fFSM = 0;
95 if(fRegion != nullptr) fFSM = fRegion->GetFa << 99 if(fRegion) fFSM = fRegion->GetFastSimulationManager();
96 return fFSM; 100 return fFSM;
97 } 101 }
98 102
99 // ******************************************* 103 // ********************************************************************
>> 104 // AddDaughter
>> 105 // ********************************************************************
>> 106 //
>> 107 inline
>> 108 void G4LogicalVolume::AddDaughter(G4VPhysicalVolume* pNewDaughter)
>> 109 {
>> 110 if( !fDaughters.empty() && fDaughters[0]->IsReplicated() )
>> 111 {
>> 112 G4cerr << "ERROR - Attempt to place a volume in a mother volume" << G4endl
>> 113 << " already containing a replicated volume." << G4endl
>> 114 << " A volume can either contain several placements" << G4endl
>> 115 << " or a unique replica or parameterised volume !" << G4endl
>> 116 << " Mother logical volume: " << GetName() << G4endl
>> 117 << " Placing volume: " << pNewDaughter->GetName() << G4endl;
>> 118 G4Exception("G4LogicalVolume::AddDaughter()", "InvalidSetup", FatalException,
>> 119 "Replica or parameterised volume must be the only daughter !");
>> 120 }
>> 121
>> 122 fMass = 0.;
>> 123 fDaughters.push_back(pNewDaughter);
>> 124
>> 125 G4LogicalVolume* pDaughterLogical = pNewDaughter->GetLogicalVolume();
>> 126
>> 127 // Propagate the Field Manager, if the daughter has no field Manager.
>> 128 //
>> 129 G4FieldManager* pDaughterFieldManager = pDaughterLogical->GetFieldManager();
>> 130
>> 131 if( pDaughterFieldManager == 0 )
>> 132 {
>> 133 pDaughterLogical->SetFieldManager(fFieldManager, true);
>> 134 }
>> 135 if (fRegion)
>> 136 {
>> 137 PropagateRegion();
>> 138 fRegion->RegionModified(true);
>> 139 }
>> 140 }
>> 141
>> 142 // ********************************************************************
100 // IsDaughter 143 // IsDaughter
101 // ******************************************* 144 // ********************************************************************
102 // 145 //
103 inline 146 inline
104 G4bool G4LogicalVolume::IsDaughter(const G4VPh 147 G4bool G4LogicalVolume::IsDaughter(const G4VPhysicalVolume* p) const
105 { 148 {
106 for (const auto & daughter : fDaughters) << 149 G4PhysicalVolumeList::const_iterator i;
>> 150 for ( i=fDaughters.begin(); i!=fDaughters.end(); ++i )
107 { 151 {
108 if (*daughter==*p) return true; << 152 if (**i==*p) return true;
109 } 153 }
110 return false; 154 return false;
111 } 155 }
112 156
113 // ******************************************* 157 // ********************************************************************
114 // CharacteriseDaughters << 158 // RemoveDaughter
115 // ******************************************* 159 // ********************************************************************
116 // 160 //
117 inline 161 inline
118 EVolume G4LogicalVolume::CharacteriseDaughters << 162 void G4LogicalVolume::RemoveDaughter(const G4VPhysicalVolume* p)
119 { 163 {
120 return fDaughtersVolumeType; << 164 G4PhysicalVolumeList::iterator i;
>> 165 for ( i=fDaughters.begin(); i!=fDaughters.end(); ++i )
>> 166 {
>> 167 if (**i==*p)
>> 168 {
>> 169 fDaughters.erase(i);
>> 170 break;
>> 171 }
>> 172 }
>> 173 if (fRegion)
>> 174 {
>> 175 fRegion->RegionModified(true);
>> 176 }
>> 177 fMass = 0.;
121 } 178 }
122 179
123 // ******************************************* 180 // ********************************************************************
124 // DeduceDaughtersType << 181 // ClearDaughters
125 // ******************************************* 182 // ********************************************************************
126 // 183 //
127 inline 184 inline
128 EVolume G4LogicalVolume::DeduceDaughtersType() << 185 void G4LogicalVolume::ClearDaughters()
129 { 186 {
130 EVolume type= kNormal; << 187 G4PhysicalVolumeList::iterator i;
131 G4VPhysicalVolume* pVol; << 188 for ( i=fDaughters.begin(); i!=fDaughters.end(); ++i )
132 << 189 {
133 if ( GetNoDaughters() >= 1 ) << 190 fDaughters.erase(i);
>> 191 }
>> 192 if (fRegion)
134 { 193 {
135 pVol = GetDaughter(0); << 194 fRegion->RegionModified(true);
136 type = pVol->VolumeType(); <<
137 } 195 }
138 return type; << 196 fMass = 0.;
139 } 197 }
140 198
141 // ******************************************* 199 // ********************************************************************
142 // GetMasterSolid << 200 // GetSolid
143 // ******************************************* 201 // ********************************************************************
144 // 202 //
145 inline 203 inline
146 G4VSolid* G4LogicalVolume::GetMasterSolid() co << 204 G4VSolid* G4LogicalVolume::GetSolid() const
147 { 205 {
148 return fSolid; 206 return fSolid;
149 } 207 }
150 208
151 // ******************************************* 209 // ********************************************************************
152 // GetMasterSensitiveDetector << 210 // SetSolid
>> 211 // ********************************************************************
>> 212 //
>> 213 inline
>> 214 void G4LogicalVolume::SetSolid(G4VSolid *pSolid)
>> 215 {
>> 216 assert(pSolid != 0);
>> 217 fSolid=pSolid;
>> 218 fMass = 0.;
>> 219 }
>> 220
>> 221 // ********************************************************************
>> 222 // GetMaterial
>> 223 // ********************************************************************
>> 224 //
>> 225 inline
>> 226 G4Material* G4LogicalVolume::GetMaterial() const
>> 227 {
>> 228 return fMaterial;
>> 229 }
>> 230
>> 231 // ********************************************************************
>> 232 // SetMaterial
>> 233 // ********************************************************************
>> 234 //
>> 235 inline
>> 236 void G4LogicalVolume::SetMaterial(G4Material *pMaterial)
>> 237 {
>> 238 fMaterial=pMaterial;
>> 239 fMass = 0.;
>> 240 }
>> 241
>> 242 // ********************************************************************
>> 243 // UpdateMaterial
153 // ******************************************* 244 // ********************************************************************
154 // 245 //
155 inline 246 inline
156 G4VSensitiveDetector* G4LogicalVolume::GetMast << 247 void G4LogicalVolume::UpdateMaterial(G4Material *pMaterial)
>> 248 {
>> 249 fMaterial=pMaterial;
>> 250 if(fRegion) { fCutsCouple = fRegion->FindCouple(pMaterial); }
>> 251 fMass = 0.;
>> 252 }
>> 253
>> 254 // ********************************************************************
>> 255 // GetSensitiveDetector
>> 256 // ********************************************************************
>> 257 //
>> 258 inline
>> 259 G4VSensitiveDetector* G4LogicalVolume::GetSensitiveDetector() const
157 { 260 {
158 return fSensitiveDetector; 261 return fSensitiveDetector;
159 } 262 }
160 263
161 // ******************************************* 264 // ********************************************************************
>> 265 // SetSensitiveDetector
>> 266 // ********************************************************************
>> 267 //
>> 268 inline
>> 269 void G4LogicalVolume::SetSensitiveDetector(G4VSensitiveDetector* pSDetector)
>> 270 {
>> 271 fSensitiveDetector = pSDetector;
>> 272 }
>> 273
>> 274 // ********************************************************************
162 // GetUserLimits 275 // GetUserLimits
163 // ******************************************* 276 // ********************************************************************
164 // 277 //
165 inline 278 inline
166 G4UserLimits* G4LogicalVolume::GetUserLimits() 279 G4UserLimits* G4LogicalVolume::GetUserLimits() const
167 { 280 {
168 if(fUserLimits != nullptr) return fUserLimit << 281 if(fUserLimits) return fUserLimits;
169 if(fRegion != nullptr) return fRegion->GetUs << 282 if(fRegion) return fRegion->GetUserLimits();
170 return nullptr; << 283 return 0;
171 } 284 }
172 285
173 // ******************************************* 286 // ********************************************************************
174 // SetUserLimits 287 // SetUserLimits
175 // ******************************************* 288 // ********************************************************************
176 // 289 //
177 inline 290 inline
178 void G4LogicalVolume::SetUserLimits(G4UserLimi 291 void G4LogicalVolume::SetUserLimits(G4UserLimits* pULimits)
179 { 292 {
180 fUserLimits = pULimits; 293 fUserLimits = pULimits;
181 } 294 }
182 295
183 // ******************************************* 296 // ********************************************************************
184 // GetVoxelHeader 297 // GetVoxelHeader
185 // ******************************************* 298 // ********************************************************************
186 // 299 //
187 inline 300 inline
188 G4SmartVoxelHeader* G4LogicalVolume::GetVoxelH 301 G4SmartVoxelHeader* G4LogicalVolume::GetVoxelHeader() const
189 { 302 {
190 return fVoxel; 303 return fVoxel;
191 } 304 }
192 305
193 // ******************************************* 306 // ********************************************************************
194 // SetVoxelHeader 307 // SetVoxelHeader
195 // ******************************************* 308 // ********************************************************************
196 // 309 //
197 inline 310 inline
198 void G4LogicalVolume::SetVoxelHeader(G4SmartVo 311 void G4LogicalVolume::SetVoxelHeader(G4SmartVoxelHeader* pVoxel)
199 { 312 {
200 fVoxel = pVoxel; 313 fVoxel = pVoxel;
201 } 314 }
202 315
203 // ******************************************* 316 // ********************************************************************
204 // GetSmartless 317 // GetSmartless
205 // ******************************************* 318 // ********************************************************************
206 // 319 //
207 inline 320 inline
208 G4double G4LogicalVolume::GetSmartless() const 321 G4double G4LogicalVolume::GetSmartless() const
209 { 322 {
210 return fSmartless; 323 return fSmartless;
211 } 324 }
212 325
213 // ******************************************* 326 // ********************************************************************
214 // SetSmartless 327 // SetSmartless
215 // ******************************************* 328 // ********************************************************************
216 // 329 //
217 inline 330 inline
218 void G4LogicalVolume::SetSmartless(G4double sm << 331 void G4LogicalVolume::SetSmartless(G4double s)
219 { 332 {
220 fSmartless = smt; << 333 fSmartless = s;
221 } 334 }
222 335
223 // ******************************************* 336 // ********************************************************************
224 // IsToOptimise 337 // IsToOptimise
225 // ******************************************* 338 // ********************************************************************
226 // 339 //
227 inline 340 inline
228 G4bool G4LogicalVolume::IsToOptimise() const 341 G4bool G4LogicalVolume::IsToOptimise() const
229 { 342 {
230 return fOptimise; 343 return fOptimise;
231 } 344 }
232 345
233 // ******************************************* 346 // ********************************************************************
234 // SetOptimisation 347 // SetOptimisation
235 // ******************************************* 348 // ********************************************************************
236 // 349 //
237 inline 350 inline
238 void G4LogicalVolume::SetOptimisation(G4bool o 351 void G4LogicalVolume::SetOptimisation(G4bool optim)
239 { 352 {
240 fOptimise = optim; 353 fOptimise = optim;
241 } 354 }
242 355
243 // ******************************************* 356 // ********************************************************************
244 // IsRootRegion 357 // IsRootRegion
245 // ******************************************* 358 // ********************************************************************
246 // 359 //
247 inline 360 inline
248 G4bool G4LogicalVolume::IsRootRegion() const 361 G4bool G4LogicalVolume::IsRootRegion() const
249 { 362 {
250 return fRootRegion; 363 return fRootRegion;
251 } 364 }
252 365
253 // ******************************************* 366 // ********************************************************************
254 // SetRegionRootFlag 367 // SetRegionRootFlag
255 // ******************************************* 368 // ********************************************************************
256 // 369 //
257 inline 370 inline
258 void G4LogicalVolume::SetRegionRootFlag(G4bool 371 void G4LogicalVolume::SetRegionRootFlag(G4bool rreg)
259 { 372 {
260 fRootRegion = rreg; 373 fRootRegion = rreg;
261 } 374 }
262 375
263 // ******************************************* 376 // ********************************************************************
264 // IsRegion 377 // IsRegion
265 // ******************************************* 378 // ********************************************************************
266 // 379 //
267 inline 380 inline
268 G4bool G4LogicalVolume::IsRegion() const 381 G4bool G4LogicalVolume::IsRegion() const
269 { 382 {
270 G4bool reg = false; 383 G4bool reg = false;
271 if (fRegion != nullptr) reg = true; << 384 if (fRegion) reg = true;
272 return reg; 385 return reg;
273 } 386 }
274 387
275 // ******************************************* 388 // ********************************************************************
276 // SetRegion 389 // SetRegion
277 // ******************************************* 390 // ********************************************************************
278 // 391 //
279 inline 392 inline
280 void G4LogicalVolume::SetRegion(G4Region* reg) 393 void G4LogicalVolume::SetRegion(G4Region* reg)
281 { 394 {
282 fRegion = reg; 395 fRegion = reg;
283 } 396 }
284 397
285 // ******************************************* 398 // ********************************************************************
286 // GetRegion 399 // GetRegion
287 // ******************************************* 400 // ********************************************************************
288 // 401 //
289 inline 402 inline
290 G4Region* G4LogicalVolume::GetRegion() const 403 G4Region* G4LogicalVolume::GetRegion() const
291 { 404 {
292 return fRegion; 405 return fRegion;
293 } 406 }
294 407
295 // ******************************************* 408 // ********************************************************************
296 // PropagateRegion 409 // PropagateRegion
297 // ******************************************* 410 // ********************************************************************
298 // 411 //
299 inline 412 inline
300 void G4LogicalVolume::PropagateRegion() 413 void G4LogicalVolume::PropagateRegion()
301 { 414 {
302 fRegion->ScanVolumeTree(this, true); 415 fRegion->ScanVolumeTree(this, true);
303 } 416 }
304 417
305 // ******************************************* 418 // ********************************************************************
>> 419 // GetMaterialCutsCouple
>> 420 // ********************************************************************
>> 421 //
>> 422 inline
>> 423 const G4MaterialCutsCouple* G4LogicalVolume::GetMaterialCutsCouple() const
>> 424 {
>> 425 return fCutsCouple;
>> 426 }
>> 427
>> 428 // ********************************************************************
>> 429 // SetMaterialCutsCouple
>> 430 // ********************************************************************
>> 431 //
>> 432 inline
>> 433 void G4LogicalVolume::SetMaterialCutsCouple(G4MaterialCutsCouple* cuts)
>> 434 {
>> 435 fCutsCouple = cuts;
>> 436 }
>> 437
>> 438 // ********************************************************************
306 // Lock 439 // Lock
307 // ******************************************* 440 // ********************************************************************
308 // 441 //
309 inline 442 inline
310 void G4LogicalVolume::Lock() 443 void G4LogicalVolume::Lock()
311 { 444 {
312 fLock = true; 445 fLock = true;
313 } 446 }
314 447
315 // ******************************************* 448 // ********************************************************************
316 // Operator == 449 // Operator ==
317 // ******************************************* 450 // ********************************************************************
318 // 451 //
319 inline 452 inline
320 G4bool G4LogicalVolume::operator == ( const G4 453 G4bool G4LogicalVolume::operator == ( const G4LogicalVolume& lv) const
321 { 454 {
322 return this==&lv; << 455 return (this==&lv) ? true : false;
323 } 456 }
324 457
325 // ******************************************* 458 // ********************************************************************
326 // SetBiasWeight << 459 // GetVisAttributes
327 // ******************************************* 460 // ********************************************************************
328 // 461 //
329 inline 462 inline
330 void G4LogicalVolume::SetBiasWeight(G4double w << 463 const G4VisAttributes* G4LogicalVolume::GetVisAttributes () const
331 { 464 {
332 fBiasWeight = weight; << 465 return fVisAttributes;
333 } 466 }
334 467
335 // ******************************************* 468 // ********************************************************************
336 // GetBiasWeight << 469 // SetVisAttributes
337 // ******************************************* 470 // ********************************************************************
338 // 471 //
339 inline 472 inline
340 G4double G4LogicalVolume::GetBiasWeight() cons << 473 void G4LogicalVolume::SetVisAttributes (const G4VisAttributes* pVA)
341 { 474 {
342 return fBiasWeight; << 475 fVisAttributes = pVA;
343 } 476 }
344 477
345 // ******************************************* 478 // ********************************************************************
346 // GetVisAttributes << 479 // SetBiasWeight
347 // ******************************************* 480 // ********************************************************************
348 // 481 //
349 inline 482 inline
350 const G4VisAttributes* G4LogicalVolume::GetVis << 483 void G4LogicalVolume::SetBiasWeight(G4double w)
>> 484 {
>> 485 fBiasWeight = w;
>> 486 }
>> 487
>> 488 // ********************************************************************
>> 489 // GetBiasWeight
>> 490 // ********************************************************************
>> 491 //
>> 492 inline
>> 493 G4double G4LogicalVolume::GetBiasWeight() const
351 { 494 {
352 return fVisAttributes.get(); << 495 return fBiasWeight;
353 } 496 }
>> 497
354 498