CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
GlobalHitsProducer.cc
Go to the documentation of this file.
1 
8 #include "Validation/GlobalHits/interface/GlobalHitsProducer.h"
9 #include "Geometry/Records/interface/CaloGeometryRecord.h"
10 
11 GlobalHitsProducer::GlobalHitsProducer(const edm::ParameterSet& iPSet) :
12  fName(""), verbosity(0), frequency(0), vtxunit(0), label(""),
13  getAllProvenances(false), printProvenanceInfo(false), nRawGenPart(0),
14  G4VtxSrc_Token_( consumes<edm::SimVertexContainer>((iPSet.getParameter<edm::InputTag>("G4VtxSrc"))) ),
15  G4TrkSrc_Token_( consumes<edm::SimTrackContainer>(iPSet.getParameter<edm::InputTag>("G4TrkSrc")) ),
16  //ECalEBSrc_(""), ECalEESrc_(""), ECalESSrc_(""), HCalSrc_(""),
17  //PxlBrlLowSrc_(""), PxlBrlHighSrc_(""), PxlFwdLowSrc_(""),
18  //PxlFwdHighSrc_(""), SiTIBLowSrc_(""), SiTIBHighSrc_(""),
19  //SiTOBLowSrc_(""), SiTOBHighSrc_(""), SiTIDLowSrc_(""),
20  //SiTIDHighSrc_(""), SiTECLowSrc_(""), SiTECHighSrc_(""),
21  //MuonDtSrc_(""), MuonCscSrc_(""), MuonRpcSrc_(""),
22  count(0)
23 {
24  std::string MsgLoggerCat = "GlobalHitsProducer_GlobalHitsProducer";
25 
26  // get information from parameter set
27  fName = iPSet.getUntrackedParameter<std::string>("Name");
28  verbosity = iPSet.getUntrackedParameter<int>("Verbosity");
29  frequency = iPSet.getUntrackedParameter<int>("Frequency");
30  vtxunit = iPSet.getUntrackedParameter<int>("VtxUnit");
31  label = iPSet.getParameter<std::string>("Label");
32  edm::ParameterSet m_Prov =
33  iPSet.getParameter<edm::ParameterSet>("ProvenanceLookup");
34  getAllProvenances =
35  m_Prov.getUntrackedParameter<bool>("GetAllProvenances");
36  printProvenanceInfo =
37  m_Prov.getUntrackedParameter<bool>("PrintProvenanceInfo");
38 
39  //get Labels to use to extract information
40  PxlBrlLowSrc_ = iPSet.getParameter<edm::InputTag>("PxlBrlLowSrc");
41  PxlBrlHighSrc_ = iPSet.getParameter<edm::InputTag>("PxlBrlHighSrc");
42  PxlFwdLowSrc_ = iPSet.getParameter<edm::InputTag>("PxlFwdLowSrc");
43  PxlFwdHighSrc_ = iPSet.getParameter<edm::InputTag>("PxlFwdHighSrc");
44 
45  SiTIBLowSrc_ = iPSet.getParameter<edm::InputTag>("SiTIBLowSrc");
46  SiTIBHighSrc_ = iPSet.getParameter<edm::InputTag>("SiTIBHighSrc");
47  SiTOBLowSrc_ = iPSet.getParameter<edm::InputTag>("SiTOBLowSrc");
48  SiTOBHighSrc_ = iPSet.getParameter<edm::InputTag>("SiTOBHighSrc");
49  SiTIDLowSrc_ = iPSet.getParameter<edm::InputTag>("SiTIDLowSrc");
50  SiTIDHighSrc_ = iPSet.getParameter<edm::InputTag>("SiTIDHighSrc");
51  SiTECLowSrc_ = iPSet.getParameter<edm::InputTag>("SiTECLowSrc");
52  SiTECHighSrc_ = iPSet.getParameter<edm::InputTag>("SiTECHighSrc");
53 
54  MuonCscSrc_ = iPSet.getParameter<edm::InputTag>("MuonCscSrc");
55  MuonDtSrc_ = iPSet.getParameter<edm::InputTag>("MuonDtSrc");
56  MuonRpcSrc_ = iPSet.getParameter<edm::InputTag>("MuonRpcSrc");
57 
58  ECalEBSrc_ = iPSet.getParameter<edm::InputTag>("ECalEBSrc");
59  ECalEESrc_ = iPSet.getParameter<edm::InputTag>("ECalEESrc");
60  ECalESSrc_ = iPSet.getParameter<edm::InputTag>("ECalESSrc");
61 
62  HCalSrc_ = iPSet.getParameter<edm::InputTag>("HCalSrc");
63 
64  // fix for consumes
65  PxlBrlLowSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("PxlBrlLowSrc"));
66  PxlBrlHighSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("PxlBrlHighSrc"));
67  PxlFwdLowSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("PxlFwdLowSrc"));
68  PxlFwdHighSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("PxlFwdHighSrc"));
69 
70  SiTIBLowSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTIBLowSrc"));
71  SiTIBHighSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTIBHighSrc"));
72  SiTOBLowSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTOBLowSrc"));
73  SiTOBHighSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTOBHighSrc"));
74  SiTIDLowSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTIDLowSrc"));
75  SiTIDHighSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTIDHighSrc"));
76  SiTECLowSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTECLowSrc"));
77  SiTECHighSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("SiTECHighSrc"));
78 
79  MuonCscSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("MuonCscSrc"));
80  MuonDtSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("MuonDtSrc"));
81  MuonRpcSrc_Token_ = consumes<edm::PSimHitContainer>(iPSet.getParameter<edm::InputTag>("MuonRpcSrc"));
82 
83  ECalEBSrc_Token_ = consumes<edm::PCaloHitContainer>(iPSet.getParameter<edm::InputTag>("ECalEBSrc"));
84  ECalEESrc_Token_ = consumes<edm::PCaloHitContainer>(iPSet.getParameter<edm::InputTag>("ECalEESrc"));
85  ECalESSrc_Token_ = consumes<edm::PCaloHitContainer>(iPSet.getParameter<edm::InputTag>("ECalESSrc"));
86  HCalSrc_Token_ = consumes<edm::PCaloHitContainer>(iPSet.getParameter<edm::InputTag>("HCalSrc"));
87 
88  // use value of first digit to determine default output level (inclusive)
89  // 0 is none, 1 is basic, 2 is fill output, 3 is gather output
90  verbosity %= 10;
91 
92  // create persistent object
93  produces<PGlobalSimHit>(label);
94 
95  // print out Parameter Set information being used
96  if (verbosity >= 0) {
97  edm::LogInfo(MsgLoggerCat)
98  << "\n===============================\n"
99  << "Initialized as EDProducer with parameter values:\n"
100  << " Name = " << fName << "\n"
101  << " Verbosity = " << verbosity << "\n"
102  << " Frequency = " << frequency << "\n"
103  << " VtxUnit = " << vtxunit << "\n"
104  << " Label = " << label << "\n"
105  << " GetProv = " << getAllProvenances << "\n"
106  << " PrintProv = " << printProvenanceInfo << "\n"
107  << " PxlBrlLowSrc = " << PxlBrlLowSrc_.label()
108  << ":" << PxlBrlLowSrc_.instance() << "\n"
109  << " PxlBrlHighSrc = " << PxlBrlHighSrc_.label()
110  << ":" << PxlBrlHighSrc_.instance() << "\n"
111  << " PxlFwdLowSrc = " << PxlFwdLowSrc_.label()
112  << ":" << PxlBrlLowSrc_.instance() << "\n"
113  << " PxlFwdHighSrc = " << PxlFwdHighSrc_.label()
114  << ":" << PxlBrlHighSrc_.instance() << "\n"
115  << " SiTIBLowSrc = " << SiTIBLowSrc_.label()
116  << ":" << SiTIBLowSrc_.instance() << "\n"
117  << " SiTIBHighSrc = " << SiTIBHighSrc_.label()
118  << ":" << SiTIBHighSrc_.instance() << "\n"
119  << " SiTOBLowSrc = " << SiTOBLowSrc_.label()
120  << ":" << SiTOBLowSrc_.instance() << "\n"
121  << " SiTOBHighSrc = " << SiTOBHighSrc_.label()
122  << ":" << SiTOBHighSrc_.instance() << "\n"
123  << " SiTIDLowSrc = " << SiTIDLowSrc_.label()
124  << ":" << SiTIDLowSrc_.instance() << "\n"
125  << " SiTIDHighSrc = " << SiTIDHighSrc_.label()
126  << ":" << SiTIDHighSrc_.instance() << "\n"
127  << " SiTECLowSrc = " << SiTECLowSrc_.label()
128  << ":" << SiTECLowSrc_.instance() << "\n"
129  << " SiTECHighSrc = " << SiTECHighSrc_.label()
130  << ":" << SiTECHighSrc_.instance() << "\n"
131  << " MuonCscSrc = " << MuonCscSrc_.label()
132  << ":" << MuonCscSrc_.instance() << "\n"
133  << " MuonDtSrc = " << MuonDtSrc_.label()
134  << ":" << MuonDtSrc_.instance() << "\n"
135  << " MuonRpcSrc = " << MuonRpcSrc_.label()
136  << ":" << MuonRpcSrc_.instance() << "\n"
137  << " ECalEBSrc = " << ECalEBSrc_.label()
138  << ":" << ECalEBSrc_.instance() << "\n"
139  << " ECalEESrc = " << ECalEESrc_.label()
140  << ":" << ECalEESrc_.instance() << "\n"
141  << " ECalESSrc = " << ECalESSrc_.label()
142  << ":" << ECalESSrc_.instance() << "\n"
143  << " HCalSrc = " << HCalSrc_.label()
144  << ":" << HCalSrc_.instance() << "\n"
145  << "===============================\n";
146  }
147 
148  // migrated here from beginJob
149  clear();
150 
151 }
152 
153 GlobalHitsProducer::~GlobalHitsProducer()
154 {
155 }
156 
157 void GlobalHitsProducer::beginJob( void )
158 {
159  return;
160 }
161 
162 void GlobalHitsProducer::endJob()
163 {
164  std::string MsgLoggerCat = "GlobalHitsProducer_endJob";
165  if (verbosity >= 0)
166  edm::LogInfo(MsgLoggerCat)
167  << "Terminating having processed " << count << " events.";
168  return;
169 }
170 
171 void GlobalHitsProducer::produce(edm::Event& iEvent,
172  const edm::EventSetup& iSetup)
173 {
174  std::string MsgLoggerCat = "GlobalHitsProducer_produce";
175 
176  // keep track of number of events processed
177  ++count;
178 
179  // get event id information
180  int nrun = iEvent.id().run();
181  int nevt = iEvent.id().event();
182 
183  if (verbosity > 0) {
184  edm::LogInfo(MsgLoggerCat)
185  << "Processing run " << nrun << ", event " << nevt
186  << " (" << count << " events total)";
187  } else if (verbosity == 0) {
188  if (nevt%frequency == 0 || nevt == 1) {
189  edm::LogInfo(MsgLoggerCat)
190  << "Processing run " << nrun << ", event " << nevt
191  << " (" << count << " events total)";
192  }
193  }
194 
195  // clear event holders
196  clear();
197 
198  // look at information available in the event
199  if (getAllProvenances) {
200 
201  std::vector<const edm::Provenance*> AllProv;
202  iEvent.getAllProvenance(AllProv);
203 
204  if (verbosity >= 0)
205  edm::LogInfo(MsgLoggerCat)
206  << "Number of Provenances = " << AllProv.size();
207 
208  if (printProvenanceInfo && (verbosity >= 0)) {
209  TString eventout("\nProvenance info:\n");
210 
211  for (unsigned int i = 0; i < AllProv.size(); ++i) {
212  eventout += "\n ******************************";
213  eventout += "\n Module : ";
214  eventout += AllProv[i]->moduleLabel();
215  eventout += "\n ProductID : ";
216  eventout += AllProv[i]->productID().id();
217  eventout += "\n ClassName : ";
218  eventout += AllProv[i]->className();
219  eventout += "\n InstanceName : ";
220  eventout += AllProv[i]->productInstanceName();
221  eventout += "\n BranchName : ";
222  eventout += AllProv[i]->branchName();
223  }
224  eventout += "\n ******************************\n";
225  edm::LogInfo(MsgLoggerCat) << eventout << "\n";
226  printProvenanceInfo = false;
227  }
228  getAllProvenances = false;
229  }
230 
231  // call fill functions
232  //gather G4MC information from event
233  fillG4MC(iEvent);
234  // gather Tracker information from event
235  fillTrk(iEvent,iSetup);
236  // gather muon information from event
237  fillMuon(iEvent, iSetup);
238  // gather Ecal information from event
239  fillECal(iEvent, iSetup);
240  // gather Hcal information from event
241  fillHCal(iEvent, iSetup);
242 
243  if (verbosity > 0)
244  edm::LogInfo (MsgLoggerCat)
245  << "Done gathering data from event.";
246 
247  // produce object to put into event
248  std::auto_ptr<PGlobalSimHit> pOut(new PGlobalSimHit);
249 
250  if (verbosity > 2)
251  edm::LogInfo (MsgLoggerCat)
252  << "Saving event contents:";
253 
254  // call store functions
255  // store G4MC information in product
256  storeG4MC(*pOut);
257  // store Tracker information in produce
258  storeTrk(*pOut);
259  // store Muon information in produce
260  storeMuon(*pOut);
261  // store ECal information in produce
262  storeECal(*pOut);
263  // store HCal information in produce
264  storeHCal(*pOut);
265 
266  // store information in event
267  iEvent.put(pOut,label);
268 
269  return;
270 }
271 
272 //==================fill and store functions================================
273 void GlobalHitsProducer::fillG4MC(edm::Event& iEvent)
274 {
275 
276  std::string MsgLoggerCat = "GlobalHitsProducer_fillG4MC";
277 
278  TString eventout;
279  if (verbosity > 0)
280  eventout = "\nGathering info:";
281 
283  // get MC information
285  edm::Handle<edm::HepMCProduct> HepMCEvt;
286  std::vector<edm::Handle<edm::HepMCProduct> > AllHepMCEvt;
287  iEvent.getManyByType(AllHepMCEvt);
288 
289  // loop through products and extract VtxSmearing if available. Any of them
290  // should have the information needed
291  for (unsigned int i = 0; i < AllHepMCEvt.size(); ++i) {
292  HepMCEvt = AllHepMCEvt[i];
293  if ((HepMCEvt.provenance()->product()).moduleLabel() == "VtxSmeared")
294  break;
295  }
296 
297  if (!HepMCEvt.isValid()) {
298  edm::LogWarning(MsgLoggerCat)
299  << "Unable to find HepMCProduct in event!";
300  return;
301  } else {
302  eventout += "\n Using HepMCProduct: ";
303  eventout += (HepMCEvt.provenance()->product()).moduleLabel();
304  }
305  const HepMC::GenEvent* MCEvt = HepMCEvt->GetEvent();
306  nRawGenPart = MCEvt->particles_size();
307 
308  if (verbosity > 1) {
309  eventout += "\n Number of Raw Particles collected:......... ";
310  eventout += nRawGenPart;
311  }
312 
314  // get G4Vertex information
316  // convert unit stored in SimVertex to mm
317  float unit = 0.;
318  if (vtxunit == 0) unit = 1.; // already in mm
319  if (vtxunit == 1) unit = 10.; // stored in cm, convert to mm
320 
321  edm::Handle<edm::SimVertexContainer> G4VtxContainer;
322  iEvent.getByToken(G4VtxSrc_Token_, G4VtxContainer);
323  if (!G4VtxContainer.isValid()) {
324  edm::LogWarning(MsgLoggerCat)
325  << "Unable to find SimVertex in event!";
326  return;
327  }
328  int i = 0;
329  edm::SimVertexContainer::const_iterator itVtx;
330  for (itVtx = G4VtxContainer->begin(); itVtx != G4VtxContainer->end();
331  ++itVtx) {
332 
333  ++i;
334 
335  const math::XYZTLorentzVector G4Vtx1(itVtx->position().x(),
336  itVtx->position().y(),
337  itVtx->position().z(),
338  itVtx->position().e());
339  double G4Vtx[4];
340  G4Vtx1.GetCoordinates(G4Vtx);
341 
342  G4VtxX.push_back((G4Vtx[0]*unit)/micrometer);
343  G4VtxY.push_back((G4Vtx[1]*unit)/micrometer);
344  G4VtxZ.push_back((G4Vtx[2]*unit)/millimeter);
345  }
346 
347  if (verbosity > 1) {
348  eventout += "\n Number of G4Vertices collected:............ ";
349  eventout += i;
350  }
351 
353  // get G4Track information
355  edm::Handle<edm::SimTrackContainer> G4TrkContainer;
356  iEvent.getByToken(G4TrkSrc_Token_, G4TrkContainer);
357 
358  if (!G4TrkContainer.isValid()) {
359  edm::LogWarning(MsgLoggerCat)
360  << "Unable to find SimTrack in event!";
361  return;
362  }
363  i = 0;
364  edm::SimTrackContainer::const_iterator itTrk;
365  for (itTrk = G4TrkContainer->begin(); itTrk != G4TrkContainer->end();
366  ++itTrk) {
367 
368  ++i;
369 
370  const math::XYZTLorentzVector G4Trk1(itTrk->momentum().x(),
371  itTrk->momentum().y(),
372  itTrk->momentum().z(),
373  itTrk->momentum().e());
374  double G4Trk[4];
375  G4Trk1.GetCoordinates(G4Trk);
376 
377  G4TrkPt.push_back(sqrt(G4Trk[0]*G4Trk[0]+G4Trk[1]*G4Trk[1])); //GeV
378  G4TrkE.push_back(G4Trk[3]); //GeV
379  }
380 
381  if (verbosity > 1) {
382  eventout += "\n Number of G4Tracks collected:.............. ";
383  eventout += i;
384  }
385 
386  if (verbosity > 0)
387  edm::LogInfo(MsgLoggerCat) << eventout << "\n";
388 
389  return;
390 }
391 
392 void GlobalHitsProducer::storeG4MC(PGlobalSimHit& product)
393 {
394  std::string MsgLoggerCat = "GlobalHitsProducer_storeG4MC";
395 
396  if (verbosity > 2) {
397  TString eventout("\n nRawGenPart = ");
398  eventout += nRawGenPart;
399  eventout += "\n nG4Vtx = ";
400  eventout += G4VtxX.size();
401  for (unsigned int i = 0; i < G4VtxX.size(); ++i) {
402  eventout += "\n (x,y,z) = (";
403  eventout += G4VtxX[i];
404  eventout += ", ";
405  eventout += G4VtxY[i];
406  eventout += ", ";
407  eventout += G4VtxZ[i];
408  eventout += ")";
409  }
410  eventout += "\n nG4Trk = ";
411  eventout += G4TrkPt.size();
412  for (unsigned int i = 0; i < G4TrkPt.size(); ++i) {
413  eventout += "\n (pt,e) = (";
414  eventout += G4TrkPt[i];
415  eventout += ", ";
416  eventout += G4TrkE[i];
417  eventout += ")";
418  }
419  edm::LogInfo(MsgLoggerCat) << eventout << "\n";
420  } // end verbose output
421 
422  product.putRawGenPart(nRawGenPart);
423  product.putG4Vtx(G4VtxX, G4VtxY, G4VtxZ);
424  product.putG4Trk(G4TrkPt, G4TrkE);
425 
426  return;
427 }
428 
429 void GlobalHitsProducer::fillTrk(edm::Event& iEvent,
430  const edm::EventSetup& iSetup)
431 {
432  std::string MsgLoggerCat = "GlobalHitsProducer_fillTrk";
433 
434  TString eventout;
435  if (verbosity > 0)
436  eventout = "\nGathering info:";
437 
438  // access the tracker geometry
439  edm::ESHandle<TrackerGeometry> theTrackerGeometry;
440  iSetup.get<TrackerDigiGeometryRecord>().get(theTrackerGeometry);
441  if (!theTrackerGeometry.isValid()) {
442  edm::LogWarning(MsgLoggerCat)
443  << "Unable to find TrackerDigiGeometryRecord in event!";
444  return;
445  }
446  const TrackerGeometry& theTracker(*theTrackerGeometry);
447 
448  // iterator to access containers
449  edm::PSimHitContainer::const_iterator itHit;
450 
452  // get Pixel Barrel information
454  edm::PSimHitContainer thePxlBrlHits;
455  // extract low container
456  edm::Handle<edm::PSimHitContainer> PxlBrlLowContainer;
457  iEvent.getByToken(PxlBrlLowSrc_Token_,PxlBrlLowContainer);
458  if (!PxlBrlLowContainer.isValid()) {
459  edm::LogWarning(MsgLoggerCat)
460  << "Unable to find TrackerHitsPixelBarrelLowTof in event!";
461  return;
462  }
463  // extract high container
464  edm::Handle<edm::PSimHitContainer> PxlBrlHighContainer;
465  iEvent.getByToken(PxlBrlHighSrc_Token_,PxlBrlHighContainer);
466  if (!PxlBrlHighContainer.isValid()) {
467  edm::LogWarning(MsgLoggerCat)
468  << "Unable to find TrackerHitsPixelBarrelHighTof in event!";
469  return;
470  }
471  // place both containers into new container
472  thePxlBrlHits.insert(thePxlBrlHits.end(),PxlBrlLowContainer->begin(),
473  PxlBrlLowContainer->end());
474  thePxlBrlHits.insert(thePxlBrlHits.end(),PxlBrlHighContainer->begin(),
475  PxlBrlHighContainer->end());
476 
477  // cycle through new container
478  int i = 0, j = 0;
479  for (itHit = thePxlBrlHits.begin(); itHit != thePxlBrlHits.end(); ++itHit) {
480 
481  ++i;
482 
483  // create a DetId from the detUnitId
484  DetId theDetUnitId(itHit->detUnitId());
485  int detector = theDetUnitId.det();
486  int subdetector = theDetUnitId.subdetId();
487 
488  // check that expected detector is returned
489  if ((detector == dTrk) && (subdetector == sdPxlBrl)) {
490 
491  // get the GeomDetUnit from the geometry using theDetUnitID
492  const GeomDetUnit *theDet = theTracker.idToDetUnit(theDetUnitId);
493 
494  if (!theDet) {
495  edm::LogWarning(MsgLoggerCat)
496  << "Unable to get GeomDetUnit from PxlBrlHits for Hit " << i;
497  continue;
498  }
499 
500  ++j;
501 
502  // get the Surface of the hit (knows how to go from local <-> global)
503  const BoundPlane& bSurface = theDet->surface();
504 
505  // gather necessary information
506  PxlBrlToF.push_back(itHit->tof());
507  PxlBrlR.push_back(bSurface.toGlobal(itHit->localPosition()).perp());
508  PxlBrlPhi.push_back(bSurface.toGlobal(itHit->localPosition()).phi());
509  PxlBrlEta.push_back(bSurface.toGlobal(itHit->localPosition()).eta());
510 
511  } else {
512  edm::LogWarning(MsgLoggerCat)
513  << "PxlBrl PSimHit " << i
514  << " is expected to be (det,subdet) = ("
515  << dTrk << "," << sdPxlBrl
516  << "); value returned is: ("
517  << detector << "," << subdetector << ")";
518  continue;
519  } // end detector type check
520  } // end loop through PxlBrl Hits
521 
522  if (verbosity > 1) {
523  eventout += "\n Number of Pixel Barrel Hits collected:..... ";
524  eventout += j;
525  }
526 
528  // get Pixel Forward information
530  edm::PSimHitContainer thePxlFwdHits;
531  // extract low container
532  edm::Handle<edm::PSimHitContainer> PxlFwdLowContainer;
533  iEvent.getByToken(PxlFwdLowSrc_Token_,PxlFwdLowContainer);
534  if (!PxlFwdLowContainer.isValid()) {
535  edm::LogWarning(MsgLoggerCat)
536  << "Unable to find TrackerHitsPixelEndcapLowTof in event!";
537  return;
538  }
539  // extract high container
540  edm::Handle<edm::PSimHitContainer> PxlFwdHighContainer;
541  iEvent.getByToken(PxlFwdHighSrc_Token_,PxlFwdHighContainer);
542  if (!PxlFwdHighContainer.isValid()) {
543  edm::LogWarning("GlobalHitsProducer_fillTrk")
544  << "Unable to find TrackerHitsPixelEndcapHighTof in event!";
545  return;
546  }
547  // place both containers into new container
548  thePxlFwdHits.insert(thePxlFwdHits.end(),PxlFwdLowContainer->begin(),
549  PxlFwdLowContainer->end());
550  thePxlFwdHits.insert(thePxlFwdHits.end(),PxlFwdHighContainer->begin(),
551  PxlFwdHighContainer->end());
552 
553  // cycle through new container
554  i = 0; j = 0;
555  for (itHit = thePxlFwdHits.begin(); itHit != thePxlFwdHits.end(); ++itHit) {
556 
557  ++i;
558 
559  // create a DetId from the detUnitId
560  DetId theDetUnitId(itHit->detUnitId());
561  int detector = theDetUnitId.det();
562  int subdetector = theDetUnitId.subdetId();
563 
564  // check that expected detector is returned
565  if ((detector == dTrk) && (subdetector == sdPxlFwd)) {
566 
567  // get the GeomDetUnit from the geometry using theDetUnitID
568  const GeomDetUnit *theDet = theTracker.idToDetUnit(theDetUnitId);
569 
570  if (!theDet) {
571  edm::LogWarning(MsgLoggerCat)
572  << "Unable to get GeomDetUnit from PxlFwdHits for Hit " << i;;
573  continue;
574  }
575 
576  ++j;
577 
578  // get the Surface of the hit (knows how to go from local <-> global)
579  const BoundPlane& bSurface = theDet->surface();
580 
581  // gather necessary information
582  PxlFwdToF.push_back(itHit->tof());
583  PxlFwdZ.push_back(bSurface.toGlobal(itHit->localPosition()).z());
584  PxlFwdPhi.push_back(bSurface.toGlobal(itHit->localPosition()).phi());
585  PxlFwdEta.push_back(bSurface.toGlobal(itHit->localPosition()).eta());
586  } else {
587  edm::LogWarning(MsgLoggerCat)
588  << "PxlFwd PSimHit " << i
589  << " is expected to be (det,subdet) = ("
590  << dTrk << "," << sdPxlFwd
591  << "); value returned is: ("
592  << detector << "," << subdetector << ")";
593  continue;
594  } // end detector type check
595  } // end loop through PxlFwd Hits
596 
597  if (verbosity > 1) {
598  eventout += "\n Number of Pixel Forward Hits collected:.... ";
599  eventout += j;
600  }
601 
603  // get Silicon Barrel information
605  edm::PSimHitContainer theSiBrlHits;
606  // extract TIB low container
607  edm::Handle<edm::PSimHitContainer> SiTIBLowContainer;
608  iEvent.getByToken(SiTIBLowSrc_Token_,SiTIBLowContainer);
609  if (!SiTIBLowContainer.isValid()) {
610  edm::LogWarning(MsgLoggerCat)
611  << "Unable to find TrackerHitsTIBLowTof in event!";
612  return;
613  }
614  // extract TIB high container
615  edm::Handle<edm::PSimHitContainer> SiTIBHighContainer;
616  iEvent.getByToken(SiTIBHighSrc_Token_,SiTIBHighContainer);
617  if (!SiTIBHighContainer.isValid()) {
618  edm::LogWarning(MsgLoggerCat)
619  << "Unable to find TrackerHitsTIBHighTof in event!";
620  return;
621  }
622  // extract TOB low container
623  edm::Handle<edm::PSimHitContainer> SiTOBLowContainer;
624  iEvent.getByToken(SiTOBLowSrc_Token_,SiTOBLowContainer);
625  if (!SiTOBLowContainer.isValid()) {
626  edm::LogWarning(MsgLoggerCat)
627  << "Unable to find TrackerHitsTOBLowTof in event!";
628  return;
629  }
630  // extract TOB high container
631  edm::Handle<edm::PSimHitContainer> SiTOBHighContainer;
632  iEvent.getByToken(SiTOBHighSrc_Token_,SiTOBHighContainer);
633  if (!SiTOBHighContainer.isValid()) {
634  edm::LogWarning(MsgLoggerCat)
635  << "Unable to find TrackerHitsTOBHighTof in event!";
636  return;
637  }
638  // place all containers into new container
639  theSiBrlHits.insert(theSiBrlHits.end(),SiTIBLowContainer->begin(),
640  SiTIBLowContainer->end());
641  theSiBrlHits.insert(theSiBrlHits.end(),SiTIBHighContainer->begin(),
642  SiTIBHighContainer->end());
643  theSiBrlHits.insert(theSiBrlHits.end(),SiTOBLowContainer->begin(),
644  SiTOBLowContainer->end());
645  theSiBrlHits.insert(theSiBrlHits.end(),SiTOBHighContainer->begin(),
646  SiTOBHighContainer->end());
647 
648  // cycle through new container
649  i = 0; j = 0;
650  for (itHit = theSiBrlHits.begin(); itHit != theSiBrlHits.end(); ++itHit) {
651 
652  ++i;
653 
654  // create a DetId from the detUnitId
655  DetId theDetUnitId(itHit->detUnitId());
656  int detector = theDetUnitId.det();
657  int subdetector = theDetUnitId.subdetId();
658 
659  // check that expected detector is returned
660  if ((detector == dTrk) &&
661  ((subdetector == sdSiTIB) ||
662  (subdetector == sdSiTOB))) {
663 
664  // get the GeomDetUnit from the geometry using theDetUnitID
665  const GeomDetUnit *theDet = theTracker.idToDetUnit(theDetUnitId);
666 
667  if (!theDet) {
668  edm::LogWarning(MsgLoggerCat)
669  << "Unable to get GeomDetUnit from SiBrlHits for Hit " << i;
670  continue;
671  }
672 
673  ++j;
674 
675  // get the Surface of the hit (knows how to go from local <-> global)
676  const BoundPlane& bSurface = theDet->surface();
677 
678  // gather necessary information
679  SiBrlToF.push_back(itHit->tof());
680  SiBrlR.push_back(bSurface.toGlobal(itHit->localPosition()).perp());
681  SiBrlPhi.push_back(bSurface.toGlobal(itHit->localPosition()).phi());
682  SiBrlEta.push_back(bSurface.toGlobal(itHit->localPosition()).eta());
683  } else {
684  edm::LogWarning(MsgLoggerCat)
685  << "SiBrl PSimHit " << i
686  << " is expected to be (det,subdet) = ("
687  << dTrk << "," << sdSiTIB
688  << " || " << sdSiTOB << "); value returned is: ("
689  << detector << "," << subdetector << ")";
690  continue;
691  } // end detector type check
692  } // end loop through SiBrl Hits
693 
694  if (verbosity > 1) {
695  eventout += "\n Number of Silicon Barrel Hits collected:... ";
696  eventout += j;
697  }
698 
700  // get Silicon Forward information
702  edm::PSimHitContainer theSiFwdHits;
703  // extract TID low container
704  edm::Handle<edm::PSimHitContainer> SiTIDLowContainer;
705  iEvent.getByToken(SiTIDLowSrc_Token_,SiTIDLowContainer);
706  if (!SiTIDLowContainer.isValid()) {
707  edm::LogWarning(MsgLoggerCat)
708  << "Unable to find TrackerHitsTIDLowTof in event!";
709  return;
710  }
711  // extract TID high container
712  edm::Handle<edm::PSimHitContainer> SiTIDHighContainer;
713  iEvent.getByToken(SiTIDHighSrc_Token_,SiTIDHighContainer);
714  if (!SiTIDHighContainer.isValid()) {
715  edm::LogWarning("GlobalHitsProducer_fillTrk")
716  << "Unable to find TrackerHitsTIDHighTof in event!";
717  return;
718  }
719  // extract TEC low container
720  edm::Handle<edm::PSimHitContainer> SiTECLowContainer;
721  iEvent.getByToken(SiTECLowSrc_Token_,SiTECLowContainer);
722  if (!SiTECLowContainer.isValid()) {
723  edm::LogWarning(MsgLoggerCat)
724  << "Unable to find TrackerHitsTECLowTof in event!";
725  return;
726  }
727  // extract TEC high container
728  edm::Handle<edm::PSimHitContainer> SiTECHighContainer;
729  iEvent.getByToken(SiTECHighSrc_Token_,SiTECHighContainer);
730  if (!SiTECHighContainer.isValid()) {
731  edm::LogWarning(MsgLoggerCat)
732  << "Unable to find TrackerHitsTECHighTof in event!";
733  return;
734  }
735  // place all containers into new container
736  theSiFwdHits.insert(theSiFwdHits.end(),SiTIDLowContainer->begin(),
737  SiTIDLowContainer->end());
738  theSiFwdHits.insert(theSiFwdHits.end(),SiTIDHighContainer->begin(),
739  SiTIDHighContainer->end());
740  theSiFwdHits.insert(theSiFwdHits.end(),SiTECLowContainer->begin(),
741  SiTECLowContainer->end());
742  theSiFwdHits.insert(theSiFwdHits.end(),SiTECHighContainer->begin(),
743  SiTECHighContainer->end());
744 
745  // cycle through container
746  i = 0; j = 0;
747  for (itHit = theSiFwdHits.begin(); itHit != theSiFwdHits.end(); ++itHit) {
748 
749  ++i;
750 
751  // create a DetId from the detUnitId
752  DetId theDetUnitId(itHit->detUnitId());
753  int detector = theDetUnitId.det();
754  int subdetector = theDetUnitId.subdetId();
755 
756  // check that expected detector is returned
757  if ((detector == dTrk) &&
758  ((subdetector == sdSiTID) ||
759  (subdetector == sdSiTEC))) {
760 
761  // get the GeomDetUnit from the geometry using theDetUnitID
762  const GeomDetUnit *theDet = theTracker.idToDetUnit(theDetUnitId);
763 
764  if (!theDet) {
765  edm::LogWarning(MsgLoggerCat)
766  << "Unable to get GeomDetUnit from SiFwdHits Hit " << i;
767  return;
768  }
769 
770  ++j;
771 
772  // get the Surface of the hit (knows how to go from local <-> global)
773  const BoundPlane& bSurface = theDet->surface();
774 
775  // gather necessary information
776  SiFwdToF.push_back(itHit->tof());
777  SiFwdZ.push_back(bSurface.toGlobal(itHit->localPosition()).z());
778  SiFwdPhi.push_back(bSurface.toGlobal(itHit->localPosition()).phi());
779  SiFwdEta.push_back(bSurface.toGlobal(itHit->localPosition()).eta());
780  } else {
781  edm::LogWarning(MsgLoggerCat)
782  << "SiFwd PSimHit " << i
783  << " is expected to be (det,subdet) = ("
784  << dTrk << "," << sdSiTOB
785  << " || " << sdSiTEC << "); value returned is: ("
786  << detector << "," << subdetector << ")";
787  continue;
788  } // end check detector type
789  } // end loop through SiFwd Hits
790 
791  if (verbosity > 1) {
792  eventout += "\n Number of Silicon Forward Hits collected:.. ";
793  eventout += j;
794  }
795 
796  if (verbosity > 0)
797  edm::LogInfo(MsgLoggerCat) << eventout << "\n";
798 
799  return;
800  }
801 
802 void GlobalHitsProducer::storeTrk(PGlobalSimHit& product)
803 {
804  std::string MsgLoggerCat = "GlobalHitsProducer_storeTrk";
805 
806  if (verbosity > 2) {
807  TString eventout("\n nPxlBrlHits = ");
808  eventout += PxlBrlToF.size();
809  for (unsigned int i = 0; i < PxlBrlToF.size(); ++i) {
810  eventout += "\n (tof,r,phi,eta) = (";
811  eventout += PxlBrlToF[i];
812  eventout += ", ";
813  eventout += PxlBrlR[i];
814  eventout += ", ";
815  eventout += PxlBrlPhi[i];
816  eventout += ", ";
817  eventout += PxlBrlEta[i];
818  eventout += ")";
819  } // end PxlBrl output
820  eventout += "\n nPxlFwdHits = ";
821  eventout += PxlFwdToF.size();
822  for (unsigned int i = 0; i < PxlFwdToF.size(); ++i) {
823  eventout += "\n (tof,z,phi,eta) = (";
824  eventout += PxlFwdToF[i];
825  eventout += ", ";
826  eventout += PxlFwdZ[i];
827  eventout += ", ";
828  eventout += PxlFwdPhi[i];
829  eventout += ", ";
830  eventout += PxlFwdEta[i];
831  eventout += ")";
832  } // end PxlFwd output
833  eventout += "\n nSiBrlHits = ";
834  eventout += SiBrlToF.size();
835  for (unsigned int i = 0; i < SiBrlToF.size(); ++i) {
836  eventout += "\n (tof,r,phi,eta) = (";
837  eventout += SiBrlToF[i];
838  eventout += ", ";
839  eventout += SiBrlR[i];
840  eventout += ", ";
841  eventout += SiBrlPhi[i];
842  eventout += ", ";
843  eventout += SiBrlEta[i];
844  eventout += ")";
845  } // end SiBrl output
846  eventout += "\n nSiFwdHits = ";
847  eventout += SiFwdToF.size();
848  for (unsigned int i = 0; i < SiFwdToF.size(); ++i) {
849  eventout += "\n (tof,z,phi,eta) = (";
850  eventout += SiFwdToF[i];
851  eventout += ", ";
852  eventout += SiFwdZ[i];
853  eventout += ", ";
854  eventout += SiFwdPhi[i];
855  eventout += ", ";
856  eventout += SiFwdEta[i];
857  eventout += ")";
858  } // end SiFwd output
859  edm::LogInfo(MsgLoggerCat) << eventout << "\n";
860  } // end verbose output
861 
862  product.putPxlBrlHits(PxlBrlToF,PxlBrlR,PxlBrlPhi,PxlBrlEta);
863  product.putPxlFwdHits(PxlFwdToF,PxlFwdZ,PxlFwdPhi,PxlFwdEta);
864  product.putSiBrlHits(SiBrlToF,SiBrlR,SiBrlPhi,SiBrlEta);
865  product.putSiFwdHits(SiFwdToF,SiFwdZ,SiFwdPhi,SiFwdEta);
866 
867  return;
868 }
869 
870 void GlobalHitsProducer::fillMuon(edm::Event& iEvent,
871  const edm::EventSetup& iSetup)
872 {
873  std::string MsgLoggerCat = "GlobalHitsProducer_fillMuon";
874 
875  TString eventout;
876  if (verbosity > 0)
877  eventout = "\nGathering info:";
878 
879  // iterator to access containers
880  edm::PSimHitContainer::const_iterator itHit;
881 
882  //int i = 0, j = 0;
884  // access the CSC Muon
886  // access the CSC Muon geometry
887  edm::ESHandle<CSCGeometry> theCSCGeometry;
888  iSetup.get<MuonGeometryRecord>().get(theCSCGeometry);
889  if (!theCSCGeometry.isValid()) {
890  edm::LogWarning(MsgLoggerCat)
891  << "Unable to find MuonGeometryRecord for the CSCGeometry in event!";
892  return;
893  }
894  const CSCGeometry& theCSCMuon(*theCSCGeometry);
895 
896  // get Muon CSC information
897  edm::Handle<edm::PSimHitContainer> MuonCSCContainer;
898  iEvent.getByToken(MuonCscSrc_Token_,MuonCSCContainer);
899  if (!MuonCSCContainer.isValid()) {
900  edm::LogWarning(MsgLoggerCat)
901  << "Unable to find MuonCSCHits in event!";
902  return;
903  }
904 
905  // cycle through container
906  int i = 0, j = 0;
907  for (itHit = MuonCSCContainer->begin(); itHit != MuonCSCContainer->end();
908  ++itHit) {
909 
910  ++i;
911 
912  // create a DetId from the detUnitId
913  DetId theDetUnitId(itHit->detUnitId());
914  int detector = theDetUnitId.det();
915  int subdetector = theDetUnitId.subdetId();
916 
917  // check that expected detector is returned
918  if ((detector == dMuon) &&
919  (subdetector == sdMuonCSC)) {
920 
921  // get the GeomDetUnit from the geometry using theDetUnitID
922  const GeomDetUnit *theDet = theCSCMuon.idToDetUnit(theDetUnitId);
923 
924  if (!theDet) {
925  edm::LogWarning(MsgLoggerCat)
926  << "Unable to get GeomDetUnit from theCSCMuon for hit " << i;
927  continue;
928  }
929 
930  ++j;
931 
932  // get the Surface of the hit (knows how to go from local <-> global)
933  const BoundPlane& bSurface = theDet->surface();
934 
935  // gather necessary information
936  MuonCscToF.push_back(itHit->tof());
937  MuonCscZ.push_back(bSurface.toGlobal(itHit->localPosition()).z());
938  MuonCscPhi.push_back(bSurface.toGlobal(itHit->localPosition()).phi());
939  MuonCscEta.push_back(bSurface.toGlobal(itHit->localPosition()).eta());
940  } else {
941  edm::LogWarning(MsgLoggerCat)
942  << "MuonCsc PSimHit " << i
943  << " is expected to be (det,subdet) = ("
944  << dMuon << "," << sdMuonCSC
945  << "); value returned is: ("
946  << detector << "," << subdetector << ")";
947  continue;
948  } // end detector type check
949  } // end loop through CSC Hits
950 
951  if (verbosity > 1) {
952  eventout += "\n Number of CSC muon Hits collected:......... ";
953  eventout += j;
954  }
955 
956  //i = 0, j = 0;
958  // access the DT Muon
960  // access the DT Muon geometry
961  edm::ESHandle<DTGeometry> theDTGeometry;
962  iSetup.get<MuonGeometryRecord>().get(theDTGeometry);
963  if (!theDTGeometry.isValid()) {
964  edm::LogWarning(MsgLoggerCat)
965  << "Unable to find MuonGeometryRecord for the DTGeometry in event!";
966  return;
967  }
968  const DTGeometry& theDTMuon(*theDTGeometry);
969 
970  // get Muon DT information
971  edm::Handle<edm::PSimHitContainer> MuonDtContainer;
972  iEvent.getByToken(MuonDtSrc_Token_,MuonDtContainer);
973  if (!MuonDtContainer.isValid()) {
974  edm::LogWarning(MsgLoggerCat)
975  << "Unable to find MuonDTHits in event!";
976  return;
977  }
978 
979  // cycle through container
980  i = 0, j = 0;
981  for (itHit = MuonDtContainer->begin(); itHit != MuonDtContainer->end();
982  ++itHit) {
983 
984  ++i;
985 
986  // create a DetId from the detUnitId
987  DetId theDetUnitId(itHit->detUnitId());
988  int detector = theDetUnitId.det();
989  int subdetector = theDetUnitId.subdetId();
990 
991  // check that expected detector is returned
992  if ((detector == dMuon) &&
993  (subdetector == sdMuonDT)) {
994 
995  // CSC uses wires and layers rather than the full detID
996  // get the wireId
997  DTWireId wireId(itHit->detUnitId());
998 
999  // get the DTLayer from the geometry using the wireID
1000  const DTLayer *theDet = theDTMuon.layer(wireId.layerId());
1001 
1002  if (!theDet) {
1003  edm::LogWarning(MsgLoggerCat)
1004  << "Unable to get GeomDetUnit from theDtMuon for hit " << i;
1005  continue;
1006  }
1007 
1008  ++j;
1009 
1010  // get the Surface of the hit (knows how to go from local <-> global)
1011  const BoundPlane& bSurface = theDet->surface();
1012 
1013  // gather necessary information
1014  MuonDtToF.push_back(itHit->tof());
1015  MuonDtR.push_back(bSurface.toGlobal(itHit->localPosition()).perp());
1016  MuonDtPhi.push_back(bSurface.toGlobal(itHit->localPosition()).phi());
1017  MuonDtEta.push_back(bSurface.toGlobal(itHit->localPosition()).eta());
1018  } else {
1019  edm::LogWarning(MsgLoggerCat)
1020  << "MuonDt PSimHit " << i
1021  << " is expected to be (det,subdet) = ("
1022  << dMuon << "," << sdMuonDT
1023  << "); value returned is: ("
1024  << detector << "," << subdetector << ")";
1025  continue;
1026  } // end detector type check
1027  } // end loop through DT Hits
1028 
1029  if (verbosity > 1) {
1030  eventout += "\n Number of DT muon Hits collected:.......... ";
1031  eventout += j;
1032  }
1033 
1034  //i = 0, j = 0;
1035  //int RPCBrl = 0, RPCFwd = 0;
1037  // access the RPC Muon
1039  // access the RPC Muon geometry
1040  edm::ESHandle<RPCGeometry> theRPCGeometry;
1041  iSetup.get<MuonGeometryRecord>().get(theRPCGeometry);
1042  if (!theRPCGeometry.isValid()) {
1043  edm::LogWarning(MsgLoggerCat)
1044  << "Unable to find MuonGeometryRecord for the RPCGeometry in event!";
1045  return;
1046  }
1047  const RPCGeometry& theRPCMuon(*theRPCGeometry);
1048 
1049  // get Muon RPC information
1050  edm::Handle<edm::PSimHitContainer> MuonRPCContainer;
1051  iEvent.getByToken(MuonRpcSrc_Token_,MuonRPCContainer);
1052  if (!MuonRPCContainer.isValid()) {
1053  edm::LogWarning(MsgLoggerCat)
1054  << "Unable to find MuonRPCHits in event!";
1055  return;
1056  }
1057 
1058  // cycle through container
1059  i = 0, j = 0;
1060  int RPCBrl =0, RPCFwd = 0;
1061  for (itHit = MuonRPCContainer->begin(); itHit != MuonRPCContainer->end();
1062  ++itHit) {
1063 
1064  ++i;
1065 
1066  // create a DetID from the detUnitId
1067  DetId theDetUnitId(itHit->detUnitId());
1068  int detector = theDetUnitId.det();
1069  int subdetector = theDetUnitId.subdetId();
1070 
1071  // check that expected detector is returned
1072  if ((detector == dMuon) &&
1073  (subdetector == sdMuonRPC)) {
1074 
1075  // get an RPCDetID from the detUnitID
1076  RPCDetId RPCId(itHit->detUnitId());
1077 
1078  // find the region of the RPC hit
1079  int region = RPCId.region();
1080 
1081  // get the GeomDetUnit from the geometry using the RPCDetId
1082  const GeomDetUnit *theDet = theRPCMuon.idToDetUnit(theDetUnitId);
1083 
1084  if (!theDet) {
1085  edm::LogWarning(MsgLoggerCat)
1086  << "Unable to get GeomDetUnit from theRPCMuon for hit " << i;
1087  continue;
1088  }
1089 
1090  ++j;
1091 
1092  // get the Surface of the hit (knows how to go from local <-> global)
1093  const BoundPlane& bSurface = theDet->surface();
1094 
1095  // gather necessary information
1096  if ((region == sdMuonRPCRgnFwdp) || (region == sdMuonRPCRgnFwdn)) {
1097  ++RPCFwd;
1098 
1099  MuonRpcFwdToF.push_back(itHit->tof());
1100  MuonRpcFwdZ.push_back(bSurface.toGlobal(itHit->localPosition()).z());
1101  MuonRpcFwdPhi.
1102  push_back(bSurface.toGlobal(itHit->localPosition()).phi());
1103  MuonRpcFwdEta.
1104  push_back(bSurface.toGlobal(itHit->localPosition()).eta());
1105  } else if (region == sdMuonRPCRgnBrl) {
1106  ++RPCBrl;
1107 
1108  MuonRpcBrlToF.push_back(itHit->tof());
1109  MuonRpcBrlR.
1110  push_back(bSurface.toGlobal(itHit->localPosition()).perp());
1111  MuonRpcBrlPhi.
1112  push_back(bSurface.toGlobal(itHit->localPosition()).phi());
1113  MuonRpcBrlEta.
1114  push_back(bSurface.toGlobal(itHit->localPosition()).eta());
1115  } else {
1116  edm::LogWarning(MsgLoggerCat)
1117  << "Invalid region for RPC Muon hit" << i;
1118  continue;
1119  } // end check of region
1120  } else {
1121  edm::LogWarning(MsgLoggerCat)
1122  << "MuonRpc PSimHit " << i
1123  << " is expected to be (det,subdet) = ("
1124  << dMuon << "," << sdMuonRPC
1125  << "); value returned is: ("
1126  << detector << "," << subdetector << ")";
1127  continue;
1128  } // end detector type check
1129  } // end loop through RPC Hits
1130 
1131  if (verbosity > 1) {
1132  eventout += "\n Number of RPC muon Hits collected:......... ";
1133  eventout += j;
1134  eventout += "\n RPC Barrel muon Hits:............ ";
1135  eventout += RPCBrl;
1136  eventout += "\n RPC Forward muon Hits:........... ";
1137  eventout += RPCFwd;
1138  }
1139 
1140  if (verbosity > 0)
1141  edm::LogInfo(MsgLoggerCat) << eventout << "\n";
1142 
1143  return;
1144 }
1145 
1146 void GlobalHitsProducer::storeMuon(PGlobalSimHit& product)
1147 {
1148  std::string MsgLoggerCat = "GlobalHitsProducer_storeMuon";
1149 
1150  if (verbosity > 2) {
1151  TString eventout("\n nMuonCSCHits = ");
1152  eventout += MuonCscToF.size();
1153  for (unsigned int i = 0; i < MuonCscToF.size(); ++i) {
1154  eventout += "\n (tof,z,phi,eta) = (";
1155  eventout += MuonCscToF[i];
1156  eventout += ", ";
1157  eventout += MuonCscZ[i];
1158  eventout += ", ";
1159  eventout += MuonCscPhi[i];
1160  eventout += ", ";
1161  eventout += MuonCscEta[i];
1162  eventout += ")";
1163  } // end MuonCsc output
1164  eventout += "\n nMuonDtHits = ";
1165  eventout += MuonDtToF.size();
1166  for (unsigned int i = 0; i < MuonDtToF.size(); ++i) {
1167  eventout += "\n (tof,r,phi,eta) = (";
1168  eventout += MuonDtToF[i];
1169  eventout += ", ";
1170  eventout += MuonDtR[i];
1171  eventout += ", ";
1172  eventout += MuonDtPhi[i];
1173  eventout += ", ";
1174  eventout += MuonDtEta[i];
1175  eventout += ")";
1176  } // end MuonDt output
1177  eventout += "\n nMuonRpcBrlHits = ";
1178  eventout += MuonRpcBrlToF.size();
1179  for (unsigned int i = 0; i < MuonRpcBrlToF.size(); ++i) {
1180  eventout += "\n (tof,r,phi,eta) = (";
1181  eventout += MuonRpcBrlToF[i];
1182  eventout += ", ";
1183  eventout += MuonRpcBrlR[i];
1184  eventout += ", ";
1185  eventout += MuonRpcBrlPhi[i];
1186  eventout += ", ";
1187  eventout += MuonRpcBrlEta[i];
1188  eventout += ")";
1189  } // end MuonRpcBrl output
1190  eventout += "\n nMuonRpcFwdHits = ";
1191  eventout += MuonRpcFwdToF.size();
1192  for (unsigned int i = 0; i < MuonRpcFwdToF.size(); ++i) {
1193  eventout += "\n (tof,z,phi,eta) = (";
1194  eventout += MuonRpcFwdToF[i];
1195  eventout += ", ";
1196  eventout += MuonRpcFwdZ[i];
1197  eventout += ", ";
1198  eventout += MuonRpcFwdPhi[i];
1199  eventout += ", ";
1200  eventout += MuonRpcFwdEta[i];
1201  eventout += ")";
1202  } // end MuonRpcFwd output
1203  edm::LogInfo(MsgLoggerCat) << eventout << "\n";
1204  } // end verbose output
1205 
1206  product.putMuonCscHits(MuonCscToF,MuonCscZ,MuonCscPhi,MuonCscEta);
1207  product.putMuonDtHits(MuonDtToF,MuonDtR,MuonDtPhi,MuonDtEta);
1208  product.putMuonRpcBrlHits(MuonRpcBrlToF,MuonRpcBrlR,MuonRpcBrlPhi,
1209  MuonRpcBrlEta);
1210  product.putMuonRpcFwdHits(MuonRpcFwdToF,MuonRpcFwdZ,MuonRpcFwdPhi,
1211  MuonRpcFwdEta);
1212 
1213  return;
1214 }
1215 
1216 void GlobalHitsProducer::fillECal(edm::Event& iEvent,
1217  const edm::EventSetup& iSetup)
1218 {
1219  std::string MsgLoggerCat = "GlobalHitsProducer_fillECal";
1220 
1221  TString eventout;
1222  if (verbosity > 0)
1223  eventout = "\nGathering info:";
1224 
1225  // access the calorimeter geometry
1226  edm::ESHandle<CaloGeometry> theCaloGeometry;
1227  iSetup.get<CaloGeometryRecord>().get(theCaloGeometry);
1228  if (!theCaloGeometry.isValid()) {
1229  edm::LogWarning(MsgLoggerCat)
1230  << "Unable to find CaloGeometryRecord in event!";
1231  return;
1232  }
1233  const CaloGeometry& theCalo(*theCaloGeometry);
1234 
1235  // iterator to access containers
1236  edm::PCaloHitContainer::const_iterator itHit;
1237 
1239  // get ECal information
1241  edm::PCaloHitContainer theECalHits;
1242  // extract EB container
1243  edm::Handle<edm::PCaloHitContainer> EBContainer;
1244  iEvent.getByToken(ECalEBSrc_Token_,EBContainer);
1245  if (!EBContainer.isValid()) {
1246  edm::LogWarning(MsgLoggerCat)
1247  << "Unable to find EcalHitsEB in event!";
1248  return;
1249  }
1250  // extract EE container
1251  edm::Handle<edm::PCaloHitContainer> EEContainer;
1252  iEvent.getByToken(ECalEESrc_Token_,EEContainer);
1253  if (!EEContainer.isValid()) {
1254  edm::LogWarning(MsgLoggerCat)
1255  << "Unable to find EcalHitsEE in event!";
1256  return;
1257  }
1258  // place both containers into new container
1259  theECalHits.insert(theECalHits.end(),EBContainer->begin(),
1260  EBContainer->end());
1261  theECalHits.insert(theECalHits.end(),EEContainer->begin(),
1262  EEContainer->end());
1263 
1264  // cycle through new container
1265  int i = 0, j = 0;
1266  for (itHit = theECalHits.begin(); itHit != theECalHits.end(); ++itHit) {
1267 
1268  ++i;
1269 
1270  // create a DetId from the detUnitId
1271  DetId theDetUnitId(itHit->id());
1272  int detector = theDetUnitId.det();
1273  int subdetector = theDetUnitId.subdetId();
1274 
1275  // check that expected detector is returned
1276  if ((detector == dEcal) &&
1277  ((subdetector == sdEcalBrl) ||
1278  (subdetector == sdEcalFwd))) {
1279 
1280  // get the Cell geometry
1281  const CaloCellGeometry *theDet = theCalo.
1282  getSubdetectorGeometry(theDetUnitId)->getGeometry(theDetUnitId);
1283 
1284  if (!theDet) {
1285  edm::LogWarning(MsgLoggerCat)
1286  << "Unable to get CaloCellGeometry from ECalHits for Hit " << i;
1287  continue;
1288  }
1289 
1290  ++j;
1291 
1292  // get the global position of the cell
1293  const GlobalPoint& globalposition = theDet->getPosition();
1294 
1295  // gather necessary information
1296  ECalE.push_back(itHit->energy());
1297  ECalToF.push_back(itHit->time());
1298  ECalPhi.push_back(globalposition.phi());
1299  ECalEta.push_back(globalposition.eta());
1300 
1301  } else {
1302  edm::LogWarning(MsgLoggerCat)
1303  << "ECal PCaloHit " << i
1304  << " is expected to be (det,subdet) = ("
1305  << dEcal << "," << sdEcalBrl
1306  << " || " << sdEcalFwd << "); value returned is: ("
1307  << detector << "," << subdetector << ")";
1308  continue;
1309  } // end detector type check
1310  } // end loop through ECal Hits
1311 
1312  if (verbosity > 1) {
1313  eventout += "\n Number of ECal Hits collected:............. ";
1314  eventout += j;
1315  }
1316 
1318  // Get Preshower information
1320  // extract PreShower container
1321  edm::Handle<edm::PCaloHitContainer> PreShContainer;
1322  iEvent.getByToken(ECalESSrc_Token_,PreShContainer);
1323  if (!PreShContainer.isValid()) {
1324  edm::LogWarning(MsgLoggerCat)
1325  << "Unable to find EcalHitsES in event!";
1326  return;
1327  }
1328 
1329  // cycle through container
1330  i = 0, j = 0;
1331  for (itHit = PreShContainer->begin();
1332  itHit != PreShContainer->end(); ++itHit) {
1333 
1334  ++i;
1335 
1336  // create a DetId from the detUnitId
1337  DetId theDetUnitId(itHit->id());
1338  int detector = theDetUnitId.det();
1339  int subdetector = theDetUnitId.subdetId();
1340 
1341  // check that expected detector is returned
1342  if ((detector == dEcal) &&
1343  (subdetector == sdEcalPS)) {
1344 
1345  // get the Cell geometry
1346  const CaloCellGeometry *theDet = theCalo.
1347  getSubdetectorGeometry(theDetUnitId)->getGeometry(theDetUnitId);
1348 
1349  if (!theDet) {
1350  edm::LogWarning(MsgLoggerCat)
1351  << "Unable to get CaloCellGeometry from PreShContainer for Hit "
1352  << i;
1353  continue;
1354  }
1355 
1356  ++j;
1357 
1358  // get the global position of the cell
1359  const GlobalPoint& globalposition = theDet->getPosition();
1360 
1361  // gather necessary information
1362  PreShE.push_back(itHit->energy());
1363  PreShToF.push_back(itHit->time());
1364  PreShPhi.push_back(globalposition.phi());
1365  PreShEta.push_back(globalposition.eta());
1366 
1367  } else {
1368  edm::LogWarning(MsgLoggerCat)
1369  << "PreSh PCaloHit " << i
1370  << " is expected to be (det,subdet) = ("
1371  << dEcal << "," << sdEcalPS
1372  << "); value returned is: ("
1373  << detector << "," << subdetector << ")";
1374  continue;
1375  } // end detector type check
1376  } // end loop through PreShower Hits
1377 
1378  if (verbosity > 1) {
1379  eventout += "\n Number of PreSh Hits collected:............ ";
1380  eventout += j;
1381  }
1382 
1383  if (verbosity > 0)
1384  edm::LogInfo(MsgLoggerCat) << eventout << "\n";
1385 
1386  return;
1387 }
1388 
1389 void GlobalHitsProducer::storeECal(PGlobalSimHit& product)
1390 {
1391  std::string MsgLoggerCat = "GlobalHitsProducer_storeECal";
1392 
1393  if (verbosity > 2) {
1394  TString eventout("\n nECalHits = ");
1395  eventout += ECalE.size();
1396  for (unsigned int i = 0; i < ECalE.size(); ++i) {
1397  eventout += "\n (e,tof,phi,eta) = (";
1398  eventout += ECalE[i];
1399  eventout += ", ";
1400  eventout += ECalToF[i];
1401  eventout += ", ";
1402  eventout += ECalPhi[i];
1403  eventout += ", ";
1404  eventout += ECalEta[i];
1405  eventout += ")";
1406  } // end ECal output
1407  eventout += "\n nPreShHits = ";
1408  eventout += PreShE.size();
1409  for (unsigned int i = 0; i < PreShE.size(); ++i) {
1410  eventout += "\n (e,tof,phi,eta) = (";
1411  eventout += PreShE[i];
1412  eventout += ", ";
1413  eventout += PreShToF[i];
1414  eventout += ", ";
1415  eventout += PreShPhi[i];
1416  eventout += ", ";
1417  eventout += PreShEta[i];
1418  eventout += ")";
1419  } // end PreShower output
1420  edm::LogInfo(MsgLoggerCat) << eventout << "\n";
1421  } // end verbose output
1422 
1423  product.putECalHits(ECalE,ECalToF,ECalPhi,ECalEta);
1424  product.putPreShHits(PreShE,PreShToF,PreShPhi,PreShEta);
1425 
1426  return;
1427 }
1428 
1429 void GlobalHitsProducer::fillHCal(edm::Event& iEvent,
1430  const edm::EventSetup& iSetup)
1431 {
1432  std::string MsgLoggerCat = "GlobalHitsProducer_fillHCal";
1433 
1434  TString eventout;
1435  if (verbosity > 0)
1436  eventout = "\nGathering info:";
1437 
1438  // access the calorimeter geometry
1439  edm::ESHandle<CaloGeometry> theCaloGeometry;
1440  iSetup.get<CaloGeometryRecord>().get(theCaloGeometry);
1441  if (!theCaloGeometry.isValid()) {
1442  edm::LogWarning(MsgLoggerCat)
1443  << "Unable to find CaloGeometryRecord in event!";
1444  return;
1445  }
1446  const CaloGeometry& theCalo(*theCaloGeometry);
1447 
1448  // iterator to access containers
1449  edm::PCaloHitContainer::const_iterator itHit;
1450 
1452  // get HCal information
1454  // extract HCal container
1455  edm::Handle<edm::PCaloHitContainer> HCalContainer;
1456  iEvent.getByToken(HCalSrc_Token_,HCalContainer);
1457  if (!HCalContainer.isValid()) {
1458  edm::LogWarning(MsgLoggerCat)
1459  << "Unable to find HCalHits in event!";
1460  return;
1461  }
1462 
1463  // cycle through container
1464  int i = 0, j = 0;
1465  for (itHit = HCalContainer->begin();
1466  itHit != HCalContainer->end(); ++itHit) {
1467 
1468  ++i;
1469 
1470  // create a DetId from the detUnitId
1471  DetId theDetUnitId(itHit->id());
1472  int detector = theDetUnitId.det();
1473  int subdetector = theDetUnitId.subdetId();
1474 
1475  // check that expected detector is returned
1476  if ((detector == dHcal) &&
1477  ((subdetector == sdHcalBrl) ||
1478  (subdetector == sdHcalEC) ||
1479  (subdetector == sdHcalOut) ||
1480  (subdetector == sdHcalFwd))) {
1481 
1482  // get the Cell geometry
1483  const CaloCellGeometry *theDet = theCalo.
1484  getSubdetectorGeometry(theDetUnitId)->getGeometry(theDetUnitId);
1485 
1486  if (!theDet) {
1487  edm::LogWarning(MsgLoggerCat)
1488  << "Unable to get CaloCellGeometry from HCalContainer for Hit " << i;
1489  continue;
1490  }
1491 
1492  ++j;
1493 
1494  // get the global position of the cell
1495  const GlobalPoint& globalposition = theDet->getPosition();
1496 
1497  // gather necessary information
1498  HCalE.push_back(itHit->energy());
1499  HCalToF.push_back(itHit->time());
1500  HCalPhi.push_back(globalposition.phi());
1501  HCalEta.push_back(globalposition.eta());
1502 
1503  } else {
1504  edm::LogWarning(MsgLoggerCat)
1505  << "HCal PCaloHit " << i
1506  << " is expected to be (det,subdet) = ("
1507  << dHcal << "," << sdHcalBrl
1508  << " || " << sdHcalEC << " || " << sdHcalOut << " || " << sdHcalFwd
1509  << "); value returned is: ("
1510  << detector << "," << subdetector << ")";
1511  continue;
1512  } // end detector type check
1513  } // end loop through HCal Hits
1514 
1515  if (verbosity > 1) {
1516  eventout += "\n Number of HCal Hits collected:............. ";
1517  eventout += j;
1518  }
1519 
1520  if (verbosity > 0)
1521  edm::LogInfo(MsgLoggerCat) << eventout << "\n";
1522 
1523  return;
1524 }
1525 
1526 void GlobalHitsProducer::storeHCal(PGlobalSimHit& product)
1527 {
1528  std::string MsgLoggerCat = "GlobalHitsProducer_storeHCal";
1529 
1530  if (verbosity > 2) {
1531  TString eventout("\n nHCalHits = ");
1532  eventout += HCalE.size();
1533  for (unsigned int i = 0; i < HCalE.size(); ++i) {
1534  eventout += "\n (e,tof,phi,eta) = (";
1535  eventout += HCalE[i];
1536  eventout += ", ";
1537  eventout += HCalToF[i];
1538  eventout += ", ";
1539  eventout += HCalPhi[i];
1540  eventout += ", ";
1541  eventout += HCalEta[i];
1542  eventout += ")";
1543  } // end HCal output
1544  edm::LogInfo(MsgLoggerCat) << eventout << "\n";
1545  } // end verbose output
1546 
1547  product.putHCalHits(HCalE,HCalToF,HCalPhi,HCalEta);
1548 
1549  return;
1550 }
1551 
1552 void GlobalHitsProducer::clear()
1553 {
1554  std::string MsgLoggerCat = "GlobalHitsProducer_clear";
1555 
1556  if (verbosity > 0)
1557  edm::LogInfo(MsgLoggerCat)
1558  << "Clearing event holders";
1559 
1560  // reset G4MC info
1561  nRawGenPart = 0;
1562  G4VtxX.clear();
1563  G4VtxY.clear();
1564  G4VtxZ.clear();
1565  G4TrkPt.clear();
1566  G4TrkE.clear();
1567 
1568  // reset electromagnetic info
1569  // reset ECal info
1570  ECalE.clear();
1571  ECalToF.clear();
1572  ECalPhi.clear();
1573  ECalEta.clear();
1574  // reset Preshower info
1575  PreShE.clear();
1576  PreShToF.clear();
1577  PreShPhi.clear();
1578  PreShEta.clear();
1579 
1580  // reset hadronic info
1581  // reset HCal info
1582  HCalE.clear();
1583  HCalToF.clear();
1584  HCalPhi.clear();
1585  HCalEta.clear();
1586 
1587  // reset tracker info
1588  // reset Pixel info
1589  PxlBrlToF.clear();
1590  PxlBrlR.clear();
1591  PxlBrlPhi.clear();
1592  PxlBrlEta.clear();
1593  PxlFwdToF.clear();
1594  PxlFwdZ.clear();
1595  PxlFwdPhi.clear();
1596  PxlFwdEta.clear();
1597  // reset strip info
1598  SiBrlToF.clear();
1599  SiBrlR.clear();
1600  SiBrlPhi.clear();
1601  SiBrlEta.clear();
1602  SiFwdToF.clear();
1603  SiFwdZ.clear();
1604  SiFwdPhi.clear();
1605  SiFwdEta.clear();
1606 
1607  // reset muon info
1608  // reset muon DT info
1609  MuonDtToF.clear();
1610  MuonDtR.clear();
1611  MuonDtPhi.clear();
1612  MuonDtEta.clear();
1613  // reset muon CSC info
1614  MuonCscToF.clear();
1615  MuonCscZ.clear();
1616  MuonCscPhi.clear();
1617  MuonCscEta.clear();
1618  // rest muon RPC info
1619  MuonRpcBrlToF.clear();
1620  MuonRpcBrlR.clear();
1621  MuonRpcBrlPhi.clear();
1622  MuonRpcBrlEta.clear();
1623  MuonRpcFwdToF.clear();
1624  MuonRpcFwdZ.clear();
1625  MuonRpcFwdPhi.clear();
1626  MuonRpcFwdEta.clear();
1627 
1628  return;
1629 }
GlobalHitsProducer::ECalESSrc_
edm::InputTag ECalESSrc_
Definition: GlobalHitsProducer.h:147
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:42
GlobalHitsProducer::storeECal
void storeECal(PGlobalSimHit &)
Definition: GlobalHitsProducer.cc:1389
GlobalHitsProducer::SiTIBHighSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > SiTIBHighSrc_Token_
Definition: GlobalHitsProducer.h:196
edm::Event::getManyByType
void getManyByType(std::vector< Handle< PROD > > &results) const
Definition: Event.h:424
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:44
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
edm::Event::getAllProvenance
void getAllProvenance(std::vector< Provenance const * > &provenances) const
Definition: Event.cc:86
i
int i
Definition: DBlmapReader.cc:9
dEcal
static const int dEcal
Definition: GlobalDigisAnalyzer.h:273
dTrk
static const int dTrk
Definition: GlobalDigisAnalyzer.h:257
GlobalHitsProducer::SiTECHighSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > SiTECHighSrc_Token_
Definition: GlobalHitsProducer.h:202
edm::PCaloHitContainer
std::vector< PCaloHit > PCaloHitContainer
Definition: PCaloHitContainer.h:8
sdMuonRPCRgnFwdn
static const int sdMuonRPCRgnFwdn
Definition: GlobalDigisAnalyzer.h:271
GlobalHitsProducer::produce
virtual void produce(edm::Event &, const edm::EventSetup &) override
Definition: GlobalHitsProducer.cc:171
sdHcalOut
static const int sdHcalOut
Definition: GlobalDigisAnalyzer.h:284
GlobalHitsProducer::ECalEESrc_
edm::InputTag ECalEESrc_
Definition: GlobalHitsProducer.h:138
PGlobalSimHit::putHCalHits
void putHCalHits(const std::vector< float > &e, const std::vector< float > &tof, const std::vector< float > &phi, const std::vector< float > &eta)
Definition: PValidationFormats.cc:79
sdMuonDT
static const int sdMuonDT
Definition: GlobalDigisAnalyzer.h:266
GlobalHitsProducer::MuonRpcSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > MuonRpcSrc_Token_
Definition: GlobalHitsProducer.h:229
GlobalHitsProducer::MuonCscSrc_
edm::InputTag MuonCscSrc_
Definition: GlobalHitsProducer.h:217
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434
PGlobalSimHit::putMuonRpcFwdHits
void putMuonRpcFwdHits(const std::vector< float > &tof, const std::vector< float > &z, const std::vector< float > &phi, const std::vector< float > &eta)
Definition: PValidationFormats.cc:196
GlobalHitsProducer::MuonDtSrc_
edm::InputTag MuonDtSrc_
Definition: GlobalHitsProducer.h:210
PGlobalSimHit::putMuonRpcBrlHits
void putMuonRpcBrlHits(const std::vector< float > &tof, const std::vector< float > &r, const std::vector< float > &phi, const std::vector< float > &eta)
Definition: PValidationFormats.cc:213
sdSiTID
static const int sdSiTID
Definition: GlobalDigisAnalyzer.h:261
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
GlobalHitsProducer::SiTOBHighSrc_
edm::InputTag SiTOBHighSrc_
Definition: GlobalHitsProducer.h:190
GlobalHitsProducer::PxlBrlHighSrc_
edm::InputTag PxlBrlHighSrc_
Definition: GlobalHitsProducer.h:170
GlobalHitsProducer::PxlBrlLowSrc_
edm::InputTag PxlBrlLowSrc_
Definition: GlobalHitsProducer.h:169
sdEcalPS
static const int sdEcalPS
Definition: GlobalDigisAnalyzer.h:276
PGlobalSimHit::putMuonCscHits
void putMuonCscHits(const std::vector< float > &tof, const std::vector< float > &z, const std::vector< float > &phi, const std::vector< float > &eta)
Definition: PValidationFormats.cc:162
PGlobalSimHit::putECalHits
void putECalHits(const std::vector< float > &e, const std::vector< float > &tof, const std::vector< float > &phi, const std::vector< float > &eta)
Definition: PValidationFormats.cc:47
GlobalHitsProducer::storeMuon
void storeMuon(PGlobalSimHit &)
Definition: GlobalHitsProducer.cc:1146
GlobalHitsProducer::PxlFwdLowSrc_
edm::InputTag PxlFwdLowSrc_
Definition: GlobalHitsProducer.h:171
sdMuonRPC
static const int sdMuonRPC
Definition: GlobalDigisAnalyzer.h:268
GlobalHitsProducer::storeG4MC
void storeG4MC(PGlobalSimHit &)
Definition: GlobalHitsProducer.cc:392
eta
T eta() const
Definition: Basic3DVectorLD.h:177
GlobalHitsProducer::PxlBrlHighSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > PxlBrlHighSrc_Token_
Definition: GlobalHitsProducer.h:174
GlobalHitsProducer::SiTIDLowSrc_
edm::InputTag SiTIDLowSrc_
Definition: GlobalHitsProducer.h:191
GlobalHitsProducer::storeTrk
void storeTrk(PGlobalSimHit &)
Definition: GlobalHitsProducer.cc:802
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:35
detailsBasic3DVector::z
float float float z
Definition: extBasic3DVector.h:15
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
DTGeometry::layer
const DTLayer * layer(DTLayerId id) const
Return a layer given its id.
Definition: DTGeometry.cc:110
GlobalHitsProducer::PxlFwdHighSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > PxlFwdHighSrc_Token_
Definition: GlobalHitsProducer.h:176
GlobalHitsProducer::SiTECHighSrc_
edm::InputTag SiTECHighSrc_
Definition: GlobalHitsProducer.h:194
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
GlobalHitsProducer::SiTIBLowSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > SiTIBLowSrc_Token_
Definition: GlobalHitsProducer.h:195
GlobalHitsProducer::ECalEBSrc_
edm::InputTag ECalEBSrc_
Definition: GlobalHitsProducer.h:137
GlobalHitsProducer::fillMuon
void fillMuon(edm::Event &, const edm::EventSetup &)
Definition: GlobalHitsProducer.cc:870
PGlobalSimHit::putG4Trk
void putG4Trk(const std::vector< float > &pt, const std::vector< float > &e)
Definition: PValidationFormats.cc:35
GlobalHitsProducer::HCalSrc_Token_
edm::EDGetTokenT< edm::PCaloHitContainer > HCalSrc_Token_
Definition: GlobalHitsProducer.h:157
PFRecoTauDiscriminationAgainstMuon2_cfi.verbosity
tuple verbosity
Definition: PFRecoTauDiscriminationAgainstMuon2_cfi.py:38
GlobalHitsProducer::MuonCscSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > MuonCscSrc_Token_
Definition: GlobalHitsProducer.h:218
iEvent
int iEvent
Definition: GenABIO.cc:230
GlobalHitsProducer::GlobalHitsProducer
GlobalHitsProducer(const edm::ParameterSet &)
Definition: GlobalHitsProducer.cc:11
sdSiTIB
static const int sdSiTIB
Definition: GlobalDigisAnalyzer.h:260
sdPxlBrl
static const int sdPxlBrl
Definition: GlobalDigisAnalyzer.h:258
edm::Event::put
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:116
csvLumiCalc.unit
string unit
Definition: csvLumiCalc.py:46
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
PGlobalSimHit::putSiBrlHits
void putSiBrlHits(const std::vector< float > &tof, const std::vector< float > &r, const std::vector< float > &phi, const std::vector< float > &eta)
Definition: PValidationFormats.cc:146
sdMuonCSC
static const int sdMuonCSC
Definition: GlobalDigisAnalyzer.h:267
PGlobalSimHit::putPxlFwdHits
void putPxlFwdHits(const std::vector< float > &tof, const std::vector< float > &z, const std::vector< float > &phi, const std::vector< float > &eta)
Definition: PValidationFormats.cc:95
sdEcalFwd
static const int sdEcalFwd
Definition: GlobalDigisAnalyzer.h:275
GlobalHitsProducer::fillTrk
void fillTrk(edm::Event &, const edm::EventSetup &)
Definition: GlobalHitsProducer.cc:429
j
int j
Definition: DBlmapReader.cc:9
GlobalHitsProducer::SiTIBLowSrc_
edm::InputTag SiTIBLowSrc_
Definition: GlobalHitsProducer.h:187
GlobalHitsProducer::SiTIDLowSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > SiTIDLowSrc_Token_
Definition: GlobalHitsProducer.h:199
GlobalHitsProducer::SiTECLowSrc_
edm::InputTag SiTECLowSrc_
Definition: GlobalHitsProducer.h:193
sdSiTOB
static const int sdSiTOB
Definition: GlobalDigisAnalyzer.h:262
sdEcalBrl
static const int sdEcalBrl
Definition: GlobalDigisAnalyzer.h:274
edm::Provenance::product
BranchDescription const & product() const
Definition: Provenance.h:42
sdMuonRPCRgnFwdp
static const int sdMuonRPCRgnFwdp
Definition: GlobalDigisAnalyzer.h:270
sdMuonRPCRgnBrl
static const int sdMuonRPCRgnBrl
Definition: GlobalDigisAnalyzer.h:269
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
GlobalHitsProducer::SiTIDHighSrc_
edm::InputTag SiTIDHighSrc_
Definition: GlobalHitsProducer.h:192
GlobalHitsProducer::SiTOBLowSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > SiTOBLowSrc_Token_
Definition: GlobalHitsProducer.h:197
PGlobalSimHit::putSiFwdHits
void putSiFwdHits(const std::vector< float > &tof, const std::vector< float > &z, const std::vector< float > &phi, const std::vector< float > &eta)
Definition: PValidationFormats.cc:129
GlobalHitsProducer::G4TrkSrc_Token_
edm::EDGetTokenT< edm::SimTrackContainer > G4TrkSrc_Token_
Definition: GlobalHitsProducer.h:129
GlobalHitsProducer::PxlFwdHighSrc_
edm::InputTag PxlFwdHighSrc_
Definition: GlobalHitsProducer.h:172
GlobalHitsProducer::ECalESSrc_Token_
edm::EDGetTokenT< edm::PCaloHitContainer > ECalESSrc_Token_
Definition: GlobalHitsProducer.h:148
DetId
Definition: DetId.h:18
dHcal
static const int dHcal
Definition: GlobalDigisAnalyzer.h:280
GlobalHitsProducer::SiTOBHighSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > SiTOBHighSrc_Token_
Definition: GlobalHitsProducer.h:198
GlobalHitsProducer::beginJob
virtual void beginJob(void)
Definition: GlobalHitsProducer.cc:157
PGlobalSimHit::putPreShHits
void putPreShHits(const std::vector< float > &e, const std::vector< float > &tof, const std::vector< float > &phi, const std::vector< float > &eta)
Definition: PValidationFormats.cc:63
GlobalHitsProducer::fillECal
void fillECal(edm::Event &, const edm::EventSetup &)
Definition: GlobalHitsProducer.cc:1216
GlobalHitsProducer::fillHCal
void fillHCal(edm::Event &, const edm::EventSetup &)
Definition: GlobalHitsProducer.cc:1429
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::SimVertexContainer
std::vector< SimVertex > SimVertexContainer
Definition: SimVertexContainer.h:13
GlobalHitsProducer::SiTIDHighSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > SiTIDHighSrc_Token_
Definition: GlobalHitsProducer.h:200
sdHcalFwd
static const int sdHcalFwd
Definition: GlobalDigisAnalyzer.h:285
sdHcalBrl
static const int sdHcalBrl
Definition: GlobalDigisAnalyzer.h:282
GlobalHitsProducer::storeHCal
void storeHCal(PGlobalSimHit &)
Definition: GlobalHitsProducer.cc:1526
GlobalHitsProducer::SiTOBLowSrc_
edm::InputTag SiTOBLowSrc_
Definition: GlobalHitsProducer.h:189
GlobalHitsProducer::PxlFwdLowSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > PxlFwdLowSrc_Token_
Definition: GlobalHitsProducer.h:175
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:42
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:76
TrackerGeometry::idToDetUnit
virtual const GeomDetUnit * idToDetUnit(DetId) const
Return the pointer to the GeomDetUnit corresponding to a given DetId.
Definition: TrackerGeometry.cc:126
GlobalHitsProducer::MuonDtSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > MuonDtSrc_Token_
Definition: GlobalHitsProducer.h:211
GlobalHitsProducer::PxlBrlLowSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > PxlBrlLowSrc_Token_
Definition: GlobalHitsProducer.h:173
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
GlobalHitsProducer::G4VtxSrc_Token_
edm::EDGetTokenT< edm::SimVertexContainer > G4VtxSrc_Token_
Definition: GlobalHitsProducer.h:128
sdSiTEC
static const int sdSiTEC
Definition: GlobalDigisAnalyzer.h:263
perp
T perp() const
Magnitude of transverse component.
Definition: Basic3DVectorLD.h:153
dMuon
static const int dMuon
Definition: GlobalDigisAnalyzer.h:265
GlobalHitsProducer::ECalEBSrc_Token_
edm::EDGetTokenT< edm::PCaloHitContainer > ECalEBSrc_Token_
Definition: GlobalHitsProducer.h:139
GlobalHitsProducer::ECalEESrc_Token_
edm::EDGetTokenT< edm::PCaloHitContainer > ECalEESrc_Token_
Definition: GlobalHitsProducer.h:140
PGlobalSimHit::putPxlBrlHits
void putPxlBrlHits(const std::vector< float > &tof, const std::vector< float > &r, const std::vector< float > &phi, const std::vector< float > &eta)
Definition: PValidationFormats.cc:112
CaloGeometry::getGeometry
const CaloCellGeometry * getGeometry(const DetId &id) const
Get the cell geometry of a given detector id.
Definition: CaloGeometry.cc:76
PGlobalSimHit::putRawGenPart
void putRawGenPart(int n)
Definition: PValidationFormats.cc:15
PGlobalSimHit::putG4Vtx
void putG4Vtx(const std::vector< float > &x, const std::vector< float > &y, const std::vector< float > &z)
Definition: PValidationFormats.cc:21
GlobalHitsProducer::fillG4MC
void fillG4MC(edm::Event &)
Definition: GlobalHitsProducer.cc:273
PGlobalSimHit::putMuonDtHits
void putMuonDtHits(const std::vector< float > &tof, const std::vector< float > &r, const std::vector< float > &phi, const std::vector< float > &eta)
Definition: PValidationFormats.cc:179
CSCGeometry::idToDetUnit
virtual const GeomDetUnit * idToDetUnit(DetId) const override
Return the pointer to the GeomDetUnit corresponding to a given DetId.
Definition: CSCGeometry.cc:93
edm::PSimHitContainer
std::vector< PSimHit > PSimHitContainer
Definition: PSimHitContainer.h:11
GlobalHitsProducer::SiTECLowSrc_Token_
edm::EDGetTokenT< edm::PSimHitContainer > SiTECLowSrc_Token_
Definition: GlobalHitsProducer.h:201
sdPxlFwd
static const int sdPxlFwd
Definition: GlobalDigisAnalyzer.h:259
edm::false
volatile std::atomic< bool > shutdown_flag false
Definition: UnixSignalHandlers.cc:22
edm::ESHandleBase::isValid
bool isValid() const
Definition: ESHandle.h:37
DetId::det
Detector det() const
get the detector field from this detid
Definition: DetId.h:35
CaloCellGeometry::getPosition
const GlobalPoint & getPosition() const
Returns the position of reference for this cell.
Definition: CaloCellGeometry.h:77
GlobalHitsProducer::SiTIBHighSrc_
edm::InputTag SiTIBHighSrc_
Definition: GlobalHitsProducer.h:188
edm::SimTrackContainer
std::vector< SimTrack > SimTrackContainer
Definition: SimTrackContainer.h:12
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:43
GlobalHitsProducer::MuonRpcSrc_
edm::InputTag MuonRpcSrc_
Definition: GlobalHitsProducer.h:228
sdHcalEC
static const int sdHcalEC
Definition: GlobalDigisAnalyzer.h:283
edm::HandleBase::provenance
Provenance const * provenance() const
Definition: HandleBase.h:85
RPCGeometry::idToDetUnit
virtual const GeomDetUnit * idToDetUnit(DetId) const override
Return the pointer to the GeomDetUnit corresponding to a given DetId.
Definition: RPCGeometry.cc:48
RPCDetId::region
int region() const
Region id: 0 for Barrel, +/-1 For +/- Endcap.
Definition: RPCDetId.h:63
phi
Definition: DDAxes.h:10
Generated for CMSSW Reference Manual by doxygen 1.8.5