CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes | Friends
RunManagerMT Class Reference

#include <RunManagerMT.h>

Public Member Functions

const SensitiveDetectorCatalogcatalog () const
 
void Connect (RunAction *)
 
const std::vector< std::string > & G4Commands () const
 
void initG4 (const DDCompactView *, const cms::DDCompactView *, const HepPDT::ParticleDataTable *)
 
void initializeUserActions ()
 
PhysicsListphysicsListForWorker () const
 
 RunManagerMT (edm::ParameterSet const &)
 
void stopG4 ()
 
const DDDWorldworld () 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< PhysicsListm_physicsList
 
const std::string m_PhysicsTablesDir
 
edm::ParameterSet m_pPhysics
 
edm::ParameterSet m_pRunAction
 
SimActivityRegistry m_registry
 
bool m_RestorePhysicsTables
 
std::unique_ptr< SimRunInterfacem_runInterface
 
bool m_runTerminated
 
G4StateManager * m_stateManager
 
bool m_StorePhysicsTables
 
CustomUIsessionm_UIsession
 
RunActionm_userRunAction
 
std::unique_ptr< DDDWorldm_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 66 of file RunManagerMT.cc.

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

67  : m_managerInitialized(false),
68  m_runTerminated(false),
69  m_PhysicsTablesDir(p.getUntrackedParameter<std::string>("PhysicsTablesDirectory", "")),
70  m_StorePhysicsTables(p.getUntrackedParameter<bool>("StorePhysicsTables", false)),
71  m_RestorePhysicsTables(p.getUntrackedParameter<bool>("RestorePhysicsTables", false)),
72  m_pPhysics(p.getParameter<edm::ParameterSet>("Physics")),
73  m_pRunAction(p.getParameter<edm::ParameterSet>("RunAction")),
74  m_g4overlap(p.getUntrackedParameter<edm::ParameterSet>("G4CheckOverlap")),
75  m_G4Commands(p.getParameter<std::vector<std::string> >("G4Commands")),
76  m_p(p) {
77  m_currentRun = nullptr;
79  m_physicsList.reset(nullptr);
80  m_world.reset(nullptr);
81 
82  m_runInterface.reset(nullptr);
83  m_userRunAction = nullptr;
84  m_currentRun = nullptr;
85 
86  m_kernel = new G4MTRunManagerKernel();
87  m_stateManager = G4StateManager::GetStateManager();
88  double th = p.getParameter<double>("ThresholdForGeometryExceptions") * CLHEP::GeV;
89  m_stateManager->SetExceptionHandler(new ExceptionHandler(th));
90  m_check = p.getUntrackedParameter<bool>("CheckGeometry", false);
91 }
const std::string m_PhysicsTablesDir
Definition: RunManagerMT.h:95
bool m_managerInitialized
Definition: RunManagerMT.h:87
const double GeV
Definition: MathUtil.h:16
G4MTRunManagerKernel * m_kernel
Definition: RunManagerMT.h:83
edm::ParameterSet m_pRunAction
Definition: RunManagerMT.h:100
bool m_runTerminated
Definition: RunManagerMT.h:88
G4StateManager * m_stateManager
Definition: RunManagerMT.h:91
std::vector< std::string > m_G4Commands
Definition: RunManagerMT.h:102
G4Run * m_currentRun
Definition: RunManagerMT.h:90
edm::ParameterSet m_p
Definition: RunManagerMT.h:103
bool m_StorePhysicsTables
Definition: RunManagerMT.h:96
std::unique_ptr< PhysicsList > m_physicsList
Definition: RunManagerMT.h:86
std::unique_ptr< DDDWorld > m_world
Definition: RunManagerMT.h:105
edm::ParameterSet m_pPhysics
Definition: RunManagerMT.h:99
bool m_RestorePhysicsTables
Definition: RunManagerMT.h:97
std::unique_ptr< SimRunInterface > m_runInterface
Definition: RunManagerMT.h:93
edm::ParameterSet m_g4overlap
Definition: RunManagerMT.h:101
RunAction * m_userRunAction
Definition: RunManagerMT.h:89
CustomUIsession * m_UIsession
Definition: RunManagerMT.h:85
RunManagerMT::~RunManagerMT ( )

Definition at line 93 of file RunManagerMT.cc.

References m_UIsession.

