#include <RunManagerMT.h>

Public Member Functions
const SensitiveDetectorCatalog &	catalog () const

void	Connect (RunAction *)

const std::vector< std::string > &	G4Commands () const

void	initG4 (const DDCompactView , const cms::DDCompactView , const HepPDT::ParticleDataTable *)

void	initializeUserActions ()

PhysicsList *	physicsListForWorker () const

	RunManagerMT (edm::ParameterSet const &)

void	stopG4 ()

const DDDWorld &	world () const

	~RunManagerMT ()

Private Member Functions
void	terminateRun ()

Private Attributes
SensitiveDetectorCatalog	m_catalog

bool	m_check

G4Run *	m_currentRun

std::vector< std::string >	m_G4Commands

edm::ParameterSet	m_g4overlap

G4MTRunManagerKernel *	m_kernel

bool	m_managerInitialized

edm::ParameterSet	m_p

std::unique_ptr< PhysicsList >	m_physicsList

const std::string	m_PhysicsTablesDir

edm::ParameterSet	m_pPhysics

edm::ParameterSet	m_pRunAction

SimActivityRegistry	m_registry

bool	m_RestorePhysicsTables

std::unique_ptr< SimRunInterface >	m_runInterface

bool	m_runTerminated

G4StateManager *	m_stateManager

bool	m_StorePhysicsTables

CustomUIsession *	m_UIsession

RunAction *	m_userRunAction

std::unique_ptr< DDDWorld >	m_world

Friends
class	RunManagerMTWorker

Detailed Description

Definition at line 51 of file RunManagerMT.h.

Constructor & Destructor Documentation

RunManagerMT::RunManagerMT ( edm::ParameterSet const & p )

explicit

Definition at line 67 of file RunManagerMT.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), m_check, m_currentRun, m_kernel, m_physicsList, m_runInterface, m_stateManager, m_UIsession, m_userRunAction, and m_world.

     : m_managerInitialized(false),
       m_runTerminated(false),
       m_PhysicsTablesDir(p.getUntrackedParameter<std::string>("PhysicsTablesDirectory", "")),
       m_StorePhysicsTables(p.getUntrackedParameter<bool>("StorePhysicsTables", false)),
       m_RestorePhysicsTables(p.getUntrackedParameter<bool>("RestorePhysicsTables", false)),
       m_pPhysics(p.getParameter<edm::ParameterSet>("Physics")),
       m_pRunAction(p.getParameter<edm::ParameterSet>("RunAction")),
       m_g4overlap(p.getUntrackedParameter<edm::ParameterSet>("G4CheckOverlap")),
       m_G4Commands(p.getParameter<std::vector<std::string> >("G4Commands")),
       m_p(p) {
   m_currentRun = nullptr;
   m_UIsession = new CustomUIsession();
   m_physicsList.reset(nullptr);
   m_world.reset(nullptr);
 
   m_runInterface.reset(nullptr);
   m_userRunAction = nullptr;
   m_currentRun = nullptr;
 
   m_kernel = new G4MTRunManagerKernel();
   m_stateManager = G4StateManager::GetStateManager();
   double th = p.getParameter<double>("ThresholdForGeometryExceptions") * CLHEP::GeV;
   bool tr = p.getParameter<bool>("TraceExceptions");
   m_stateManager->SetExceptionHandler(new ExceptionHandler(th, tr));
   m_check = p.getUntrackedParameter<bool>("CheckGeometry", false);
 }

RunManagerMT::~RunManagerMT ( )

Definition at line 95 of file RunManagerMT.cc.

References m_UIsession.

95 { delete m_UIsession; }

RunManagerMT::m_UIsession

CustomUIsession * m_UIsession

Definition: RunManagerMT.h:85

Member Function Documentation

const SensitiveDetectorCatalog& RunManagerMT::catalog ( ) const

inline

Definition at line 72 of file RunManagerMT.h.

References m_catalog.

Referenced by RunManagerMTWorker::initializeG4().

