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.8 2002/05/17 17:59:47 gcosmo Exp $ >> 25 // GEANT4 tag $Name: geant4-05-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 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 fDaughters.push_back(pNewDaughter); >> 106 >> 107 // Propagates the mother's FastSimulationManager pointer. If we are in >> 108 // the World logical volume, propagates only if pointer != 0. >> 109 // (one could also parameterise all the World volume for a particle type!) >> 110 >> 111 G4VPhysicalVolume* myPhysical = pNewDaughter->GetMother(); >> 112 G4VPhysicalVolume* pMotherPhys; >> 113 G4LogicalVolume* pDaughterLogical; >> 114 >> 115 if( myPhysical !=0 ) 107 { 116 { 108 if (*daughter==*p) return true; << 117 pMotherPhys = myPhysical->GetMother(); >> 118 pDaughterLogical = pNewDaughter->GetLogicalVolume(); >> 119 >> 120 if( (pMotherPhys!=0) || ( pMotherPhys==0) && (fFastSimulationManager!=0) ) >> 121 { >> 122 if( pDaughterLogical->GetFastSimulationManager()!=fFastSimulationManager ) >> 123 pDaughterLogical->SetFastSimulationManager(fFastSimulationManager,false); >> 124 } >> 125 else >> 126 { >> 127 pDaughterLogical->SetFastSimulationManager(0, false); >> 128 } >> 129 } >> 130 >> 131 // Propagate the Field Manager, if the daughter has no field Manager. >> 132 // >> 133 pDaughterLogical = pNewDaughter->GetLogicalVolume(); >> 134 G4FieldManager* pDaughterFieldManager = pDaughterLogical->GetFieldManager(); >> 135 >> 136 if( pDaughterFieldManager == 0 ) >> 137 { >> 138 pDaughterLogical->SetFieldManager(fFieldManager, true); 109 } 139 } 110 return false; << 111 } 140 } 112 141 113 // ******************************************* 142 // ******************************************************************** 114 // CharacteriseDaughters << 143 // IsDaughter 115 // ******************************************* 144 // ******************************************************************** 116 // 145 // 117 inline 146 inline 118 EVolume G4LogicalVolume::CharacteriseDaughters << 147 G4bool G4LogicalVolume::IsDaughter(const G4VPhysicalVolume* p) const 119 { 148 { 120 return fDaughtersVolumeType; << 149 G4PhysicalVolumeList::const_iterator i; >> 150 for ( i=fDaughters.begin(); i!=fDaughters.end(); ++i ) >> 151 { >> 152 if (**i==*p) return true; >> 153 } >> 154 return false; 121 } 155 } 122 156 123 // ******************************************* 157 // ******************************************************************** 124 // DeduceDaughtersType << 158 // RemoveDaughter 125 // ******************************************* 159 // ******************************************************************** 126 // 160 // 127 inline 161 inline 128 EVolume G4LogicalVolume::DeduceDaughtersType() << 162 void G4LogicalVolume::RemoveDaughter(const G4VPhysicalVolume* p) 129 { 163 { 130 EVolume type= kNormal; << 164 G4PhysicalVolumeList::iterator i; 131 G4VPhysicalVolume* pVol; << 165 for ( i=fDaughters.begin(); i!=fDaughters.end(); ++i ) 132 << 133 if ( GetNoDaughters() >= 1 ) << 134 { 166 { 135 pVol = GetDaughter(0); << 167 if (**i==*p) 136 type = pVol->VolumeType(); << 168 { >> 169 fDaughters.erase(i); >> 170 break; >> 171 } 137 } 172 } 138 return type; << 139 } 173 } 140 174 141 // ******************************************* 175 // ******************************************************************** 142 // GetMasterSolid << 176 // GetSolid 143 // ******************************************* 177 // ******************************************************************** 144 // 178 // 145 inline 179 inline 146 G4VSolid* G4LogicalVolume::GetMasterSolid() co << 180 G4VSolid* G4LogicalVolume::GetSolid() const 147 { 181 { 148 return fSolid; 182 return fSolid; 149 } 183 } 150 184 151 // ******************************************* 185 // ******************************************************************** 152 // GetMasterSensitiveDetector << 186 // SetSolid 153 // ******************************************* 187 // ******************************************************************** 154 // 188 // 155 inline 189 inline 156 G4VSensitiveDetector* G4LogicalVolume::GetMast << 190 void G4LogicalVolume::SetSolid(G4VSolid *pSolid) 157 { 191 { 158 return fSensitiveDetector; << 192 assert(pSolid != 0); >> 193 fSolid=pSolid; 159 } 194 } 160 195 161 // ******************************************* 196 // ******************************************************************** 162 // GetUserLimits << 197 // GetMaterial 163 // ******************************************* 198 // ******************************************************************** 164 // 199 // 165 inline 200 inline 166 G4UserLimits* G4LogicalVolume::GetUserLimits() << 201 G4Material* G4LogicalVolume::GetMaterial() const 167 { 202 { 168 if(fUserLimits != nullptr) return fUserLimit << 203 return fMaterial; 169 if(fRegion != nullptr) return fRegion->GetUs << 170 return nullptr; << 171 } 204 } 172 205 173 // ******************************************* 206 // ******************************************************************** 174 // SetUserLimits << 207 // SetMaterial 175 // ******************************************* 208 // ******************************************************************** 176 // 209 // 177 inline 210 inline 178 void G4LogicalVolume::SetUserLimits(G4UserLimi << 211 void G4LogicalVolume::SetMaterial(G4Material *pMaterial) 179 { 212 { 180 fUserLimits = pULimits; << 213 fMaterial=pMaterial; 181 } 214 } 182 215 183 // ******************************************* 216 // ******************************************************************** 184 // GetVoxelHeader << 217 // GetSensitiveDetector 185 // ******************************************* 218 // ******************************************************************** 186 // 219 // 187 inline 220 inline 188 G4SmartVoxelHeader* G4LogicalVolume::GetVoxelH << 221 G4VSensitiveDetector* G4LogicalVolume::GetSensitiveDetector() const 189 { 222 { 190 return fVoxel; << 223 return fSensitiveDetector; 191 } 224 } 192 225 193 // ******************************************* 226 // ******************************************************************** 194 // SetVoxelHeader << 227 // SetSensitiveDetector 195 // ******************************************* 228 // ******************************************************************** 196 // 229 // 197 inline 230 inline 198 void G4LogicalVolume::SetVoxelHeader(G4SmartVo << 231 void G4LogicalVolume::SetSensitiveDetector(G4VSensitiveDetector* pSDetector) 199 { 232 { 200 fVoxel = pVoxel; << 233 fSensitiveDetector = pSDetector; 201 } 234 } 202 235 203 // ******************************************* 236 // ******************************************************************** 204 // GetSmartless << 237 // GetUserLimits 205 // ******************************************* 238 // ******************************************************************** 206 // 239 // 207 inline 240 inline 208 G4double G4LogicalVolume::GetSmartless() const << 241 G4UserLimits* G4LogicalVolume::GetUserLimits() const 209 { 242 { 210 return fSmartless; << 243 return fUserLimits; 211 } 244 } 212 245 213 // ******************************************* 246 // ******************************************************************** 214 // SetSmartless << 247 // SetUserLimits 215 // ******************************************* 248 // ******************************************************************** 216 // 249 // 217 inline 250 inline 218 void G4LogicalVolume::SetSmartless(G4double sm << 251 void G4LogicalVolume::SetUserLimits(G4UserLimits* pULimits) 219 { 252 { 220 fSmartless = smt; << 253 fUserLimits = pULimits; 221 } 254 } 222 255 223 // ******************************************* 256 // ******************************************************************** 224 // IsToOptimise << 257 // GetVoxelHeader 225 // ******************************************* 258 // ******************************************************************** 226 // 259 // 227 inline 260 inline 228 G4bool G4LogicalVolume::IsToOptimise() const << 261 G4SmartVoxelHeader* G4LogicalVolume::GetVoxelHeader() const 229 { 262 { 230 return fOptimise; << 263 return fVoxel; 231 } 264 } 232 265 233 // ******************************************* 266 // ******************************************************************** 234 // SetOptimisation << 267 // SetVoxelHeader 235 // ******************************************* 268 // ******************************************************************** 236 // 269 // 237 inline 270 inline 238 void G4LogicalVolume::SetOptimisation(G4bool o << 271 void G4LogicalVolume::SetVoxelHeader(G4SmartVoxelHeader* pVoxel) 239 { 272 { 240 fOptimise = optim; << 273 fVoxel = pVoxel; 241 } 274 } 242 275 243 // ******************************************* 276 // ******************************************************************** 244 // IsRootRegion << 277 // GetSmartless 245 // ******************************************* 278 // ******************************************************************** 246 // 279 // 247 inline 280 inline 248 G4bool G4LogicalVolume::IsRootRegion() const << 281 G4double G4LogicalVolume::GetSmartless() const 249 { 282 { 250 return fRootRegion; << 283 return fSmartless; 251 } 284 } 252 285 253 // ******************************************* 286 // ******************************************************************** 254 // SetRegionRootFlag << 287 // SetSmartless 255 // ******************************************* 288 // ******************************************************************** 256 // 289 // 257 inline 290 inline 258 void G4LogicalVolume::SetRegionRootFlag(G4bool << 291 void G4LogicalVolume::SetSmartless(G4double s) 259 { 292 { 260 fRootRegion = rreg; << 293 fSmartless = s; 261 } 294 } 262 295 263 // ******************************************* 296 // ******************************************************************** 264 // IsRegion << 297 // IsToOptimise 265 // ******************************************* 298 // ******************************************************************** 266 // 299 // 267 inline 300 inline 268 G4bool G4LogicalVolume::IsRegion() const << 301 G4bool G4LogicalVolume::IsToOptimise() const 269 { 302 { 270 G4bool reg = false; << 303 return fOptimise; 271 if (fRegion != nullptr) reg = true; << 272 return reg; << 273 } 304 } 274 305 275 // ******************************************* 306 // ******************************************************************** 276 // SetRegion << 307 // SetOptimisation 277 // ******************************************* 308 // ******************************************************************** 278 // 309 // 279 inline 310 inline 280 void G4LogicalVolume::SetRegion(G4Region* reg) << 311 void G4LogicalVolume::SetOptimisation(G4bool optim) 281 { 312 { 282 fRegion = reg; << 313 fOptimise = optim; 283 } 314 } 284 315 285 // ******************************************* 316 // ******************************************************************** 286 // GetRegion << 317 // Operator == 287 // ******************************************* 318 // ******************************************************************** 288 // 319 // 289 inline 320 inline 290 G4Region* G4LogicalVolume::GetRegion() const << 321 G4bool G4LogicalVolume::operator == ( const G4LogicalVolume& lv) const 291 { 322 { 292 return fRegion; << 323 return (this==&lv) ? true : false; 293 } 324 } 294 325 295 // ******************************************* 326 // ******************************************************************** 296 // PropagateRegion << 327 // GetVisAttributes 297 // ******************************************* 328 // ******************************************************************** 298 // 329 // 299 inline 330 inline 300 void G4LogicalVolume::PropagateRegion() << 331 const G4VisAttributes* G4LogicalVolume::GetVisAttributes () const 301 { 332 { 302 fRegion->ScanVolumeTree(this, true); << 333 return fVisAttributes; 303 } 334 } 304 335 305 // ******************************************* 336 // ******************************************************************** 306 // Lock << 337 // SetVisAttributes 307 // ******************************************* 338 // ******************************************************************** 308 // 339 // 309 inline 340 inline 310 void G4LogicalVolume::Lock() << 341 void G4LogicalVolume::SetVisAttributes (const G4VisAttributes* pVA) 311 { 342 { 312 fLock = true; << 343 fVisAttributes = pVA; 313 } 344 } 314 345 315 // ******************************************* << 316 // Operator == << 317 // ******************************************* << 318 // << 319 inline 346 inline 320 G4bool G4LogicalVolume::operator == ( const G4 << 347 void G4LogicalVolume::SetVisAttributes (const G4VisAttributes& VA) 321 { 348 { 322 return this==&lv; << 349 fVisAttributes = &VA; 323 } 350 } 324 351 325 // ******************************************* 352 // ******************************************************************** 326 // SetBiasWeight << 353 // BecomeEnvelopeForFastSimulation 327 // ******************************************* 354 // ******************************************************************** 328 // 355 // 329 inline 356 inline 330 void G4LogicalVolume::SetBiasWeight(G4double w << 357 void >> 358 G4LogicalVolume::BecomeEnvelopeForFastSimulation(G4FastSimulationManager* pPA) 331 { 359 { 332 fBiasWeight = weight; << 360 SetFastSimulationManager(pPA, true); 333 } 361 } 334 362 335 // ******************************************* 363 // ******************************************************************** 336 // GetBiasWeight << 364 // SetBiasWeight 337 // ******************************************* 365 // ******************************************************************** 338 // 366 // 339 inline 367 inline 340 G4double G4LogicalVolume::GetBiasWeight() cons << 368 void G4LogicalVolume::SetBiasWeight(G4double w) 341 { 369 { 342 return fBiasWeight; << 370 fBiasWeight = w; 343 } 371 } 344 372 345 // ******************************************* 373 // ******************************************************************** 346 // GetVisAttributes << 374 // GetBiasWeight 347 // ******************************************* 375 // ******************************************************************** 348 // 376 // 349 inline 377 inline 350 const G4VisAttributes* G4LogicalVolume::GetVis << 378 G4double G4LogicalVolume::GetBiasWeight() const 351 { 379 { 352 return fVisAttributes.get(); << 380 return fBiasWeight; 353 } 381 } 354 382