#include <G4WorkerRunManager.hh>

Inheritance diagram for G4WorkerRunManager:

Public Member Functions
	G4WorkerRunManager ()

	~G4WorkerRunManager () override

void	InitializeGeometry () override

void	RunInitialization () override

void	DoEventLoop (G4int n_event, const char *macroFile=nullptr, G4int n_select=-1) override

void	ProcessOneEvent (G4int i_event) override

G4Event *	GenerateEvent (G4int i_event) override

void	RunTermination () override

void	TerminateEventLoop () override

virtual void	DoWork ()

void	SetWorkerThread (G4WorkerThread *wc)

void	SetUserInitialization (G4VUserPhysicsList *userInit) override

void	SetUserInitialization (G4VUserDetectorConstruction *userInit) override

void	SetUserInitialization (G4VUserActionInitialization *userInit) override

void	SetUserInitialization (G4UserWorkerInitialization *userInit) override

void	SetUserInitialization (G4UserWorkerThreadInitialization *userInit) override

void	SetUserAction (G4UserRunAction *userAction) override

void	SetUserAction (G4VUserPrimaryGeneratorAction *userAction) override

void	SetUserAction (G4UserEventAction *userAction) override

void	SetUserAction (G4UserStackingAction *userAction) override

void	SetUserAction (G4UserTrackingAction *userAction) override

void	SetUserAction (G4UserSteppingAction *userAction) override

void	RestoreRndmEachEvent (G4bool flag) override

Public Member Functions inherited from G4RunManager
	G4RunManager ()

virtual	~G4RunManager ()

	G4RunManager (const G4RunManager &)=delete

G4RunManager &	operator= (const G4RunManager &)=delete

virtual void	BeamOn (G4int n_event, const char *macroFile=nullptr, G4int n_select=-1)

virtual void	Initialize ()

virtual void	DefineWorldVolume (G4VPhysicalVolume *worldVol, G4bool topologyIsChanged=true)

virtual void	AbortRun (G4bool softAbort=false)

virtual void	AbortEvent ()

virtual void	InitializePhysics ()

virtual G4bool	ConfirmBeamOnCondition ()

virtual void	InitializeEventLoop (G4int n_event, const char *macroFile=nullptr, G4int n_select=-1)

virtual void	TerminateOneEvent ()

virtual void	AnalyzeEvent (G4Event *anEvent)

virtual void	SetNumberOfThreads (G4int)

virtual G4int	GetNumberOfThreads () const

void	DumpRegion (const G4String &rname) const

void	DumpRegion (G4Region *region=nullptr) const

void	GeometryHasBeenModified (G4bool prop=true)

void	ReinitializeGeometry (G4bool destroyFirst=false, G4bool prop=true)

void	PhysicsHasBeenModified ()

void	CutOffHasBeenModified ()

void	ReOptimizeMotherOf (G4VPhysicalVolume *)

void	ReOptimize (G4LogicalVolume *)

void	SetGeometryToBeOptimized (G4bool vl)

G4bool	GetGeometryToBeOptimized ()

void	GeometryDirectlyUpdated (G4bool val=true)

virtual void	RestoreRandomNumberStatus (const G4String &fileN)

const G4VUserDetectorConstruction *	GetUserDetectorConstruction () const

const G4VUserPhysicsList *	GetUserPhysicsList () const

const G4VUserActionInitialization *	GetUserActionInitialization () const

G4VUserActionInitialization *	GetNonConstUserActionInitialization () const

const G4UserWorkerInitialization *	GetUserWorkerInitialization () const

const G4UserWorkerThreadInitialization *	GetUserWorkerThreadInitialization () const

const G4UserRunAction *	GetUserRunAction () const

const G4VUserPrimaryGeneratorAction *	GetUserPrimaryGeneratorAction () const

const G4UserEventAction *	GetUserEventAction () const

const G4UserStackingAction *	GetUserStackingAction () const

const G4UserTrackingAction *	GetUserTrackingAction () const

const G4UserSteppingAction *	GetUserSteppingAction () const

void	SetNumberOfAdditionalWaitingStacks (G4int iAdd)

void	SetDefaultClassification (G4TrackStatus ts, G4ClassificationOfNewTrack val, G4ExceptionSeverity es=G4ExceptionSeverity::IgnoreTheIssue)

void	SetDefaultClassification (const G4ParticleDefinition *pd, G4ClassificationOfNewTrack val, G4ExceptionSeverity es=G4ExceptionSeverity::IgnoreTheIssue)

const G4String &	GetVersionString () const

void	SetPrimaryTransformer (G4PrimaryTransformer *pt)

void	StoreRandomNumberStatusToG4Event (G4int vl)

G4int	GetFlagRandomNumberStatusToG4Event () const

void	SetRandomNumberStore (G4bool flag)

G4bool	GetRandomNumberStore () const

void	SetRandomNumberStoreDir (const G4String &dir)

const G4String &	GetRandomNumberStoreDir () const

const G4String &	GetRandomNumberStatusForThisRun () const

const G4String &	GetRandomNumberStatusForThisEvent () const

void	SetRandomNumberStorePerEvent (G4bool flag)

G4bool	GetRandomNumberStorePerEvent () const

void	SetVerboseLevel (G4int vl)

G4int	GetVerboseLevel () const

G4int	GetPrintProgress ()

void	SetPrintProgress (G4int i)

void	SetNumberOfEventsToBeStored (G4int val)

const G4Run *	GetCurrentRun () const

G4Run *	GetNonConstCurrentRun () const

const G4Event *	GetCurrentEvent () const

const G4Event *	GetPreviousEvent (G4int i) const

void	SetRunIDCounter (G4int i)

G4int	GetNumberOfParallelWorld () const

void	SetNumberOfEventsToBeProcessed (G4int val)

G4int	GetNumberOfEventsToBeProcessed () const

G4int	GetNumberOfSelectEvents () const

const G4String &	GetSelectMacro () const

void	SetDCtable (G4DCtable *DCtbl)

RMType	GetRunManagerType () const

virtual void	RegisterSubEventType (G4int, G4int)

