GGSRngStateAction.cpp

Go to the documentation of this file.

/*
 * GGSRngStateAction.cpp
 *
 *  Created on: 4 Dec 2020
 *      Author: Nicola Mori
 */

// GGS headers
#include "montecarlo/useractions/GGSRngStateAction.h"

#include "montecarlo/mccore/GGSRunManager.h"
#include "montecarlo/pluginmanagers/GGSMCPluginMacros.h"

// Geant4 headers
#include "G4Event.hh"
#include "G4Run.hh"

// ROOT headers
#include "TFile.h"
#include "TTree.h"

RegisterUA(GGSRngStateAction);

GGSRngStateAction::GGSRngStateAction() : _outRootFile(nullptr), _outTree(nullptr) {}

void GGSRngStateAction::BeginOfEventAction(const G4Event *) {

  GGSRunManagerExtensions &rmExt = dynamic_cast<GGSRunManagerExtensions &>(*(G4RunManager::GetRunManager()));
  _rngState = rmExt.GetRandomStateAtBeginOfEvent();
}

void GGSRngStateAction::BeginOfRunAction(const G4Run *run) {
  _outRootFile = GGSRootFileService::GetInstance().GetFileForThisRun("", run);
  _outTree = GGSRootFileService::GetInstance().GetDefaultTree(_outRootFile);
  _outTree->SetTitle(TString(_outTree->GetTitle()) + "RngState ");
  _outTree->Branch("rngState", &_rngState);

  if (GGSRootFileService::GetInstance().GetOutputMode() == GGSRootFileService::OutputMode::MULTIFILE ||
      !(G4Threading::IsMultithreadedApplication())) {
    std::string rngName = G4Random::getTheEngine()->name();
    _outRootFile->WriteObjectAny(&rngName, TClass::GetClass("std::string"), "rngName");
  }
}

void GGSRngStateAction::EndOfRunAction(const G4Run *run) {
  GGSRootFileService::GetInstance().CloseFileForThisRun("", run);
  _outRootFile = nullptr;
  _outTree = nullptr;
}