72 { return m_catalog; }

RunManagerMT::m_catalog

SensitiveDetectorCatalog m_catalog

Definition: RunManagerMT.h:107

void RunManagerMT::Connect ( RunAction * runAction )

Definition at line 266 of file RunManagerMT.cc.

References SimActivityRegistry::beginOfRunSignal_, sim_act::Signaler< T >::connect(), SimActivityRegistry::endOfRunSignal_, RunAction::m_beginOfRunSignal, RunAction::m_endOfRunSignal, and m_registry.

Referenced by SimRunInterface::Connect(), and initializeUserActions().

                                                {
   runAction->m_beginOfRunSignal.connect(m_registry.beginOfRunSignal_);
   runAction->m_endOfRunSignal.connect(m_registry.endOfRunSignal_);
 }

const std::vector<std::string>& RunManagerMT::G4Commands ( ) const

inline

Definition at line 74 of file RunManagerMT.h.

References m_G4Commands.

Referenced by RunManagerMTWorker::initializeG4().

74 { return m_G4Commands; }

RunManagerMT::m_G4Commands

std::vector< std::string > m_G4Commands

Definition: RunManagerMT.h:102

void RunManagerMT::initG4	(	const DDCompactView *	pDD,
		const cms::DDCompactView *	pDD4hep,
		const HepPDT::ParticleDataTable *	fPDGTable
	)

Definition at line 97 of file RunManagerMT.cc.

