G4AtlasSvc.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#include "G4AtlasSvc.h"
#include "GaudiKernel/ServiceHandle.h"
 
// header files from Geant4
#include "G4VUserPhysicsList.hh"
#include "G4VModularPhysicsList.hh"
#include "G4ParallelWorldPhysics.hh"
#include "G4RunManager.hh"
 
G4AtlasSvc::G4AtlasSvc( const std::string& name, ISvcLocator* pSvcLocator )
  : base_class(name,pSvcLocator)
{
}
 
StatusCode G4AtlasSvc::initialize(){
  // go through all tools and retrieve them
  //  This fires initialize() for each of those tools
 
  ATH_MSG_DEBUG( "this is G4AtlasSvc::initialize() " );
  auto* rm = G4RunManager::GetRunManager();
  if (!rm) {
    ATH_MSG_ERROR("Run manager retrieval has failed");
    return StatusCode::FAILURE;
  }
  rm->Initialize();     // Initialization differs slightly in multi-threading.
  // TODO: add more details about why this is here.
  if (!m_isMT && rm->ConfirmBeamOnCondition()) {
    rm->RunInitialization();
  }
 
  ATH_CHECK(m_detGeoSvc.retrieve());
 
  ATH_CHECK(m_physicsListSvc.retrieve());
  ATH_CHECK(m_userLimitsSvc.retrieve());
 
  if (m_activateParallelGeometries) {
    G4VModularPhysicsList* thePhysicsList=dynamic_cast<G4VModularPhysicsList*>(m_physicsListSvc->GetPhysicsList());
    if (!thePhysicsList) {
      ATH_MSG_FATAL("Failed dynamic_cast!! this is not a G4VModularPhysicsList!");
      return StatusCode::FAILURE;
    }
#if G4VERSION_NUMBER >= 1010
      std::vector<std::string>& parallelWorldNames=m_detGeoSvc->GetParallelWorldNames();
      for (auto& it: parallelWorldNames) {
        thePhysicsList->RegisterPhysics(new G4ParallelWorldPhysics(it,true));
      }
#endif
  }
 
  return StatusCode::SUCCESS;
}