CMS 3D CMS Logo
Phase2ITValidateRecHit.cc
Go to the documentation of this file.
1 // Package: Phase2ITValidateRecHit
2 // Class: Phase2ITValidateRecHit
3 //
7 //
8 // Author: Shubhi Parolia, Suvankar Roy Chowdhury
9 // Date: June 2020
10 //
11 // system include files
12 #include <memory>
13 #include <map>
14 #include <vector>
15 #include <algorithm>
16 #include "FWCore/MessageLogger/interface/MessageLogger.h"
17 #include "FWCore/Framework/interface/ESWatcher.h"
18 #include "FWCore/Framework/interface/Frameworkfwd.h"
19 #include "FWCore/Framework/interface/Event.h"
20 #include "FWCore/Framework/interface/MakerMacros.h"
21 #include "FWCore/Framework/interface/ESHandle.h"
22 #include "FWCore/ParameterSet/interface/ParameterSet.h"
23 #include "FWCore/ServiceRegistry/interface/Service.h"
24 #include "FWCore/Utilities/interface/InputTag.h"
25 #include "DataFormats/Common/interface/Handle.h"
26 #include "DataFormats/Common/interface/DetSetVector.h"
27 #include "DataFormats/DetId/interface/DetId.h"
28 #include "DataFormats/GeometrySurface/interface/LocalError.h"
29 #include "DataFormats/GeometryVector/interface/LocalPoint.h"
30 #include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
31 #include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHitCollection.h"
32 #include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
33 #include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
34 #include "Geometry/CommonDetUnit/interface/GeomDet.h"
35 #include "Geometry/CommonDetUnit/interface/TrackerGeomDet.h"
36 #include "Geometry/CommonDetUnit/interface/PixelGeomDetUnit.h"
37 #include "Geometry/CommonDetUnit/interface/PixelGeomDetType.h"
38 #include "Geometry/Records/interface/TrackerTopologyRcd.h"
39 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
40 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
41 //--- for SimHit association
42 #include "SimDataFormats/Track/interface/SimTrackContainer.h"
43 #include "SimDataFormats/TrackingHit/interface/PSimHit.h"
44 #include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
45 #include "SimTracker/TrackerHitAssociation/interface/TrackerHitAssociator.h"
46 //DQM
47 #include "DQMServices/Core/interface/DQMEDAnalyzer.h"
48 #include "DQMServices/Core/interface/DQMStore.h"
49 #include "DQMServices/Core/interface/MonitorElement.h"
50 #include "Validation/SiTrackerPhase2V/interface/TrackerPhase2ValidationUtil.h"
51 #include "DQM/SiTrackerPhase2/interface/TrackerPhase2DQMUtil.h"
52 
53 class Phase2ITValidateRecHit : public DQMEDAnalyzer {
54 public:
55  explicit Phase2ITValidateRecHit(const edm::ParameterSet&);
56  ~Phase2ITValidateRecHit() override;
57  void bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iRun, edm::EventSetup const& iSetup) override;
58  void analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) override;
59  void dqmBeginRun(const edm::Run&, const edm::EventSetup&) override;
60  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
61 
62 private:
63  void fillITHistos(const edm::Event& iEvent,
64  const TrackerHitAssociator& associateRecHit,
65  const std::vector<edm::Handle<edm::PSimHitContainer>>& simHits,
66  const std::map<unsigned int, SimTrack>& selectedSimTrackMap);
67 
68  void bookLayerHistos(DQMStore::IBooker& ibooker, unsigned int det_id, std::string& subdir);
69 
70  edm::ParameterSet config_;
71  TrackerHitAssociator::Config trackerHitAssociatorConfig_;
72  const double simtrackminpt_;
73  std::string geomType_;
74  const edm::EDGetTokenT<SiPixelRecHitCollection> tokenRecHitsIT_;
75  const edm::EDGetTokenT<edm::SimTrackContainer> simTracksToken_;
76  std::vector<edm::EDGetTokenT<edm::PSimHitContainer>> simHitTokens_;
77  const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> geomToken_;
78  const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> topoToken_;
79  const TrackerGeometry* tkGeom_ = nullptr;
80  const TrackerTopology* tTopo_ = nullptr;
81  struct RecHitME {
82  MonitorElement* deltaX = nullptr;
83  MonitorElement* deltaY = nullptr;
84  MonitorElement* pullX = nullptr;
85  MonitorElement* pullY = nullptr;
86  MonitorElement* deltaX_eta = nullptr;
87  MonitorElement* deltaY_eta = nullptr;
88  MonitorElement* pullX_eta = nullptr;
89  MonitorElement* pullY_eta = nullptr;
90  //For rechits matched to primary simhits
91  MonitorElement* numberRecHitsprimary = nullptr;
92  MonitorElement* pullX_primary;
93  MonitorElement* pullY_primary;
94  MonitorElement* deltaX_primary;
95  MonitorElement* deltaY_primary;
96  };
97  std::map<std::string, RecHitME> layerMEs_;
98 };
99 
100 Phase2ITValidateRecHit::Phase2ITValidateRecHit(const edm::ParameterSet& iConfig)
101  : config_(iConfig),
102  trackerHitAssociatorConfig_(iConfig, consumesCollector()),
103  simtrackminpt_(iConfig.getParameter<double>("SimTrackMinPt")),
104  tokenRecHitsIT_(consumes<SiPixelRecHitCollection>(iConfig.getParameter<edm::InputTag>("rechitsSrc"))),
105  simTracksToken_(consumes<edm::SimTrackContainer>(iConfig.getParameter<edm::InputTag>("simTracksSrc"))),
106  geomToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>()),
107  topoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd, edm::Transition::BeginRun>()) {
108  edm::LogInfo("Phase2ITValidateRecHit") << ">>> Construct Phase2ITValidateRecHit ";
109  for (const auto& itName : config_.getParameter<std::vector<std::string>>("ROUList")) {
110  simHitTokens_.push_back(consumes<std::vector<PSimHit>>(edm::InputTag("g4SimHits", itName)));
111  }
112 }
113 //
114 Phase2ITValidateRecHit::~Phase2ITValidateRecHit() {
115  edm::LogInfo("Phase2ITValidateRecHit") << ">>> Destroy Phase2ITValidateRecHit ";
116 }
117 
118 void Phase2ITValidateRecHit::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
119  std::vector<edm::Handle<edm::PSimHitContainer>> simHits;
120  for (const auto& itoken : simHitTokens_) {
121  edm::Handle<edm::PSimHitContainer> simHitHandle;
122  iEvent.getByToken(itoken, simHitHandle);
123  if (!simHitHandle.isValid())
124  continue;
125  simHits.emplace_back(simHitHandle);
126  }
127  // Get the SimTracks and push them in a map of id, SimTrack
128  edm::Handle<edm::SimTrackContainer> simTracks;
129  iEvent.getByToken(simTracksToken_, simTracks);
130 
131  std::map<unsigned int, SimTrack> selectedSimTrackMap;
132  for (const auto& simTrackIt : *simTracks) {
133  if (simTrackIt.momentum().pt() > simtrackminpt_) {
134  selectedSimTrackMap.emplace(simTrackIt.trackId(), simTrackIt);
135  }
136  }
137  TrackerHitAssociator associateRecHit(iEvent, trackerHitAssociatorConfig_);
138  fillITHistos(iEvent, associateRecHit, simHits, selectedSimTrackMap);
139 }
140 
141 void Phase2ITValidateRecHit::fillITHistos(const edm::Event& iEvent,
142  const TrackerHitAssociator& associateRecHit,
143  const std::vector<edm::Handle<edm::PSimHitContainer>>& simHits,
144  const std::map<unsigned int, SimTrack>& selectedSimTrackMap) {
145  // Get the RecHits
146  edm::Handle<SiPixelRecHitCollection> rechits;
147  iEvent.getByToken(tokenRecHitsIT_, rechits);
148  if (!rechits.isValid())
149  return;
150  std::map<std::string, unsigned int> nrechitLayerMap_primary;
151  // Loop over modules
152  for (const auto& DSViter : *rechits) {
153  // Get the detector unit's id
154  unsigned int rawid(DSViter.detId());
155  DetId detId(rawid);
156  // Get the geomdet
157  const GeomDetUnit* geomDetunit(tkGeom_->idToDetUnit(detId));
158  if (!geomDetunit)
159  continue;
160  // determine the detector we are in
161  std::string key = phase2tkutil::getITHistoId(detId.rawId(), tTopo_);
162  if (nrechitLayerMap_primary.find(key) == nrechitLayerMap_primary.end()) {
163  nrechitLayerMap_primary.emplace(key, DSViter.size());
164  } else {
165  nrechitLayerMap_primary[key] += DSViter.size();
166  }
167  //loop over rechits for a single detId
168  for (const auto& rechit : DSViter) {
169  //GetSimHits
170  const std::vector<SimHitIdpr>& matchedId = associateRecHit.associateHitId(rechit);
171  const PSimHit* simhitClosest = nullptr;
172  float minx = 10000;
173  LocalPoint lp = rechit.localPosition();
174  for (const auto& simHitCol : simHits) {
175  for (const auto& simhitIt : *simHitCol) {
176  if (detId.rawId() != simhitIt.detUnitId())
177  continue;
178  for (const auto& mId : matchedId) {
179  if (simhitIt.trackId() == mId.first) {
180  if (!simhitClosest || abs(simhitIt.localPosition().x() - lp.x()) < minx) {
181  minx = abs(simhitIt.localPosition().x() - lp.x());
182  simhitClosest = &simhitIt;
183  }
184  }
185  }
186  } //end loop over PSimhitcontainers
187  } //end loop over simHits
188  if (!simhitClosest)
189  continue;
190  auto simTrackIt(selectedSimTrackMap.find(simhitClosest->trackId()));
191  bool isPrimary = false;
192  //check if simhit is primary
193  if (simTrackIt != selectedSimTrackMap.end())
194  isPrimary = phase2tkutil::isPrimary(simTrackIt->second, simhitClosest);
195  Local3DPoint simlp(simhitClosest->localPosition());
196  const LocalError& lperr = rechit.localPositionError();
197  double dx = phase2tkutil::cmtomicron * (lp.x() - simlp.x());
198  double dy = phase2tkutil::cmtomicron * (lp.y() - simlp.y());
199  double pullx = 999.;
200  double pully = 999.;
201  if (lperr.xx())
202  pullx = (lp.x() - simlp.x()) / std::sqrt(lperr.xx());
203  if (lperr.yy())
204  pully = (lp.y() - simlp.y()) / std::sqrt(lperr.yy());
205  float eta = geomDetunit->surface().toGlobal(lp).eta();
206  layerMEs_[key].deltaX->Fill(dx);
207  layerMEs_[key].deltaY->Fill(dy);
208  layerMEs_[key].pullX->Fill(pullx);
209  layerMEs_[key].pullY->Fill(pully);
210  layerMEs_[key].deltaX_eta->Fill(eta, dx);
211  layerMEs_[key].deltaY_eta->Fill(eta, dy);
212  layerMEs_[key].pullX_eta->Fill(eta, pullx);
213  layerMEs_[key].pullY_eta->Fill(eta, pully);
214  if (isPrimary) {
215  layerMEs_[key].deltaX_primary->Fill(dx);
216  layerMEs_[key].deltaY_primary->Fill(dy);
217  layerMEs_[key].pullX_primary->Fill(pullx);
218  layerMEs_[key].pullY_primary->Fill(pully);
219  } else
220  nrechitLayerMap_primary[key]--;
221  } //end loop over rechits of a detId
222  } //End loop over DetSetVector
223 
224  //fill nRecHit counter per layer
225  for (const auto& lme : nrechitLayerMap_primary) {
226  layerMEs_[lme.first].numberRecHitsprimary->Fill(nrechitLayerMap_primary[lme.first]);
227  }
228 }
229 
230 void Phase2ITValidateRecHit::dqmBeginRun(const edm::Run& iRun, const edm::EventSetup& iSetup) {
231  edm::ESHandle<TrackerGeometry> geomHandle = iSetup.getHandle(geomToken_);
232  tkGeom_ = &(*geomHandle);
233  edm::ESHandle<TrackerTopology> tTopoHandle = iSetup.getHandle(topoToken_);
234  tTopo_ = tTopoHandle.product();
235 }
236 //
237 // -- Book Histograms
238 //
239 void Phase2ITValidateRecHit::bookHistograms(DQMStore::IBooker& ibooker,
240  edm::Run const& iRun,
241  edm::EventSetup const& iSetup) {
242  std::string top_folder = config_.getParameter<std::string>("TopFolderName");
243  edm::LogInfo("Phase2ITValidateRecHit") << " Booking Histograms in : " << top_folder;
244  edm::ESWatcher<TrackerDigiGeometryRecord> theTkDigiGeomWatcher;
245  if (theTkDigiGeomWatcher.check(iSetup)) {
246  for (auto const& det_u : tkGeom_->detUnits()) {
247  //Always check TrackerNumberingBuilder before changing this part
248  if (!(det_u->subDetector() == GeomDetEnumerators::SubDetector::P2PXB ||
249  det_u->subDetector() == GeomDetEnumerators::SubDetector::P2PXEC))
250  continue;
251  unsigned int detId_raw = det_u->geographicalId().rawId();
252  bookLayerHistos(ibooker, detId_raw, top_folder);
253  }
254  }
255 }
256 //
257 void Phase2ITValidateRecHit::bookLayerHistos(DQMStore::IBooker& ibooker, unsigned int det_id, std::string& subdir) {
258  ibooker.cd();
259  std::string key = phase2tkutil::getITHistoId(det_id, tTopo_);
260  if (key.empty())
261  return;
262  if (layerMEs_.find(key) == layerMEs_.end()) {
263  ibooker.cd();
264  RecHitME local_histos;
265  ibooker.setCurrentFolder(subdir + "/" + key);
266  edm::LogInfo("Phase2ITValidateRecHit") << " Booking Histograms in : " << (subdir + "/" + key);
267 
268  local_histos.deltaX = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaX"), ibooker);
269 
270  local_histos.deltaY = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaY"), ibooker);
271 
272  local_histos.pullX = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("PullX"), ibooker);
273 
274  local_histos.pullY = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("PullY"), ibooker);
275 
276  local_histos.deltaX_eta =
277  phase2tkutil::bookProfile1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaX_eta"), ibooker);
278 
279  local_histos.deltaY_eta =
280  phase2tkutil::bookProfile1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaY_eta"), ibooker);
281 
282  local_histos.pullX_eta =
283  phase2tkutil::bookProfile1DFromPSet(config_.getParameter<edm::ParameterSet>("PullX_eta"), ibooker);
284 
285  local_histos.pullY_eta =
286  phase2tkutil::bookProfile1DFromPSet(config_.getParameter<edm::ParameterSet>("PullY_eta"), ibooker);
287  ibooker.setCurrentFolder(subdir + "/" + key + "/PrimarySimHits");
288  //all histos for Primary particles
289  local_histos.numberRecHitsprimary =
290  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("nRecHits_primary"), ibooker);
291 
292  local_histos.deltaX_primary =
293  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaX_primary"), ibooker);
294 
295  local_histos.deltaY_primary =
296  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaY_primary"), ibooker);
297 
298  local_histos.pullX_primary =
299  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("PullX_primary"), ibooker);
300 
301  local_histos.pullY_primary =
302  phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("PullY_primary"), ibooker);
303 
304  layerMEs_.emplace(key, local_histos);
305  }
306 }
307 void Phase2ITValidateRecHit::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
308  // rechitValidIT
309  edm::ParameterSetDescription desc;
310  {
311  edm::ParameterSetDescription psd0;
312  psd0.add<std::string>("name", "Delta_X");
313  psd0.add<std::string>("title", "Delta_X;RecHit resolution X coordinate [#mum]");
314  psd0.add<bool>("switch", true);
315  psd0.add<double>("xmin", -100.0);
316  psd0.add<double>("xmax", 100.0);
317  psd0.add<int>("NxBins", 100);
318  desc.add<edm::ParameterSetDescription>("DeltaX", psd0);
319  }
320  {
321  edm::ParameterSetDescription psd0;
322  psd0.add<std::string>("name", "Delta_Y");
323  psd0.add<std::string>("title", "Delta_Y;RecHit resolution Y coordinate [#mum];");
324  psd0.add<bool>("switch", true);
325  psd0.add<double>("xmin", -100.0);
326  psd0.add<double>("xmax", 100.0);
327  psd0.add<int>("NxBins", 100);
328  desc.add<edm::ParameterSetDescription>("DeltaY", psd0);
329  }
330  {
331  edm::ParameterSetDescription psd0;
332  psd0.add<std::string>("name", "Pull_X");
333  psd0.add<std::string>("title", "Pull_X;pull x;");
334  psd0.add<double>("xmin", -4.0);
335  psd0.add<bool>("switch", true);
336  psd0.add<double>("xmax", 4.0);
337  psd0.add<int>("NxBins", 100);
338  desc.add<edm::ParameterSetDescription>("PullX", psd0);
339  }
340  {
341  edm::ParameterSetDescription psd0;
342  psd0.add<std::string>("name", "Pull_Y");
343  psd0.add<std::string>("title", "Pull_Y;pull y;");
344  psd0.add<double>("xmin", -4.0);
345  psd0.add<bool>("switch", true);
346  psd0.add<double>("xmax", 4.0);
347  psd0.add<int>("NxBins", 100);
348  desc.add<edm::ParameterSetDescription>("PullY", psd0);
349  }
350  {
351  edm::ParameterSetDescription psd0;
352  psd0.add<std::string>("name", "Delta_X_vs_Eta");
353  psd0.add<std::string>("title", "Delta_X_vs_Eta;#eta;#Delta x [#mum]");
354  psd0.add<double>("ymin", -100.0);
355  psd0.add<double>("ymax", 100.0);
356  psd0.add<int>("NxBins", 82);
357  psd0.add<bool>("switch", true);
358  psd0.add<double>("xmax", 4.1);
359  psd0.add<double>("xmin", -4.1);
360  desc.add<edm::ParameterSetDescription>("DeltaX_eta", psd0);
361  }
362  {
363  edm::ParameterSetDescription psd0;
364  psd0.add<std::string>("name", "Delta_Y_vs_Eta");
365  psd0.add<std::string>("title", "Delta_Y_vs_Eta;#eta;#Delta y [#mum]");
366  psd0.add<double>("ymin", -100.0);
367  psd0.add<double>("ymax", 100.0);
368  psd0.add<int>("NxBins", 82);
369  psd0.add<bool>("switch", true);
370  psd0.add<double>("xmax", 4.1);
371  psd0.add<double>("xmin", -4.1);
372  desc.add<edm::ParameterSetDescription>("DeltaY_eta", psd0);
373  }
374  {
375  edm::ParameterSetDescription psd0;
376  psd0.add<std::string>("name", "Pull_X_vs_Eta");
377  psd0.add<std::string>("title", "Pull_X_vs_Eta;#eta;pull x");
378  psd0.add<double>("ymax", 4.0);
379  psd0.add<int>("NxBins", 82);
380  psd0.add<bool>("switch", true);
381  psd0.add<double>("xmax", 4.1);
382  psd0.add<double>("xmin", -4.1);
383  psd0.add<double>("ymin", -4.0);
384  desc.add<edm::ParameterSetDescription>("PullX_eta", psd0);
385  }
386  {
387  edm::ParameterSetDescription psd0;
388  psd0.add<std::string>("name", "Pull_Y_vs_Eta");
389  psd0.add<std::string>("title", "Pull_Y_vs_Eta;#eta;pull y");
390  psd0.add<double>("ymax", 4.0);
391  psd0.add<int>("NxBins", 82);
392  psd0.add<bool>("switch", true);
393  psd0.add<double>("xmax", 4.1);
394  psd0.add<double>("xmin", -4.1);
395  psd0.add<double>("ymin", -4.0);
396  desc.add<edm::ParameterSetDescription>("PullY_eta", psd0);
397  }
398  //simhits primary
399  {
400  edm::ParameterSetDescription psd0;
401  psd0.add<std::string>("name", "Number_RecHits_matched_PrimarySimTrack");
402  psd0.add<std::string>("title", "Number of RecHits matched to primary SimTrack;;");
403  psd0.add<double>("xmin", 0.0);
404  psd0.add<bool>("switch", true);
405  psd0.add<double>("xmax", 0.0);
406  psd0.add<int>("NxBins", 100);
407  desc.add<edm::ParameterSetDescription>("nRecHits_primary", psd0);
408  }
409  {
410  edm::ParameterSetDescription psd0;
411  psd0.add<std::string>("name", "Delta_X_SimHitPrimary");
412  psd0.add<std::string>("title", "Delta_X_SimHitPrimary;#delta x [#mum];");
413  psd0.add<double>("xmin", -100.0);
414  psd0.add<bool>("switch", true);
415  psd0.add<double>("xmax", 100.0);
416  psd0.add<int>("NxBins", 100);
417  desc.add<edm::ParameterSetDescription>("DeltaX_primary", psd0);
418  }
419  {
420  edm::ParameterSetDescription psd0;
421  psd0.add<std::string>("name", "Delta_Y_SimHitPrimary");
422  psd0.add<std::string>("title", "Delta_Y_SimHitPrimary;#Delta y [#mum];");
423  psd0.add<double>("xmin", -100.0);
424  psd0.add<bool>("switch", true);
425  psd0.add<double>("xmax", 100.0);
426  psd0.add<int>("NxBins", 100);
427  desc.add<edm::ParameterSetDescription>("DeltaY_primary", psd0);
428  }
429  {
430  edm::ParameterSetDescription psd0;
431  psd0.add<std::string>("name", "Pull_X_SimHitPrimary");
432  psd0.add<std::string>("title", "Pull_X_SimHitPrimary;pull x;");
433  psd0.add<double>("ymax", 4.0);
434  psd0.add<int>("NxBins", 82);
435  psd0.add<bool>("switch", true);
436  psd0.add<double>("xmax", 4.1);
437  psd0.add<double>("xmin", -4.1);
438  psd0.add<double>("ymin", -4.0);
439  desc.add<edm::ParameterSetDescription>("PullX_primary", psd0);
440  }
441  {
442  edm::ParameterSetDescription psd0;
443  psd0.add<std::string>("name", "Pull_Y_SimHitPrimary");
444  psd0.add<std::string>("title", "Pull_Y_SimHitPrimary;pull y;");
445  psd0.add<double>("ymax", 4.0);
446  psd0.add<int>("NxBins", 82);
447  psd0.add<bool>("switch", true);
448  psd0.add<double>("xmax", 4.1);
449  psd0.add<double>("xmin", -4.1);
450  psd0.add<double>("ymin", -4.0);
451  desc.add<edm::ParameterSetDescription>("PullY_primary", psd0);
452  }
453  //to be used in TrackerHitAssociator
454  desc.add<bool>("associatePixel", true);
455  desc.add<bool>("associateStrip", false);
456  desc.add<bool>("usePhase2Tracker", true);
457  desc.add<bool>("associateRecoTracks", false);
458  desc.add<bool>("associateHitbySimTrack", true);
459  desc.add<edm::InputTag>("pixelSimLinkSrc", edm::InputTag("simSiPixelDigis", "Pixel"));
460  desc.add<std::vector<std::string>>("ROUList",
461  {
462  "TrackerHitsPixelBarrelLowTof",
463  "TrackerHitsPixelBarrelHighTof",
464  "TrackerHitsPixelEndcapLowTof",
465  "TrackerHitsPixelEndcapHighTof",
466  });
467  //
468  desc.add<edm::InputTag>("simTracksSrc", edm::InputTag("g4SimHits"));
469  desc.add<edm::InputTag>("SimVertexSource", edm::InputTag("g4SimHits"));
470  desc.add<double>("SimTrackMinPt", 2.0);
471  desc.add<edm::InputTag>("rechitsSrc", edm::InputTag("siPixelRecHits"));
472  desc.add<std::string>("TopFolderName", "TrackerPhase2ITRecHitV");
473  desc.add<bool>("Verbosity", false);
474  descriptions.add("Phase2ITValidateRecHit", desc);
475 }
476 //define this as a plug-in
477 DEFINE_FWK_MODULE(Phase2ITValidateRecHit);
edm::ESWatcher::check
bool check(const edm::EventSetup &iSetup)
Definition: ESWatcher.h:52
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
Phase2ITValidateRecHit::RecHitME::deltaY
MonitorElement * deltaY
Definition: Phase2ITValidateRecHit.cc:83
dqm::impl::MonitorElement
Definition: MonitorElement.h:98
phase2tkutil::book1DFromPSet
MonitorElement * book1DFromPSet(const edm::ParameterSet &hpars, DQMStore::IBooker &ibooker)
Definition: TrackerPhase2DQMUtil.cc:50
Handle.h
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
PixelSubdetector.h
Phase2ITValidateRecHit::RecHitME::deltaX_eta
MonitorElement * deltaX_eta
Definition: Phase2ITValidateRecHit.cc:86
MessageLogger.h
TrackerGeometry.h
GeomDet
Definition: GeomDet.h:27
edm::ESWatcher< TrackerDigiGeometryRecord >
ESHandle.h
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
Phase2ITValidateRecHit
Definition: Phase2ITValidateRecHit.cc:53
edm::Run
Definition: Run.h:45
edm::EDGetTokenT
Definition: EDGetToken.h:33
edm
HLT enums.
Definition: AlignableModifier.h:19
Phase2ITValidateRecHit::tokenRecHitsIT_
const edm::EDGetTokenT< SiPixelRecHitCollection > tokenRecHitsIT_
Definition: Phase2ITValidateRecHit.cc:74
TrackerTopology
Definition: TrackerTopology.h:16
Phase2ITValidateRecHit::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: Phase2ITValidateRecHit.cc:307
PSimHitContainer.h
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89281
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
Phase2ITValidateRecHit::RecHitME::pullY_primary
MonitorElement * pullY_primary
Definition: Phase2ITValidateRecHit.cc:93
TrackerPhase2ValidationUtil.h
FastTrackerRecHitCombiner_cfi.simHits
simHits
Definition: FastTrackerRecHitCombiner_cfi.py:5
TrackCandidateProducer_cfi.simTracks
simTracks
Definition: TrackCandidateProducer_cfi.py:15
DQMStore.h
Phase2ITValidateRecHit::fillITHistos
void fillITHistos(const edm::Event &iEvent, const TrackerHitAssociator &associateRecHit, const std::vector< edm::Handle< edm::PSimHitContainer >> &simHits, const std::map< unsigned int, SimTrack > &selectedSimTrackMap)
Definition: Phase2ITValidateRecHit.cc:141
Phase2ITValidateRecHit::config_
edm::ParameterSet config_
Definition: Phase2ITValidateRecHit.cc:70
Phase2ITValidateRecHit::bookHistograms
void bookHistograms(DQMStore::IBooker &ibooker, edm::Run const &iRun, edm::EventSetup const &iSetup) override
Definition: Phase2ITValidateRecHit.cc:239
TrackerHitAssociator.h
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
TrackerHitAssociator::Config
Definition: TrackerHitAssociator.h:57
edm::Handle< edm::PSimHitContainer >
Phase2ITValidateRecHit::tTopo_
const TrackerTopology * tTopo_
Definition: Phase2ITValidateRecHit.cc:80
Phase2ITValidateRecHit::RecHitME
Definition: Phase2ITValidateRecHit.cc:81
TrackerGeometry::idToDetUnit
const TrackerGeomDet * idToDetUnit(DetId) const override
Return the pointer to the GeomDetUnit corresponding to a given DetId.
Definition: TrackerGeometry.cc:183
Phase2ITValidateRecHit::RecHitME::pullY
MonitorElement * pullY
Definition: Phase2ITValidateRecHit.cc:85
DetId
Definition: DetId.h:17
Phase2ITValidateRecHit::layerMEs_
std::map< std::string, RecHitME > layerMEs_
Definition: Phase2ITValidateRecHit.cc:97
MakerMacros.h
TrackerTopology.h
Phase2ITValidateRecHit::RecHitME::deltaY_primary
MonitorElement * deltaY_primary
Definition: Phase2ITValidateRecHit.cc:95
Phase2ITValidateRecHit::RecHitME::pullX_primary
MonitorElement * pullX_primary
Definition: Phase2ITValidateRecHit.cc:92
PSimHit.h
TrackerTopologyRcd.h
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
Phase2ITValidateRecHit::RecHitME::deltaX
MonitorElement * deltaX
Definition: Phase2ITValidateRecHit.cc:82
LocalError.h
TrackerGeometry::detUnits
const DetContainer & detUnits() const override
Returm a vector of all GeomDet.
Definition: TrackerGeometry.h:61
Service.h
PVValHelper::eta
Definition: PVValidationHelpers.h:70
TrackerDigiGeometryRecord
Definition: TrackerDigiGeometryRecord.h:15
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
Phase2ITValidateRecHit::RecHitME::pullY_eta
MonitorElement * pullY_eta
Definition: Phase2ITValidateRecHit.cc:89
TrackerPhase2DQMUtil.h
edm::ESHandle< TrackerGeometry >
PSimHit::localPosition
Local3DPoint localPosition() const
Definition: PSimHit.h:52
dqm::implementation::NavigatorBase::cd
virtual void cd()
Definition: DQMStore.cc:29
Phase2ITValidateRecHit::tkGeom_
const TrackerGeometry * tkGeom_
Definition: Phase2ITValidateRecHit.cc:79
Phase2ITValidateRecHit::RecHitME::numberRecHitsprimary
MonitorElement * numberRecHitsprimary
Definition: Phase2ITValidateRecHit.cc:91
HI_PhotonSkim_cff.rechits
rechits
Definition: HI_PhotonSkim_cff.py:76
Point3DBase< float, LocalTag >
PixelGeomDetType.h
Phase2ITValidateRecHit::Phase2ITValidateRecHit
Phase2ITValidateRecHit(const edm::ParameterSet &)
Definition: Phase2ITValidateRecHit.cc:100
DQMEDAnalyzer.h
edm::ConfigurationDescriptions
Definition: ConfigurationDescriptions.h:28
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
DQMEDAnalyzer
Definition: DQMEDAnalyzer.py:1
Phase2ITValidateRecHit::simHitTokens_
std::vector< edm::EDGetTokenT< edm::PSimHitContainer > > simHitTokens_
Definition: Phase2ITValidateRecHit.cc:76
TrackerDigiGeometryRecord.h
SiPixelRecHitCollection.h
Phase2ITValidateRecHit::RecHitME::deltaY_eta
MonitorElement * deltaY_eta
Definition: Phase2ITValidateRecHit.cc:87
MonitorElement.h
edm::ParameterSet
Definition: ParameterSet.h:47
edm::Transition
Transition
Definition: Transition.h:12
Event.h
Phase2ITValidateRecHit::geomToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
Definition: Phase2ITValidateRecHit.cc:77
TrackerGeomDet.h
LocalError
Definition: LocalError.h:12
phase2tkutil::isPrimary
bool isPrimary(const SimTrack &simTrk, const PSimHit *simHit)
Definition: TrackerPhase2ValidationUtil.cc:2
Phase2ITValidateRecHit::RecHitME::deltaX_primary
MonitorElement * deltaX_primary
Definition: Phase2ITValidateRecHit.cc:94
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
Phase2ITValidateRecHit::dqmBeginRun
void dqmBeginRun(const edm::Run &, const edm::EventSetup &) override
Definition: Phase2ITValidateRecHit.cc:230
iEvent
int iEvent
Definition: GenABIO.cc:224
edm::EventSetup::getHandle
ESHandle< T > getHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:155
Phase2ITValidateRecHit::simtrackminpt_
const double simtrackminpt_
Definition: Phase2ITValidateRecHit.cc:72
trackerHitRTTI::vector
Definition: trackerHitRTTI.h:21
Phase2ITValidateRecHit::~Phase2ITValidateRecHit
~Phase2ITValidateRecHit() override
Definition: Phase2ITValidateRecHit.cc:114
Phase2ITValidateRecHit::bookLayerHistos
void bookLayerHistos(DQMStore::IBooker &ibooker, unsigned int det_id, std::string &subdir)
Definition: Phase2ITValidateRecHit.cc:257
PVValHelper::dy
Definition: PVValidationHelpers.h:50
edm::EventSetup
Definition: EventSetup.h:58
phase2tkutil::cmtomicron
static constexpr float cmtomicron
Definition: TrackerPhase2ValidationUtil.h:9
DetSetVector.h
edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >
InputTag.h
TrackerHitAssociator
Definition: TrackerHitAssociator.h:55
GeomDet.h
edmNew::DetSetVector
Definition: DetSetNew.h:13
submitPVResolutionJobs.desc
string desc
Definition: submitPVResolutionJobs.py:251
phase2tkutil::getITHistoId
std::string getITHistoId(uint32_t det_id, const TrackerTopology *tTopo)
Definition: TrackerPhase2DQMUtil.cc:2
Phase2ITValidateRecHit::trackerHitAssociatorConfig_
TrackerHitAssociator::Config trackerHitAssociatorConfig_
Definition: Phase2ITValidateRecHit.cc:71
phase2tkutil::bookProfile1DFromPSet
MonitorElement * bookProfile1DFromPSet(const edm::ParameterSet &hpars, DQMStore::IBooker &ibooker)
Definition: TrackerPhase2DQMUtil.cc:77
edm::SimTrackContainer
std::vector< SimTrack > SimTrackContainer
Definition: SimTrackContainer.h:12
Phase2ITValidateRecHit::geomType_
std::string geomType_
Definition: Phase2ITValidateRecHit.cc:73
DetId.h
Frameworkfwd.h
ESWatcher.h
LocalPoint.h
PixelGeomDetUnit.h
Phase2ITValidateRecHit::analyze
void analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup) override
Definition: Phase2ITValidateRecHit.cc:118
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
dqm::implementation::IBooker
Definition: DQMStore.h:43
PSimHit::trackId
unsigned int trackId() const
Definition: PSimHit.h:106
Phase2ITValidateRecHit::RecHitME::pullX
MonitorElement * pullX
Definition: Phase2ITValidateRecHit.cc:84
GeomDetEnumerators::P2PXB
Definition: GeomDetEnumerators.h:27
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
Phase2ITValidateRecHit::simTracksToken_
const edm::EDGetTokenT< edm::SimTrackContainer > simTracksToken_
Definition: Phase2ITValidateRecHit.cc:75
TrackerTopologyRcd
Definition: TrackerTopologyRcd.h:10
ParameterSet.h
Phase2ITValidateRecHit::RecHitME::pullX_eta
MonitorElement * pullX_eta
Definition: Phase2ITValidateRecHit.cc:88
PSimHit
Definition: PSimHit.h:15
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
DeDxTools::esConsumes
ESGetTokenH3DDVariant esConsumes(std::string const &Reccord, edm::ConsumesCollector &)
Definition: DeDxTools.cc:283
edm::Event
Definition: Event.h:73
crabWrapper.key
key
Definition: crabWrapper.py:19
TrackerHitAssociator::associateHitId
std::vector< SimHitIdpr > associateHitId(const TrackingRecHit &thit) const
Definition: TrackerHitAssociator.cc:332
SimTrackContainer.h
edm::InputTag
Definition: InputTag.h:15
PVValHelper::dx
Definition: PVValidationHelpers.h:49
Phase2ITValidateRecHit::topoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > topoToken_
Definition: Phase2ITValidateRecHit.cc:78
GeomDetEnumerators::P2PXEC
Definition: GeomDetEnumerators.h:28
TrackerGeometry
Definition: TrackerGeometry.h:14
Generated for CMSSW Reference Manual by doxygen 1.8.16