References RunAction::BeginOfRunAction(), trackerTree::check(), mps_setup::cmd, mps_check::command, beamerCreator::create(), cuts, SimActivityRegistry::dddWorldSignal_, DeadROC_duringRun::dir, Exception, get, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), recoMuon::in, initializeUserActions(), m_catalog, m_check, m_currentRun, m_G4Commands, m_g4overlap, m_kernel, m_managerInitialized, m_p, m_physicsList, m_PhysicsTablesDir, m_pPhysics, m_registry, m_RestorePhysicsTables, m_stateManager, m_StorePhysicsTables, m_UIsession, m_userRunAction, m_world, SiStripPI::max, gpuVertexFinder::nn, AlCaHLTBitMon_QueryRunRegistry::string, and world().

                                                                     {
   if (m_managerInitialized) {
     edm::LogWarning("SimG4CoreApplication") << "RunManagerMT::initG4 was already done - exit";
     return;
   }
   bool geoFromDD4hep = m_p.getParameter<bool>("g4GeometryDD4hepSource");
   bool cuts = m_pPhysics.getParameter<bool>("CutsPerRegion");
   bool protonCut = m_pPhysics.getParameter<bool>("CutsOnProton");
   int verb = m_pPhysics.getUntrackedParameter<int>("Verbosity", 0);
   int stepverb = m_p.getUntrackedParameter<int>("SteppingVerbosity", 0);
   edm::LogVerbatim("SimG4CoreApplication")
       << "RunManagerMT: start initialising of geometry DD4hep: " << geoFromDD4hep << "\n"
       << "              cutsPerRegion: " << cuts << " cutForProton: " << protonCut << "\n"
       << "              G4 verbosity: " << verb;
 
   G4Timer timer;
   timer.Start();
 
   G4UImanager::GetUIpointer()->SetCoutDestination(m_UIsession);
   G4UImanager::GetUIpointer()->SetMasterUIManager(true);
 
   m_world = std::make_unique<DDDWorld>(pDD, pDD4hep, m_catalog, verb, cuts, protonCut);
   G4VPhysicalVolume* world = m_world.get()->GetWorldVolume();
 
   m_kernel->SetVerboseLevel(verb);
   edm::LogVerbatim("SimG4CoreApplication")
       << "RunManagerMT: Define cuts: " << cuts << " Geant4 run manager verbosity: " << verb;
 
   const G4RegionStore* regStore = G4RegionStore::GetInstance();
   const G4PhysicalVolumeStore* pvs = G4PhysicalVolumeStore::GetInstance();
   const G4LogicalVolumeStore* lvs = G4LogicalVolumeStore::GetInstance();
   unsigned int numPV = pvs->size();
   unsigned int numLV = lvs->size();
   unsigned int nn = regStore->size();
   edm::LogVerbatim("SimG4CoreApplication")
       << "RunManagerMT: " << numPV << " physical volumes; " << numLV << " logical volumes; " << nn << " regions.";
 
   if (m_check) {
     m_kernel->SetVerboseLevel(2);
   }
   m_kernel->DefineWorldVolume(world, true);
   m_registry.dddWorldSignal_(m_world.get());
   m_stateManager->SetNewState(G4State_PreInit);
 
   // Create physics list
   edm::LogVerbatim("SimG4CoreApplication") << "RunManagerMT: create PhysicsList";
 
   std::unique_ptr<PhysicsListMakerBase> physicsMaker(
       PhysicsListFactory::get()->create(m_pPhysics.getParameter<std::string>("type")));
   if (physicsMaker.get() == nullptr) {
     throw cms::Exception("Configuration") << "Unable to find the Physics list requested";
   }
   m_physicsList = physicsMaker->make(m_pPhysics, m_registry);
 
   PhysicsList* phys = m_physicsList.get();
   if (phys == nullptr) {
     throw cms::Exception("Configuration") << "Physics list construction failed!";
   }
   if (stepverb > 0) {
     verb = std::max(verb, 1);
   }
   G4HadronicParameters::Instance()->SetVerboseLevel(verb);
   G4EmParameters::Instance()->SetVerbose(verb);
   G4EmParameters::Instance()->SetWorkerVerbose(std::max(verb - 1, 0));
 
   // exotic particle physics
   double monopoleMass = m_pPhysics.getUntrackedParameter<double>("MonopoleMass", 0);
   if (monopoleMass > 0.0) {
     phys->RegisterPhysics(new CMSMonopolePhysics(fPDGTable, m_pPhysics));
   }
   bool exotica = m_pPhysics.getUntrackedParameter<bool>("ExoticaTransport", false);
   if (exotica) {
     CMSExoticaPhysics exo(phys, m_pPhysics);
   }
 
   // adding GFlash, Russian Roulette for eletrons and gamma,
   // step limiters on top of any Physics Lists
   phys->RegisterPhysics(new ParametrisedEMPhysics("EMoptions", m_pPhysics));
 
   if (m_RestorePhysicsTables) {
     m_physicsList->SetPhysicsTableRetrieved(m_PhysicsTablesDir);
   }
   edm::LogVerbatim("SimG4CoreApplication") << "RunManagerMT: start initialisation of PhysicsList for master";
 
   m_physicsList->SetDefaultCutValue(m_pPhysics.getParameter<double>("DefaultCutValue") * CLHEP::cm);
   m_physicsList->SetCutsWithDefault();
   m_kernel->SetPhysics(phys);
 
   edm::LogVerbatim("SimG4CoreApplication") << "RunManagerMT: PhysicsList and cuts are defined";
 
   // Geant4 UI commands before initialisation of physics
   if (!m_G4Commands.empty()) {
     edm::LogVerbatim("SimG4CoreApplication") << "RunManagerMT: Requested UI commands: ";
     for (const std::string& command : m_G4Commands) {
       edm::LogVerbatim("SimG4CoreApplication") << "    " << command;
       G4UImanager::GetUIpointer()->ApplyCommand(command);
     }
   }
 
   m_stateManager->SetNewState(G4State_Init);
   edm::LogVerbatim("SimG4CoreApplication") << "RunManagerMT: G4State is Init";
   m_kernel->InitializePhysics();
   m_kernel->SetUpDecayChannels();
 
   if (m_kernel->RunInitialization()) {
     m_managerInitialized = true;
   } else {
     throw cms::Exception("LogicError") << "G4RunManagerKernel initialization failed!";
   }
 
   if (m_StorePhysicsTables) {
     std::ostringstream dir;
     dir << m_PhysicsTablesDir << '\0';
     std::string cmd = std::string("/control/shell mkdir -p ") + m_PhysicsTablesDir;
     if (!std::ifstream(dir.str().c_str(), std::ios::in))
       G4UImanager::GetUIpointer()->ApplyCommand(cmd);
     m_physicsList->StorePhysicsTable(m_PhysicsTablesDir);
   }
   // Appload nuclear level data up to Z=84
   G4NuclearLevelData::GetInstance()->UploadNuclearLevelData(84);
 
   if (verb > 1) {
     m_physicsList->DumpCutValuesTable();
   }
   edm::LogVerbatim("SimG4CoreApplication")
       << "RunManagerMT: Physics is initilized, now initialise user actions, verb=" << verb;
 
   initializeUserActions();
 
   // geometry dump
   auto writeFile = m_p.getUntrackedParameter<std::string>("FileNameGDML", "");
   if (!writeFile.empty()) {
     G4GDMLParser gdml;
     gdml.SetRegionExport(true);
     gdml.SetEnergyCutsExport(true);
     gdml.Write(writeFile, m_world->GetWorldVolume(), true);
   }
 
   // G4Region dump file name
   auto regionFile = m_p.getUntrackedParameter<std::string>("FileNameRegions", "");
 
   // Geometry checks
   if (m_check || !regionFile.empty()) {
     CMSG4CheckOverlap check(m_g4overlap, regionFile, m_UIsession, world);
   }
 
   m_stateManager->SetNewState(G4State_PreInit);
   G4HadronicParameters::Instance()->SetVerboseLevel(std::max(verb - 1, 0));
 
   // If the Geant4 particle table is needed, decomment the lines below
   //
   //G4ParticleTable::GetParticleTable()->DumpTable("ALL");
   //
   m_stateManager->SetNewState(G4State_GeomClosed);
   m_currentRun = new G4Run();
   m_userRunAction->BeginOfRunAction(m_currentRun);
   timer.Stop();
   G4cout.precision(4);
   G4cout << "RunManagerMT: initG4 done " << timer << G4endl;
 }