virtual void	MergeTrajectories (const G4SubEvent , const G4Event )

virtual void	UpdateScoringForSubEvent (const G4SubEvent , const G4Event )

virtual const G4SubEvent *	GetSubEvent (G4int, G4bool &, G4long &, G4long &, G4long &, G4bool)

virtual void	SubEventFinished (const G4SubEvent , const G4Event )

virtual G4int	GetSubEventType () const

virtual void	SetSubEventType (G4int)

virtual std::size_t	GetMaxNTrack () const

virtual void	TrajectoriesToBeMerged (G4bool)

virtual void	ReportEventDeletion (const G4Event *evt)

void	ResetNavigatorAtInitialization (G4bool val=true)

Static Public Member Functions
static G4WorkerRunManager *	GetWorkerRunManager ()

static G4WorkerRunManagerKernel *	GetWorkerRunManagerKernel ()

Static Public Member Functions inherited from G4RunManager
static G4RunManager *	GetRunManager ()

static G4bool	IfGeometryHasBeenDestroyed ()

Protected Member Functions
void	ConstructScoringWorlds () override

void	StoreRNGStatus (const G4String &filenamePrefix) override

void	rndmSaveThisRun () override

void	rndmSaveThisEvent () override

virtual void	MergePartialResults (G4bool mergeEvents=true)

virtual void	SetupDefaultRNGEngine ()

Protected Member Functions inherited from G4RunManager
	G4RunManager (RMType rmType)

virtual void	StackPreviousEvent (G4Event *anEvent)

virtual void	CleanUpPreviousEvents ()

virtual void	CleanUpUnnecessaryEvents (G4int keepNEvents)

void	UpdateScoring (const G4Event *evt=nullptr)

virtual void	DeleteUserInitializations ()

Protected Attributes
G4WorkerThread *	workerContext = nullptr

G4bool	eventLoopOnGoing = false

G4bool	runIsSeeded = false

G4int	nevModulo = -1

G4int	currEvID = -1

G4int	luxury = -1

G4SeedsQueue	seedsQueue

G4bool	readStatusFromFile = false

Protected Attributes inherited from G4RunManager
G4RunManagerKernel *	kernel = nullptr

G4EventManager *	eventManager = nullptr

G4VUserDetectorConstruction *	userDetector = nullptr

G4VUserPhysicsList *	physicsList = nullptr

G4VUserActionInitialization *	userActionInitialization = nullptr

G4UserWorkerInitialization *	userWorkerInitialization = nullptr

G4UserWorkerThreadInitialization *	userWorkerThreadInitialization = nullptr

G4UserRunAction *	userRunAction = nullptr

G4VUserPrimaryGeneratorAction *	userPrimaryGeneratorAction = nullptr

G4UserEventAction *	userEventAction = nullptr

G4UserStackingAction *	userStackingAction = nullptr

G4UserTrackingAction *	userTrackingAction = nullptr

G4UserSteppingAction *	userSteppingAction = nullptr

G4bool	geometryInitialized = false

G4bool	physicsInitialized = false

G4bool	runAborted = false

G4bool	initializedAtLeastOnce = false

G4bool	geometryToBeOptimized = true

G4int	runIDCounter = 0

G4int	verboseLevel = 0

G4int	printModulo = -1

G4Timer *	timer = nullptr

G4DCtable *	DCtable = nullptr

G4Run *	currentRun = nullptr

G4Event *	currentEvent = nullptr

std::list< G4Event * > *	previousEvents = nullptr

G4int	n_perviousEventsToBeStored = 0

G4int	numberOfEventToBeProcessed = 0

G4bool	storeRandomNumberStatus = false

G4int	storeRandomNumberStatusToG4Event = 0

G4String	randomNumberStatusDir = "./"

G4String	randomNumberStatusForThisRun = ""

G4String	randomNumberStatusForThisEvent = ""

G4bool	rngStatusEventsFlag = false

G4VPhysicalVolume *	currentWorld = nullptr

G4int	nParallelWorlds = 0

G4String	msgText = " "

G4int	n_select_msg = -1

G4int	numberOfEventProcessed = 0

G4String	selectMacro = ""

G4bool	fakeRun = false

G4bool	isScoreNtupleWriter = false

G4bool	geometryDirectlyUpdated = false

RMType	runManagerType

Additional Inherited Members
Public Types inherited from G4RunManager
enum	RMType { sequentialRM , masterRM , workerRM , subEventMasterRM , subEventWorkerRM }

Static Protected Attributes inherited from G4RunManager
static G4RUN_DLL G4bool	fGeometryHasBeenDestroyed = false

Detailed Description

Definition at line 48 of file G4WorkerRunManager.hh.

Constructor & Destructor Documentation

◆ G4WorkerRunManager()

G4WorkerRunManager::G4WorkerRunManager ( )

Definition at line 76 of file G4WorkerRunManager.cc.

                                       : G4RunManager(workerRM)
{
  // This constructor should never be called in non-multithreaded mode
 
#ifndef G4MULTITHREADED
  G4ExceptionDescription msg;
  msg << "Geant4 code is compiled without multi-threading support "
         "(-DG4MULTITHREADED "
         "is set to off).";
  msg << " This type of RunManager can only be used in mult-threaded "
         "applications.";
  G4Exception("G4WorkerRunManager::G4WorkerRunManager()", "Run0103", FatalException, msg);
#endif
  // G4ParticleTable::GetParticleTable()->WorkerG4ParticleTable();
     // WorkerG4ParticleTable() would be performed twice,
     // as it is called already from G4ParticleTable::GetParticleTable()
     // which is called beforehand by other Geant4 classes
  G4ScoringManager* masterScM = G4MTRunManager::GetMasterScoringManager();
  if (masterScM != nullptr) G4ScoringManager::GetScoringManager();  // TLS instance for a worker
 
  // Properly initialise luxury level for Ranlux* engines...
  //
  if (dynamic_cast<const CLHEP::Ranlux64Engine*>(G4Random::getTheEngine()) != nullptr) {
    const auto theEngine = dynamic_cast<const CLHEP::Ranlux64Engine*>(G4Random::getTheEngine());
    luxury = theEngine->getLuxury();
  }
  else if (dynamic_cast<const CLHEP::RanluxEngine*>(G4Random::getTheEngine()) != nullptr) {
    const auto theEngine = dynamic_cast<const CLHEP::RanluxEngine*>(G4Random::getTheEngine());
    luxury = theEngine->getLuxury();
  }
 
  G4UImanager::GetUIpointer()->SetIgnoreCmdNotFound(true);
 
#ifdef G4MULTITHREADED
  G4VVisManager* pVVis = G4VVisManager::GetConcreteInstance();
  if (pVVis != nullptr) {
    pVVis->SetUpForAThread();
    visIsSetUp = true;
  }
#endif
}