93 { delete 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; }
SensitiveDetectorCatalog m_catalog
Definition: RunManagerMT.h:107
void RunManagerMT::Connect ( RunAction runAction)

Definition at line 264 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().

264  {
267 }
SimActivityRegistry::EndOfRunSignal m_endOfRunSignal
Definition: RunAction.h:24
BeginOfRunSignal beginOfRunSignal_
EndOfRunSignal endOfRunSignal_
void connect(Observer< const T * > *iObs)
does not take ownership of memory
Definition: Signaler.h:55
SimActivityRegistry::BeginOfRunSignal m_beginOfRunSignal
Definition: RunAction.h:23
SimActivityRegistry m_registry
Definition: RunManagerMT.h:106
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; }
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 95 of file RunManagerMT.cc.

References RunAction::BeginOfRunAction(), trackerTree::check(), mps_setup::cmd, mps_check::command, edm::errors::Configuration, beamerCreator::create(), cuts, SimActivityRegistry::dddWorldSignal_, DeadROC_duringRun::dir, Exception, get, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), recoMuon::in, initializeUserActions(), edm::errors::LogicError, 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().

97  {
99  edm::LogWarning("SimG4CoreApplication") << "RunManagerMT::initG4 was already done - exit";
100  return;
101  }
102  bool geoFromDD4hep = m_p.getParameter<bool>("g4GeometryDD4hepSource");
103  bool cuts = m_pPhysics.getParameter<bool>("CutsPerRegion");
104  bool protonCut = m_pPhysics.getParameter<bool>("CutsOnProton");
105  int verb = m_pPhysics.getUntrackedParameter<int>("Verbosity", 0);
106  int stepverb = m_p.getUntrackedParameter<int>("SteppingVerbosity", 0);
107  edm::LogVerbatim("SimG4CoreApplication")
108  << "RunManagerMT: start initialising of geometry DD4hep: " << geoFromDD4hep << "\n"
109  << " cutsPerRegion: " << cuts << " cutForProton: " << protonCut << "\n"
110  << " G4 verbosity: " << verb;
111 
112  G4Timer timer;
113  timer.Start();
114 
115  G4UImanager::GetUIpointer()->SetCoutDestination(m_UIsession);
116  G4UImanager::GetUIpointer()->SetMasterUIManager(true);
117 
118  m_world = std::make_unique<DDDWorld>(pDD, pDD4hep, m_catalog, verb, cuts, protonCut);
119  G4VPhysicalVolume* world = m_world.get()->GetWorldVolume();
120 
121  m_kernel->SetVerboseLevel(verb);
122  edm::LogVerbatim("SimG4CoreApplication")
123  << "RunManagerMT: Define cuts: " << cuts << " Geant4 run manager verbosity: " << verb;
124 
125  const G4RegionStore* regStore = G4RegionStore::GetInstance();
126  const G4PhysicalVolumeStore* pvs = G4PhysicalVolumeStore::GetInstance();
127  const G4LogicalVolumeStore* lvs = G4LogicalVolumeStore::GetInstance();
128  unsigned int numPV = pvs->size();
129  unsigned int numLV = lvs->size();
130  unsigned int nn = regStore->size();
131  edm::LogVerbatim("SimG4CoreApplication")
132  << "RunManagerMT: " << numPV << " physical volumes; " << numLV << " logical volumes; " << nn << " regions.";
133 
134  if (m_check) {
135  m_kernel->SetVerboseLevel(2);
136  }
137  m_kernel->DefineWorldVolume(world, true);
139  m_stateManager->SetNewState(G4State_PreInit);
140 
141  // Create physics list
142  edm::LogVerbatim("SimG4CoreApplication") << "RunManagerMT: create PhysicsList";
143 
144  std::unique_ptr<PhysicsListMakerBase> physicsMaker(
146  if (physicsMaker.get() == nullptr) {
147  throw edm::Exception(edm::errors::Configuration) << "Unable to find the Physics list requested";
148  }
149  m_physicsList = physicsMaker->make(m_pPhysics, m_registry);
150 
151  PhysicsList* phys = m_physicsList.get();
152  if (phys == nullptr) {
153  throw edm::Exception(edm::errors::Configuration, "Physics list construction failed!");
154  }
155  if (stepverb > 0) {
156  verb = std::max(verb, 1);
157  }
158  G4HadronicParameters::Instance()->SetVerboseLevel(verb);
159  G4EmParameters::Instance()->SetVerbose(verb);
160  G4EmParameters::Instance()->SetWorkerVerbose(std::max(verb - 1, 0));
161 
162  // exotic particle physics
163  double monopoleMass = m_pPhysics.getUntrackedParameter<double>("MonopoleMass", 0);
164  if (monopoleMass > 0.0) {
165  phys->RegisterPhysics(new CMSMonopolePhysics(fPDGTable, m_pPhysics));
166  }
167  bool exotica = m_pPhysics.getUntrackedParameter<bool>("ExoticaTransport", false);
168  if (exotica) {
169  CMSExoticaPhysics exo(phys, m_pPhysics);
170  }
171 
172  // adding GFlash, Russian Roulette for eletrons and gamma,
173  // step limiters on top of any Physics Lists
174  phys->RegisterPhysics(new ParametrisedEMPhysics("EMoptions", m_pPhysics));
175 
177  m_physicsList->SetPhysicsTableRetrieved(m_PhysicsTablesDir);
178  }
179  edm::LogVerbatim("SimG4CoreApplication") << "RunManagerMT: start initialisation of PhysicsList for master";
180 
181  m_physicsList->SetDefaultCutValue(m_pPhysics.getParameter<double>("DefaultCutValue") * CLHEP::cm);
182  m_physicsList->SetCutsWithDefault();
183  m_kernel->SetPhysics(phys);
184 
185  edm::LogVerbatim("SimG4CoreApplication") << "RunManagerMT: PhysicsList and cuts are defined";
186 
187  // Geant4 UI commands before initialisation of physics
188  if (!m_G4Commands.empty()) {
189  edm::LogVerbatim("SimG4CoreApplication") << "RunManagerMT: Requested UI commands: ";
190  for (const std::string& command : m_G4Commands) {
191  edm::LogVerbatim("SimG4CoreApplication") << " " << command;
192  G4UImanager::GetUIpointer()->ApplyCommand(command);
193  }
194  }
195 
196  m_stateManager->SetNewState(G4State_Init);
197  edm::LogVerbatim("SimG4CoreApplication") << "RunManagerMT: G4State is Init";
198  m_kernel->InitializePhysics();
199  m_kernel->SetUpDecayChannels();
200 
201  if (m_kernel->RunInitialization()) {
202  m_managerInitialized = true;
203  } else {
204  throw edm::Exception(edm::errors::LogicError, "G4RunManagerKernel initialization failed!");
205  }
206 
207  if (m_StorePhysicsTables) {
208  std::ostringstream dir;
209  dir << m_PhysicsTablesDir << '\0';
210  std::string cmd = std::string("/control/shell mkdir -p ") + m_PhysicsTablesDir;
211  if (!std::ifstream(dir.str().c_str(), std::ios::in))
212  G4UImanager::GetUIpointer()->ApplyCommand(cmd);
213  m_physicsList->StorePhysicsTable(m_PhysicsTablesDir);
214  }
215  // Appload nuclear level data up to Z=84
216  G4NuclearLevelData::GetInstance()->UploadNuclearLevelData(84);
217 
218  if (verb > 1) {
219  m_physicsList->DumpCutValuesTable();
220  }
221  edm::LogVerbatim("SimG4CoreApplication")
222  << "RunManagerMT: Physics is initilized, now initialise user actions, verb=" << verb;
223 
225 
226  // geometry dump
227  auto writeFile = m_p.getUntrackedParameter<std::string>("FileNameGDML", "");
228  if (!writeFile.empty()) {
229  G4GDMLParser gdml;
230  gdml.SetRegionExport(true);
231  gdml.SetEnergyCutsExport(true);
232  gdml.Write(writeFile, m_world->GetWorldVolume(), true);
233  }
234 
235  // G4Region dump file name
236  auto regionFile = m_p.getUntrackedParameter<std::string>("FileNameRegions", "");
237 
238  // Geometry checks
239  if (m_check || !regionFile.empty()) {
240  CMSG4CheckOverlap check(m_g4overlap, regionFile, m_UIsession, world);
241  }
242 
243  m_stateManager->SetNewState(G4State_PreInit);
244  G4HadronicParameters::Instance()->SetVerboseLevel(std::max(verb - 1, 0));
245 
246  // If the Geant4 particle table is needed, decomment the lines below
247  //
248  //G4ParticleTable::GetParticleTable()->DumpTable("ALL");
249  //
250  m_stateManager->SetNewState(G4State_GeomClosed);
251  m_currentRun = new G4Run();
253  timer.Stop();
254  G4cout.precision(4);
255  G4cout << "RunManagerMT: initG4 done " << timer << G4endl;
256 }
const std::string m_PhysicsTablesDir
Definition: RunManagerMT.h:95
Log< level::Info, true > LogVerbatim
bool m_managerInitialized
Definition: RunManagerMT.h:87
T getUntrackedParameter(std::string const &, T const &) const
void BeginOfRunAction(const G4Run *aRun) override
Definition: RunAction.cc:17
G4MTRunManagerKernel * m_kernel
Definition: RunManagerMT.h:83
SensitiveDetectorCatalog m_catalog
Definition: RunManagerMT.h:107
G4StateManager * m_stateManager
Definition: RunManagerMT.h:91
std::vector< std::string > m_G4Commands
Definition: RunManagerMT.h:102
G4Run * m_currentRun
Definition: RunManagerMT.h:90
void initializeUserActions()
const DDDWorld & world() const
Definition: RunManagerMT.h:70
TkSoA const *__restrict__ CAHitNtupletGeneratorKernelsGPU::QualityCuts cuts
edm::ParameterSet m_p
Definition: RunManagerMT.h:103
bool m_StorePhysicsTables
Definition: RunManagerMT.h:96
std::unique_ptr< PhysicsList > m_physicsList
Definition: RunManagerMT.h:86
std::unique_ptr< DDDWorld > m_world
Definition: RunManagerMT.h:105
edm::ParameterSet m_pPhysics
Definition: RunManagerMT.h:99
DDDWorldSignal dddWorldSignal_
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
bool m_RestorePhysicsTables
Definition: RunManagerMT.h:97
list command
Definition: mps_check.py:25
list cmd
Definition: mps_setup.py:244
int32_t *__restrict__ nn
edm::ParameterSet m_g4overlap
Definition: RunManagerMT.h:101
#define get
Log< level::Warning, false > LogWarning
RunAction * m_userRunAction
Definition: RunManagerMT.h:89
CustomUIsession * m_UIsession
Definition: RunManagerMT.h:85
SimActivityRegistry m_registry
Definition: RunManagerMT.h:106
void RunManagerMT::initializeUserActions ( )