void RunManagerMT::initializeUserActions ( )

Definition at line 260 of file RunManagerMT.cc.

References Connect(), m_pRunAction, m_runInterface, and m_userRunAction.

Referenced by initG4().

                                          {
   m_runInterface = std::make_unique<SimRunInterface>(this, true);
   m_userRunAction = new RunAction(m_pRunAction, m_runInterface.get(), true);
   Connect(m_userRunAction);
 }

PhysicsList* RunManagerMT::physicsListForWorker ( ) const

inline

Definition at line 78 of file RunManagerMT.h.

References m_physicsList.

Referenced by RunManagerMTWorker::initializeG4().

78 { return m_physicsList.get(); }

RunManagerMT::m_physicsList

std::unique_ptr< PhysicsList > m_physicsList

Definition: RunManagerMT.h:86

void RunManagerMT::stopG4 ( )

Definition at line 271 of file RunManagerMT.cc.

References m_runTerminated, m_stateManager, and terminateRun().

                           {
   G4GeometryManager::GetInstance()->OpenGeometry();
   m_stateManager->SetNewState(G4State_Quit);
   if (!m_runTerminated) {
     terminateRun();
   }
 }

void RunManagerMT::terminateRun ( )

private

Definition at line 279 of file RunManagerMT.cc.

References RunAction::EndOfRunAction(), m_currentRun, m_kernel, m_runTerminated, and m_userRunAction.

Referenced by stopG4().

                                 {
   if (nullptr != m_userRunAction) {
     m_userRunAction->EndOfRunAction(m_currentRun);
     delete m_userRunAction;
     m_userRunAction = nullptr;
   }
   if (!m_runTerminated) {
     m_kernel->RunTermination();
   }
   m_runTerminated = true;
   edm::LogVerbatim("SimG4CoreApplication") << "RunManagerMT:: terminateRun done";
 }