Referenced by GetWorkerRunManager().

◆ ~G4WorkerRunManager()

G4WorkerRunManager::~G4WorkerRunManager ( )

override

Definition at line 119 of file G4WorkerRunManager.cc.

{
  CleanUpPreviousEvents();
  // Put these pointers to zero: owned by master thread
  // If not to zero, the base class destructor will attempt to
  // delete them
  userDetector = nullptr;
  userWorkerInitialization = nullptr;
  userWorkerThreadInitialization = nullptr;
  userActionInitialization = nullptr;
  physicsList->TerminateWorker();
  physicsList = nullptr;
  if (verboseLevel > 1) G4cout << "Destroying WorkerRunManager (" << this << ")" << G4endl;
}

Member Function Documentation

◆ ConstructScoringWorlds()

void G4WorkerRunManager::ConstructScoringWorlds ( )

overrideprotectedvirtual

Reimplemented from G4RunManager.

Definition at line 482 of file G4WorkerRunManager.cc.

{
  using MeshShape = G4VScoringMesh::MeshShape;
 
  // Return if unnecessary
  G4ScoringManager* ScM = G4ScoringManager::GetScoringManagerIfExist();
  if (ScM == nullptr) return;
  auto nPar = (G4int)ScM->GetNumberOfMesh();
  if (nPar < 1) return;
 
  // Update thread-local G4TransportationManager of all the world volumes
  kernel->WorkerUpdateWorldVolume();
 
  G4ScoringManager* masterScM = G4MTRunManager::GetMasterScoringManager();
 
  auto particleIterator = G4ParticleTable::GetParticleTable()->GetIterator();
 
  for (G4int iw = 0; iw < nPar; ++iw) {
    G4VScoringMesh* mesh = ScM->GetMesh(iw);
    if (fGeometryHasBeenDestroyed) mesh->GeometryHasBeenDestroyed();
    G4VPhysicalVolume* pWorld = nullptr;
    if (mesh->GetShape() != MeshShape::realWorldLogVol) {
      pWorld =
        G4TransportationManager::GetTransportationManager()->IsWorldExisting(ScM->GetWorldName(iw));
      if (pWorld == nullptr) {
        G4ExceptionDescription ed;
        ed << "Mesh name <" << ScM->GetWorldName(iw) << "> is not found in the master thread.";
        G4Exception("G4WorkerRunManager::ConstructScoringWorlds()", "RUN79001", FatalException, ed);
      }
    }
    if ((mesh->GetMeshElementLogical()) == nullptr) {
      G4AutoLock l(&ConstructScoringWorldsMutex);
      G4VScoringMesh* masterMesh = masterScM->GetMesh(iw);
      mesh->SetMeshElementLogical(masterMesh->GetMeshElementLogical());
      l.unlock();
 
      if (mesh->GetShape() != MeshShape::realWorldLogVol) {
        G4ParallelWorldProcess* theParallelWorldProcess = mesh->GetParallelWorldProcess();
        if (theParallelWorldProcess != nullptr) {
          theParallelWorldProcess->SetParallelWorld(ScM->GetWorldName(iw));
        }
        else {
          theParallelWorldProcess = new G4ParallelWorldProcess(ScM->GetWorldName(iw));
          mesh->SetParallelWorldProcess(theParallelWorldProcess);
          theParallelWorldProcess->SetParallelWorld(ScM->GetWorldName(iw));
 
          particleIterator->reset();
          while ((*particleIterator)()) {
            G4ParticleDefinition* particle = particleIterator->value();
            G4ProcessManager* pmanager = particle->GetProcessManager();
            if (pmanager != nullptr) {
              pmanager->AddProcess(theParallelWorldProcess);
              if (theParallelWorldProcess->IsAtRestRequired(particle)) {
                pmanager->SetProcessOrdering(theParallelWorldProcess, idxAtRest, 9900);
              }
              pmanager->SetProcessOrderingToSecond(theParallelWorldProcess, idxAlongStep);
              pmanager->SetProcessOrdering(theParallelWorldProcess, idxPostStep, 9900);
            }  // if(pmanager)
          }  // while
        }
        theParallelWorldProcess->SetLayeredMaterialFlag(mesh->LayeredMassFlg());
      }
    }
    mesh->WorkerConstruct(pWorld);
  }
}

Referenced by G4WorkerSubEvtRunManager::DoWork(), and G4WorkerTaskRunManager::DoWork().

◆ DoEventLoop()

void G4WorkerRunManager::DoEventLoop	(	G4int	n_event,
		const char *	macroFile = nullptr,
		G4int	n_select = -1 )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager, and G4WorkerTaskRunManager.

Definition at line 252 of file G4WorkerRunManager.cc.

{
  if (userPrimaryGeneratorAction == nullptr) {
    G4Exception("G4RunManager::GenerateEvent()", "Run0032", FatalException,
                "G4VUserPrimaryGeneratorAction is not defined!");
  }
 
  // This is the same as in the sequential case, just the for-loop indexes are
  // different
  InitializeEventLoop(n_event, macroFile, n_select);
 
  // Reset random number seeds queue
  while (!seedsQueue.empty()) {
    seedsQueue.pop();
  }
  // for each run, worker should receive at least one set of random number
  // seeds.
  runIsSeeded = false;
 
  // Event loop
  eventLoopOnGoing = true;
  ///////    G4int i_event = workerContext->GetThreadId();
  G4int i_event = -1;
  nevModulo = -1;
  currEvID = -1;
 
  while (eventLoopOnGoing) {
    ProcessOneEvent(i_event);
    if (eventLoopOnGoing) {
      TerminateOneEvent();
      if (runAborted) {
        eventLoopOnGoing = false;
      }
    }
  }
 
  TerminateEventLoop();
}