Definition at line 258 of file RunManagerMT.cc.

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

Referenced by initG4().

258  {
259  m_runInterface = std::make_unique<SimRunInterface>(this, true);
261  Connect(m_userRunAction);
262 }
void Connect(RunAction *)
edm::ParameterSet m_pRunAction
Definition: RunManagerMT.h:100
std::unique_ptr< SimRunInterface > m_runInterface
Definition: RunManagerMT.h:93
RunAction * m_userRunAction
Definition: RunManagerMT.h:89
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(); }
std::unique_ptr< PhysicsList > m_physicsList
Definition: RunManagerMT.h:86
void RunManagerMT::stopG4 ( )

Definition at line 269 of file RunManagerMT.cc.

References m_runTerminated, m_stateManager, and terminateRun().

269  {
270  G4GeometryManager::GetInstance()->OpenGeometry();
271  m_stateManager->SetNewState(G4State_Quit);
272  if (!m_runTerminated) {
273  terminateRun();
274  }
275 }
void terminateRun()
bool m_runTerminated
Definition: RunManagerMT.h:88
G4StateManager * m_stateManager
Definition: RunManagerMT.h:91
void RunManagerMT::terminateRun ( )
private

Definition at line 277 of file RunManagerMT.cc.

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

Referenced by stopG4().

