| 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.19 2004/11/16 13:43:12 gcosmo Exp $
>> 25 // GEANT4 tag $Name: geant4-07-00-cand-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 <<
133 if ( GetNoDaughters() >= 1 ) <<
134 { 199 {
135 pVol = GetDaughter(0); << 200 fDaughters.erase(i);
136 type = pVol->VolumeType(); <<
137 } 201 }
138 return type; << 202 if (fRegion)
>> 203 {
>> 204 fRegion->RegionModified(true);
>> 205 }
>> 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
>> 221 // ********************************************************************
>> 222 //
>> 223 inline
>> 224 void G4LogicalVolume::SetSolid(G4VSolid *pSolid)
>> 225 {
>> 226 assert(pSolid != 0);
>> 227 fSolid=pSolid;
>> 228 fMass = 0.;
>> 229 }
>> 230
>> 231 // ********************************************************************
>> 232 // GetMaterial
153 // ******************************************* 233 // ********************************************************************
154 // 234 //
155 inline 235 inline
156 G4VSensitiveDetector* G4LogicalVolume::GetMast << 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 return fUserLimits;
169 if(fRegion != nullptr) return fRegion->GetUs <<
170 return nullptr; <<
171 } 292 }
172 293
173 // ******************************************* 294 // ********************************************************************
174 // SetUserLimits 295 // SetUserLimits
175 // ******************************************* 296 // ********************************************************************
176 // 297 //
177 inline 298 inline
178 void G4LogicalVolume::SetUserLimits(G4UserLimi 299 void G4LogicalVolume::SetUserLimits(G4UserLimits* pULimits)
179 { 300 {
180 fUserLimits = pULimits; 301 fUserLimits = pULimits;
181 } 302 }
182 303
183 // ******************************************* 304 // ********************************************************************
184 // GetVoxelHeader 305 // GetVoxelHeader
185 // ******************************************* 306 // ********************************************************************
186 // 307 //
187 inline 308 inline
188 G4SmartVoxelHeader* G4LogicalVolume::GetVoxelH 309 G4SmartVoxelHeader* G4LogicalVolume::GetVoxelHeader() const
189 { 310 {
190 return fVoxel; 311 return fVoxel;
191 } 312 }
192 313
193 // ******************************************* 314 // ********************************************************************
194 // SetVoxelHeader 315 // SetVoxelHeader
195 // ******************************************* 316 // ********************************************************************
196 // 317 //
197 inline 318 inline
198 void G4LogicalVolume::SetVoxelHeader(G4SmartVo 319 void G4LogicalVolume::SetVoxelHeader(G4SmartVoxelHeader* pVoxel)
199 { 320 {
200 fVoxel = pVoxel; 321 fVoxel = pVoxel;
201 } 322 }
202 323
203 // ******************************************* 324 // ********************************************************************
204 // GetSmartless 325 // GetSmartless
205 // ******************************************* 326 // ********************************************************************
206 // 327 //
207 inline 328 inline
208 G4double G4LogicalVolume::GetSmartless() const 329 G4double G4LogicalVolume::GetSmartless() const
209 { 330 {
210 return fSmartless; 331 return fSmartless;
211 } 332 }
212 333
213 // ******************************************* 334 // ********************************************************************
214 // SetSmartless 335 // SetSmartless
215 // ******************************************* 336 // ********************************************************************
216 // 337 //
217 inline 338 inline
218 void G4LogicalVolume::SetSmartless(G4double sm << 339 void G4LogicalVolume::SetSmartless(G4double s)
219 { 340 {
220 fSmartless = smt; << 341 fSmartless = s;
221 } 342 }
222 343
223 // ******************************************* 344 // ********************************************************************
224 // IsToOptimise 345 // IsToOptimise
225 // ******************************************* 346 // ********************************************************************
226 // 347 //
227 inline 348 inline
228 G4bool G4LogicalVolume::IsToOptimise() const 349 G4bool G4LogicalVolume::IsToOptimise() const
229 { 350 {
230 return fOptimise; 351 return fOptimise;
231 } 352 }
232 353
233 // ******************************************* 354 // ********************************************************************
234 // SetOptimisation 355 // SetOptimisation
235 // ******************************************* 356 // ********************************************************************
236 // 357 //
237 inline 358 inline
238 void G4LogicalVolume::SetOptimisation(G4bool o 359 void G4LogicalVolume::SetOptimisation(G4bool optim)
239 { 360 {
240 fOptimise = optim; 361 fOptimise = optim;
241 } 362 }
242 363
243 // ******************************************* 364 // ********************************************************************
244 // IsRootRegion 365 // IsRootRegion
245 // ******************************************* 366 // ********************************************************************
246 // 367 //
247 inline 368 inline
248 G4bool G4LogicalVolume::IsRootRegion() const 369 G4bool G4LogicalVolume::IsRootRegion() const
249 { 370 {
250 return fRootRegion; 371 return fRootRegion;
251 } 372 }
252 373
253 // ******************************************* 374 // ********************************************************************
254 // SetRegionRootFlag 375 // SetRegionRootFlag
255 // ******************************************* 376 // ********************************************************************
256 // 377 //
257 inline 378 inline
258 void G4LogicalVolume::SetRegionRootFlag(G4bool 379 void G4LogicalVolume::SetRegionRootFlag(G4bool rreg)
259 { 380 {
260 fRootRegion = rreg; 381 fRootRegion = rreg;
261 } 382 }
262 383
263 // ******************************************* 384 // ********************************************************************
264 // IsRegion 385 // IsRegion
265 // ******************************************* 386 // ********************************************************************
266 // 387 //
267 inline 388 inline
268 G4bool G4LogicalVolume::IsRegion() const 389 G4bool G4LogicalVolume::IsRegion() const
269 { 390 {
270 G4bool reg = false; 391 G4bool reg = false;
271 if (fRegion != nullptr) reg = true; << 392 if (fRegion) reg = true;
272 return reg; 393 return reg;
273 } 394 }
274 395
275 // ******************************************* 396 // ********************************************************************
276 // SetRegion 397 // SetRegion
277 // ******************************************* 398 // ********************************************************************
278 // 399 //
279 inline 400 inline
280 void G4LogicalVolume::SetRegion(G4Region* reg) 401 void G4LogicalVolume::SetRegion(G4Region* reg)
281 { 402 {
282 fRegion = reg; 403 fRegion = reg;
283 } 404 }
284 405
285 // ******************************************* 406 // ********************************************************************
286 // GetRegion 407 // GetRegion
287 // ******************************************* 408 // ********************************************************************
288 // 409 //
289 inline 410 inline
290 G4Region* G4LogicalVolume::GetRegion() const 411 G4Region* G4LogicalVolume::GetRegion() const
291 { 412 {
292 return fRegion; 413 return fRegion;
293 } 414 }
294 415
295 // ******************************************* 416 // ********************************************************************
296 // PropagateRegion 417 // PropagateRegion
297 // ******************************************* 418 // ********************************************************************
298 // 419 //
299 inline 420 inline
300 void G4LogicalVolume::PropagateRegion() 421 void G4LogicalVolume::PropagateRegion()
301 { 422 {
302 fRegion->ScanVolumeTree(this, true); 423 fRegion->ScanVolumeTree(this, true);
303 } 424 }
304 425
305 // ******************************************* 426 // ********************************************************************
306 // Lock << 427 // GetMaterialCutsCouple
307 // ******************************************* 428 // ********************************************************************
308 // 429 //
309 inline 430 inline
310 void G4LogicalVolume::Lock() << 431 const G4MaterialCutsCouple* G4LogicalVolume::GetMaterialCutsCouple() const
311 { 432 {
312 fLock = true; << 433 return fCutsCouple;
>> 434 }
>> 435
>> 436 // ********************************************************************
>> 437 // SetMaterialCutsCouple
>> 438 // ********************************************************************
>> 439 //
>> 440 inline
>> 441 void G4LogicalVolume::SetMaterialCutsCouple(G4MaterialCutsCouple* cuts)
>> 442 {
>> 443 fCutsCouple = cuts;
313 } 444 }
314 445
315 // ******************************************* 446 // ********************************************************************
316 // Operator == 447 // Operator ==
317 // ******************************************* 448 // ********************************************************************
318 // 449 //
319 inline 450 inline
320 G4bool G4LogicalVolume::operator == ( const G4 451 G4bool G4LogicalVolume::operator == ( const G4LogicalVolume& lv) const
321 { 452 {
322 return this==&lv; << 453 return (this==&lv) ? true : false;
323 } 454 }
324 455
325 // ******************************************* 456 // ********************************************************************
326 // SetBiasWeight << 457 // GetVisAttributes
327 // ******************************************* 458 // ********************************************************************
328 // 459 //
329 inline 460 inline
330 void G4LogicalVolume::SetBiasWeight(G4double w << 461 const G4VisAttributes* G4LogicalVolume::GetVisAttributes () const
331 { 462 {
332 fBiasWeight = weight; << 463 return fVisAttributes;
333 } 464 }
334 465
335 // ******************************************* 466 // ********************************************************************
336 // GetBiasWeight << 467 // SetVisAttributes
337 // ******************************************* 468 // ********************************************************************
338 // 469 //
339 inline 470 inline
340 G4double G4LogicalVolume::GetBiasWeight() cons << 471 void G4LogicalVolume::SetVisAttributes (const G4VisAttributes* pVA)
341 { 472 {
342 return fBiasWeight; << 473 fVisAttributes = pVA;
>> 474 }
>> 475
>> 476 inline
>> 477 void G4LogicalVolume::SetVisAttributes (const G4VisAttributes& VA)
>> 478 {
>> 479 fVisAttributes = &VA;
343 } 480 }
344 481
345 // ******************************************* 482 // ********************************************************************
346 // GetVisAttributes << 483 // BecomeEnvelopeForFastSimulation
347 // ******************************************* 484 // ********************************************************************
348 // 485 //
349 inline 486 inline
350 const G4VisAttributes* G4LogicalVolume::GetVis << 487 void
>> 488 G4LogicalVolume::BecomeEnvelopeForFastSimulation(G4FastSimulationManager* pPA)
>> 489 {
>> 490 SetFastSimulationManager(pPA, true);
>> 491 }
>> 492
>> 493 // ********************************************************************
>> 494 // SetBiasWeight
>> 495 // ********************************************************************
>> 496 //
>> 497 inline
>> 498 void G4LogicalVolume::SetBiasWeight(G4double w)
351 { 499 {
352 return fVisAttributes.get(); << 500 fBiasWeight = w;
>> 501 }
>> 502
>> 503 // ********************************************************************
>> 504 // GetBiasWeight
>> 505 // ********************************************************************
>> 506 //
>> 507 inline
>> 508 G4double G4LogicalVolume::GetBiasWeight() const
>> 509 {
>> 510 return fBiasWeight;
353 } 511 }
354 512