CMS 3D CMS Logo

SiPixelRecHitsValid.cc
Go to the documentation of this file.
1 // SiPixelRecHitsValid.cc
2 // Description: see SiPixelRecHitsValid.h
3 // Author: Jason Shaev, JHU
4 // Created 6/7/06
5 //
6 // G. Giurgiu, JHU (ggiurgiu@pha.jhu.edu)
7 // added pull distributions (12/27/06)
8 //--------------------------------
9 
11 
23 
25 
32 
40 
43 
44 #include <cmath>
45 
47  : tTopoEsToken(esConsumes()),
48  tGeomEsToken(esConsumes()),
49  trackerHitAssociatorConfig_(ps, consumesCollector()),
50  siPixelRecHitCollectionToken_(consumes<SiPixelRecHitCollection>(ps.getParameter<edm::InputTag>("src"))) {}
51 
53 
55  ibooker.setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Pixel/clustBPIX");
56 
57  Char_t histo[200];
58 
59  // ---------------------------------------------------------------
60  // All histograms that depend on plaquette number have 7 indexes.
61  // The first 4 (0-3) correspond to Panel 1 plaquettes 1-4.
62  // The last 3 (4-6) correspond to Panel 2 plaquettes 1-3.
63  // ---------------------------------------------------------------
64 
65  //Cluster y-size by module number for barrel
66  for (int i = 0; i < 8; i++) {
67  sprintf(histo, "Clust_y_size_Module%d", i + 1);
68  clustYSizeModule[i] = ibooker.book1D(histo, "Cluster y-size by Module", 20, 0.5, 20.5);
69  } // end for
70 
71  //Cluster x-size by layer for barrel
72  for (int i = 0; i < 3; i++) {
73  sprintf(histo, "Clust_x_size_Layer%d", i + 1);
74  clustXSizeLayer[i] = ibooker.book1D(histo, "Cluster x-size by Layer", 20, 0.5, 20.5);
75  } // end for
76 
77  //Cluster charge by module for 3 layers of barrel
78  for (int i = 0; i < 8; i++) {
79  //Cluster charge by module for Layer1
80  sprintf(histo, "Clust_charge_Layer1_Module%d", i + 1);
81  clustChargeLayer1Modules[i] = ibooker.book1D(histo, "Cluster charge Layer 1 by Module", 50, 0., 200000.);
82 
83  //Cluster charge by module for Layer2
84  sprintf(histo, "Clust_charge_Layer2_Module%d", i + 1);
85  clustChargeLayer2Modules[i] = ibooker.book1D(histo, "Cluster charge Layer 2 by Module", 50, 0., 200000.);
86 
87  //Cluster charge by module for Layer3
88  sprintf(histo, "Clust_charge_Layer3_Module%d", i + 1);
89  clustChargeLayer3Modules[i] = ibooker.book1D(histo, "Cluster charge Layer 3 by Module", 50, 0., 200000.);
90  } // end for
91 
92  ibooker.setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Pixel/clustFPIX");
93  //Cluster x-size, y-size, and charge by plaquette for Disks in Forward
94  for (int i = 0; i < 7; i++) {
95  //Cluster x-size for Disk1 by Plaquette
96  sprintf(histo, "Clust_x_size_Disk1_Plaquette%d", i + 1);
97  clustXSizeDisk1Plaquettes[i] = ibooker.book1D(histo, "Cluster X-size for Disk1 by Plaquette", 20, 0.5, 20.5);
98 
99  //Cluster x-size for Disk2 by Plaquette
100  sprintf(histo, "Clust_x_size_Disk2_Plaquette%d", i + 1);
101  clustXSizeDisk2Plaquettes[i] = ibooker.book1D(histo, "Cluster X-size for Disk2 by Plaquette", 20, 0.5, 20.5);
102 
103  //Cluster y-size for Disk1 by Plaquette
104  sprintf(histo, "Clust_y_size_Disk1_Plaquette%d", i + 1);
105  clustYSizeDisk1Plaquettes[i] = ibooker.book1D(histo, "Cluster Y-size for Disk1 by Plaquette", 20, 0.5, 20.5);
106 
107  //Cluster y-size for Disk2 by Plaquette
108  sprintf(histo, "Clust_y_size_Disk2_Plaquette%d", i + 1);
109  clustYSizeDisk2Plaquettes[i] = ibooker.book1D(histo, "Cluster Y-size for Disk2 by Plaquette", 20, 0.5, 20.5);
110 
111  //Cluster charge for Disk1 by Plaquette
112  sprintf(histo, "Clust_charge_Disk1_Plaquette%d", i + 1);
113  clustChargeDisk1Plaquettes[i] = ibooker.book1D(histo, "Cluster charge for Disk1 by Plaquette", 50, 0., 200000.);
114 
115  //Cluster charge for Disk2 by Plaquette
116  sprintf(histo, "Clust_charge_Disk2_Plaquette%d", i + 1);
117  clustChargeDisk2Plaquettes[i] = ibooker.book1D(histo, "Cluster charge for Disk2 by Plaquette", 50, 0., 200000.);
118  } // end for
119 
120  ibooker.setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Pixel/recHitBPIX");
121  //RecHit Bunch crossing all barrel hits
122  recHitBunchB = ibooker.book1D("RecHit_Bunch_Barrel", "RecHit Bunch Crossing, Barrel", 20, -10., 10.);
123 
124  //RecHit Event, in-time bunch, all barrel hits
125  recHitEventB = ibooker.book1D("RecHit_Event_Barrel", "RecHit Event (in-time bunch), Barrel", 100, 0., 100.);
126 
127  //RecHit X Resolution all barrel hits
128  recHitXResAllB = ibooker.book1D("RecHit_xres_b_All", "RecHit X Res All Modules in Barrel", 100, -200., 200.);
129 
130  //RecHit Y Resolution all barrel hits
131  recHitYResAllB = ibooker.book1D("RecHit_yres_b_All", "RecHit Y Res All Modules in Barrel", 100, -200., 200.);
132 
133  //RecHit X distribution for full modules for barrel
134  recHitXFullModules = ibooker.book1D("RecHit_x_FullModules", "RecHit X distribution for full modules", 100, -2., 2.);
135 
136  //RecHit X distribution for half modules for barrel
137  recHitXHalfModules = ibooker.book1D("RecHit_x_HalfModules", "RecHit X distribution for half modules", 100, -1., 1.);
138 
139  //RecHit Y distribution all modules for barrel
140  recHitYAllModules = ibooker.book1D("RecHit_y_AllModules", "RecHit Y distribution for all modules", 100, -4., 4.);
141 
142  //RecHit X resolution for flipped and unflipped ladders by layer for barrel
143  for (int i = 0; i < 3; i++) {
144  //RecHit no. of matched simHits all ladders by layer
145  sprintf(histo, "RecHit_NsimHit_Layer%d", i + 1);
146  recHitNsimHitLayer[i] = ibooker.book1D(histo, "RecHit Number of simHits by Layer", 30, 0., 30.);
147 
148  //RecHit X resolution for flipped ladders by layer
149  sprintf(histo, "RecHit_XRes_FlippedLadder_Layer%d", i + 1);
150  recHitXResFlippedLadderLayers[i] = ibooker.book1D(histo, "RecHit XRes Flipped Ladders by Layer", 100, -200., 200.);
151 
152  //RecHit X resolution for unflipped ladders by layer
153  sprintf(histo, "RecHit_XRes_UnFlippedLadder_Layer%d", i + 1);
155  ibooker.book1D(histo, "RecHit XRes NonFlipped Ladders by Layer", 100, -200., 200.);
156  } // end for
157 
158  //RecHit Y resolutions for layers by module for barrel
159  for (int i = 0; i < 8; i++) {
160  //Rec Hit Y resolution by module for Layer1
161  sprintf(histo, "RecHit_YRes_Layer1_Module%d", i + 1);
162  recHitYResLayer1Modules[i] = ibooker.book1D(histo, "RecHit YRes Layer1 by module", 100, -200., 200.);
163 
164  //RecHit Y resolution by module for Layer2
165  sprintf(histo, "RecHit_YRes_Layer2_Module%d", i + 1);
166  recHitYResLayer2Modules[i] = ibooker.book1D(histo, "RecHit YRes Layer2 by module", 100, -200., 200.);
167 
168  //RecHit Y resolution by module for Layer3
169  sprintf(histo, "RecHit_YRes_Layer3_Module%d", i + 1);
170  recHitYResLayer3Modules[i] = ibooker.book1D(histo, "RecHit YRes Layer3 by module", 100, -200., 200.);
171  } // end for
172 
173  ibooker.setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Pixel/recHitFPIX");
174  //RecHit Bunch crossing all plaquettes
175  recHitBunchF = ibooker.book1D("RecHit_Bunch_Forward", "RecHit Bunch Crossing, Forward", 20, -10., 10.);
176 
177  //RecHit Event, in-time bunch, all plaquettes
178  recHitEventF = ibooker.book1D("RecHit_Event_Forward", "RecHit Event (in-time bunch), Forward", 100, 0., 100.);
179 
180  //RecHit No. of simHits, by disk
181  recHitNsimHitDisk1 = ibooker.book1D("RecHit_NsimHit_Disk1", "RecHit Number of simHits, Disk1", 30, 0., 30.);
182  recHitNsimHitDisk2 = ibooker.book1D("RecHit_NsimHit_Disk2", "RecHit Number of simHits, Disk2", 30, 0., 30.);
183 
184  //RecHit X resolution all plaquettes
185  recHitXResAllF = ibooker.book1D("RecHit_xres_f_All", "RecHit X Res All in Forward", 100, -200., 200.);
186 
187  //RecHit Y resolution all plaquettes
188  recHitYResAllF = ibooker.book1D("RecHit_yres_f_All", "RecHit Y Res All in Forward", 100, -200., 200.);
189 
190  //RecHit X distribution for plaquette with x-size 1 in forward
192  ibooker.book1D("RecHit_x_Plaquette_xsize1", "RecHit X Distribution for plaquette x-size1", 100, -2., 2.);
193 
194  //RecHit X distribution for plaquette with x-size 2 in forward
196  ibooker.book1D("RecHit_x_Plaquette_xsize2", "RecHit X Distribution for plaquette x-size2", 100, -2., 2.);
197 
198  //RecHit Y distribution for plaquette with y-size 2 in forward
200  ibooker.book1D("RecHit_y_Plaquette_ysize2", "RecHit Y Distribution for plaquette y-size2", 100, -4., 4.);
201 
202  //RecHit Y distribution for plaquette with y-size 3 in forward
204  ibooker.book1D("RecHit_y_Plaquette_ysize3", "RecHit Y Distribution for plaquette y-size3", 100, -4., 4.);
205 
206  //RecHit Y distribution for plaquette with y-size 4 in forward
208  ibooker.book1D("RecHit_y_Plaquette_ysize4", "RecHit Y Distribution for plaquette y-size4", 100, -4., 4.);
209 
210  //RecHit Y distribution for plaquette with y-size 5 in forward
212  ibooker.book1D("RecHit_y_Plaquette_ysize5", "RecHit Y Distribution for plaquette y-size5", 100, -4., 4.);
213 
214  //X and Y resolutions for both disks by plaquette in forward
215  for (int i = 0; i < 7; i++) {
216  //X resolution for Disk1 by plaquette
217  sprintf(histo, "RecHit_XRes_Disk1_Plaquette%d", i + 1);
218  recHitXResDisk1Plaquettes[i] = ibooker.book1D(histo, "RecHit XRes Disk1 by plaquette", 100, -200., 200.);
219  //X resolution for Disk2 by plaquette
220  sprintf(histo, "RecHit_XRes_Disk2_Plaquette%d", i + 1);
221  recHitXResDisk2Plaquettes[i] = ibooker.book1D(histo, "RecHit XRes Disk2 by plaquette", 100, -200., 200.);
222 
223  //Y resolution for Disk1 by plaquette
224  sprintf(histo, "RecHit_YRes_Disk1_Plaquette%d", i + 1);
225  recHitYResDisk1Plaquettes[i] = ibooker.book1D(histo, "RecHit YRes Disk1 by plaquette", 100, -200., 200.);
226  //Y resolution for Disk2 by plaquette
227  sprintf(histo, "RecHit_YRes_Disk2_Plaquette%d", i + 1);
228  recHitYResDisk2Plaquettes[i] = ibooker.book1D(histo, "RecHit YRes Disk2 by plaquette", 100, -200., 200.);
229  }
230 
231  ibooker.setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Pixel/recHitPullsBPIX");
232  recHitXPullAllB = ibooker.book1D("RecHit_xres_b_All", "RecHit X Pull All Modules in Barrel", 100, -10.0, 10.0);
233  recHitYPullAllB = ibooker.book1D("RecHit_yres_b_All", "RecHit Y Pull All Modules in Barrel", 100, -10.0, 10.0);
234 
235  for (int i = 0; i < 3; i++) {
236  sprintf(histo, "RecHit_XPull_FlippedLadder_Layer%d", i + 1);
238  ibooker.book1D(histo, "RecHit XPull Flipped Ladders by Layer", 100, -10.0, 10.0);
239 
240  sprintf(histo, "RecHit_XPull_UnFlippedLadder_Layer%d", i + 1);
242  ibooker.book1D(histo, "RecHit XPull NonFlipped Ladders by Layer", 100, -10.0, 10.0);
243  }
244 
245  for (int i = 0; i < 8; i++) {
246  sprintf(histo, "RecHit_YPull_Layer1_Module%d", i + 1);
247  recHitYPullLayer1Modules[i] = ibooker.book1D(histo, "RecHit YPull Layer1 by module", 100, -10.0, 10.0);
248 
249  sprintf(histo, "RecHit_YPull_Layer2_Module%d", i + 1);
250  recHitYPullLayer2Modules[i] = ibooker.book1D(histo, "RecHit YPull Layer2 by module", 100, -10.0, 10.0);
251 
252  sprintf(histo, "RecHit_YPull_Layer3_Module%d", i + 1);
253  recHitYPullLayer3Modules[i] = ibooker.book1D(histo, "RecHit YPull Layer3 by module", 100, -10.0, 10.0);
254  }
255 
256  ibooker.setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Pixel/recHitPullsFPIX");
257  recHitXPullAllF = ibooker.book1D("RecHit_XPull_f_All", "RecHit X Pull All in Forward", 100, -10.0, 10.0);
258 
259  recHitYPullAllF = ibooker.book1D("RecHit_YPull_f_All", "RecHit Y Pull All in Forward", 100, -10.0, 10.0);
260 
261  for (int i = 0; i < 7; i++) {
262  sprintf(histo, "RecHit_XPull_Disk1_Plaquette%d", i + 1);
263  recHitXPullDisk1Plaquettes[i] = ibooker.book1D(histo, "RecHit XPull Disk1 by plaquette", 100, -10.0, 10.0);
264  sprintf(histo, "RecHit_XPull_Disk2_Plaquette%d", i + 1);
265  recHitXPullDisk2Plaquettes[i] = ibooker.book1D(histo, "RecHit XPull Disk2 by plaquette", 100, -10.0, 10.0);
266 
267  sprintf(histo, "RecHit_YPull_Disk1_Plaquette%d", i + 1);
268  recHitYPullDisk1Plaquettes[i] = ibooker.book1D(histo, "RecHit YPull Disk1 by plaquette", 100, -10.0, 10.0);
269 
270  sprintf(histo, "RecHit_YPull_Disk2_Plaquette%d", i + 1);
271  recHitYPullDisk2Plaquettes[i] = ibooker.book1D(histo, "RecHit YPull Disk2 by plaquette", 100, -10.0, 10.0);
272  }
273 }
274 
276  //Retrieve tracker topology from geometry
277  const TrackerTopology* tTopo = &es.getData(tTopoEsToken);
278 
279  edm::LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();
280  if (e.id().event() % 1000 == 0)
281  std::cout << " Run = " << e.id().run() << " Event = " << e.id().event() << std::endl;
282 
283  //Get RecHits
285  e.getByToken(siPixelRecHitCollectionToken_, recHitColl);
286 
287  //Get event setup
288  const TrackerGeometry* theTracker = &es.getData(tGeomEsToken);
289 
291 
292  //iterate over detunits
293  //for (TrackerGeometry::DetContainer::const_iterator it = geom->dets().begin(); it != geom->dets().end(); it++) {
294  for (const auto& it : theTracker->dets()) {
295  DetId detId = it->geographicalId();
296  unsigned int subid = detId.subdetId();
297 
298  if (!((subid == 1) || (subid == 2)))
299  continue;
300 
301  const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*>(theTracker->idToDet(detId));
302 
303  SiPixelRecHitCollection::const_iterator pixeldet = recHitColl->find(detId);
304  if (pixeldet == recHitColl->end())
305  continue;
306  SiPixelRecHitCollection::DetSet pixelrechitRange = *pixeldet;
307  SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorBegin = pixelrechitRange.begin();
308  SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorEnd = pixelrechitRange.end();
309  SiPixelRecHitCollection::DetSet::const_iterator pixeliter = pixelrechitRangeIteratorBegin;
310  std::vector<PSimHit> matched;
311 
312  //----Loop over rechits for this detId
313  for (; pixeliter != pixelrechitRangeIteratorEnd; pixeliter++) {
314  matched.clear();
315  matched = associate.associateHit(*pixeliter);
316 
317  if (!matched.empty()) {
318  float closest = 9999.9;
319  std::vector<PSimHit>::const_iterator closestit = matched.begin();
320  LocalPoint lp = pixeliter->localPosition();
321  float rechit_x = lp.x();
322  float rechit_y = lp.y();
323 
324  //loop over sim hits and fill closet
325  for (std::vector<PSimHit>::const_iterator m = matched.begin(); m < matched.end(); m++) {
326  float sim_x1 = (*m).entryPoint().x();
327  float sim_x2 = (*m).exitPoint().x();
328  float sim_xpos = 0.5 * (sim_x1 + sim_x2);
329 
330  float sim_y1 = (*m).entryPoint().y();
331  float sim_y2 = (*m).exitPoint().y();
332  float sim_ypos = 0.5 * (sim_y1 + sim_y2);
333 
334  float x_res = fabs(sim_xpos - rechit_x);
335  float y_res = fabs(sim_ypos - rechit_y);
336 
337  float dist = sqrt(x_res * x_res + y_res * y_res);
338 
339  if (dist < closest) {
340  closest = dist;
341  closestit = m;
342  }
343  } // end sim hit loop
344 
345  if (subid == 1) { //<----------barrel
346  fillBarrel(*pixeliter, *closestit, detId, theGeomDet, tTopo);
347  } // end barrel
348  if (subid == 2) { // <-------forward
349  fillForward(*pixeliter, *closestit, detId, theGeomDet, tTopo);
350  }
351 
352  } // end matched emtpy
353 
354  int NsimHit = matched.size();
355  if (subid == 1) { //<----------barrel
356  for (unsigned int i = 0; i < 3; i++)
357  if (tTopo->pxbLayer(detId) == i + 1)
358  recHitNsimHitLayer[i]->Fill(NsimHit);
359  } // end barrel
360  if (subid == 2) { // <-------forward
361  if (tTopo->pxfDisk(detId) == 1)
362  recHitNsimHitDisk1->Fill(NsimHit);
363  else
364  recHitNsimHitDisk2->Fill(NsimHit);
365  }
366  } // <-----end rechit loop
367  } // <------ end detunit loop
368 }
369 
371  const PSimHit& simHit,
372  DetId detId,
373  const PixelGeomDetUnit* theGeomDet,
374  const TrackerTopology* tTopo) {
375  const float cmtomicron = 10000.0;
376 
377  int bunch = simHit.eventId().bunchCrossing();
378  int event = simHit.eventId().event();
379 
380  recHitBunchB->Fill(bunch);
381  if (bunch == 0)
383 
384  LocalPoint lp = recHit.localPosition();
385  float lp_y = lp.y();
386  float lp_x = lp.x();
387 
388  LocalError lerr = recHit.localPositionError();
389  float lerr_x = sqrt(lerr.xx());
390  float lerr_y = sqrt(lerr.yy());
391 
392  recHitYAllModules->Fill(lp_y);
393 
394  float sim_x1 = simHit.entryPoint().x();
395  float sim_x2 = simHit.exitPoint().x();
396  float sim_xpos = 0.5 * (sim_x1 + sim_x2);
397  float res_x = (lp.x() - sim_xpos) * cmtomicron;
398 
399  recHitXResAllB->Fill(res_x);
400 
401  float sim_y1 = simHit.entryPoint().y();
402  float sim_y2 = simHit.exitPoint().y();
403  float sim_ypos = 0.5 * (sim_y1 + sim_y2);
404  float res_y = (lp.y() - sim_ypos) * cmtomicron;
405 
406  recHitYResAllB->Fill(res_y);
407 
408  float pull_x = (lp_x - sim_xpos) / lerr_x;
409  float pull_y = (lp_y - sim_ypos) / lerr_y;
410 
411  recHitXPullAllB->Fill(pull_x);
412  recHitYPullAllB->Fill(pull_y);
413 
414  int rows = theGeomDet->specificTopology().nrows();
415 
416  if (rows == 160) {
417  recHitXFullModules->Fill(lp_x);
418  } else if (rows == 80) {
419  recHitXHalfModules->Fill(lp_x);
420  }
421 
422  float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0., 0., 0.)).perp();
423  float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0., 0., 1.)).perp();
424 
425  if (tmp2 < tmp1) { // flipped
426  for (unsigned int i = 0; i < 3; i++) {
427  if (tTopo->pxbLayer(detId) == i + 1) {
430  }
431  }
432  } else {
433  for (unsigned int i = 0; i < 3; i++) {
434  if (tTopo->pxbLayer(detId) == i + 1) {
437  }
438  }
439  }
440 
441  //get cluster
442  SiPixelRecHit::ClusterRef const& clust = recHit.cluster();
443 
444  // fill module dependent info
445  for (unsigned int i = 0; i < 8; i++) {
446  if (tTopo->pxbModule(detId) == i + 1) {
447  int sizeY = (*clust).sizeY();
449 
450  if (tTopo->pxbLayer(detId) == 1) {
451  float charge = (*clust).charge();
453  recHitYResLayer1Modules[i]->Fill(res_y);
454  recHitYPullLayer1Modules[i]->Fill(pull_y);
455  } else if (tTopo->pxbLayer(detId) == 2) {
456  float charge = (*clust).charge();
458  recHitYResLayer2Modules[i]->Fill(res_y);
459  recHitYPullLayer2Modules[i]->Fill(pull_y);
460  } else if (tTopo->pxbLayer(detId) == 3) {
461  float charge = (*clust).charge();
463  recHitYResLayer3Modules[i]->Fill(res_y);
464  recHitYPullLayer3Modules[i]->Fill(pull_y);
465  }
466  }
467  }
468  int sizeX = (*clust).sizeX();
469  if (tTopo->pxbLayer(detId) == 1)
471  if (tTopo->pxbLayer(detId) == 2)
473  if (tTopo->pxbLayer(detId) == 3)
475 }
476 
478  const PSimHit& simHit,
479  DetId detId,
480  const PixelGeomDetUnit* theGeomDet,
481  const TrackerTopology* tTopo) {
482  int rows = theGeomDet->specificTopology().nrows();
483  int cols = theGeomDet->specificTopology().ncolumns();
484 
485  const float cmtomicron = 10000.0;
486 
487  int bunch = simHit.eventId().bunchCrossing();
488  int event = simHit.eventId().event();
489 
490  recHitBunchF->Fill(bunch);
491  if (bunch == 0)
493 
494  LocalPoint lp = recHit.localPosition();
495  float lp_x = lp.x();
496  float lp_y = lp.y();
497 
498  LocalError lerr = recHit.localPositionError();
499  float lerr_x = sqrt(lerr.xx());
500  float lerr_y = sqrt(lerr.yy());
501 
502  float sim_x1 = simHit.entryPoint().x();
503  float sim_x2 = simHit.exitPoint().x();
504  float sim_xpos = 0.5 * (sim_x1 + sim_x2);
505 
506  float sim_y1 = simHit.entryPoint().y();
507  float sim_y2 = simHit.exitPoint().y();
508  float sim_ypos = 0.5 * (sim_y1 + sim_y2);
509 
510  float pull_x = (lp_x - sim_xpos) / lerr_x;
511  float pull_y = (lp_y - sim_ypos) / lerr_y;
512 
513  if (rows == 80) {
515  } else if (rows == 160) {
517  }
518 
519  if (cols == 104) {
521  } else if (cols == 156) {
523  } else if (cols == 208) {
525  } else if (cols == 260) {
527  }
528 
529  float res_x = (lp.x() - sim_xpos) * cmtomicron;
530 
531  recHitXResAllF->Fill(res_x);
532  recHitXPullAllF->Fill(pull_x);
533 
534  float res_y = (lp.y() - sim_ypos) * cmtomicron;
535 
536  recHitYPullAllF->Fill(pull_y);
537 
538  // get cluster
539  SiPixelRecHit::ClusterRef const& clust = recHit.cluster();
540 
541  // fill plaquette dependent info
542  for (unsigned int i = 0; i < 7; i++) {
543  if (tTopo->pxfModule(detId) == i + 1) {
544  if (tTopo->pxfDisk(detId) == 1) {
545  int sizeX = (*clust).sizeX();
547 
548  int sizeY = (*clust).sizeY();
550 
551  float charge = (*clust).charge();
553 
556 
559  } else {
560  int sizeX = (*clust).sizeX();
562 
563  int sizeY = (*clust).sizeY();
565 
566  float charge = (*clust).charge();
568 
571 
574 
575  } // end else
576  } // end if module
577  else if (tTopo->pxfPanel(detId) == 2 && (tTopo->pxfModule(detId) + 4) == i + 1) {
578  if (tTopo->pxfDisk(detId) == 1) {
579  int sizeX = (*clust).sizeX();
581 
582  int sizeY = (*clust).sizeY();
584 
585  float charge = (*clust).charge();
587 
590 
593  } else {
594  int sizeX = (*clust).sizeX();
596 
597  int sizeY = (*clust).sizeY();
599 
600  float charge = (*clust).charge();
602 
605 
608 
609  } // end else
610  } // end else
611  } // end for
612 }
SiPixelRecHitsValid::bookHistograms
void bookHistograms(DQMStore::IBooker &ibooker, const edm::Run &run, const edm::EventSetup &es) override
Definition: SiPixelRecHitsValid.cc:54
muonTagProbeFilters_cff.matched
matched
Definition: muonTagProbeFilters_cff.py:62
TrackerGeometry::idToDet
const TrackerGeomDet * idToDet(DetId) const override
Definition: TrackerGeometry.cc:193
SiPixelRecHitsValid::clustYSizeDisk1Plaquettes
MonitorElement * clustYSizeDisk1Plaquettes[7]
Definition: SiPixelRecHitsValid.h:60
SiPixelRecHitsValid::recHitXPullDisk1Plaquettes
MonitorElement * recHitXPullDisk1Plaquettes[7]
Definition: SiPixelRecHitsValid.h:113
Handle.h
SiPixelRecHitsValid::recHitXPullAllB
MonitorElement * recHitXPullAllB
Definition: SiPixelRecHitsValid.h:100
mps_fire.i
i
Definition: mps_fire.py:428
PixelSubdetector.h
MessageLogger.h
SiPixelRecHitsValid::recHitBunchB
MonitorElement * recHitBunchB
Definition: SiPixelRecHitsValid.h:76
TrackerGeometry.h
SiPixelRecHitsValid::clustChargeLayer2Modules
MonitorElement * clustChargeLayer2Modules[8]
Definition: SiPixelRecHitsValid.h:54
PixelTopology.h
ESHandle.h
SiPixelRecHitsValid::fillBarrel
void fillBarrel(const SiPixelRecHit &, const PSimHit &, DetId, const PixelGeomDetUnit *, const TrackerTopology *tTopo)
Definition: SiPixelRecHitsValid.cc:370
SiPixelRecHitsValid::recHitYResLayer2Modules
MonitorElement * recHitYResLayer2Modules[8]
Definition: SiPixelRecHitsValid.h:74
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
edm::Run
Definition: Run.h:45
SiPixelRecHitsValid::tTopoEsToken
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > tTopoEsToken
Definition: SiPixelRecHitsValid.h:47
edm
HLT enums.
Definition: AlignableModifier.h:19
TrackerTopology
Definition: TrackerTopology.h:16
gather_cfg.cout
cout
Definition: gather_cfg.py:144
PSimHitContainer.h
SiPixelCluster.h
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89301
TrackerHitAssociator::associateHit
std::vector< PSimHit > associateHit(const TrackingRecHit &thit) const
Definition: TrackerHitAssociator.cc:212
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
SiPixelRecHitsValid::recHitXPullDisk2Plaquettes
MonitorElement * recHitXPullDisk2Plaquettes[7]
Definition: SiPixelRecHitsValid.h:114
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:197
GeomDetType.h
timingPdfMaker.histo
histo
Definition: timingPdfMaker.py:278
DQMStore.h
SiPixelRecHitsValid::tGeomEsToken
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > tGeomEsToken
Definition: SiPixelRecHitsValid.h:48
SiPixelRecHitsValid::recHitXResAllB
MonitorElement * recHitXResAllB
Definition: SiPixelRecHitsValid.h:66
SiPixelRecHitsValid::clustYSizeDisk2Plaquettes
MonitorElement * clustYSizeDisk2Plaquettes[7]
Definition: SiPixelRecHitsValid.h:61
SiPixelRecHitsValid::recHitYPullLayer2Modules
MonitorElement * recHitYPullLayer2Modules[8]
Definition: SiPixelRecHitsValid.h:106
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
SiPixelRecHitsValid::recHitXPullAllF
MonitorElement * recHitXPullAllF
Definition: SiPixelRecHitsValid.h:110
edm::Handle
Definition: AssociativeIterator.h:50
SiPixelRecHitsValid::clustChargeDisk1Plaquettes
MonitorElement * clustChargeDisk1Plaquettes[7]
Definition: SiPixelRecHitsValid.h:62
SiPixelRecHitsValid::clustChargeDisk2Plaquettes
MonitorElement * clustChargeDisk2Plaquettes[7]
Definition: SiPixelRecHitsValid.h:63
rpcPointValidation_cfi.recHit
recHit
Definition: rpcPointValidation_cfi.py:7
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:23
SiPixelRecHitsValid::recHitXPullFlippedLadderLayers
MonitorElement * recHitXPullFlippedLadderLayers[3]
Definition: SiPixelRecHitsValid.h:103
edm::Ref
Definition: AssociativeIterator.h:58
SiPixelRecHitsValid::recHitYAllModules
MonitorElement * recHitYAllModules
Definition: SiPixelRecHitsValid.h:70
TrackerTopology::pxbLayer
unsigned int pxbLayer(const DetId &id) const
Definition: TrackerTopology.h:144
TrackerTopology::pxfPanel
unsigned int pxfPanel(const DetId &id) const
Definition: TrackerTopology.h:450
SiPixelRecHitsValid::recHitYPullLayer1Modules
MonitorElement * recHitYPullLayer1Modules[8]
Definition: SiPixelRecHitsValid.h:105
DetId
Definition: DetId.h:17
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
edmNew::DetSet::end
iterator end()
Definition: DetSetNew.h:56
MakerMacros.h
TrackerTopology.h
edmNew::DetSetVector::find
const_iterator find(id_type i, bool update=false) const
Definition: DetSetVectorNew.h:533
SiPixelRecHitsValid::recHitYPullLayer3Modules
MonitorElement * recHitYPullLayer3Modules[8]
Definition: SiPixelRecHitsValid.h:107
PSimHit.h
SiPixelRecHitsValid::recHitXResFlippedLadderLayers
MonitorElement * recHitXResFlippedLadderLayers[3]
Definition: SiPixelRecHitsValid.h:71
SiPixelRecHitsValid::recHitNsimHitDisk1
MonitorElement * recHitNsimHitDisk1
Definition: SiPixelRecHitsValid.h:95
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
LocalError::xx
float xx() const
Definition: LocalError.h:22
SiPixelRecHitsValid::recHitYResLayer3Modules
MonitorElement * recHitYResLayer3Modules[8]
Definition: SiPixelRecHitsValid.h:75
Service.h
SiPixelRecHitsValid::SiPixelRecHitsValid
SiPixelRecHitsValid(const edm::ParameterSet &conf)
Definition: SiPixelRecHitsValid.cc:46
rpcPointValidation_cfi.simHit
simHit
Definition: rpcPointValidation_cfi.py:24
visualization-live-secondInstance_cfg.m
m
Definition: visualization-live-secondInstance_cfg.py:78
SiPixelRecHitsValid::fillForward
void fillForward(const SiPixelRecHit &, const PSimHit &, DetId, const PixelGeomDetUnit *, const TrackerTopology *tTopo)
Definition: SiPixelRecHitsValid.cc:477
PixelGeomDetUnit
Definition: PixelGeomDetUnit.h:15
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:290
Surface::toGlobal
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:79
SiPixelRecHitsValid::recHitXPlaquetteSize1
MonitorElement * recHitXPlaquetteSize1
Definition: SiPixelRecHitsValid.h:83
edmNew::DetSet
Definition: DetSetNew.h:22
PixelTopology::ncolumns
virtual int ncolumns() const =0
RecoTauValidation_cfi.sizeY
sizeY
Definition: RecoTauValidation_cfi.py:290
SiPixelRecHitsValid::recHitYPlaquetteSize3
MonitorElement * recHitYPlaquetteSize3
Definition: SiPixelRecHitsValid.h:86
SiPixelRecHitsValid::recHitEventF
MonitorElement * recHitEventF
Definition: SiPixelRecHitsValid.h:94
SiPixelRecHitsValid::recHitYPlaquetteSize4
MonitorElement * recHitYPlaquetteSize4
Definition: SiPixelRecHitsValid.h:87
Point3DBase< float, LocalTag >
SiPixelRecHit.h
SiPixelRecHitsValid::recHitNsimHitDisk2
MonitorElement * recHitNsimHitDisk2
Definition: SiPixelRecHitsValid.h:96
PixelGeomDetType.h
SiPixelRecHitsValid::clustXSizeDisk2Plaquettes
MonitorElement * clustXSizeDisk2Plaquettes[7]
Definition: SiPixelRecHitsValid.h:59
SiPixelRecHitsValid::recHitYResAllB
MonitorElement * recHitYResAllB
Definition: SiPixelRecHitsValid.h:67
DetId::subdetId
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum)
Definition: DetId.h:48
SiPixelRecHitsValid::recHitYPullDisk2Plaquettes
MonitorElement * recHitYPullDisk2Plaquettes[7]
Definition: SiPixelRecHitsValid.h:116
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
SiPixelRecHitsValid::trackerHitAssociatorConfig_
TrackerHitAssociator::Config trackerHitAssociatorConfig_
Definition: SiPixelRecHitsValid.h:118
SiPixelRecHitsValid::recHitYResLayer1Modules
MonitorElement * recHitYResLayer1Modules[8]
Definition: SiPixelRecHitsValid.h:73
SiPixelRecHitsValid::clustYSizeModule
MonitorElement * clustYSizeModule[8]
Definition: SiPixelRecHitsValid.h:51
edm::ParameterSet
Definition: ParameterSet.h:47
SiPixelRecHitsValid::clustXSizeDisk1Plaquettes
MonitorElement * clustXSizeDisk1Plaquettes[7]
Definition: SiPixelRecHitsValid.h:58
RecoTauValidation_cfi.sizeX
sizeX
Definition: RecoTauValidation_cfi.py:289
Event.h
LocalError
Definition: LocalError.h:12
SiPixelRecHitsValid::recHitXResAllF
MonitorElement * recHitXResAllF
Definition: SiPixelRecHitsValid.h:81
PixelGeomDetUnit::specificTopology
virtual const PixelTopology & specificTopology() const
Returns a reference to the pixel proxy topology.
Definition: PixelGeomDetUnit.cc:17
SiPixelRecHitsValid::recHitYPullAllB
MonitorElement * recHitYPullAllB
Definition: SiPixelRecHitsValid.h:101
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
SiPixelRecHitsValid::recHitXPullNonFlippedLadderLayers
MonitorElement * recHitXPullNonFlippedLadderLayers[3]
Definition: SiPixelRecHitsValid.h:104
TrackerTopology::pxfModule
unsigned int pxfModule(const DetId &id) const
Definition: TrackerTopology.h:163
TrackerTopology::pxfDisk
unsigned int pxfDisk(const DetId &id) const
Definition: TrackerTopology.h:446
SiPixelRecHitsValid::recHitXResDisk1Plaquettes
MonitorElement * recHitXResDisk1Plaquettes[7]
Definition: SiPixelRecHitsValid.h:89
edmNew::DetSet::begin
iterator begin()
Definition: DetSetNew.h:54
SiPixelRecHitsValid::recHitBunchF
MonitorElement * recHitBunchF
Definition: SiPixelRecHitsValid.h:93
edm::EventSetup
Definition: EventSetup.h:58
TrackerGeometry::dets
const DetContainer & dets() const override
Returm a vector of all GeomDet (including all GeomDetUnits)
Definition: TrackerGeometry.h:62
DetSetVector.h
SiPixelRecHitsValid.h
GeometricDet.h
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:1053
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:127
SiPixelRecHitsValid::recHitYPlaquetteSize5
MonitorElement * recHitYPlaquetteSize5
Definition: SiPixelRecHitsValid.h:88
SiPixelRecHitsValid::clustXSizeLayer
MonitorElement * clustXSizeLayer[3]
Definition: SiPixelRecHitsValid.h:52
SiPixelRecHitsValid::recHitYResDisk2Plaquettes
MonitorElement * recHitYResDisk2Plaquettes[7]
Definition: SiPixelRecHitsValid.h:92
TrackerHitAssociator
Definition: TrackerHitAssociator.h:55
GeomDet.h
edmNew::DetSetVector
Definition: DetSetNew.h:13
SiPixelRecHitsValid::recHitEventB
MonitorElement * recHitEventB
Definition: SiPixelRecHitsValid.h:77
SiPixelRecHitsValid::recHitXPlaquetteSize2
MonitorElement * recHitXPlaquetteSize2
Definition: SiPixelRecHitsValid.h:84
writedatasetfile.run
run
Definition: writedatasetfile.py:27
Ref.h
SiPixelRecHitsValid::recHitYPlaquetteSize2
MonitorElement * recHitYPlaquetteSize2
Definition: SiPixelRecHitsValid.h:85
SiPixelRecHitsValid::recHitXFullModules
MonitorElement * recHitXFullModules
Definition: SiPixelRecHitsValid.h:68
SiPixelRecHitsValid::~SiPixelRecHitsValid
~SiPixelRecHitsValid() override
Definition: SiPixelRecHitsValid.cc:52
SiPixelRecHitsValid::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c) override
Definition: SiPixelRecHitsValid.cc:275
TrackerTopology::pxbModule
unsigned int pxbModule(const DetId &id) const
Definition: TrackerTopology.h:160
DetId.h
SiPixelRecHitsValid::recHitXHalfModules
MonitorElement * recHitXHalfModules
Definition: SiPixelRecHitsValid.h:69
LocalPoint.h
edmNew::DetSetVector::end
const_iterator end(bool update=false) const
Definition: DetSetVectorNew.h:550
SiPixelRecHitsValid::recHitYPullAllF
MonitorElement * recHitYPullAllF
Definition: SiPixelRecHitsValid.h:111
PixelGeomDetUnit.h
EventSetup.h
dqm::implementation::IBooker
Definition: DQMStore.h:43
SiPixelRecHitsValid::recHitYPullDisk1Plaquettes
MonitorElement * recHitYPullDisk1Plaquettes[7]
Definition: SiPixelRecHitsValid.h:115
SiPixelRecHitsValid::recHitNsimHitLayer
MonitorElement * recHitNsimHitLayer[3]
Definition: SiPixelRecHitsValid.h:78
PSimHit
Definition: PSimHit.h:15
SiPixelRecHitsValid::clustChargeLayer1Modules
MonitorElement * clustChargeLayer1Modules[8]
Definition: SiPixelRecHitsValid.h:53
OwnVector.h
SiPixelRecHitsValid
Definition: SiPixelRecHitsValid.h:30
event
Definition: event.py:1
DeDxTools::esConsumes
ESGetTokenH3DDVariant esConsumes(std::string const &Reccord, edm::ConsumesCollector &)
Definition: DeDxTools.cc:283
edm::Event
Definition: Event.h:73
PixelTopology::nrows
virtual int nrows() const =0
postprocess-scan-build.rows
rows
Definition: postprocess-scan-build.py:23
SiPixelRecHitsValid::recHitXResNonFlippedLadderLayers
MonitorElement * recHitXResNonFlippedLadderLayers[3]
Definition: SiPixelRecHitsValid.h:72
SiPixelRecHitsValid::recHitYResDisk1Plaquettes
MonitorElement * recHitYResDisk1Plaquettes[7]
Definition: SiPixelRecHitsValid.h:91
GlobalPoint.h
LocalError::yy
float yy() const
Definition: LocalError.h:24
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:69
SiPixelRecHitsValid::clustChargeLayer3Modules
MonitorElement * clustChargeLayer3Modules[8]
Definition: SiPixelRecHitsValid.h:55
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
SiPixelRecHitsValid::siPixelRecHitCollectionToken_
edm::EDGetTokenT< SiPixelRecHitCollection > siPixelRecHitCollectionToken_
Definition: SiPixelRecHitsValid.h:119
Local3DPoint
Point3DBase< float, LocalTag > Local3DPoint
Definition: LocalPoint.h:9
cmtomicron
const float cmtomicron
Definition: PixelResolutionHistograms.cc:54
SiPixelRecHitsValid::recHitXResDisk2Plaquettes
MonitorElement * recHitXResDisk2Plaquettes[7]
Definition: SiPixelRecHitsValid.h:90
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
TrackerGeometry
Definition: TrackerGeometry.h:14
SiPixelRecHitsValid::recHitYResAllF
MonitorElement * recHitYResAllF
Definition: SiPixelRecHitsValid.h:82