Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // G4ImportanceConfigurator 27 // 28 // Author: Michael Dressel, CERN 29 // ---------------------------------------------------------------------- 30 31 #include "G4ImportanceConfigurator.hh" 32 33 #include "G4ImportanceProcess.hh" 34 #include "G4ProcessPlacer.hh" 35 #include "G4ImportanceAlgorithm.hh" 36 37 #include "G4TransportationManager.hh" 38 #include "G4AutoLock.hh" 39 40 #ifdef G4MULTITHREADED 41 namespace 42 { 43 G4Mutex BiasConfigMutex = G4MUTEX_INITIALIZER; 44 } 45 #endif 46 47 G4ImportanceConfigurator:: 48 G4ImportanceConfigurator(const G4VPhysicalVolume* worldvolume, 49 const G4String& particlename, 50 G4VIStore& istore, 51 const G4VImportanceAlgorithm* ialg, G4bool para) 52 : fWorld(worldvolume), 53 fWorldName(worldvolume->GetName()), 54 fPlacer(particlename), 55 fIStore(istore), 56 fDeleteIalg( ( ! ialg) ), 57 fIalgorithm(( (fDeleteIalg == true) ? new G4ImportanceAlgorithm : ialg )), 58 paraflag(para) 59 {;} 60 61 G4ImportanceConfigurator:: 62 G4ImportanceConfigurator(const G4String& worldvolumeName, 63 const G4String& particlename, 64 G4VIStore& istore, 65 const G4VImportanceAlgorithm* ialg, G4bool para) 66 : fWorldName(worldvolumeName), 67 fPlacer(particlename), 68 fIStore(istore), 69 fDeleteIalg( ( ! ialg) ), 70 fIalgorithm(( (fDeleteIalg == true) ? new G4ImportanceAlgorithm : ialg )), 71 paraflag(para) 72 { 73 fWorld = G4TransportationManager::GetTransportationManager()->GetNavigatorForTracking()->GetWorldVolume(); 74 if(paraflag) fWorld = G4TransportationManager::GetTransportationManager()->GetParallelWorld(fWorldName); 75 } 76 77 G4ImportanceConfigurator::~G4ImportanceConfigurator() 78 { 79 if (fImportanceProcess) 80 { 81 fPlacer.RemoveProcess(fImportanceProcess); 82 delete fImportanceProcess; 83 } 84 if (fDeleteIalg) 85 { 86 delete fIalgorithm; 87 } 88 } 89 90 void 91 G4ImportanceConfigurator::Configure(G4VSamplerConfigurator* preConf) 92 { 93 G4cout << "G4ImportanceConfigurator:: entering importance configure, paraflag " << paraflag << G4endl; 94 const G4VTrackTerminator *terminator = 0; 95 if (preConf) 96 { 97 terminator = preConf->GetTrackTerminator(); 98 } 99 100 #ifdef G4MULTITHREADED 101 G4AutoLock l(&BiasConfigMutex); 102 #endif 103 fImportanceProcess = 104 new G4ImportanceProcess(*fIalgorithm, 105 fIStore, 106 terminator,"ImportanceProcess",paraflag); 107 108 if (!fImportanceProcess) 109 { 110 G4Exception("G4ImportanceConfigurator::Configure()", 111 "FatalError", FatalException, 112 "Failed allocation of G4ImportanceProcess !"); 113 } 114 115 if(paraflag) fImportanceProcess->SetParallelWorld(fWorld->GetName()); 116 #ifdef G4MULTITHREADED 117 l.unlock(); 118 #endif 119 fPlacer.AddProcessAsSecondDoIt(fImportanceProcess); 120 } 121 122 const G4VTrackTerminator* 123 G4ImportanceConfigurator::GetTrackTerminator() const 124 { 125 return fImportanceProcess; 126 } 127 128 void G4ImportanceConfigurator::SetWorldName(const G4String& name) 129 { 130 G4cout << " G4ImportanceConfigurator:: setting world name: " << name << G4endl; 131 fWorldName = name; 132 } 133