277  {
278  if (nullptr != m_userRunAction) {
280  delete m_userRunAction;
281  m_userRunAction = nullptr;
282  }
283  if (!m_runTerminated) {
284  m_kernel->RunTermination();
285  }
286  m_runTerminated = true;
287  edm::LogVerbatim("SimG4CoreApplication") << "RunManagerMT:: terminateRun done";
288 }
Log< level::Info, true > LogVerbatim
G4MTRunManagerKernel * m_kernel
Definition: RunManagerMT.h:83
bool m_runTerminated
Definition: RunManagerMT.h:88
G4Run * m_currentRun
Definition: RunManagerMT.h:90
RunAction * m_userRunAction
Definition: RunManagerMT.h:89
void EndOfRunAction(const G4Run *aRun) override
Definition: RunAction.cc:26
const DDDWorld& RunManagerMT::world ( ) const
inline

Definition at line 70 of file RunManagerMT.h.

References m_world.

Referenced by initG4(), and RunManagerMTWorker::initializeG4().

70 { return *m_world; }
std::unique_ptr< DDDWorld > m_world
Definition: RunManagerMT.h:105

Friends And Related Function Documentation

friend class RunManagerMTWorker
friend

Definition at line 52 of file RunManagerMT.h.

Member Data Documentation

SensitiveDetectorCatalog RunManagerMT::m_catalog
private

