67 #ifdef theParticleIterator
68 #undef theParticleIterator
74 {
return fRunManagerKernel; }
77 : physicsList(0),currentWorld(0),
78 geometryInitialized(false),physicsInitialized(false),
79 geometryToBeOptimized(true),
80 physicsNeedsToBeReBuilt(true),verboseLevel(0),
81 numberOfParallelWorld(0),geometryNeedsToBeClosed(true),
82 numberOfStaticAllocators(0)
85 InvalidOperationDetection();
88 if(allocList) numberOfStaticAllocators = allocList->
Size();
92 G4Exception(
"G4RunManagerKernel::G4RunManagerKernel()",
"Run0001",
93 FatalException,
"More than one G4RunManagerKernel is constructed.");
95 fRunManagerKernel =
this;
102 ED<<
"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"<<
G4endl;
103 ED<<
" G4RunManagerKernel fatal exception"<<
G4endl;
104 ED<<
" -- Following particles have already been registered"<<
G4endl;
105 ED<<
" before G4RunManagerKernel is instantiated."<<
G4endl;
106 for(
int i=0;i<particleTable->
entries();i++)
108 ED<<
"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"<<
G4endl;
109 G4Exception(
"G4RunManagerKernel::G4RunManagerKernel()",
"Run0002",
129 vs = vs.substr(1,vs.size()-2);
130 versionString =
" Geant4 version ";
132 versionString +=
" ";
133 versionString += G4Date;
135 <<
"*************************************************************" <<
G4endl
136 << versionString <<
G4endl
137 <<
" Copyright : Geant4 Collaboration" <<
G4endl
138 <<
" Reference : NIM A 506 (2003), 250-303" <<
G4endl
139 <<
" WWW : http://cern.ch/geant4" <<
G4endl
140 <<
"*************************************************************" <<
G4endl
145 : physicsList(0),currentWorld(0),
146 geometryInitialized(false),physicsInitialized(false),
147 geometryToBeOptimized(true),
148 physicsNeedsToBeReBuilt(true),verboseLevel(0),
149 numberOfParallelWorld(0),geometryNeedsToBeClosed(true),
150 numberOfStaticAllocators(0)
153 #ifndef G4MULTITHREADED
155 msg<<
"Geant4 code is compiled without multi-threading support (-DG4MULTITHREADED is set to off).";
156 msg<<
" This type of RunManagerKernel can only be used in mult-threaded applications.";
161 InvalidOperationDetection();
165 if(fRunManagerKernel)
167 G4Exception(
"G4RunManagerKernel::G4RunManagerKernel()",
"Run0001",
168 FatalException,
"More than one G4RunManagerKernel is constructed.");
170 fRunManagerKernel =
this;
193 msgx<<
" This type of RunManagerKernel can only be used in mult-threaded applications.";
203 vs = vs.substr(1,vs.size()-2);
207 versionString =
" Geant4 version ";
209 versionString +=
" ";
210 versionString += G4Date;
212 <<
"*************************************************************" <<
G4endl
213 << versionString <<
G4endl
214 <<
" << in Multi-threaded mode >> " <<
G4endl
215 <<
" Copyright : Geant4 Collaboration" <<
G4endl
216 <<
" Reference : NIM A 506 (2003), 250-303" <<
G4endl
217 <<
" WWW : http://cern.ch/geant4" <<
G4endl
218 <<
"*************************************************************" <<
G4endl
222 versionString =
" Local thread RunManagerKernel version ";
225 <<
"^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^" << G4endl
226 << versionString << G4endl
227 <<
"^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^" << G4endl
244 "Default world region should have a unique logical volume.");
246 std::vector<G4LogicalVolume*>::iterator lvItr
250 <<
"Obsolete world logical volume is removed from the default region." <<
G4endl;
261 if(verboseLevel>0)
G4cout <<
"G4 kernel has come to Quit state." <<
G4endl;
273 if(verboseLevel>1)
G4cout <<
"G4SDManager deleted." <<
G4endl;
276 if(verboseLevel>1)
G4cout <<
"EventManager deleted." <<
G4endl;
279 if(verboseLevel>1)
G4cout <<
"Units table cleared." <<
G4endl;
283 allocList->
Destroy(numberOfStaticAllocators,verboseLevel);
285 if(verboseLevel>1)
G4cout <<
"G4Allocator objects are deleted." <<
G4endl;
291 G4cout <<
"Thread-local UImanager is to be deleted." << G4endl
292 <<
"There should not be any thread-local G4cout/G4cerr hereafter."
296 if(pUImanager)
delete pUImanager;
297 if(verboseLevel>1)
G4cout <<
"UImanager deleted." <<
G4endl;
299 delete pStateManager;
300 if(verboseLevel>1)
G4cout <<
"StateManager deleted." <<
G4endl;
301 delete defaultExceptionHandler;
302 if(verboseLevel>0)
G4cout <<
"RunManagerKernel is deleted. Good bye :)" <<
G4endl;
303 fRunManagerKernel = 0;
311 G4MTRunManager::masterWorlds_t::iterator itrMW = masterWorlds.begin();
312 for(;itrMW!=masterWorlds.end();itrMW++)
330 G4Exception(
"G4RunManagerKernel::DefineWorldVolume",
331 "DefineWorldVolumeAtIncorrectState",
333 "Geant4 kernel is not PreInit or Idle state : Method ignored.");
334 if(verboseLevel>1)
G4cerr <<
"Current application state is "
339 currentWorld = worldVol;
343 G4MTRunManager::masterWorlds_t::iterator itrMW = masterWorlds.begin();
344 for(;itrMW!=masterWorlds.end();itrMW++)
346 if((*itrMW).first == 0)
348 if((*itrMW).second != currentWorld)
350 G4Exception(
"G4RunManagerKernel::WorkerDefineWorldVolume",
"RUN3091",
370 geometryInitialized =
true;
382 G4Exception(
"G4RunManagerKernel::DefineWorldVolume",
385 "Geant4 kernel is not PreInit or Idle state : Method ignored.");
395 ED <<
"The world volume has a user-defined region <"
398 ED <<
"World would have a default region assigned by RunManagerKernel."
408 currentWorld = worldVol;
416 <<
" is registered to the default region." <<
G4endl;
430 geometryInitialized =
true;
443 G4cout <<
"List of instantiated particles ============================================" <<
G4endl;
445 for(
G4int i=0;i<nPtcl;i++)
504 G4Exception(
"G4RunManagerKernel::InitializePhysics",
506 "Geant4 kernel is not PreInit or Idle state : Method ignored.");
512 G4Exception(
"G4RunManagerKernel::InitializePhysics",
514 "G4VUserPhysicsList is not defined");
518 if(verboseLevel>1)
G4cout <<
"physicsList->Construct() start." <<
G4endl;
522 if(verboseLevel>1)
G4cout <<
"physicsList->CheckParticleList() start." <<
G4endl;
526 if(verboseLevel>1)
G4cout <<
"physicsList->setCut() start." <<
G4endl;
532 static G4bool createIsomerOnlyOnce =
false;
535 if(!createIsomerOnlyOnce)
537 createIsomerOnlyOnce =
true;
554 physicsInitialized =
true;
564 if(!geometryInitialized)
566 G4Exception(
"G4RunManagerKernel::RunInitialization",
569 "Geometry has not yet initialized : method ignored.");
573 if(!physicsInitialized)
575 G4Exception(
"G4RunManagerKernel::RunInitialization",
578 "Physics has not yet initialized : method ignored.");
584 G4Exception(
"G4RunManagerKernel::RunInitialization",
587 "Geant4 kernel not in Idle state : method ignored.");
632 if(verboseLevel>1)
G4cout <<
"Start closing geometry." <<
G4endl;
635 geomManager->
CloseGeometry(geometryToBeOptimized, verboseLevel>1);
649 "Geant4 kernel not in Idle state : method ignored.");
665 || physicsNeedsToBeReBuilt)
669 physicsNeedsToBeReBuilt =
false;
681 std::vector<G4VPhysicalVolume*>::iterator wItr;
693 for(
size_t iw=0;iw<nWorlds;iw++)
697 if(*wItr==currentWorld)
712 <<
"> does not have specific production cuts," <<
G4endl
713 <<
"even though it appears in the current tracking world." <<
G4endl;
714 G4cerr <<
"Default cuts are used for this region." <<
G4endl;
721 ->GetDefaultProductionCuts());
731 for(
size_t iw=0;iw<nWorlds;iw++)
734 if(*wItr!=currentWorld)
770 G4cout <<
" -- appears in <"
774 {
G4cout <<
" -- is not associated to any world."; }
777 {
G4cout <<
" This region is in the mass world." <<
G4endl; }
779 {
G4cout <<
" This region is in the parallel world." <<
G4endl; }
781 G4cout <<
" Root logical volume(s) : ";
784 for(
size_t j=0;j<nRootLV;j++)
785 {
G4cout << (*lvItr)->GetName() <<
" "; lvItr++; }
793 G4cout <<
" Materials : ";
796 for(
size_t iMate=0;iMate<nMaterial;iMate++)
798 G4cout << (*mItr)->GetName() <<
" ";
806 <<
"> does not have specific production cuts." <<
G4endl;
807 G4cerr <<
"Default cuts are used for this region." <<
G4endl;
813 G4cout <<
" Production cuts : "
828 void G4RunManagerKernel::CheckRegularGeometry()
831 for(G4LogicalVolumeStore::iterator pos=store->begin(); pos!=store->end(); pos++)
833 if((*pos)&&((*pos)->GetNoDaughters()==1))
835 if((*pos)->GetDaughter(0)->IsRegularStructure())
849 G4bool G4RunManagerKernel::ConfirmCoupledTransportation()
853 theParticleIterator->
reset();
854 while((*theParticleIterator)())
869 void G4RunManagerKernel::SetScoreSplitter()
877 if( ! InitSplitter ) {
880 theParticleIterator->
reset();
881 while( (*theParticleIterator)() )
891 G4cout <<
"G4RunManagerKernel -- G4ScoreSplittingProcess is appended to all particles." <<
G4endl;
900 theParticleIterator->
reset();
902 while((*theParticleIterator)())
909 for (
G4int idx = 0 ; idx<procs.
size() ; ++idx)
911 const G4VProcess* masterP = procs[idx]->GetMasterProcess();
916 procs[idx]->SetMasterProcess(const_cast<G4VProcess*>(procs[idx]));
virtual void GeometryHasChanged()=0
static G4Geantino * GeantinoDefinition()
G4VPhysicalVolume * IsWorldExisting(const G4String &worldName)
G4int GetParticleDefinitionID() const
G4Region * defaultRegionForParallelWorld
void InitializeLightIons()
G4ProductionCuts * GetProductionCuts() const
std::vector< G4Material * >::const_iterator GetMaterialIterator() const
std::map< G4int, G4VPhysicalVolume * > masterWorlds_t
virtual void SetupShadowProcess() const
G4Region * GetRegion(const G4String &name, G4bool verbose=true) const
void DumpCutValuesTableIfRequested()
const G4String & GetName() const
void SetPhysics(G4VUserPhysicsList *uPhys)
void SetParticleDefinitionID(G4int id=-1)
void UpdateMaterialList(G4VPhysicalVolume *currentWorld=0)
void DumpTable(const G4String &particle_name="ALL")
std::ostringstream G4ExceptionDescription
void AddRootLogicalVolume(G4LogicalVolume *lv)
void UsedInParallelGeometry(G4bool val=true)
static G4VVisManager * GetConcreteInstance()
static G4AllocatorList * GetAllocatorListIfExist()
std::vector< G4VPhysicalVolume * >::iterator GetWorldsIterator()
G4double GetProductionCut(G4int index) const
G4VUserRegionInformation * GetUserInformation() const
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
static masterWorlds_t & GetMasterWorlds()
void SetReadiness(G4bool val=true)
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
G4Region * GetRegion() const
G4ParticleDefinition * GetGenericIon() const
G4ProcessManager * GetProcessManager() const
static G4RunManagerKernel * GetRunManagerKernel()
#define G4MUTEX_INITIALIZER
G4bool RunInitialization(G4bool fakeRun=false)
const G4String & GetParticleName() const
void WorkerDefineWorldVolume(G4VPhysicalVolume *worldVol, G4bool topologyIsChanged=true)
G4bool IsGeneralIon() const
static void ConstructParticle()
static G4UImanager * GetUIpointer()
static G4RegionStore * GetInstance()
static G4StateManager * GetStateManager()
void SetRegion(G4Region *reg)
G4VPhysicalVolume * GetWorldPhysical() const
G4ParticleDefinition * GetParticle(G4int index) const
G4IonTable * GetIonTable() const
void DumpCutValuesTable(G4int flag=1)
G4GLOB_DLL std::ostream G4cout
G4bool SetNewState(G4ApplicationState requestedState)
void PhysicsTableUpdated()
void reset(G4bool ifSkipIon=true)
const G4String & GetName() const
void DumpRegion(const G4String &rname) const
static G4MTRunManager * GetMasterRunManager()
void UsedInMassGeometry(G4bool val=true)
static G4LogicalVolumeStore * GetInstance()
G4ApplicationState GetCurrentState() const
RMKType runManagerKernelType
static G4GeometryManager * GetInstance()
const G4String & GetProcessName() const
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
void SetupDefaultRegion()
static G4TransportationManager * GetTransportationManager()
static void ClearUnitsTable()
static G4ProductionCutsTable * GetProductionCutsTable()
G4bool BelongsTo(G4VPhysicalVolume *thePhys) const
void UseCoupledTransportation(G4bool vl=true)
G4LogicalVolume * GetLogicalVolume() const
void SetWorld(G4VPhysicalVolume *wp)
static G4ParticleTable * GetParticleTable()
G4PrimaryTransformer * GetPrimaryTransformer() const
G4bool IsInMassGeometry() const
G4FastSimulationManager * GetFastSimulationManager() const
size_t GetNoWorlds() const
void Destroy(G4int nStat=0, G4int verboseLevel=0)
void SetProductionCuts(G4ProductionCuts *cut)
void OpenGeometry(G4VPhysicalVolume *vol=0)
void UpdateCoupleTable(G4VPhysicalVolume *currentWorld)
virtual void ConstructParticle()=0
G4bool RegisterWorld(G4VPhysicalVolume *aWorld)
size_t GetNumberOfMaterials() const
void SetWorldForTracking(G4VPhysicalVolume *theWorld)
void BuildPhysicsTables(G4bool fakeRun)
G4ProcessVector * GetAlongStepProcessVector(G4ProcessVectorTypeIndex typ=typeGPIL) const
size_t GetNumberOfRootVolumes() const
virtual ~G4RunManagerKernel()
void RemoveRootLogicalVolume(G4LogicalVolume *lv, G4bool scan=true)
void WorkerUpdateWorldVolume()
G4bool IsInParallelGeometry() const
G4PTblDicIterator * GetIterator() const
G4UserLimits * GetUserLimits() const
G4UserSteppingAction * GetRegionalSteppingAction() const
G4bool CloseGeometry(G4bool pOptimise=true, G4bool verbose=false, G4VPhysicalVolume *vol=0)
#define theParticleIterator
G4ProcessVector * GetProcessList() const
G4GLOB_DLL std::ostream G4cerr
void DefineWorldVolume(G4VPhysicalVolume *worldVol, G4bool topologyIsChanged=true)
G4String GetStateString(G4ApplicationState aState) const
static G4SDManager * GetSDMpointerIfExist()
std::vector< G4LogicalVolume * >::iterator GetRootLogicalVolumeIterator()
G4bool geometryNeedsToBeClosed