◆ DoWork()

void G4WorkerRunManager::DoWork ( )

virtual

Reimplemented in G4WorkerSubEvtRunManager, and G4WorkerTaskRunManager.

Definition at line 706 of file G4WorkerRunManager.cc.

{
  G4MTRunManager* mrm = G4MTRunManager::GetMasterRunManager();
  G4MTRunManager::WorkerActionRequest nextAction = mrm->ThisWorkerWaitForNextAction();
  while (nextAction != G4MTRunManager::WorkerActionRequest::ENDWORKER) {
    if (nextAction == G4MTRunManager::WorkerActionRequest::NEXTITERATION)  // start the next
                                                                           // run
    {
      // The following code deals with changing materials between runs
      static G4ThreadLocal G4bool skipInitialization = true;
      if (skipInitialization) {
        // re-initialization is not necessary for the first run
        skipInitialization = false;
      }
      else {
        //        ReinitializeGeometry();
        workerContext->UpdateGeometryAndPhysicsVectorFromMaster();
      }
 
      // Execute UI commands stored in the master UI manager
      std::vector<G4String> cmds = mrm->GetCommandStack();
      G4UImanager* uimgr = G4UImanager::GetUIpointer();  // TLS instance
      for (const auto& cmd : cmds) {
        uimgr->ApplyCommand(cmd);
      }
      // Start this run
      G4int numevents = mrm->GetNumberOfEventsToBeProcessed();
      G4String macroFile = mrm->GetSelectMacro();
      G4int numSelect = mrm->GetNumberOfSelectEvents();
      if (macroFile.empty() || macroFile == " ") {
        this->BeamOn(numevents);
      }
      else {
        this->BeamOn(numevents, macroFile, numSelect);
      }
    }
    else if (nextAction == G4MTRunManager::WorkerActionRequest::PROCESSUI) {
      std::vector<G4String> cmds = mrm->GetCommandStack();
      G4UImanager* uimgr = G4UImanager::GetUIpointer();  // TLS instance
      for (const auto& cmd : cmds) {
        uimgr->ApplyCommand(cmd);
      }
      mrm->ThisWorkerProcessCommandsStackDone();
    }
    else {
      G4ExceptionDescription d;
      d << "Cannot continue, this worker has been requested an unknown action: "
        << static_cast<std::underlying_type<G4MTRunManager::WorkerActionRequest>::type>(nextAction);
      G4Exception("G4WorkerRunManager::DoWork", "Run0104", FatalException, d);
    }
 
    // Now wait for master thread to signal new action to be performed
    nextAction = mrm->ThisWorkerWaitForNextAction();
  }  // No more actions to perform
 
  return;
}

Referenced by G4MTRunManagerKernel::StartThread().

◆ GenerateEvent()

G4Event * G4WorkerRunManager::GenerateEvent ( G4int i_event )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager, and G4WorkerTaskRunManager.

Definition at line 305 of file G4WorkerRunManager.cc.

{
  auto anEvent = new G4Event(i_event);
  G4long s1 = 0;
  G4long s2 = 0;
  G4long s3 = 0;
  G4bool eventHasToBeSeeded = true;
  if (G4MTRunManager::SeedOncePerCommunication() == 1 && runIsSeeded) {
    eventHasToBeSeeded = false;
  }
 
  if (i_event < 0) {
    G4int nevM = G4MTRunManager::GetMasterRunManager()->GetEventModulo();
    if (nevM == 1) {
      eventLoopOnGoing = G4MTRunManager::GetMasterRunManager()->SetUpAnEvent(anEvent, s1, s2, s3,
                                                                             eventHasToBeSeeded);
      runIsSeeded = true;
    }
    else {
      if (nevModulo <= 0) {
        G4int nevToDo = G4MTRunManager::GetMasterRunManager()->SetUpNEvents(anEvent, &seedsQueue,
                                                                            eventHasToBeSeeded);
        if (nevToDo == 0) {
          eventLoopOnGoing = false;
        }
        else {
          currEvID = anEvent->GetEventID();
          nevModulo = nevToDo - 1;
        }
      }
      else {
        if (G4MTRunManager::SeedOncePerCommunication() > 0) eventHasToBeSeeded = false;
        anEvent->SetEventID(++currEvID);
        --nevModulo;
      }
      if (eventLoopOnGoing && eventHasToBeSeeded) {
        s1 = seedsQueue.front();
        seedsQueue.pop();
        s2 = seedsQueue.front();
        seedsQueue.pop();
      }
    }
 
    if (!eventLoopOnGoing) {
      delete anEvent;
      return nullptr;
    }
  }
  else if (eventHasToBeSeeded) {
    // Need to reseed random number generator
    G4RNGHelper* helper = G4RNGHelper::GetInstance();
    s1 = helper->GetSeed(i_event * 2);
    s2 = helper->GetSeed(i_event * 2 + 1);
  }
 
  if (eventHasToBeSeeded) {
    G4long seeds[3] = {s1, s2, 0};
    G4Random::setTheSeeds(seeds, luxury);
    runIsSeeded = true;
  }
 
  // Read from file seed.
  // Andrea Dotti 4 November 2015
  // This is required for strong-reproducibility, in MT mode we have that each
  // thread produces, for each event a status file, we want to do that.
  // Search a random file with the format run{%d}evt{%d}.rndm
 
  // This is the filename base constructed from run and event
  const auto filename = [&] {
    std::ostringstream os;
    os << "run" << currentRun->GetRunID() << "evt" << anEvent->GetEventID();
    return os.str();
  };
 
  G4bool RNGstatusReadFromFile = false;
  if (readStatusFromFile) {
    // Build full path of RNG status file for this event
    std::ostringstream os;
    os << filename() << ".rndm";
    const G4String& randomStatusFile = os.str();
    std::ifstream ifile(randomStatusFile.c_str());
    if (ifile) {  // File valid and readable
      RNGstatusReadFromFile = true;
      G4Random::restoreEngineStatus(randomStatusFile.c_str());
    }
  }
 
  if (storeRandomNumberStatusToG4Event == 1 || storeRandomNumberStatusToG4Event == 3) {
    std::ostringstream oss;
    G4Random::saveFullState(oss);
    randomNumberStatusForThisEvent = oss.str();
    anEvent->SetRandomNumberStatus(randomNumberStatusForThisEvent);
  }
 
  if (storeRandomNumberStatus && !RNGstatusReadFromFile)
  {  // If reading from file, avoid to rewrite the same
    G4String fileN = "currentEvent";
    if (rngStatusEventsFlag) {
      fileN = filename();
    }
    StoreRNGStatus(fileN);
  }
 
  if (printModulo > 0 && anEvent->GetEventID() % printModulo == 0) {
    G4cout << "--> Event " << anEvent->GetEventID() << " starts";
    if (eventHasToBeSeeded) {
      G4cout << " with initial seeds (" << s1 << "," << s2 << ")";
    }
    G4cout << "." << G4endl;
  }
  userPrimaryGeneratorAction->GeneratePrimaries(anEvent);
  return anEvent;
}

