GGSRunManager.cpp

Go to the documentation of this file.

/*
 * GGSRunManager.cpp
 *
 *  Created on: 19/ott/2013
 *      Author: mori
 */

#include "montecarlo/mccore/GGSRunManager.h"

#include "G4Run.hh"
#include "G4TiMemory.hh"
#include "G4UImanager.hh"

void GGSRunManager::DoEventLoop(G4int n_event, const char *macroFile, G4int n_select) {
  // Implementation copy-pasted from G4 10.06.p02 with no modification, to just allow for binding the call to
  // ProcessOneEvent to the override defined in this class.
  TIMEMORY_AUTO_TIMER("");
  InitializeEventLoop(n_event, macroFile, n_select);

  // Event loop
  for (G4int i_event = 0; i_event < n_event; i_event++) {
    ProcessOneEvent(i_event);
    TerminateOneEvent();
    if (runAborted)
      break;
  }

  // For G4MTRunManager, TerminateEventLoop() is invoked after all threads are finished.
  if (runManagerType == sequentialRM)
    TerminateEventLoop();
}

void GGSRunManager::ProcessOneEvent(G4int i_event) {
  TIMEMORY_AUTO_TIMER("");
  currentEvent = _SetupRandomNumbers(i_event);
  BeginOfEventProcessing();
  userPrimaryGeneratorAction->GeneratePrimaries(currentEvent);
  DumpStatus(i_event);
  bool resim = true;
  while (resim) {
    BeginOfEventSimulation();
    eventManager->ProcessOneEvent(currentEvent);
    resim = IsKilledAndToBeSimulatedAgain(currentEvent);
    if (resim) {
      G4Event *newEvent = new G4Event(currentEvent->GetEventID());
      delete currentEvent;
      currentEvent = newEvent;
      userPrimaryGeneratorAction->GeneratePrimaries(newEvent);
    }
  }
  AnalyzeEvent(currentEvent);
  UpdateScoring();
  if (i_event < n_select_msg)
    G4UImanager::GetUIpointer()->ApplyCommand(msgText);
}

G4Event *GGSRunManager::_SetupRandomNumbers(G4int i_event) {
  // Implementation copy-pasted from G4RunManager::GenerateEvent in G4 10.06.p02, with these modifications:
  // 1. don't generate of the event at the end.
  // 2. offload the handling of seeds files to GGSRunManagerExtensions::HandleEventRngStateFile

  TIMEMORY_AUTO_TIMER("");
  if (!userPrimaryGeneratorAction) {
    G4Exception("G4RunManager::GenerateEvent()", "Run0032", FatalException,
                "G4VUserPrimaryGeneratorAction is not defined!");
    return 0;
  }

  HandleEventRngStateFile(currentRun->GetRunID(), i_event);

  G4Event *anEvent = new G4Event(i_event);

  if (storeRandomNumberStatusToG4Event == 1 || storeRandomNumberStatusToG4Event == 3) {
    std::ostringstream oss;
    G4Random::saveFullState(oss);
    randomNumberStatusForThisEvent = oss.str();
    anEvent->SetRandomNumberStatus(randomNumberStatusForThisEvent);
  }

  if (printModulo > 0 && anEvent->GetEventID() % printModulo == 0) {
    G4cout << "--> Event " << anEvent->GetEventID() << " starts." << G4endl;
  }
  return anEvent;
}