58 #include "HepPDT/ParticleDataTable.hh" 61 #include "G4GeometryManager.hh" 62 #include "G4StateManager.hh" 63 #include "G4ApplicationState.hh" 64 #include "G4RunManagerKernel.hh" 65 #include "G4UImanager.hh" 67 #include "G4EventManager.hh" 70 #include "G4TransportationManager.hh" 71 #include "G4ParticleTable.hh" 73 #include "G4FieldManager.hh" 75 #include "G4GDMLParser.hh" 76 #include "G4SystemOfUnits.hh" 88 std::vector<std::shared_ptr<SimWatcher> >& oWatchers,
89 std::vector<std::shared_ptr<SimProducer> >& oProds
95 vector<ParameterSet> watchers = iP.
getParameter<vector<ParameterSet> >(
"Watchers");
97 for(vector<ParameterSet>::iterator itWatcher = watchers.begin();
98 itWatcher != watchers.end();
100 std::shared_ptr<SimWatcherMakerBase> maker(
103 if(maker.get()==
nullptr) {
105 <<
"Unable to find the requested Watcher";
108 std::shared_ptr<SimWatcher> watcherTemp;
109 std::shared_ptr<SimProducer> producerTemp;
110 maker->make(*itWatcher,iReg,watcherTemp,producerTemp);
111 oWatchers.push_back(watcherTemp);
113 oProds.push_back(producerTemp);
122 m_nonBeam(p.getParameter<
bool>(
"NonBeamEvent")),
124 m_managerInitialized(
false),
127 m_pUseMagneticField(p.getParameter<
bool>(
"UseMagneticField")),
129 m_PhysicsTablesDir(p.getParameter<
std::
string>(
"PhysicsTablesDirectory")),
130 m_StorePhysicsTables(p.getParameter<
bool>(
"StorePhysicsTables")),
131 m_RestorePhysicsTables(p.getParameter<
bool>(
"RestorePhysicsTables")),
132 m_EvtMgrVerbosity(p.getUntrackedParameter<
int>(
"G4EventManagerVerbosity",0)),
138 m_pStackingAction(p.getParameter<
edm::
ParameterSet>(
"StackingAction")),
139 m_pTrackingAction(p.getParameter<
edm::
ParameterSet>(
"TrackingAction")),
140 m_pSteppingAction(p.getParameter<
edm::
ParameterSet>(
"SteppingAction")),
142 m_G4Commands(p.getParameter<
std::vector<
std::
string> >(
"G4Commands")),
146 m_kernel =
new G4RunManagerKernel();
147 G4StateManager::GetStateManager()->SetExceptionHandler(
new ExceptionHandler());
169 std::vector<edm::ParameterSet> watchers
170 = p.
getParameter<std::vector<edm::ParameterSet> >(
"Watchers");
181 G4StateManager::GetStateManager()->SetNewState(G4State_Quit);
182 G4GeometryManager::GetInstance()->OpenGeometry();
193 <<
"[SimG4Core RunManager]\n" 194 <<
"The Geometry configuration is changed during the job execution\n" 195 <<
"this is not allowed, the geometry must stay unchanged\n";
201 <<
"[SimG4Core RunManager]\n" 202 <<
"The MagneticField configuration is changed during the job execution\n" 203 <<
"this is not allowed, the MagneticField must stay unchanged\n";
227 G4TransportationManager * tM =
228 G4TransportationManager::GetTransportationManager();
229 tM->SetFieldManager(fieldManager);
230 fieldBuilder.build( fieldManager, tM->GetPropagatorInField());
243 std::pair< std::vector<SensitiveTkDetector*>,
244 std::vector<SensitiveCaloDetector*> > sensDets =
251 <<
" RunManager: Sensitive Detector " 252 <<
"building finished; found " 254 <<
" Tk type Producers, and " 256 <<
" Calo type producers ";
264 std::unique_ptr<PhysicsListMakerBase>
267 if (physicsMaker.get()==
nullptr) {
269 <<
"Unable to find the Physics list requested";
276 <<
"Physics list construction failed!";
281 if(monopoleMass > 0.0) {
297 <<
"RunManager: start initialisation of PhysicsList";
316 <<
"G4RunManagerKernel initialization failed!";
321 std::ostringstream
dir;
322 dir << tableDir <<
'\0';
325 G4UImanager::GetUIpointer()->ApplyCommand(cmd);
335 std::vector<int> ve =
m_p.
getParameter<std::vector<int> >(
"VerboseEvents");
336 std::vector<int> vn =
m_p.
getParameter<std::vector<int> >(
"VertexNumber");
337 std::vector<int> vt =
m_p.
getParameter<std::vector<int> >(
"VerboseTracks");
345 G4cout <<
"RunManager: Requested UI commands: " << G4endl;
348 G4UImanager::GetUIpointer()->ApplyCommand(
m_G4Commands[it]);
354 gdml.SetRegionExport(
true);
355 gdml.SetEnergyCutsExport(
true);
377 G4StateManager::GetStateManager()->SetNewState(G4State_Quit);
395 std::stringstream ss;
396 ss <<
" RunManager::produce(): event " << inpevt.
id().
event()
397 <<
" with no G4PrimaryVertices\n" ;
403 <<
"RunManager::produce: start Event " << inpevt.
id().
event()
412 <<
" RunManager::produce: ended Event " << inpevt.
id().
event();
425 G4Event * evt =
new G4Event(evtid);
454 m_kernel->GetEventManager()->GetTrackingManager()->GetTrack();
455 t->SetTrackStatus(fStopAndKill) ;
464 m_kernel->GetEventManager()->GetStackManager()->clear() ;
465 m_kernel->GetEventManager()->GetTrackingManager()->EventAborted() ;
467 G4StateManager* stateManager = G4StateManager::GetStateManager();
468 stateManager->SetNewState(G4State_GeomClosed);
478 G4EventManager * eventManager =
m_kernel->GetEventManager();
486 eventManager->SetUserAction(userEventAction);
491 eventManager->SetUserAction(userTrackingAction);
496 eventManager->SetUserAction(userSteppingAction);
498 eventManager->SetUserAction(
new StackingAction(userTrackingAction,
503 <<
"No generator; initialized " 504 <<
"only RunAction!";
512 G4StateManager::GetStateManager()->SetNewState(G4State_GeomClosed);
586 <<
" RunManager WARNING : " 587 <<
"error opening file <" <<
m_FieldFile <<
"> for magnetic field";
589 double rmax = 9000*mm;
590 double zmax = 16000*mm;
595 int nr = (
int)(rmax/dr);
596 int nz = 2*(
int)(zmax/dz);
603 double cosf =
cos(phi);
604 double sinf =
sin(phi);
606 double point[4] = {0.0,0.0,0.0,0.0};
607 double bfield[3] = {0.0,0.0,0.0};
609 fout << std::setprecision(6);
610 for(
int i=0;
i<=nr; ++
i) {
612 for(
int j=0; j<=nz; ++j) {
616 field->GetFieldValue(point, bfield);
617 fout <<
"R(mm)= " << r/mm <<
" phi(deg)= " << phi/degree
618 <<
" Z(mm)= " << z/mm <<
" Bz(tesla)= " << bfield[2]/tesla
619 <<
" Br(tesla)= " << (bfield[0]*cosf + bfield[1]*sinf)/tesla
620 <<
" Bphi(tesla)= " << (bfield[0]*sinf - bfield[1]*cosf)/tesla
T getParameter(std::string const &) const
EventNumber_t event() const
T getUntrackedParameter(std::string const &, T const &) const
edm::EDGetTokenT< edm::HepMCProduct > m_HepMC
virtual const math::XYZTLorentzVector * genVertex() const
SimTrackManager * GetSimTrackManager()
void BeginOfRunAction(const G4Run *aRun) override
edm::ESWatcher< IdealMagneticFieldRecord > idealMagRcdWatcher_
const G4VPhysicalVolume * GetWorldVolume() const
RunManager(edm::ParameterSet const &p, edm::ConsumesCollector &&i)
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)
def create(alignables, pedeDump, additionalData, outputFile, config)
G4StepSignal g4StepSignal_
edm::EDGetTokenT< edm::LHCTransportLinkContainer > m_LHCtr
SimActivityRegistry m_registry
virtual const HepMC::GenEvent * genEvent() const
SimActivityRegistry::G4StepSignal m_g4StepSignal
HepPDT::ParticleDataTable ParticleDataTable
bool getByToken(EDGetToken token, Handle< PROD > &result) const
std::vector< std::shared_ptr< SimProducer > > m_producers
void connect(Observer< const BeginOfJob * > *iObject)
void HepMC2G4(const HepMC::GenEvent *g, G4Event *e)
EndOfEventSignal endOfEventSignal_
Sin< T >::type sin(const T &t)
void initializeUserActions()
SimActivityRegistry::EndOfRunSignal m_endOfRunSignal
BeginOfRunSignal beginOfRunSignal_
SimActivityRegistry::EndOfEventSignal m_endOfEventSignal
static void createWatchers(const edm::ParameterSet &iP, SimActivityRegistry &iReg, std::vector< std::shared_ptr< SimWatcher > > &oWatchers, std::vector< std::shared_ptr< SimProducer > > &oProds)
edm::ESWatcher< IdealGeometryRecord > idealGeomRcdWatcher_
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)
RunAction
list of unwanted particles (gluons and quarks)
SimRunInterface * m_runInterface
Cos< T >::type cos(const T &t)
bool m_RestorePhysicsTables
void hepEvent(const HepMC::GenEvent *r)
std::vector< std::shared_ptr< SimWatcher > > m_watchers
std::unique_ptr< CMSSteppingVerbose > m_sVerbose
void setGenEvent(const HepMC::GenEvent *inpevt)
BeginOfEventSignal beginOfEventSignal_
SimActivityRegistry::EndOfTrackSignal m_endOfTrackSignal
void resetGenParticleId(edm::Event &inpevt)
EndOfRunSignal endOfRunSignal_
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
edm::ParameterSet m_g4overlap
const HepMC::GenEvent * GetEvent() const
std::pair< std::vector< SensitiveTkDetector * >, std::vector< SensitiveCaloDetector * > > create(const DDCompactView &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *, SimActivityRegistry ®) const
T const * product() const
unsigned int nGenParts() const
RunAction * m_userRunAction
DDDWorldSignal dddWorldSignal_
void nonBeamEvent2G4(const HepMC::GenEvent *g, G4Event *e)
bool check(const edm::EventSetup &iSetup)
G4Event * generateEvent(edm::Event &inpevt)
void ReportRegions(const std::string &ss)
void DumpMagneticField(const G4Field *) const
SimActivityRegistry::BeginOfRunSignal m_beginOfRunSignal
void Connect(RunAction *)
std::vector< LHCTransportLink > LHCTransportLinkContainer
PrimaryTransformer * m_primaryTransformer
std::string m_PhysicsTablesDir
void PostUserTrackingAction(const G4Track *aTrack) override
bool m_StorePhysicsTables
std::unique_ptr< CustomUIsession > m_UIsession
EndOfTrackSignal endOfTrackSignal_
std::unique_ptr< DDG4ProductionCuts > m_prodCuts
edm::ParameterSet m_pRunAction
std::unique_ptr< SimTrackManager > m_trackManager
T const * product() const
edm::ParameterSet m_pTrackingAction
void collisionPoint(const math::XYZTLorentzVectorD &v)
std::vector< std::string > m_G4Commands
SimActivityRegistry::BeginOfEventSignal m_beginOfEventSignal
std::unique_ptr< PhysicsList > m_physicsList
bool m_managerInitialized
void initG4(const edm::EventSetup &es)
T get(const Candidate &c)
void EndOfRunAction(const G4Run *aRun) override
*vegas h *****************************************************used in the default bin number in original ***version of VEGAS is ***a higher bin number might help to derive a more precise ***grade subtle point
BeginOfTrackSignal beginOfTrackSignal_
G4RunManagerKernel * m_kernel
void produce(edm::Event &inpevt, const edm::EventSetup &es)
SimActivityRegistry::BeginOfTrackSignal m_beginOfTrackSignal