Referenced by ProcessOneEvent().

◆ GetWorkerRunManager()

G4WorkerRunManager * G4WorkerRunManager::GetWorkerRunManager ( )

static

Definition at line 64 of file G4WorkerRunManager.cc.

{
  return static_cast<G4WorkerRunManager*>(G4RunManager::GetRunManager());
}

Referenced by GetWorkerRunManagerKernel(), G4RTWorkerInitialization::WorkerRunEnd(), and G4RTWorkerInitialization::WorkerRunStart().

◆ GetWorkerRunManagerKernel()

G4WorkerRunManagerKernel * G4WorkerRunManager::GetWorkerRunManagerKernel ( )

static

Definition at line 70 of file G4WorkerRunManager.cc.

{
  return static_cast<G4WorkerRunManagerKernel*>(GetWorkerRunManager()->kernel);
}

◆ InitializeGeometry()

void G4WorkerRunManager::InitializeGeometry ( )

overridevirtual

Reimplemented from G4RunManager.

Definition at line 135 of file G4WorkerRunManager.cc.

{
  if (userDetector == nullptr) {
    G4Exception("G4RunManager::InitializeGeometry", "Run0033", FatalException,
                "G4VUserDetectorConstruction is not defined!");
    return;
  }
 
  if (fGeometryHasBeenDestroyed) {
    G4TransportationManager::GetTransportationManager()->ClearParallelWorlds();
  }
 
  // Step 0: Contribute to the voxelisation of the geometry
  G4GeometryManager* geomManager = G4GeometryManager::GetInstance();
  if( geomManager->IsParallelOptimisationConfigured() ) {
    G4cout << "G4RunManager::InitializeGeometry calling GeometryManager's UndertakeOptimisation" 
           << G4endl;  // TODO - suppress / delete this in final version
    geomManager->UndertakeOptimisation();
  }
  //  A barrier must ensure that all that all threads have finished this work.
  //  Currently we rely on the (later) barrier at the end of initialisation.
  
  // Step1: Get pointer to the physiWorld (note: needs to get the "super
  // pointer, i.e. the one shared by all threads"
  G4RunManagerKernel* masterKernel = G4MTRunManager::GetMasterRunManagerKernel();
  G4VPhysicalVolume* worldVol = masterKernel->GetCurrentWorld();
  // Step2:, Call a new "WorkerDefineWorldVolume( pointer from 2-, false);
  kernel->WorkerDefineWorldVolume(worldVol, false);
  kernel->SetNumberOfParallelWorld(masterKernel->GetNumberOfParallelWorld());
  // Step3: Call user's ConstructSDandField()
  userDetector->ConstructSDandField();
  userDetector->ConstructParallelSD();
  geometryInitialized = true;
}

◆ MergePartialResults()

void G4WorkerRunManager::MergePartialResults ( G4bool mergeEvents = true )

protectedvirtual

Definition at line 420 of file G4WorkerRunManager.cc.

{
  // Merge partial results into global run
  G4MTRunManager* mtRM = G4MTRunManager::GetMasterRunManager();
  G4ScoringManager* ScM = G4ScoringManager::GetScoringManagerIfExist();
  if (ScM != nullptr) mtRM->MergeScores(ScM);
#ifdef G4VERBOSE
  if(mergeEvents && verboseLevel>3) {
    auto eventVector = currentRun->GetEventVector();
    if(eventVector!=nullptr || !(eventVector->empty())) {
      G4cout<<"G4WorkerRunManager::MergePartialResults : merging "
            <<eventVector->size()<<" events."<<G4endl;
    }
  }
#endif
  if(mergeEvents) mtRM->MergeRun(currentRun);
}

Referenced by RunTermination(), G4WorkerSubEvtRunManager::RunTermination(), and G4WorkerTaskRunManager::RunTermination().

◆ ProcessOneEvent()

void G4WorkerRunManager::ProcessOneEvent ( G4int i_event )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager, and G4WorkerTaskRunManager.

Definition at line 292 of file G4WorkerRunManager.cc.

{
  currentEvent = GenerateEvent(i_event);
  if (eventLoopOnGoing) {
    eventManager->ProcessOneEvent(currentEvent);
    AnalyzeEvent(currentEvent);
    UpdateScoring();
    if (currentEvent->GetEventID() < n_select_msg)
      G4UImanager::GetUIpointer()->ApplyCommand(msgText);
  }
}

Referenced by DoEventLoop().

◆ RestoreRndmEachEvent()

void G4WorkerRunManager::RestoreRndmEachEvent ( G4bool flag )

inlineoverridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager, and G4WorkerTaskRunManager.

Definition at line 85 of file G4WorkerRunManager.hh.

85{ readStatusFromFile = flag; }

◆ rndmSaveThisEvent()

