39 #include "HepPDT/defs.h"
40 #include "HepPDT/TableBuilder.hh"
41 #include "HepPDT/ParticleDataTable.hh"
44 #include "G4StateManager.hh"
45 #include "G4ApplicationState.hh"
46 #include "G4RunManagerKernel.hh"
47 #include "G4UImanager.hh"
49 #include "G4EventManager.hh"
52 #include "G4TransportationManager.hh"
53 #include "G4ParticleTable.hh"
68 std::vector<boost::shared_ptr<SimWatcher> >& oWatchers,
69 std::vector<boost::shared_ptr<SimProducer> >& oProds
74 vector<ParameterSet> watchers;
76 watchers = iP.
getParameter<vector<ParameterSet> >(
"Watchers");
80 for(vector<ParameterSet>::iterator itWatcher = watchers.begin();
81 itWatcher != watchers.end();
83 std::auto_ptr<SimWatcherMakerBase>
maker(
85 (itWatcher->getParameter<std::string> (
"type")) );
90 boost::shared_ptr<SimWatcher> watcherTemp;
91 boost::shared_ptr<SimProducer> producerTemp;
92 maker->make(*itWatcher,iReg,watcherTemp,producerTemp);
93 oWatchers.push_back(watcherTemp);
95 oProds.push_back(producerTemp);
117 : m_generator(0), m_nonBeam(p.getParameter<bool>(
"NonBeamEvent")),
118 m_primaryTransformer(0),
119 m_managerInitialized(
false),
123 m_pUseMagneticField(p.getParameter<bool>(
"UseMagneticField")),
124 m_currentRun(0), m_currentEvent(0), m_simEvent(0),
125 m_PhysicsTablesDir(p.getParameter<std::string>(
"PhysicsTablesDirectory")),
126 m_StorePhysicsTables(p.getParameter<bool>(
"StorePhysicsTables")),
127 m_RestorePhysicsTables(p.getParameter<bool>(
"RestorePhysicsTables")),
128 m_EvtMgrVerbosity(p.getUntrackedParameter<int>(
"G4EventManagerVerbosity",0)),
129 m_Override(p.getParameter<bool>(
"OverrideUserStackingAction")),
130 m_pField(p.getParameter<edm::
ParameterSet>(
"MagneticField")),
131 m_pGenerator(p.getParameter<edm::
ParameterSet>(
"Generator")),
132 m_pPhysics(p.getParameter<edm::
ParameterSet>(
"Physics")),
133 m_pRunAction(p.getParameter<edm::
ParameterSet>(
"RunAction")),
134 m_pEventAction(p.getParameter<edm::
ParameterSet>(
"EventAction")),
135 m_pStackingAction(p.getParameter<edm::
ParameterSet>(
"StackingAction")),
136 m_pTrackingAction(p.getParameter<edm::
ParameterSet>(
"TrackingAction")),
137 m_pSteppingAction(p.getParameter<edm::
ParameterSet>(
"SteppingAction")),
138 m_G4Commands(p.getParameter<std::vector<std::string> >(
"G4Commands")),
139 m_p(p), m_fieldBuilder(0)
142 m_kernel = G4RunManagerKernel::GetRunManagerKernel();
173 <<
"[SimG4Core RunManager]\n"
174 <<
"The Geometry configuration is changed during the job execution\n"
175 <<
"this is not allowed, the geometry must stay unchanged\n";
181 <<
"[SimG4Core RunManager]\n"
182 <<
"The MagneticField configuration is changed during the job execution\n"
183 <<
"this is not allowed, the MagneticField must stay unchanged\n";
207 G4TransportationManager * tM = G4TransportationManager::GetTransportationManager();
217 std::pair< std::vector<SensitiveTkDetector*>,
225 edm::LogInfo(
"SimG4CoreApplication") <<
" RunManager: Sensitive Detector building finished; found " <<
m_sensTkDets.size()
226 <<
" Tk type Producers, and " <<
m_sensCaloDets.size() <<
" Calo type producers ";
237 std::auto_ptr<PhysicsListMakerBase> physicsMaker(
240 if (physicsMaker.get()==0)
throw SimG4Exception(
"Unable to find the Physics list requested");
251 else throw SimG4Exception(
"G4RunManagerKernel initialization failed!");
255 std::ostringstream
dir;
256 dir << tableDir <<
'\0';
257 std::string
cmd = std::string(
"/control/shell mkdir -p ")+tableDir;
259 G4UImanager::GetUIpointer()->ApplyCommand(cmd);
270 edm::LogInfo(
"SimG4CoreApplication") <<
"RunManager:: Requests UI: "
272 G4UImanager::GetUIpointer()->ApplyCommand(
m_G4Commands[it]);
301 <<
" RunManager::produce event " << inpevt.
id().
event()
302 <<
" with no G4PrimaryVertices \n Aborting Run" ;
324 G4Event *
e =
new G4Event(inpevt.
id().
event());
335 if (!HepMCEvt.isValid())
337 throw SimG4Exception(
"Unable to find HepMCProduct(HepMC::GenEvent) in edm::Event ");
365 m_kernel->GetEventManager()->GetTrackingManager()->GetTrack();
366 t->SetTrackStatus(fStopAndKill) ;
387 m_kernel->GetEventManager()->GetStackManager()->clear() ;
388 m_kernel->GetEventManager()->GetTrackingManager()->EventAborted() ;
390 G4StateManager* stateManager = G4StateManager::GetStateManager();
391 stateManager->SetNewState(G4State_GeomClosed);
405 G4EventManager * eventManager =
m_kernel->GetEventManager();
413 eventManager->SetUserAction(userEventAction);
417 eventManager->SetUserAction(userTrackingAction);
421 eventManager->SetUserAction(userSteppingAction);
424 edm::LogInfo(
"SimG4CoreApplication") <<
" RunManager: user StackingAction overridden " ;
430 edm::LogWarning(
"SimG4CoreApplication") <<
" RunManager: WARNING : No generator; initialized only RunAction!";
443 G4StateManager::GetStateManager()->SetNewState(G4State_GeomClosed);
T getParameter(std::string const &) const
EventNumber_t event() const
edm::ParameterSet m_pGenerator
T getUntrackedParameter(std::string const &, T const &) const
std::auto_ptr< SimTrackManager > m_trackManager
std::pair< std::vector< SensitiveTkDetector * >, std::vector< SensitiveCaloDetector * > > create(const DDDWorld &w, const DDCompactView &cpv, SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *m, SimActivityRegistry ®) const
virtual const math::XYZTLorentzVector * genVertex() const
edm::ESWatcher< IdealMagneticFieldRecord > idealMagRcdWatcher_
G4UserRunAction * m_userRunAction
virtual const double eventWeight() const
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< double > > XYZTLorentzVectorD
Lorentz vector with cylindrical internal representation using pseudorapidity.
void abortRun(bool softAbort=false)
G4StepSignal g4StepSignal_
SimActivityRegistry m_registry
virtual const HepMC::GenEvent * genEvent() const
SimActivityRegistry::G4StepSignal m_g4StepSignal
HepPDT::ParticleDataTable ParticleDataTable
void HepMC2G4(const HepMC::GenEvent *g, G4Event *e)
EndOfEventSignal endOfEventSignal_
void initializeUserActions()
SimActivityRegistry::EndOfRunSignal m_endOfRunSignal
void build(G4FieldManager *fM=0, G4PropagatorInField *fP=0)
BeginOfRunSignal beginOfRunSignal_
bool getByType(Handle< PROD > &result) const
SimActivityRegistry::EndOfEventSignal m_endOfEventSignal
edm::ESWatcher< IdealGeometryRecord > idealGeomRcdWatcher_
ExceptionHandler * m_CustomExceptionHandler
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e g
std::vector< SensitiveCaloDetector * > m_sensCaloDets
std::vector< SensitiveTkDetector * > m_sensTkDets
BeginOfJobSignal beginOfJobSignal_
U second(std::pair< T, U > const &p)
std::auto_ptr< PhysicsList > m_physicsList
edm::serviceregistry::AllArgsMaker< DTDataMonitorInterface, DTDataIntegrityTask > maker
void collisionPoint(math::XYZTLorentzVectorD v)
bool m_RestorePhysicsTables
void hepEvent(const HepMC::GenEvent *r)
void setGenEvent(const HepMC::GenEvent *inpevt)
BeginOfEventSignal beginOfEventSignal_
SimActivityRegistry::EndOfTrackSignal m_endOfTrackSignal
void resetGenParticleId(edm::Event &inpevt)
EndOfRunSignal endOfRunSignal_
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
edm::ParameterSet m_pStackingAction
unsigned int nTracks() const
void connect(Observer< const T * > *iObs)
does not take ownership of memory
edm::ParameterSet m_pEventAction
unsigned int nVertices() const
edm::ParameterSet m_pField
edm::ParameterSet m_pPhysics
edm::ParameterSet m_pSteppingAction
virtual void PostUserTrackingAction(const G4Track *aTrack)
unsigned int nGenParts() const
DDDWorldSignal dddWorldSignal_
T const * product() const
void nonBeamEvent2G4(const HepMC::GenEvent *g, G4Event *e)
bool check(const edm::EventSetup &iSetup)
G4Event * generateEvent(edm::Event &inpevt)
SimActivityRegistry::BeginOfRunSignal m_beginOfRunSignal
static void createWatchers(const edm::ParameterSet &iP, SimActivityRegistry &iReg, std::vector< boost::shared_ptr< SimWatcher > > &oWatchers, std::vector< boost::shared_ptr< SimProducer > > &oProds)
PrimaryTransformer * m_primaryTransformer
sim::FieldBuilder * m_fieldBuilder
std::string m_PhysicsTablesDir
bool m_StorePhysicsTables
void connect(SimActivityRegistry &iOther)
forwards our signals to slots connected to iOther
EndOfTrackSignal endOfTrackSignal_
std::vector< boost::shared_ptr< SimProducer > > m_producers
std::vector< boost::shared_ptr< SimWatcher > > m_watchers
edm::ParameterSet m_pRunAction
RunManager(edm::ParameterSet const &p)
edm::ParameterSet m_pTrackingAction
std::vector< std::string > m_G4Commands
SimActivityRegistry::BeginOfEventSignal m_beginOfEventSignal
bool m_managerInitialized
void initG4(const edm::EventSetup &es)
T get(const Candidate &c)
BeginOfTrackSignal beginOfTrackSignal_
G4RunManagerKernel * m_kernel
void produce(edm::Event &inpevt, const edm::EventSetup &es)
SimActivityRegistry::BeginOfTrackSignal m_beginOfTrackSignal