Definition at line 107 of file RunManagerMT.h.

Referenced by catalog(), and initG4().

bool RunManagerMT::m_check
private

Definition at line 98 of file RunManagerMT.h.

Referenced by initG4(), and RunManagerMT().

G4Run* RunManagerMT::m_currentRun
private

Definition at line 90 of file RunManagerMT.h.

Referenced by initG4(), RunManagerMT(), and terminateRun().

std::vector<std::string> RunManagerMT::m_G4Commands
private

Definition at line 102 of file RunManagerMT.h.

Referenced by G4Commands(), and initG4().

edm::ParameterSet RunManagerMT::m_g4overlap
private

Definition at line 101 of file RunManagerMT.h.

Referenced by initG4().

G4MTRunManagerKernel* RunManagerMT::m_kernel
private

Definition at line 83 of file RunManagerMT.h.

Referenced by initG4(), RunManagerMT(), and terminateRun().

bool RunManagerMT::m_managerInitialized
private

Definition at line 87 of file RunManagerMT.h.

Referenced by initG4().

edm::ParameterSet RunManagerMT::m_p
private

Definition at line 103 of file RunManagerMT.h.

Referenced by initG4().

std::unique_ptr<PhysicsList> RunManagerMT::m_physicsList
private

Definition at line 86 of file RunManagerMT.h.

Referenced by initG4(), physicsListForWorker(), and RunManagerMT().

const std::string RunManagerMT::m_PhysicsTablesDir
private

Definition at line 95 of file RunManagerMT.h.

Referenced by initG4().

edm::ParameterSet RunManagerMT::m_pPhysics
private

Definition at line 99 of file RunManagerMT.h.

Referenced by initG4().

edm::ParameterSet RunManagerMT::m_pRunAction
private

Definition at line 100 of file RunManagerMT.h.

Referenced by initializeUserActions().

SimActivityRegistry RunManagerMT::m_registry
private

Definition at line 106 of file RunManagerMT.h.

Referenced by Connect(), and initG4().

bool RunManagerMT::m_RestorePhysicsTables
private

Definition at line 97 of file RunManagerMT.h.

Referenced by initG4().

std::unique_ptr<SimRunInterface> RunManagerMT::m_runInterface
private

Definition at line 93 of file RunManagerMT.h.

Referenced by initializeUserActions(), and RunManagerMT().

bool RunManagerMT::m_runTerminated
private

Definition at line 88 of file RunManagerMT.h.

Referenced by stopG4(), and terminateRun().

G4StateManager* RunManagerMT::m_stateManager
private

Definition at line 91 of file RunManagerMT.h.

Referenced by initG4(), RunManagerMT(), and stopG4().

bool RunManagerMT::m_StorePhysicsTables
private

Definition at line 96 of file RunManagerMT.h.

Referenced by initG4().

CustomUIsession* RunManagerMT::m_UIsession
private

Definition at line 85 of file RunManagerMT.h.

Referenced by initG4(), RunManagerMT(), and ~RunManagerMT().

RunAction* RunManagerMT::m_userRunAction
private

Definition at line 89 of file RunManagerMT.h.

Referenced by initG4(), initializeUserActions(), RunManagerMT(), and terminateRun().

std::unique_ptr<DDDWorld> RunManagerMT::m_world
private

Definition at line 105 of file RunManagerMT.h.

Referenced by initG4(), RunManagerMT(), and world().