void G4WorkerRunManager::rndmSaveThisEvent ( )

overrideprotectedvirtual

Reimplemented from G4RunManager.

Definition at line 672 of file G4WorkerRunManager.cc.

{
  if (currentEvent == nullptr) {
    G4cerr << "Warning from G4RunManager::rndmSaveThisEvent():"
           << " there is no currentEvent available." << G4endl << "Command ignored." << G4endl;
    return;
  }
 
  if (!storeRandomNumberStatus) {
    G4cerr << "Warning from G4RunManager::rndmSaveThisEvent():"
           << " Random number engine status is not available." << G4endl
           << "/random/setSavingFlag command must be issued "
           << "prior to the start of the run. Command ignored." << G4endl;
    return;
  }
 
  std::ostringstream oos;
  oos << "G4Worker" << workerContext->GetThreadId() << "_"
      << "currentEvent.rndm"
      << "\0";
  G4fs::path fileIn = randomNumberStatusDir + oos.str();
 
  std::ostringstream os;
  os << "run" << currentRun->GetRunID() << "evt" << currentEvent->GetEventID() << ".rndm" << '\0';
  G4fs::path fileOut = randomNumberStatusDir + os.str();
 
  if (G4CopyRandomState(fileIn, fileOut, "G4WorkerRunManager::rndmSaveThisEvent()")
      && verboseLevel > 0)
  {
    G4cout << fileIn << " is copied to " << fileOut << G4endl;
  }
}

◆ rndmSaveThisRun()

void G4WorkerRunManager::rndmSaveThisRun ( )

overrideprotectedvirtual

Reimplemented from G4RunManager.

Definition at line 642 of file G4WorkerRunManager.cc.

{
  G4int runNumber = 0;
  if (currentRun != nullptr) runNumber = currentRun->GetRunID();
  if (!storeRandomNumberStatus) {
    G4cerr << "Warning from G4RunManager::rndmSaveThisRun():"
           << " Random number status was not stored prior to this run." << G4endl
           << "/random/setSavingFlag command must be issued. "
           << "Command ignored." << G4endl;
    return;
  }
 
  std::ostringstream oos;
  oos << "G4Worker" << workerContext->GetThreadId() << "_"
      << "currentRun.rndm"
      << "\0";
  G4fs::path fileIn = randomNumberStatusDir + oos.str();
 
  std::ostringstream os;
  os << "run" << runNumber << ".rndm" << '\0';
  G4fs::path fileOut = randomNumberStatusDir + os.str();
 
  if (G4CopyRandomState(fileIn, fileOut, "G4WorkerRunManager::rndmSaveThisRun()")
      && verboseLevel > 0)
  {
    G4cout << fileIn << " is copied to " << fileOut << G4endl;
  }
}

◆ RunInitialization()

void G4WorkerRunManager::RunInitialization ( )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager, and G4WorkerTaskRunManager.

Definition at line 171 of file G4WorkerRunManager.cc.

{
#ifdef G4MULTITHREADED
  if (!visIsSetUp) {
    G4VVisManager* pVVis = G4VVisManager::GetConcreteInstance();
    if (pVVis != nullptr) {
      pVVis->SetUpForAThread();
      visIsSetUp = true;
    }
  }
#endif
 
  if (!(kernel->RunInitialization(fakeRun))) return;
 
  // Signal this thread can start event loop.
  // Note this will return only when all threads reach this point
  G4MTRunManager::GetMasterRunManager()->ThisWorkerReady();
  if (fakeRun) return;
 
  const G4UserWorkerInitialization* uwi =
    G4MTRunManager::GetMasterRunManager()->GetUserWorkerInitialization();
  CleanUpPreviousEvents();
  delete currentRun;
  currentRun = nullptr;
 
  if (fGeometryHasBeenDestroyed) G4ParallelWorldProcessStore::GetInstance()->UpdateWorlds();
  // Call a user hook: this is guaranteed all threads are "synchronized"
  if (uwi != nullptr) uwi->WorkerRunStart();
 
  if (userRunAction != nullptr) currentRun = userRunAction->GenerateRun();
  if (currentRun == nullptr) currentRun = new G4Run();
 
  currentRun->SetRunID(runIDCounter);
  currentRun->SetNumberOfEventToBeProcessed(numberOfEventToBeProcessed);
 
  currentRun->SetDCtable(DCtable);
  G4SDManager* fSDM = G4SDManager::GetSDMpointerIfExist();
  if (fSDM != nullptr) {
    currentRun->SetHCtable(fSDM->GetHCtable());
  }
 
  if (G4VScoreNtupleWriter::Instance() != nullptr) {
    auto hce = (fSDM != nullptr) ? fSDM->PrepareNewEvent() : nullptr;
    isScoreNtupleWriter = G4VScoreNtupleWriter::Instance()->Book(hce);
    delete hce;
  }
 
  std::ostringstream oss;
  G4Random::saveFullState(oss);
  randomNumberStatusForThisRun = oss.str();
  currentRun->SetRandomNumberStatus(randomNumberStatusForThisRun);
 
  for (G4int i_prev = 0; i_prev < n_perviousEventsToBeStored; ++i_prev) {
    previousEvents->push_back(nullptr);
  }
 
  if (printModulo > 0 || verboseLevel > 0) {
    G4cout << "### Run " << currentRun->GetRunID() << " starts on worker thread "
           << G4Threading::G4GetThreadId() << "." << G4endl;
  }
  if (userRunAction != nullptr) userRunAction->BeginOfRunAction(currentRun);
 
  if (isScoreNtupleWriter) {
    G4VScoreNtupleWriter::Instance()->OpenFile();
  }
 
  if (storeRandomNumberStatus) {
    G4String fileN = "currentRun";
    if (rngStatusEventsFlag) {
      std::ostringstream os;
      os << "run" << currentRun->GetRunID();
      fileN = os.str();
    }
    StoreRNGStatus(fileN);
  }
 
  runAborted = false;
  numberOfEventProcessed = 0;
}

◆ RunTermination()

void G4WorkerRunManager::RunTermination ( )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager, and G4WorkerTaskRunManager.

