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