| 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.21 2005/05/22 16:41:57 allison Exp $
>> 25 // GEANT4 tag $Name: geant4-07-01-patch-01 $
>> 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 G4cerr << "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()", "InvalidSetup", FatalException,
>> 114 "Replica or parameterised volume must be the only daughter !");
>> 115 }
>> 116
>> 117 fMass = 0.;
>> 118 fDaughters.push_back(pNewDaughter);
>> 119
>> 120 // Propagates the mother's FastSimulationManager pointer. If we are in
>> 121 // the World logical volume, propagates only if pointer != 0.
>> 122 // (one could also parameterise all the World volume for a particle type!)
>> 123
>> 124 G4LogicalVolume* myLogical = pNewDaughter->GetMotherLogical();
>> 125 G4LogicalVolume* pDaughterLogical = pNewDaughter->GetLogicalVolume();
>> 126
>> 127 if( myLogical !=0 )
>> 128 {
>> 129 if( pDaughterLogical->GetFastSimulationManager()!=fFastSimulationManager )
>> 130 pDaughterLogical->SetFastSimulationManager(fFastSimulationManager,false);
>> 131 }
>> 132 else
>> 133 {
>> 134 pDaughterLogical->SetFastSimulationManager(0, false);
>> 135 }
>> 136
>> 137 // Propagate the Field Manager, if the daughter has no field Manager.
>> 138 //
>> 139 G4FieldManager* pDaughterFieldManager = pDaughterLogical->GetFieldManager();
>> 140
>> 141 if( pDaughterFieldManager == 0 )
>> 142 {
>> 143 pDaughterLogical->SetFieldManager(fFieldManager, true);
>> 144 }
>> 145 if (fRegion)
>> 146 {
>> 147 PropagateRegion();
>> 148 fRegion->RegionModified(true);
>> 149 }
97 } 150 }
98 151
99 // ******************************************* 152 // ********************************************************************
100 // IsDaughter 153 // IsDaughter
101 // ******************************************* 154 // ********************************************************************
102 // 155 //
103 inline 156 inline
104 G4bool G4LogicalVolume::IsDaughter(const G4VPh 157 G4bool G4LogicalVolume::IsDaughter(const G4VPhysicalVolume* p) const
105 { 158 {
106 for (const auto & daughter : fDaughters) << 159 G4PhysicalVolumeList::const_iterator i;
>> 160 for ( i=fDaughters.begin(); i!=fDaughters.end(); ++i )
107 { 161 {
108 if (*daughter==*p) return true; << 162 if (**i==*p) return true;
109 } 163 }
110 return false; 164 return false;
111 } 165 }
112 166
113 // ******************************************* 167 // ********************************************************************
114 // CharacteriseDaughters << 168 // RemoveDaughter
115 // ******************************************* 169 // ********************************************************************
116 // 170 //
117 inline 171 inline
118 EVolume G4LogicalVolume::CharacteriseDaughters << 172 void G4LogicalVolume::RemoveDaughter(const G4VPhysicalVolume* p)
119 { 173 {
120 return fDaughtersVolumeType; << 174 G4PhysicalVolumeList::iterator i;
>> 175 for ( i=fDaughters.begin(); i!=fDaughters.end(); ++i )
>> 176 {
>> 177 if (**i==*p)
>> 178 {
>> 179 fDaughters.erase(i);
>> 180 break;
>> 181 }
>> 182 }
>> 183 if (fRegion)
>> 184 {
>> 185 fRegion->RegionModified(true);
>> 186 }
>> 187 fMass = 0.;
121 } 188 }
122 189
123 // ******************************************* 190 // ********************************************************************
124 // DeduceDaughtersType << 191 // ClearDaughters
125 // ******************************************* 192 // ********************************************************************
126 // 193 //
127 inline 194 inline
128 EVolume G4LogicalVolume::DeduceDaughtersType() << 195 void G4LogicalVolume::ClearDaughters()
129 { 196 {
130 EVolume type= kNormal; << 197 G4PhysicalVolumeList::iterator i;
131 G4VPhysicalVolume* pVol; << 198 for ( i=fDaughters.begin(); i!=fDaughters.end(); ++i )
132 << 199 {
133 if ( GetNoDaughters() >= 1 ) << 200 fDaughters.erase(i);
>> 201 }
>> 202 if (fRegion)
134 { 203 {
135 pVol = GetDaughter(0); << 204 fRegion->RegionModified(true);
136 type = pVol->VolumeType(); <<
137 } 205 }
138 return type; << 206 fMass = 0.;
139 } 207 }
140 208
141 // ******************************************* 209 // ********************************************************************
142 // GetMasterSolid << 210 // GetSolid
143 // ******************************************* 211 // ********************************************************************
144 // 212 //
145 inline 213 inline
146 G4VSolid* G4LogicalVolume::GetMasterSolid() co << 214 G4VSolid* G4LogicalVolume::GetSolid() const
147 { 215 {
148 return fSolid; 216 return fSolid;
149 } 217 }
150 218
151 // ******************************************* 219 // ********************************************************************
152 // GetMasterSensitiveDetector << 220 // SetSolid
153 // ******************************************* 221 // ********************************************************************
154 // 222 //
155 inline 223 inline
156 G4VSensitiveDetector* G4LogicalVolume::GetMast << 224 void G4LogicalVolume::SetSolid(G4VSolid *pSolid)
>> 225 {
>> 226 assert(pSolid != 0);
>> 227 fSolid=pSolid;
>> 228 fMass = 0.;
>> 229 }
>> 230
>> 231 // ********************************************************************
>> 232 // GetMaterial
>> 233 // ********************************************************************
>> 234 //
>> 235 inline
>> 236 G4Material* G4LogicalVolume::GetMaterial() const
>> 237 {
>> 238 return fMaterial;
>> 239 }
>> 240
>> 241 // ********************************************************************
>> 242 // SetMaterial
>> 243 // ********************************************************************
>> 244 //
>> 245 inline
>> 246 void G4LogicalVolume::SetMaterial(G4Material *pMaterial)
>> 247 {
>> 248 fMaterial=pMaterial;
>> 249 fMass = 0.;
>> 250 }
>> 251
>> 252 // ********************************************************************
>> 253 // UpdateMaterial
>> 254 // ********************************************************************
>> 255 //
>> 256 inline
>> 257 void G4LogicalVolume::UpdateMaterial(G4Material *pMaterial)
>> 258 {
>> 259 fMaterial=pMaterial;
>> 260 if(fRegion) { fCutsCouple = fRegion->FindCouple(pMaterial); }
>> 261 fMass = 0.;
>> 262 }
>> 263
>> 264 // ********************************************************************
>> 265 // GetSensitiveDetector
>> 266 // ********************************************************************
>> 267 //
>> 268 inline
>> 269 G4VSensitiveDetector* G4LogicalVolume::GetSensitiveDetector() const
157 { 270 {
158 return fSensitiveDetector; 271 return fSensitiveDetector;
159 } 272 }
160 273
161 // ******************************************* 274 // ********************************************************************
>> 275 // SetSensitiveDetector
>> 276 // ********************************************************************
>> 277 //
>> 278 inline
>> 279 void G4LogicalVolume::SetSensitiveDetector(G4VSensitiveDetector* pSDetector)
>> 280 {
>> 281 fSensitiveDetector = pSDetector;
>> 282 }
>> 283
>> 284 // ********************************************************************
162 // GetUserLimits 285 // GetUserLimits
163 // ******************************************* 286 // ********************************************************************
164 // 287 //
165 inline 288 inline
166 G4UserLimits* G4LogicalVolume::GetUserLimits() 289 G4UserLimits* G4LogicalVolume::GetUserLimits() const
167 { 290 {
168 if(fUserLimits != nullptr) return fUserLimit << 291 if(fUserLimits) return fUserLimits;
169 if(fRegion != nullptr) return fRegion->GetUs << 292 if(fRegion) return fRegion->GetUserLimits();
170 return nullptr; << 293 return 0;
171 } 294 }
172 295
173 // ******************************************* 296 // ********************************************************************
174 // SetUserLimits 297 // SetUserLimits
175 // ******************************************* 298 // ********************************************************************
176 // 299 //
177 inline 300 inline
178 void G4LogicalVolume::SetUserLimits(G4UserLimi 301 void G4LogicalVolume::SetUserLimits(G4UserLimits* pULimits)
179 { 302 {
180 fUserLimits = pULimits; 303 fUserLimits = pULimits;
181 } 304 }
182 305
183 // ******************************************* 306 // ********************************************************************
184 // GetVoxelHeader 307 // GetVoxelHeader
185 // ******************************************* 308 // ********************************************************************
186 // 309 //
187 inline 310 inline
188 G4SmartVoxelHeader* G4LogicalVolume::GetVoxelH 311 G4SmartVoxelHeader* G4LogicalVolume::GetVoxelHeader() const
189 { 312 {
190 return fVoxel; 313 return fVoxel;
191 } 314 }
192 315
193 // ******************************************* 316 // ********************************************************************
194 // SetVoxelHeader 317 // SetVoxelHeader
195 // ******************************************* 318 // ********************************************************************
196 // 319 //
197 inline 320 inline
198 void G4LogicalVolume::SetVoxelHeader(G4SmartVo 321 void G4LogicalVolume::SetVoxelHeader(G4SmartVoxelHeader* pVoxel)
199 { 322 {
200 fVoxel = pVoxel; 323 fVoxel = pVoxel;
201 } 324 }
202 325
203 // ******************************************* 326 // ********************************************************************
204 // GetSmartless 327 // GetSmartless
205 // ******************************************* 328 // ********************************************************************
206 // 329 //
207 inline 330 inline
208 G4double G4LogicalVolume::GetSmartless() const 331 G4double G4LogicalVolume::GetSmartless() const
209 { 332 {
210 return fSmartless; 333 return fSmartless;
211 } 334 }
212 335
213 // ******************************************* 336 // ********************************************************************
214 // SetSmartless 337 // SetSmartless
215 // ******************************************* 338 // ********************************************************************
216 // 339 //
217 inline 340 inline
218 void G4LogicalVolume::SetSmartless(G4double sm << 341 void G4LogicalVolume::SetSmartless(G4double s)
219 { 342 {
220 fSmartless = smt; << 343 fSmartless = s;
221 } 344 }
222 345
223 // ******************************************* 346 // ********************************************************************
224 // IsToOptimise 347 // IsToOptimise
225 // ******************************************* 348 // ********************************************************************
226 // 349 //
227 inline 350 inline
228 G4bool G4LogicalVolume::IsToOptimise() const 351 G4bool G4LogicalVolume::IsToOptimise() const
229 { 352 {
230 return fOptimise; 353 return fOptimise;
231 } 354 }
232 355
233 // ******************************************* 356 // ********************************************************************
234 // SetOptimisation 357 // SetOptimisation
235 // ******************************************* 358 // ********************************************************************
236 // 359 //
237 inline 360 inline
238 void G4LogicalVolume::SetOptimisation(G4bool o 361 void G4LogicalVolume::SetOptimisation(G4bool optim)
239 { 362 {
240 fOptimise = optim; 363 fOptimise = optim;
241 } 364 }
242 365
243 // ******************************************* 366 // ********************************************************************
244 // IsRootRegion 367 // IsRootRegion
245 // ******************************************* 368 // ********************************************************************
246 // 369 //
247 inline 370 inline
248 G4bool G4LogicalVolume::IsRootRegion() const 371 G4bool G4LogicalVolume::IsRootRegion() const
249 { 372 {
250 return fRootRegion; 373 return fRootRegion;
251 } 374 }
252 375
253 // ******************************************* 376 // ********************************************************************
254 // SetRegionRootFlag 377 // SetRegionRootFlag
255 // ******************************************* 378 // ********************************************************************
256 // 379 //
257 inline 380 inline
258 void G4LogicalVolume::SetRegionRootFlag(G4bool 381 void G4LogicalVolume::SetRegionRootFlag(G4bool rreg)
259 { 382 {
260 fRootRegion = rreg; 383 fRootRegion = rreg;
261 } 384 }
262 385
263 // ******************************************* 386 // ********************************************************************
264 // IsRegion 387 // IsRegion
265 // ******************************************* 388 // ********************************************************************
266 // 389 //
267 inline 390 inline
268 G4bool G4LogicalVolume::IsRegion() const 391 G4bool G4LogicalVolume::IsRegion() const
269 { 392 {
270 G4bool reg = false; 393 G4bool reg = false;
271 if (fRegion != nullptr) reg = true; << 394 if (fRegion) reg = true;
272 return reg; 395 return reg;
273 } 396 }
274 397
275 // ******************************************* 398 // ********************************************************************
276 // SetRegion 399 // SetRegion
277 // ******************************************* 400 // ********************************************************************
278 // 401 //
279 inline 402 inline
280 void G4LogicalVolume::SetRegion(G4Region* reg) 403 void G4LogicalVolume::SetRegion(G4Region* reg)
281 { 404 {
282 fRegion = reg; 405 fRegion = reg;
283 } 406 }
284 407
285 // ******************************************* 408 // ********************************************************************
286 // GetRegion 409 // GetRegion
287 // ******************************************* 410 // ********************************************************************
288 // 411 //
289 inline 412 inline
290 G4Region* G4LogicalVolume::GetRegion() const 413 G4Region* G4LogicalVolume::GetRegion() const
291 { 414 {
292 return fRegion; 415 return fRegion;
293 } 416 }
294 417
295 // ******************************************* 418 // ********************************************************************
296 // PropagateRegion 419 // PropagateRegion
297 // ******************************************* 420 // ********************************************************************
298 // 421 //
299 inline 422 inline
300 void G4LogicalVolume::PropagateRegion() 423 void G4LogicalVolume::PropagateRegion()
301 { 424 {
302 fRegion->ScanVolumeTree(this, true); 425 fRegion->ScanVolumeTree(this, true);
303 } 426 }
304 427
305 // ******************************************* 428 // ********************************************************************
306 // Lock << 429 // GetMaterialCutsCouple
>> 430 // ********************************************************************
>> 431 //
>> 432 inline
>> 433 const G4MaterialCutsCouple* G4LogicalVolume::GetMaterialCutsCouple() const
>> 434 {
>> 435 return fCutsCouple;
>> 436 }
>> 437
>> 438 // ********************************************************************
>> 439 // SetMaterialCutsCouple
307 // ******************************************* 440 // ********************************************************************
308 // 441 //
309 inline 442 inline
310 void G4LogicalVolume::Lock() << 443 void G4LogicalVolume::SetMaterialCutsCouple(G4MaterialCutsCouple* cuts)
311 { 444 {
312 fLock = true; << 445 fCutsCouple = cuts;
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 // BecomeEnvelopeForFastSimulation
347 // ******************************************* 480 // ********************************************************************
348 // 481 //
349 inline 482 inline
350 const G4VisAttributes* G4LogicalVolume::GetVis << 483 void
>> 484 G4LogicalVolume::BecomeEnvelopeForFastSimulation(G4FastSimulationManager* pPA)
351 { 485 {
352 return fVisAttributes.get(); << 486 SetFastSimulationManager(pPA, true);
>> 487 }
>> 488
>> 489 // ********************************************************************
>> 490 // SetBiasWeight
>> 491 // ********************************************************************
>> 492 //
>> 493 inline
>> 494 void G4LogicalVolume::SetBiasWeight(G4double w)
>> 495 {
>> 496 fBiasWeight = w;
>> 497 }
>> 498
>> 499 // ********************************************************************
>> 500 // GetBiasWeight
>> 501 // ********************************************************************
>> 502 //
>> 503 inline
>> 504 G4double G4LogicalVolume::GetBiasWeight() const
>> 505 {
>> 506 return fBiasWeight;
353 } 507 }
354 508