Definition at line 439 of file G4WorkerRunManager.cc.

{
  if (!fakeRun) {
    MergePartialResults();
 
    // Call a user hook: note this is before the next barrier
    // so threads execute this method asyncrhonouzly
    //(TerminateRun allows for synch via G4RunAction::EndOfRun)
    const G4UserWorkerInitialization* uwi =
      G4MTRunManager::GetMasterRunManager()->GetUserWorkerInitialization();
    if (uwi != nullptr) uwi->WorkerRunEnd();
  }
 
  G4RunManager::RunTermination();
  // Signal this thread has finished envent-loop.
  // Note this will return only whan all threads reach this point
  G4MTRunManager::GetMasterRunManager()->ThisWorkerEndEventLoop();
}

◆ SetupDefaultRNGEngine()

void G4WorkerRunManager::SetupDefaultRNGEngine ( )

protectedvirtual

Reimplemented in G4WorkerSubEvtRunManager, and G4WorkerTaskRunManager.

Definition at line 592 of file G4WorkerRunManager.cc.

{
  const CLHEP::HepRandomEngine* mrnge =
    G4MTRunManager::GetMasterRunManager()->getMasterRandomEngine();
  const G4UserWorkerThreadInitialization* uwti =
    G4MTRunManager::GetMasterRunManager()->GetUserWorkerThreadInitialization();
  uwti->SetupRNGEngine(mrnge);
}

◆ SetUserAction() [1/6]

void G4WorkerRunManager::SetUserAction ( G4UserEventAction * userAction )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager.

Definition at line 604 of file G4WorkerRunManager.cc.

{
  G4RunManager::SetUserAction(ua);
}

◆ SetUserAction() [2/6]

void G4WorkerRunManager::SetUserAction ( G4UserRunAction * userAction )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager.

Definition at line 585 of file G4WorkerRunManager.cc.

{
  G4RunManager::SetUserAction(userAction);
  if (userAction != nullptr) userAction->SetMaster(false);
}

Referenced by G4RTWorkerInitialization::WorkerRunEnd(), and G4RTWorkerInitialization::WorkerRunStart().

◆ SetUserAction() [3/6]

void G4WorkerRunManager::SetUserAction ( G4UserStackingAction * userAction )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager.

Definition at line 616 of file G4WorkerRunManager.cc.

{
  G4RunManager::SetUserAction(ua);
}

◆ SetUserAction() [4/6]

void G4WorkerRunManager::SetUserAction ( G4UserSteppingAction * userAction )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager.

Definition at line 628 of file G4WorkerRunManager.cc.

{
  G4RunManager::SetUserAction(ua);
}

◆ SetUserAction() [5/6]

void G4WorkerRunManager::SetUserAction ( G4UserTrackingAction * userAction )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager.

Definition at line 622 of file G4WorkerRunManager.cc.

{
  G4RunManager::SetUserAction(ua);
}

◆ SetUserAction() [6/6]

void G4WorkerRunManager::SetUserAction ( G4VUserPrimaryGeneratorAction * userAction )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager.

Definition at line 610 of file G4WorkerRunManager.cc.

{
  G4RunManager::SetUserAction(ua);
}

◆ SetUserInitialization() [1/5]

void G4WorkerRunManager::SetUserInitialization ( G4UserWorkerInitialization * userInit )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager.

Definition at line 550 of file G4WorkerRunManager.cc.

{
  G4Exception("G4RunManager::SetUserInitialization(G4UserWorkerInitialization*)", "Run0118",
              FatalException, "This method should be used only with an instance of G4MTRunManager");
}

◆ SetUserInitialization() [2/5]

void G4WorkerRunManager::SetUserInitialization ( G4UserWorkerThreadInitialization * userInit )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager.

Definition at line 557 of file G4WorkerRunManager.cc.

{
  G4Exception("G4RunManager::SetUserInitialization(G4UserWorkerThreadInitialization*)", "Run0119",
              FatalException, "This method should be used only with an instance of G4MTRunManager");
}

◆ SetUserInitialization() [3/5]

void G4WorkerRunManager::SetUserInitialization ( G4VUserActionInitialization * userInit )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager.

Definition at line 564 of file G4WorkerRunManager.cc.

{
  G4Exception("G4RunManager::SetUserInitialization(G4VUserActionInitialization*)", "Run0120",
              FatalException, "This method should be used only with an instance of G4MTRunManager");
}

◆ SetUserInitialization() [4/5]

void G4WorkerRunManager::SetUserInitialization ( G4VUserDetectorConstruction * userInit )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager.

Definition at line 571 of file G4WorkerRunManager.cc.

{
  G4Exception("G4RunManager::SetUserInitialization(G4VUserDetectorConstruction*)", "Run0121",
              FatalException, "This method should be used only with an instance of G4MTRunManager");
}

◆ SetUserInitialization() [5/5]

void G4WorkerRunManager::SetUserInitialization ( G4VUserPhysicsList * userInit )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager.

Definition at line 578 of file G4WorkerRunManager.cc.

{
  pl->InitializeWorker();
  G4RunManager::SetUserInitialization(pl);
}

Referenced by G4MTRunManagerKernel::StartThread().

◆ SetWorkerThread()

void G4WorkerRunManager::SetWorkerThread ( G4WorkerThread * wc )

inline

Definition at line 71 of file G4WorkerRunManager.hh.

71{ workerContext = wc; }

Referenced by G4MTRunManagerKernel::StartThread().

◆ StoreRNGStatus()

void G4WorkerRunManager::StoreRNGStatus ( const G4String & filenamePrefix )

overrideprotectedvirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager, and G4WorkerTaskRunManager.

Definition at line 634 of file G4WorkerRunManager.cc.

{
  std::ostringstream os;
  os << randomNumberStatusDir << "G4Worker" << workerContext->GetThreadId() << "_" << fn << ".rndm";
  G4Random::saveEngineStatus(os.str().c_str());
}

Referenced by GenerateEvent(), and RunInitialization().

◆ TerminateEventLoop()

void G4WorkerRunManager::TerminateEventLoop ( )

