| 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.11 2003/03/05 18:07:36 gcosmo Exp $
>> 25 // GEANT4 tag $Name: geant4-05-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 return fFSM; <<
97 } 96 }
98 97
99 // ******************************************* 98 // ********************************************************************
100 // IsDaughter << 99 // AddDaughter
101 // ******************************************* 100 // ********************************************************************
102 // 101 //
103 inline 102 inline
104 G4bool G4LogicalVolume::IsDaughter(const G4VPh << 103 void G4LogicalVolume::AddDaughter(G4VPhysicalVolume* pNewDaughter)
105 { 104 {
106 for (const auto & daughter : fDaughters) << 105 if( !fDaughters.empty() && fDaughters[0]->IsReplicated() )
107 { 106 {
108 if (*daughter==*p) return true; << 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);
109 } 155 }
110 return false; <<
111 } 156 }
112 157
113 // ******************************************* 158 // ********************************************************************
114 // CharacteriseDaughters << 159 // IsDaughter
115 // ******************************************* 160 // ********************************************************************
116 // 161 //
117 inline 162 inline
118 EVolume G4LogicalVolume::CharacteriseDaughters << 163 G4bool G4LogicalVolume::IsDaughter(const G4VPhysicalVolume* p) const
119 { 164 {
120 return fDaughtersVolumeType; << 165 G4PhysicalVolumeList::const_iterator i;
>> 166 for ( i=fDaughters.begin(); i!=fDaughters.end(); ++i )
>> 167 {
>> 168 if (**i==*p) return true;
>> 169 }
>> 170 return false;
121 } 171 }
122 172
123 // ******************************************* 173 // ********************************************************************
124 // DeduceDaughtersType << 174 // RemoveDaughter
125 // ******************************************* 175 // ********************************************************************
126 // 176 //
127 inline 177 inline
128 EVolume G4LogicalVolume::DeduceDaughtersType() << 178 void G4LogicalVolume::RemoveDaughter(const G4VPhysicalVolume* p)
129 { 179 {
130 EVolume type= kNormal; << 180 G4PhysicalVolumeList::iterator i;
131 G4VPhysicalVolume* pVol; << 181 for ( i=fDaughters.begin(); i!=fDaughters.end(); ++i )
132 << 182 {
133 if ( GetNoDaughters() >= 1 ) << 183 if (**i==*p)
>> 184 {
>> 185 fDaughters.erase(i);
>> 186 break;
>> 187 }
>> 188 }
>> 189 if (fRegion)
134 { 190 {
135 pVol = GetDaughter(0); << 191 fRegion->RegionModified(true);
136 type = pVol->VolumeType(); <<
137 } 192 }
138 return type; <<
139 } 193 }
140 194
141 // ******************************************* 195 // ********************************************************************
142 // GetMasterSolid << 196 // GetSolid
143 // ******************************************* 197 // ********************************************************************
144 // 198 //
145 inline 199 inline
146 G4VSolid* G4LogicalVolume::GetMasterSolid() co << 200 G4VSolid* G4LogicalVolume::GetSolid() const
147 { 201 {
148 return fSolid; 202 return fSolid;
149 } 203 }
150 204
151 // ******************************************* 205 // ********************************************************************
152 // GetMasterSensitiveDetector << 206 // SetSolid
>> 207 // ********************************************************************
>> 208 //
>> 209 inline
>> 210 void G4LogicalVolume::SetSolid(G4VSolid *pSolid)
>> 211 {
>> 212 assert(pSolid != 0);
>> 213 fSolid=pSolid;
>> 214 }
>> 215
>> 216 // ********************************************************************
>> 217 // GetMaterial
153 // ******************************************* 218 // ********************************************************************
154 // 219 //
155 inline 220 inline
156 G4VSensitiveDetector* G4LogicalVolume::GetMast << 221 G4Material* G4LogicalVolume::GetMaterial() const
>> 222 {
>> 223 return fMaterial;
>> 224 }
>> 225
>> 226 // ********************************************************************
>> 227 // SetMaterial
>> 228 // ********************************************************************
>> 229 //
>> 230 inline
>> 231 void G4LogicalVolume::SetMaterial(G4Material *pMaterial)
>> 232 {
>> 233 fMaterial=pMaterial;
>> 234 }
>> 235
>> 236 // ********************************************************************
>> 237 // GetSensitiveDetector
>> 238 // ********************************************************************
>> 239 //
>> 240 inline
>> 241 G4VSensitiveDetector* G4LogicalVolume::GetSensitiveDetector() const
157 { 242 {
158 return fSensitiveDetector; 243 return fSensitiveDetector;
159 } 244 }
160 245
161 // ******************************************* 246 // ********************************************************************
>> 247 // SetSensitiveDetector
>> 248 // ********************************************************************
>> 249 //
>> 250 inline
>> 251 void G4LogicalVolume::SetSensitiveDetector(G4VSensitiveDetector* pSDetector)
>> 252 {
>> 253 fSensitiveDetector = pSDetector;
>> 254 }
>> 255
>> 256 // ********************************************************************
162 // GetUserLimits 257 // GetUserLimits
163 // ******************************************* 258 // ********************************************************************
164 // 259 //
165 inline 260 inline
166 G4UserLimits* G4LogicalVolume::GetUserLimits() 261 G4UserLimits* G4LogicalVolume::GetUserLimits() const
167 { 262 {
168 if(fUserLimits != nullptr) return fUserLimit << 263 return fUserLimits;
169 if(fRegion != nullptr) return fRegion->GetUs <<
170 return nullptr; <<
171 } 264 }
172 265
173 // ******************************************* 266 // ********************************************************************
174 // SetUserLimits 267 // SetUserLimits
175 // ******************************************* 268 // ********************************************************************
176 // 269 //
177 inline 270 inline
178 void G4LogicalVolume::SetUserLimits(G4UserLimi 271 void G4LogicalVolume::SetUserLimits(G4UserLimits* pULimits)
179 { 272 {
180 fUserLimits = pULimits; 273 fUserLimits = pULimits;
181 } 274 }
182 275
183 // ******************************************* 276 // ********************************************************************
184 // GetVoxelHeader 277 // GetVoxelHeader
185 // ******************************************* 278 // ********************************************************************
186 // 279 //
187 inline 280 inline
188 G4SmartVoxelHeader* G4LogicalVolume::GetVoxelH 281 G4SmartVoxelHeader* G4LogicalVolume::GetVoxelHeader() const
189 { 282 {
190 return fVoxel; 283 return fVoxel;
191 } 284 }
192 285
193 // ******************************************* 286 // ********************************************************************
194 // SetVoxelHeader 287 // SetVoxelHeader
195 // ******************************************* 288 // ********************************************************************
196 // 289 //
197 inline 290 inline
198 void G4LogicalVolume::SetVoxelHeader(G4SmartVo 291 void G4LogicalVolume::SetVoxelHeader(G4SmartVoxelHeader* pVoxel)
199 { 292 {
200 fVoxel = pVoxel; 293 fVoxel = pVoxel;
201 } 294 }
202 295
203 // ******************************************* 296 // ********************************************************************
204 // GetSmartless 297 // GetSmartless
205 // ******************************************* 298 // ********************************************************************
206 // 299 //
207 inline 300 inline
208 G4double G4LogicalVolume::GetSmartless() const 301 G4double G4LogicalVolume::GetSmartless() const
209 { 302 {
210 return fSmartless; 303 return fSmartless;
211 } 304 }
212 305
213 // ******************************************* 306 // ********************************************************************
214 // SetSmartless 307 // SetSmartless
215 // ******************************************* 308 // ********************************************************************
216 // 309 //
217 inline 310 inline
218 void G4LogicalVolume::SetSmartless(G4double sm << 311 void G4LogicalVolume::SetSmartless(G4double s)
219 { 312 {
220 fSmartless = smt; << 313 fSmartless = s;
221 } 314 }
222 315
223 // ******************************************* 316 // ********************************************************************
224 // IsToOptimise 317 // IsToOptimise
225 // ******************************************* 318 // ********************************************************************
226 // 319 //
227 inline 320 inline
228 G4bool G4LogicalVolume::IsToOptimise() const 321 G4bool G4LogicalVolume::IsToOptimise() const
229 { 322 {
230 return fOptimise; 323 return fOptimise;
231 } 324 }
232 325
233 // ******************************************* 326 // ********************************************************************
234 // SetOptimisation 327 // SetOptimisation
235 // ******************************************* 328 // ********************************************************************
236 // 329 //
237 inline 330 inline
238 void G4LogicalVolume::SetOptimisation(G4bool o 331 void G4LogicalVolume::SetOptimisation(G4bool optim)
239 { 332 {
240 fOptimise = optim; 333 fOptimise = optim;
241 } 334 }
242 335
243 // ******************************************* 336 // ********************************************************************
244 // IsRootRegion 337 // IsRootRegion
245 // ******************************************* 338 // ********************************************************************
246 // 339 //
247 inline 340 inline
248 G4bool G4LogicalVolume::IsRootRegion() const 341 G4bool G4LogicalVolume::IsRootRegion() const
249 { 342 {
250 return fRootRegion; 343 return fRootRegion;
251 } 344 }
252 345
253 // ******************************************* 346 // ********************************************************************
254 // SetRegionRootFlag 347 // SetRegionRootFlag
255 // ******************************************* 348 // ********************************************************************
256 // 349 //
257 inline 350 inline
258 void G4LogicalVolume::SetRegionRootFlag(G4bool 351 void G4LogicalVolume::SetRegionRootFlag(G4bool rreg)
259 { 352 {
260 fRootRegion = rreg; 353 fRootRegion = rreg;
261 } 354 }
262 355
263 // ******************************************* 356 // ********************************************************************
264 // IsRegion 357 // IsRegion
265 // ******************************************* 358 // ********************************************************************
266 // 359 //
267 inline 360 inline
268 G4bool G4LogicalVolume::IsRegion() const 361 G4bool G4LogicalVolume::IsRegion() const
269 { 362 {
270 G4bool reg = false; 363 G4bool reg = false;
271 if (fRegion != nullptr) reg = true; << 364 if (fRegion) reg = true;
272 return reg; 365 return reg;
273 } 366 }
274 367
275 // ******************************************* 368 // ********************************************************************
276 // SetRegion 369 // SetRegion
277 // ******************************************* 370 // ********************************************************************
278 // 371 //
279 inline 372 inline
280 void G4LogicalVolume::SetRegion(G4Region* reg) 373 void G4LogicalVolume::SetRegion(G4Region* reg)
281 { 374 {
282 fRegion = reg; 375 fRegion = reg;
283 } 376 }
284 377
285 // ******************************************* 378 // ********************************************************************
286 // GetRegion 379 // GetRegion
287 // ******************************************* 380 // ********************************************************************
288 // 381 //
289 inline 382 inline
290 G4Region* G4LogicalVolume::GetRegion() const 383 G4Region* G4LogicalVolume::GetRegion() const
291 { 384 {
292 return fRegion; 385 return fRegion;
293 } 386 }
294 387
295 // ******************************************* 388 // ********************************************************************
296 // PropagateRegion 389 // PropagateRegion
297 // ******************************************* 390 // ********************************************************************
298 // 391 //
299 inline 392 inline
300 void G4LogicalVolume::PropagateRegion() 393 void G4LogicalVolume::PropagateRegion()
301 { 394 {
302 fRegion->ScanVolumeTree(this, true); 395 fRegion->ScanVolumeTree(this, true);
303 } 396 }
304 397
305 // ******************************************* 398 // ********************************************************************
306 // Lock << 399 // GetMaterialCutsCouple
307 // ******************************************* 400 // ********************************************************************
308 // 401 //
309 inline 402 inline
310 void G4LogicalVolume::Lock() << 403 const G4MaterialCutsCouple* G4LogicalVolume::GetMaterialCutsCouple() const
311 { 404 {
312 fLock = true; << 405 return fCutsCouple;
>> 406 }
>> 407
>> 408 // ********************************************************************
>> 409 // SetMaterialCutsCouple
>> 410 // ********************************************************************
>> 411 //
>> 412 inline
>> 413 void G4LogicalVolume::SetMaterialCutsCouple(G4MaterialCutsCouple* cuts)
>> 414 {
>> 415 fCutsCouple = cuts;
313 } 416 }
314 417
315 // ******************************************* 418 // ********************************************************************
316 // Operator == 419 // Operator ==
317 // ******************************************* 420 // ********************************************************************
318 // 421 //
319 inline 422 inline
320 G4bool G4LogicalVolume::operator == ( const G4 423 G4bool G4LogicalVolume::operator == ( const G4LogicalVolume& lv) const
321 { 424 {
322 return this==&lv; << 425 return (this==&lv) ? true : false;
323 } 426 }
324 427
325 // ******************************************* 428 // ********************************************************************
326 // SetBiasWeight << 429 // GetVisAttributes
327 // ******************************************* 430 // ********************************************************************
328 // 431 //
329 inline 432 inline
330 void G4LogicalVolume::SetBiasWeight(G4double w << 433 const G4VisAttributes* G4LogicalVolume::GetVisAttributes () const
331 { 434 {
332 fBiasWeight = weight; << 435 return fVisAttributes;
333 } 436 }
334 437
335 // ******************************************* 438 // ********************************************************************
336 // GetBiasWeight << 439 // SetVisAttributes
337 // ******************************************* 440 // ********************************************************************
338 // 441 //
339 inline 442 inline
340 G4double G4LogicalVolume::GetBiasWeight() cons << 443 void G4LogicalVolume::SetVisAttributes (const G4VisAttributes* pVA)
341 { 444 {
342 return fBiasWeight; << 445 fVisAttributes = pVA;
>> 446 }
>> 447
>> 448 inline
>> 449 void G4LogicalVolume::SetVisAttributes (const G4VisAttributes& VA)
>> 450 {
>> 451 fVisAttributes = &VA;
343 } 452 }
344 453
345 // ******************************************* 454 // ********************************************************************
346 // GetVisAttributes << 455 // BecomeEnvelopeForFastSimulation
347 // ******************************************* 456 // ********************************************************************
348 // 457 //
349 inline 458 inline
350 const G4VisAttributes* G4LogicalVolume::GetVis << 459 void
>> 460 G4LogicalVolume::BecomeEnvelopeForFastSimulation(G4FastSimulationManager* pPA)
351 { 461 {
352 return fVisAttributes.get(); << 462 SetFastSimulationManager(pPA, true);
>> 463 }
>> 464
>> 465 // ********************************************************************
>> 466 // SetBiasWeight
>> 467 // ********************************************************************
>> 468 //
>> 469 inline
>> 470 void G4LogicalVolume::SetBiasWeight(G4double w)
>> 471 {
>> 472 fBiasWeight = w;
>> 473 }
>> 474
>> 475 // ********************************************************************
>> 476 // GetBiasWeight
>> 477 // ********************************************************************
>> 478 //
>> 479 inline
>> 480 G4double G4LogicalVolume::GetBiasWeight() const
>> 481 {
>> 482 return fBiasWeight;
353 } 483 }
354 484