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