overridevirtual

Reimplemented from G4RunManager.

Reimplemented in G4WorkerSubEvtRunManager, and G4WorkerTaskRunManager.

Definition at line 459 of file G4WorkerRunManager.cc.

{
  if (verboseLevel > 0 && !fakeRun) {
    timer->Stop();
    G4cout << "Thread-local run terminated." << G4endl;
    G4cout << "Run Summary" << G4endl;
    if (runAborted) {
      G4cout << "  Run Aborted after " << numberOfEventProcessed << " events processed." << G4endl;
    }
    else {
      G4cout << "  Number of events processed : " << numberOfEventProcessed << G4endl;
    }
    G4cout << "  " << *timer << G4endl;
  }
}

Referenced by DoEventLoop().

Member Data Documentation

◆ currEvID

G4int G4WorkerRunManager::currEvID = -1

protected

Definition at line 106 of file G4WorkerRunManager.hh.

Referenced by DoEventLoop(), G4WorkerSubEvtRunManager::DoEventLoop(), G4WorkerTaskRunManager::DoEventLoop(), GenerateEvent(), G4WorkerSubEvtRunManager::GenerateEvent(), and G4WorkerTaskRunManager::GenerateEvent().

◆ eventLoopOnGoing

G4bool G4WorkerRunManager::eventLoopOnGoing = false

protected

Definition at line 103 of file G4WorkerRunManager.hh.

Referenced by DoEventLoop(), G4WorkerSubEvtRunManager::DoEventLoop(), G4WorkerTaskRunManager::DoEventLoop(), GenerateEvent(), G4WorkerSubEvtRunManager::GenerateEvent(), G4WorkerTaskRunManager::GenerateEvent(), ProcessOneEvent(), G4WorkerSubEvtRunManager::ProcessOneEvent(), and G4WorkerTaskRunManager::ProcessOneEvent().

◆ luxury

G4int G4WorkerRunManager::luxury = -1

protected

Definition at line 107 of file G4WorkerRunManager.hh.

Referenced by G4WorkerRunManager(), and GenerateEvent().

◆ nevModulo

G4int G4WorkerRunManager::nevModulo = -1

protected

Definition at line 105 of file G4WorkerRunManager.hh.

Referenced by DoEventLoop(), G4WorkerSubEvtRunManager::DoEventLoop(), G4WorkerTaskRunManager::DoEventLoop(), GenerateEvent(), G4WorkerSubEvtRunManager::GenerateEvent(), and G4WorkerTaskRunManager::GenerateEvent().

◆ readStatusFromFile

G4bool G4WorkerRunManager::readStatusFromFile = false

protected

Definition at line 109 of file G4WorkerRunManager.hh.

Referenced by GenerateEvent(), G4WorkerSubEvtRunManager::GenerateEvent(), G4WorkerTaskRunManager::GenerateEvent(), RestoreRndmEachEvent(), G4WorkerSubEvtRunManager::RestoreRndmEachEvent(), and G4WorkerTaskRunManager::RestoreRndmEachEvent().

◆ runIsSeeded

G4bool G4WorkerRunManager::runIsSeeded = false

protected

Definition at line 104 of file G4WorkerRunManager.hh.

Referenced by DoEventLoop(), GenerateEvent(), G4WorkerSubEvtRunManager::GenerateEvent(), G4WorkerTaskRunManager::GenerateEvent(), G4WorkerSubEvtRunManager::RunInitialization(), and G4WorkerTaskRunManager::RunInitialization().

◆ seedsQueue

G4SeedsQueue G4WorkerRunManager::seedsQueue

protected

Definition at line 108 of file G4WorkerRunManager.hh.

Referenced by DoEventLoop(), G4WorkerSubEvtRunManager::DoEventLoop(), G4WorkerTaskRunManager::DoEventLoop(), GenerateEvent(), G4WorkerSubEvtRunManager::GenerateEvent(), and G4WorkerTaskRunManager::GenerateEvent().

◆ workerContext

G4WorkerThread* G4WorkerRunManager::workerContext = nullptr

protected

Definition at line 98 of file G4WorkerRunManager.hh.

Referenced by DoWork(), G4WorkerSubEvtRunManager::DoWork(), G4WorkerTaskRunManager::DoWork(), G4WorkerTaskRunManager::GetWorkerThread(), rndmSaveThisEvent(), rndmSaveThisRun(), SetWorkerThread(), StoreRNGStatus(), G4WorkerSubEvtRunManager::StoreRNGStatus(), G4WorkerTaskRunManager::StoreRNGStatus(), G4WorkerSubEvtRunManager::TerminateEventLoop(), and G4WorkerTaskRunManager::TerminateEventLoop().

The documentation for this class was generated from the following files:

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4WorkerRunManager()

◆ ~G4WorkerRunManager()

Member Function Documentation

◆ ConstructScoringWorlds()

◆ DoEventLoop()

◆ DoWork()

◆ GenerateEvent()

◆ GetWorkerRunManager()

◆ GetWorkerRunManagerKernel()

◆ InitializeGeometry()

◆ MergePartialResults()

◆ ProcessOneEvent()

◆ RestoreRndmEachEvent()

◆ rndmSaveThisEvent()

◆ rndmSaveThisRun()

◆ RunInitialization()

◆ RunTermination()

◆ SetupDefaultRNGEngine()

◆ SetUserAction() [1/6]

◆ SetUserAction() [2/6]

◆ SetUserAction() [3/6]

◆ SetUserAction() [4/6]

◆ SetUserAction() [5/6]

◆ SetUserAction() [6/6]

◆ SetUserInitialization() [1/5]

◆ SetUserInitialization() [2/5]

◆ SetUserInitialization() [3/5]

◆ SetUserInitialization() [4/5]

◆ SetUserInitialization() [5/5]

◆ SetWorkerThread()

◆ StoreRNGStatus()

◆ TerminateEventLoop()

Member Data Documentation

◆ currEvID

◆ eventLoopOnGoing

◆ luxury

◆ nevModulo

◆ readStatusFromFile

◆ runIsSeeded

◆ seedsQueue

◆ workerContext