CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
SiPixelTrackingRecHitsValid.cc
Go to the documentation of this file.
1 
2 // File: SiPixelTrackingRecHitsValid.cc
3 // Authors: Xingtao Huang (Puerto Rico Univ.)
4 // Gavril Giurgiu (JHU)
5 // Creation Date: Oct. 2006
6 
7 #include <memory>
8 #include <string>
9 #include <iostream>
10 #include <TMath.h>
11 #include "Validation/RecoTrack/interface/SiPixelTrackingRecHitsValid.h"
12 
13 #include "DataFormats/TrajectorySeed/interface/TrajectorySeedCollection.h"
14 #include "FWCore/Framework/interface/ESHandle.h"
15 #include "DataFormats/GeometryVector/interface/GlobalPoint.h"
16 #include "DataFormats/GeometryVector/interface/GlobalVector.h"
17 #include "DataFormats/GeometryVector/interface/LocalVector.h"
18 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
19 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
20 #include "TrackingTools/Records/interface/TransientRecHitRecord.h"
21 #include "Geometry/TrackerGeometryBuilder/interface/GluedGeomDet.h"
22 
23 #include "DataFormats/TrackReco/interface/Track.h"
24 #include "DataFormats/TrackReco/interface/TrackExtra.h"
25 #include "DataFormats/DetId/interface/DetId.h"
26 #include "DataFormats/GeometryVector/interface/LocalPoint.h"
27 #include "DataFormats/GeometryVector/interface/GlobalPoint.h"
28 
29 #include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHit.h"
30 #include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
31 #include "Geometry/Records/interface/IdealGeometryRecord.h"
32 #include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
33 
34 #include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
35 
36 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
37 
38 #include "Geometry/TrackerGeometryBuilder/interface/PixelTopologyBuilder.h"
39 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
40 #include "Geometry/CommonDetUnit/interface/GeomDetType.h"
41 #include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"
42 
43 
44 #include <TTree.h>
45 #include <TFile.h>
46 
47 // End job: write and close the ntuple file
48 void SiPixelTrackingRecHitsValid::endJob()
49 {
50  if(debugNtuple_.size()!=0){
51  tfile_->Write();
52  tfile_->Close();
53  }
54 }
55 
56 void SiPixelTrackingRecHitsValid::beginJob()
57 {
58 
59  if(debugNtuple_.size()!=0){
60  tfile_ = new TFile (debugNtuple_.c_str() , "RECREATE");
61 
62  t_ = new TTree("Ntuple", "Ntuple");
63  int bufsize = 64000;
64 
65  t_->Branch("subdetId", &subdetId, "subdetId/I", bufsize);
66 
67  t_->Branch("layer" , &layer , "layer/I" , bufsize);
68  t_->Branch("ladder", &ladder, "ladder/I", bufsize);
69  t_->Branch("mod" , &mod , "mod/I" , bufsize);
70  t_->Branch("side" , &side , "side/I" , bufsize);
71  t_->Branch("disk" , &disk , "disk/I" , bufsize);
72  t_->Branch("blade" , &blade , "blade/I" , bufsize);
73  t_->Branch("panel" , &panel , "panel/I" , bufsize);
74  t_->Branch("plaq" , &plaq , "plaq/I" , bufsize);
75 
76  t_->Branch("rechitx" , &rechitx , "rechitx/F" , bufsize);
77  t_->Branch("rechity" , &rechity , "rechity/F" , bufsize);
78  t_->Branch("rechitz" , &rechitz , "rechitz/F" , bufsize);
79  t_->Branch("rechiterrx" , &rechiterrx , "rechiterrx/F" , bufsize);
80  t_->Branch("rechiterry" , &rechiterry , "rechiterry/F" , bufsize);
81  t_->Branch("rechitresx" , &rechitresx , "rechitresx/F" , bufsize);
82  t_->Branch("rechitresy" , &rechitresy , "rechitresy/F" , bufsize);
83  t_->Branch("rechitpullx", &rechitpullx, "rechitpullx/F", bufsize);
84  t_->Branch("rechitpully", &rechitpully, "rechitpully/F", bufsize);
85 
86  t_->Branch("npix" , &npix , "npix/I" , bufsize);
87  t_->Branch("nxpix" , &nxpix , "nxpix/I" , bufsize);
88  t_->Branch("nypix" , &nypix , "nypix/I" , bufsize);
89  t_->Branch("charge", &charge, "charge/F", bufsize);
90 
91  t_->Branch("alpha", &alpha, "alpha/F", bufsize);
92  t_->Branch("beta" , &beta , "beta/F" , bufsize);
93 
94  t_->Branch("phi", &phi, "phi/F", bufsize);
95  t_->Branch("eta", &eta, "eta/F", bufsize);
96 
97  t_->Branch("half" , &half , "half/I" , bufsize);
98  t_->Branch("flipped", &flipped, "flipped/I", bufsize);
99 
100  t_->Branch("simhitx", &simhitx, "simhitx/F", bufsize);
101  t_->Branch("simhity", &simhity, "simhity/F", bufsize);
102 
103  t_->Branch("nsimhit", &nsimhit, "nsimhit/I", bufsize);
104  t_->Branch("pidhit" , &pidhit , "pidhit/I" , bufsize);
105 
106  t_->Branch("evt", &evt, "evt/I", bufsize);
107  t_->Branch("run", &run, "run/I", bufsize);
108  }
109 
110 }
111 
112 SiPixelTrackingRecHitsValid::SiPixelTrackingRecHitsValid(const edm::ParameterSet& ps):conf_(ps), dbe_(0), tfile_(0), t_(0)
113 {
114  //Read config file
115  MTCCtrack_ = ps.getParameter<bool>("MTCCtrack");
116  outputFile_ = ps.getUntrackedParameter<std::string>("outputFile", "pixeltrackingrechitshisto.root");
117  siPixelRecHitCollectionToken_ = consumes<SiPixelRecHitCollection>( edm::InputTag( "siPixelRecHits" ) );
118  recoTrackCollectionToken_ = consumes<reco::TrackCollection>( edm::InputTag( ps.getUntrackedParameter<std::string>( "src" ) ) );
119  builderName_ = ps.getParameter<std::string>("TTRHBuilder");
120  checkType_ = ps.getParameter<bool>("checkType");
121  genType_ = ps.getParameter<int>("genType");
122  debugNtuple_=ps.getUntrackedParameter<std::string>("debugNtuple", "SiPixelTrackingRecHitsValid_Ntuple.root");
123 
124  // Book histograms
125  dbe_ = edm::Service<DQMStore>().operator->();
126  //dbe_->showDirStructure();
127 
128  //float math_pi = 3.14159265;
129  //float radtodeg = 180.0 / math_pi;
130 
131  // Histogram ranges (low and high)
132 
133  float xl = -1.0;
134  float xh = 1.0;
135  float errxl = 0.0;
136  float errxh = 0.003;
137  float resxl = -0.02;
138  float resxh = 0.02;
139  float pullxl = -10.0;
140  float pullxh = 10.0;
141 
142  float yl = -4.0;
143  float yh = 4.0;
144  float erryl = 0.0;
145  float erryh = 0.010;
146  float resyl = -0.04;
147  float resyh = 0.04;
148  float pullyl = -10.0;
149  float pullyh = 10.0;
150 
151  float barrel_alphal = 80.0;
152  float barrel_alphah = 100.0;
153  float barrel_betal = 10.0;
154  float barrel_betah = 170.0;
155  //float barrel_phil = -180.0;
156  //float barrel_phih = 180.0;
157  //float barrel_etal = -2.5;
158  //float barrel_etah = 2.5;
159 
160  float forward_p1_alphal = 100.0;
161  float forward_p1_alphah = 115.0;
162  float forward_p2_alphal = 65.0;
163  float forward_p2_alphah = 80.0;
164  float forward_neg_betal = 67.0;
165  float forward_neg_betah = 73.0;
166  float forward_pos_betal = 107.0;
167  float forward_pos_betah = 113.0;
168  //float forward_phil = -180.0;
169  //float forward_phih = 180.0;
170  //float forward_neg_etal = -2.5;
171  //float forward_neg_etah = -1.5;
172  //float forward_pos_etal = 1.5;
173  //float forward_pos_etah = 2.5;
174 
175  // special ranges for pulls
176  float pull_barrel_alphal = 80.0;
177  float pull_barrel_alphah = 100.0;
178  float pull_barrel_betal = 10.0;
179  float pull_barrel_betah = 170.0;
180  float pull_barrel_phil = -180.0;
181  float pull_barrel_phih = 180.0;
182  float pull_barrel_etal = -2.4;
183  float pull_barrel_etah = 2.4;
184 
185  float pull_forward_p1_alphal = 100.0;
186  float pull_forward_p1_alphah = 112.0;
187  float pull_forward_p2_alphal = 68.0;
188  float pull_forward_p2_alphah = 80.0;
189  float pull_forward_neg_betal = 68.0;
190  float pull_forward_neg_betah = 72.0;
191  float pull_forward_pos_betal = 108.0;
192  float pull_forward_pos_betah = 112.0;
193  float pull_forward_phil = -180.0;
194  float pull_forward_phih = 180.0;
195  float pull_forward_neg_etal = -2.4;
196  float pull_forward_neg_etah = -1.4;
197  float pull_forward_pos_etal = 1.5;
198  float pull_forward_pos_etah = 2.5;
199 
200  int npixl = 0;
201  int npixh = 20;
202  int nxpixl = 0;
203  int nxpixh = 10;
204  int nypixl = 0;
205  int nypixh = 20;
206 
207  float barrel_chargel = 0.0;
208  float barrel_chargeh = 250000.0;
209  float forward_chargel = 0.0;
210  float forward_chargeh = 100000.0;
211 
212  dbe_->setCurrentFolder("Tracking/TrackingRecHits/Pixel/Histograms_per_ring-layer_or_disk-plaquette");
213 
214  // Pixel barrel has 3 layers and 8 rings; book a histogram for each module given by the (layer, ring) pair
215  for (int i=0; i<3 ; i++) // loop over layers
216  {
217  Char_t chisto[100];
218 
219  sprintf(chisto, "meResxBarrelLayer_%d", i+1);
220  meResxBarrelLayer[i] = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
221  sprintf(chisto, "meResyBarrelLayer_%d", i+1);
222  meResyBarrelLayer[i] = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
223  sprintf(chisto, "mePullxBarrelLayer_%d", i+1);
224  mePullxBarrelLayer[i] = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
225  sprintf(chisto, "mePullyBarrelLayer_%d", i+1);
226  mePullyBarrelLayer[i] = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
227 
228  sprintf(chisto, "meResXvsAlphaBarrelFlippedLaddersLayer_%d", i+1);
229  meResXvsAlphaBarrelFlippedLaddersLayer[i] = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh, "");
230  sprintf(chisto, "meResYvsAlphaBarrelFlippedLaddersLayer_%d", i+1);
231  meResYvsAlphaBarrelFlippedLaddersLayer[i] = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh, "");
232  sprintf(chisto, "meResXvsBetaBarrelFlippedLaddersLayer_%d", i+1);
233  meResXvsBetaBarrelFlippedLaddersLayer[i] = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh, "");
234  sprintf(chisto, "meResYvsBetaBarrelFlippedLaddersLayer_%d", i+1);
235  meResYvsBetaBarrelFlippedLaddersLayer[i] = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh, "");
236 
237  sprintf(chisto, "meResXvsAlphaBarrelNonFlippedLaddersLayer_%d", i+1);
238  meResXvsAlphaBarrelNonFlippedLaddersLayer[i]
239  = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh, "");
240  sprintf(chisto, "meResYvsAlphaBarrelNonFlippedLaddersLayer_%d", i+1);
241  meResYvsAlphaBarrelNonFlippedLaddersLayer[i]
242  = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh, "");
243  sprintf(chisto, "meResXvsBetaBarrelNonFlippedLaddersLayer_%d", i+1);
244  meResXvsBetaBarrelNonFlippedLaddersLayer[i]
245  = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh, "");
246  sprintf(chisto, "meResYvsBetaBarrelNonFlippedLaddersLayer_%d", i+1);
247  meResYvsBetaBarrelNonFlippedLaddersLayer[i]
248  = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh, "");
249 
250  for (int j=0; j<8; j++) // loop over rings
251  {
252  sprintf(chisto, "mePosxBarrelLayerModule_%d_%d", i+1, j+1);
253  mePosxBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, xl, xh);
254  sprintf(chisto, "mePosyBarrelLayerModule_%d_%d", i+1, j+1);
255  mePosyBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, yl, yh);
256  sprintf(chisto, "meErrxBarrelLayerModule_%d_%d", i+1, j+1);
257  meErrxBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
258  sprintf(chisto, "meErryBarrelLayerModule_%d_%d", i+1, j+1);
259  meErryBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
260  sprintf(chisto, "meResxBarrelLayerModule_%d_%d", i+1, j+1);
261  meResxBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
262  sprintf(chisto, "meResyBarrelLayerModule_%d_%d", i+1, j+1);
263  meResyBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
264  sprintf(chisto, "mePullxBarrelLayerModule_%d_%d", i+1, j+1);
265  mePullxBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
266  sprintf(chisto, "mePullyBarrelLayerModule_%d_%d", i+1, j+1);
267  mePullyBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
268  sprintf(chisto, "meNpixBarrelLayerModule_%d_%d", i+1, j+1);
269  meNpixBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
270  sprintf(chisto, "meNxpixBarrelLayerModule_%d_%d", i+1, j+1);
271  meNxpixBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
272  sprintf(chisto, "meNypixBarrelLayerModule_%d_%d", i+1, j+1);
273  meNypixBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
274  sprintf(chisto, "meChargeBarrelLayerModule_%d_%d", i+1, j+1);
275  meChargeBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, barrel_chargel, barrel_chargeh);
276 
277  sprintf(chisto, "meResXvsAlphaBarrelLayerModule_%d_%d", i+1, j+1);
278  meResXvsAlphaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh, "");
279  sprintf(chisto, "meResYvsAlphaBarrelLayerModule_%d_%d", i+1, j+1);
280  meResYvsAlphaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh, "");
281  sprintf(chisto, "meResXvsBetaBarrelLayerModule_%d_%d", i+1, j+1);
282  meResXvsBetaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh, "");
283  sprintf(chisto, "meResYvsBetaBarrelLayerModule_%d_%d", i+1, j+1);
284  meResYvsBetaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh, "");
285 
286  sprintf(chisto, "mePullXvsAlphaBarrelLayerModule_%d_%d", i+1, j+1);
287  mePullXvsAlphaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah,
288  100, pullxl, pullxh, "");
289  sprintf(chisto, "mePullYvsAlphaBarrelLayerModule_%d_%d", i+1, j+1);
290  mePullYvsAlphaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah,
291  100, pullyl, pullyh, "");
292  sprintf(chisto, "mePullXvsBetaBarrelLayerModule_%d_%d", i+1, j+1);
293  mePullXvsBetaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah,
294  100, pullxl, pullxh, "");
295  sprintf(chisto, "mePullYvsBetaBarrelLayerModule_%d_%d", i+1, j+1);
296  mePullYvsBetaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah,
297  100, pullyl, pullyh, "");
298  sprintf(chisto, "mePullXvsPhiBarrelLayerModule_%d_%d", i+1, j+1);
299  mePullXvsPhiBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100,
300  pullxl, pullxh, "");
301  sprintf(chisto, "mePullYvsPhiBarrelLayerModule_%d_%d", i+1, j+1);
302  mePullYvsPhiBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih,
303  100, pullyl, pullyh, "");
304  sprintf(chisto, "mePullXvsEtaBarrelLayerModule_%d_%d", i+1, j+1);
305  mePullXvsEtaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah,
306  100, pullxl, pullxh, "");
307  sprintf(chisto, "mePullYvsEtaBarrelLayerModule_%d_%d", i+1, j+1);
308  mePullYvsEtaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah,
309  100, pullyl, pullyh, "");
310  } // for (int j=0; j<8; j++) // loop over rings
311 
312  } // for (int i=0; i<3 ; i++) // loop over layers
313 
314  // Pixel forward detector has 2 disks, 2 panels and either 3 or 4 plaquettes
315  // Panel 1 has 4 plaquettes
316  // Panel 2 has 3 plaquettes
317 
318  // Panel 1: 2 disks, 4 plaquets
319  // Panel 2: 2 disks, 3 plaquets
320  for (int i=0; i<2 ; i++) // loop over disks
321  for (int j=0; j<4; j++) // loop over plaquettes
322  {
323  Char_t chisto[100];
324 
325  sprintf(chisto, "mePosxZmPanel1DiskPlaq_%d_%d", i+1, j+1);
326  mePosxZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, xl, xh);
327  sprintf(chisto, "mePosyZmPanel1DiskPlaq_%d_%d", i+1, j+1);
328  mePosyZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, yl, yh);
329  sprintf(chisto, "meErrxZmPanel1DiskPlaq_%d_%d", i+1, j+1);
330  meErrxZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
331  sprintf(chisto, "meErryZmPanel1DiskPlaq_%d_%d", i+1, j+1);
332  meErryZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
333  sprintf(chisto, "meResxZmPanel1DiskPlaq_%d_%d", i+1, j+1);
334  meResxZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
335  sprintf(chisto, "meResyZmPanel1DiskPlaq_%d_%d", i+1, j+1);
336  meResyZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
337  sprintf(chisto, "mePullxZmPanel1DiskPlaq_%d_%d", i+1, j+1);
338  mePullxZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
339  sprintf(chisto, "mePullyZmPanel1DiskPlaq_%d_%d", i+1, j+1);
340  mePullyZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
341  sprintf(chisto, "meNpixZmPanel1DiskPlaq_%d_%d", i+1, j+1);
342  meNpixZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
343  sprintf(chisto, "meNxpixZmPanel1DiskPlaq_%d_%d", i+1, j+1);
344  meNxpixZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
345  sprintf(chisto, "meNypixZmPanel1DiskPlaq_%d_%d", i+1, j+1);
346  meNypixZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
347  sprintf(chisto, "meChargeZmPanel1DiskPlaq_%d_%d", i+1, j+1);
348  meChargeZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
349 
350  sprintf(chisto, "meResXvsAlphaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
351  meResXvsAlphaZmPanel1DiskPlaq[i][j]
352  = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh, "");
353  sprintf(chisto, "meResYvsAlphaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
354  meResYvsAlphaZmPanel1DiskPlaq[i][j]
355  = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh, "");
356  sprintf(chisto, "meResXvsBetaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
357  meResXvsBetaZmPanel1DiskPlaq[i][j]
358  = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh, "");
359  sprintf(chisto, "meResYvsBetaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
360  meResYvsBetaZmPanel1DiskPlaq[i][j]
361  = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh, "");
362 
363  sprintf(chisto, "mePullXvsAlphaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
364  mePullXvsAlphaZmPanel1DiskPlaq[i][j]
365  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl,pullxh, "");
366  sprintf(chisto, "mePullYvsAlphaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
367  mePullYvsAlphaZmPanel1DiskPlaq[i][j]
368  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl,pullyh, "");
369  sprintf(chisto, "mePullXvsBetaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
370  mePullXvsBetaZmPanel1DiskPlaq[i][j]
371  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl,pullxh, "");
372  sprintf(chisto, "mePullYvsBetaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
373  mePullYvsBetaZmPanel1DiskPlaq[i][j]
374  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl,pullyh, "");
375  sprintf(chisto, "mePullXvsPhiZmPanel1DiskPlaq_%d_%d", i+1, j+1);
376  mePullXvsPhiZmPanel1DiskPlaq[i][j]
377  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl,pullxh, "");
378  sprintf(chisto, "mePullYvsPhiZmPanel1DiskPlaq_%d_%d", i+1, j+1);
379  mePullYvsPhiZmPanel1DiskPlaq[i][j]
380  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl,pullyh, "");
381  sprintf(chisto, "mePullXvsEtaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
382  mePullXvsEtaZmPanel1DiskPlaq[i][j]
383  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl,pullxh, "");
384  sprintf(chisto, "mePullYvsEtaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
385  mePullYvsEtaZmPanel1DiskPlaq[i][j]
386  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl,pullyh, "");
387 
388  sprintf(chisto, "mePosxZpPanel1DiskPlaq_%d_%d", i+1, j+1);
389  mePosxZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, xl, xh);
390  sprintf(chisto, "mePosyZpPanel1DiskPlaq_%d_%d", i+1, j+1);
391  mePosyZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, yl, yh);
392  sprintf(chisto, "meErrxZpPanel1DiskPlaq_%d_%d", i+1, j+1);
393  meErrxZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
394  sprintf(chisto, "meErryZpPanel1DiskPlaq_%d_%d", i+1, j+1);
395  meErryZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
396  sprintf(chisto, "meResxZpPanel1DiskPlaq_%d_%d", i+1, j+1);
397  meResxZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
398  sprintf(chisto, "meResyZpPanel1DiskPlaq_%d_%d", i+1, j+1);
399  meResyZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
400  sprintf(chisto, "mePullxZpPanel1DiskPlaq_%d_%d", i+1, j+1);
401  mePullxZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
402  sprintf(chisto, "mePullyZpPanel1DiskPlaq_%d_%d", i+1, j+1);
403  mePullyZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
404  sprintf(chisto, "meNpixZpPanel1DiskPlaq_%d_%d", i+1, j+1);
405  meNpixZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
406  sprintf(chisto, "meNxpixZpPanel1DiskPlaq_%d_%d", i+1, j+1);
407  meNxpixZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
408  sprintf(chisto, "meNypixZpPanel1DiskPlaq_%d_%d", i+1, j+1);
409  meNypixZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
410  sprintf(chisto, "meChargeZpPanel1DiskPlaq_%d_%d", i+1, j+1);
411  meChargeZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
412  sprintf(chisto, "meResXvsAlphaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
413  meResXvsAlphaZpPanel1DiskPlaq[i][j]
414  = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh, "");
415  sprintf(chisto, "meResYvsAlphaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
416  meResYvsAlphaZpPanel1DiskPlaq[i][j]
417  = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh, "");
418  sprintf(chisto, "meResXvsBetaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
419  meResXvsBetaZpPanel1DiskPlaq[i][j]
420  = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh, "");
421  sprintf(chisto, "meResYvsBetaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
422  meResYvsBetaZpPanel1DiskPlaq[i][j]
423  = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh, "");
424 
425  sprintf(chisto, "mePullXvsAlphaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
426  mePullXvsAlphaZpPanel1DiskPlaq[i][j]
427  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl,pullxh, "");
428  sprintf(chisto, "mePullYvsAlphaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
429  mePullYvsAlphaZpPanel1DiskPlaq[i][j]
430  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl,pullyh, "");
431  sprintf(chisto, "mePullXvsBetaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
432  mePullXvsBetaZpPanel1DiskPlaq[i][j]
433  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl,pullxh, "");
434  sprintf(chisto, "mePullYvsBetaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
435  mePullYvsBetaZpPanel1DiskPlaq[i][j]
436  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl,pullyh, "");
437  sprintf(chisto, "mePullXvsPhiZpPanel1DiskPlaq_%d_%d", i+1, j+1);
438  mePullXvsPhiZpPanel1DiskPlaq[i][j]
439  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl,pullxh, "");
440  sprintf(chisto, "mePullYvsPhiZpPanel1DiskPlaq_%d_%d", i+1, j+1);
441  mePullYvsPhiZpPanel1DiskPlaq[i][j]
442  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl,pullyh, "");
443  sprintf(chisto, "mePullXvsEtaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
444  mePullXvsEtaZpPanel1DiskPlaq[i][j]
445  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl,pullxh, "");
446  sprintf(chisto, "mePullYvsEtaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
447  mePullYvsEtaZpPanel1DiskPlaq[i][j]
448  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl,pullyh, "");
449 
450  if ( j>2 ) continue; // panel 2 has only 3 plaquettes
451 
452  sprintf(chisto, "mePosxZmPanel2DiskPlaq_%d_%d", i+1, j+1);
453  mePosxZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, xl, xh);
454  sprintf(chisto, "mePosyZmPanel2DiskPlaq_%d_%d", i+1, j+1);
455  mePosyZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, yl, yh);
456  sprintf(chisto, "meErrxZmPanel2DiskPlaq_%d_%d", i+1, j+1);
457  meErrxZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
458  sprintf(chisto, "meErryZmPanel2DiskPlaq_%d_%d", i+1, j+1);
459  meErryZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
460  sprintf(chisto, "meResxZmPanel2DiskPlaq_%d_%d", i+1, j+1);
461  meResxZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
462  sprintf(chisto, "meResyZmPanel2DiskPlaq_%d_%d", i+1, j+1);
463  meResyZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
464  sprintf(chisto, "mePullxZmPanel2DiskPlaq_%d_%d", i+1, j+1);
465  mePullxZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
466  sprintf(chisto, "mePullyZmPanel2DiskPlaq_%d_%d", i+1, j+1);
467  mePullyZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
468  sprintf(chisto, "meNpixZmPanel2DiskPlaq_%d_%d", i+1, j+1);
469  meNpixZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
470  sprintf(chisto, "meNxpixZmPanel2DiskPlaq_%d_%d", i+1, j+1);
471  meNxpixZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
472  sprintf(chisto, "meNypixZmPanel2DiskPlaq_%d_%d", i+1, j+1);
473  meNypixZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
474  sprintf(chisto, "meChargeZmPanel2DiskPlaq_%d_%d", i+1, j+1);
475  meChargeZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
476  sprintf(chisto, "meResXvsAlphaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
477  meResXvsAlphaZmPanel2DiskPlaq[i][j]
478  = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh, "");
479  sprintf(chisto, "meResYvsAlphaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
480  meResYvsAlphaZmPanel2DiskPlaq[i][j]
481  = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh, "");
482  sprintf(chisto, "meResXvsBetaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
483  meResXvsBetaZmPanel2DiskPlaq[i][j]
484  = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh, "");
485  sprintf(chisto, "meResYvsBetaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
486  meResYvsBetaZmPanel2DiskPlaq[i][j]
487  = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh, "");
488 
489  sprintf(chisto, "mePullXvsAlphaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
490  mePullXvsAlphaZmPanel2DiskPlaq[i][j]
491  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl,pullxh, "");
492  sprintf(chisto, "mePullYvsAlphaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
493  mePullYvsAlphaZmPanel2DiskPlaq[i][j]
494  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl,pullyh, "");
495  sprintf(chisto, "mePullXvsBetaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
496  mePullXvsBetaZmPanel2DiskPlaq[i][j]
497  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl,pullxh, "");
498  sprintf(chisto, "mePullYvsBetaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
499  mePullYvsBetaZmPanel2DiskPlaq[i][j]
500  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl,pullyh, "");
501  sprintf(chisto, "mePullXvsPhiZmPanel2DiskPlaq_%d_%d", i+1, j+1);
502  mePullXvsPhiZmPanel2DiskPlaq[i][j]
503  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl,pullxh, "");
504  sprintf(chisto, "mePullYvsPhiZmPanel2DiskPlaq_%d_%d", i+1, j+1);
505  mePullYvsPhiZmPanel2DiskPlaq[i][j]
506  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl,pullyh, "");
507  sprintf(chisto, "mePullXvsEtaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
508  mePullXvsEtaZmPanel2DiskPlaq[i][j]
509  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh, "");
510  sprintf(chisto, "mePullYvsEtaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
511  mePullYvsEtaZmPanel2DiskPlaq[i][j]
512  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh, "");
513 
514  sprintf(chisto, "mePosxZpPanel2DiskPlaq_%d_%d", i+1, j+1);
515  mePosxZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, xl, xh);
516  sprintf(chisto, "mePosyZpPanel2DiskPlaq_%d_%d", i+1, j+1);
517  mePosyZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, yl, yh);
518  sprintf(chisto, "meErrxZpPanel2DiskPlaq_%d_%d", i+1, j+1);
519  meErrxZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
520  sprintf(chisto, "meErryZpPanel2DiskPlaq_%d_%d", i+1, j+1);
521  meErryZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
522  sprintf(chisto, "meResxZpPanel2DiskPlaq_%d_%d", i+1, j+1);
523  meResxZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
524  sprintf(chisto, "meResyZpPanel2DiskPlaq_%d_%d", i+1, j+1);
525  meResyZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
526  sprintf(chisto, "mePullxZpPanel2DiskPlaq_%d_%d", i+1, j+1);
527  mePullxZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
528  sprintf(chisto, "mePullyZpPanel2DiskPlaq_%d_%d", i+1, j+1);
529  mePullyZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
530  sprintf(chisto, "meNpixZpPanel2DiskPlaq_%d_%d", i+1, j+1);
531  meNpixZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
532  sprintf(chisto, "meNxpixZpPanel2DiskPlaq_%d_%d", i+1, j+1);
533  meNxpixZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
534  sprintf(chisto, "meNypixZpPanel2DiskPlaq_%d_%d", i+1, j+1);
535  meNypixZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
536  sprintf(chisto, "meChargeZpPanel2DiskPlaq_%d_%d", i+1, j+1);
537  meChargeZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
538  sprintf(chisto, "meResXvsAlphaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
539  meResXvsAlphaZpPanel2DiskPlaq[i][j]
540  = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh, "");
541  sprintf(chisto, "meResYvsAlphaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
542  meResYvsAlphaZpPanel2DiskPlaq[i][j]
543  = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh, "");
544  sprintf(chisto, "meResXvsBetaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
545  meResXvsBetaZpPanel2DiskPlaq[i][j]
546  = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh, "");
547  sprintf(chisto, "meResYvsBetaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
548  meResYvsBetaZpPanel2DiskPlaq[i][j]
549  = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh, "");
550 
551  sprintf(chisto, "mePullXvsAlphaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
552  mePullXvsAlphaZpPanel2DiskPlaq[i][j]
553  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl,pullxh, "");
554  sprintf(chisto, "mePullYvsAlphaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
555  mePullYvsAlphaZpPanel2DiskPlaq[i][j]
556  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl,pullyh, "");
557  sprintf(chisto, "mePullXvsBetaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
558  mePullXvsBetaZpPanel2DiskPlaq[i][j]
559  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl,pullxh, "");
560  sprintf(chisto, "mePullYvsBetaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
561  mePullYvsBetaZpPanel2DiskPlaq[i][j]
562  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl,pullyh, "");
563  sprintf(chisto, "mePullXvsPhiZpPanel2DiskPlaq_%d_%d", i+1, j+1);
564  mePullXvsPhiZpPanel2DiskPlaq[i][j]
565  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl,pullxh, "");
566  sprintf(chisto, "mePullYvsPhiZpPanel2DiskPlaq_%d_%d", i+1, j+1);
567  mePullYvsPhiZpPanel2DiskPlaq[i][j]
568  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl,pullyh, "");
569  sprintf(chisto, "mePullXvsEtaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
570  mePullXvsEtaZpPanel2DiskPlaq[i][j]
571  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl,pullxh, "");
572  sprintf(chisto, "mePullYvsEtaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
573  mePullYvsEtaZpPanel2DiskPlaq[i][j]
574  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl,pullyh, "");
575 
576  } // for (int j=0; j<4; j++) // loop over plaquettes
577 
578  dbe_->setCurrentFolder("Tracking/TrackingRecHits/Pixel/Histograms_all");
579 
580  Char_t chisto[100];
581  sprintf(chisto, "mePosxBarrel");
582  mePosxBarrel = dbe_->book1D(chisto, chisto, 100, xl, xh);
583  sprintf(chisto, "mePosyBarrel");
584  mePosyBarrel = dbe_->book1D(chisto, chisto, 100, yl, yh);
585  sprintf(chisto, "meErrxBarrel");
586  meErrxBarrel = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
587  sprintf(chisto, "meErryBarrel");
588  meErryBarrel = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
589  sprintf(chisto, "meResxBarrel");
590  meResxBarrel = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
591  sprintf(chisto, "meResyBarrel");
592  meResyBarrel = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
593  sprintf(chisto, "mePullxBarrel");
594  mePullxBarrel = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
595  sprintf(chisto, "mePullyBarrel");
596  mePullyBarrel = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
597  sprintf(chisto, "meNpixBarrel");
598  meNpixBarrel = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
599  sprintf(chisto, "meNxpixBarrel");
600  meNxpixBarrel = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
601  sprintf(chisto, "meNypixBarrel");
602  meNypixBarrel = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
603  sprintf(chisto, "meChargeBarrel");
604  meChargeBarrel = dbe_->book1D(chisto, chisto, 100, barrel_chargel, barrel_chargeh);
605  sprintf(chisto, "meResXvsAlphaBarrel");
606  meResXvsAlphaBarrel = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh, "");
607  sprintf(chisto, "meResYvsAlphaBarrel");
608  meResYvsAlphaBarrel = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh, "");
609  sprintf(chisto, "meResXvsBetaBarrel");
610  meResXvsBetaBarrel = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh, "");
611  sprintf(chisto, "meResYvsBetaBarrel");
612  meResYvsBetaBarrel = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh, "");
613 
614  sprintf(chisto, "mePullXvsAlphaBarrel");
615  mePullXvsAlphaBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh, "");
616  sprintf(chisto, "mePullYvsAlphaBarrel");
617  mePullYvsAlphaBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh, "");
618  sprintf(chisto, "mePullXvsBetaBarrel");
619  mePullXvsBetaBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh, "");
620  sprintf(chisto, "mePullYvsBetaBarrel");
621  mePullYvsBetaBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh, "");
622  sprintf(chisto, "mePullXvsPhiBarrel");
623  mePullXvsPhiBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh, "");
624  sprintf(chisto, "mePullYvsPhiBarrel");
625  mePullYvsPhiBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh, "");
626  sprintf(chisto, "mePullXvsEtaBarrel");
627  mePullXvsEtaBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh, "");
628  sprintf(chisto, "mePullYvsEtaBarrel");
629  mePullYvsEtaBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh, "");
630 
631  sprintf(chisto, "mePosxBarrelHalfModule");
632  mePosxBarrelHalfModule = dbe_->book1D(chisto, chisto, 100, xl, xh);
633  sprintf(chisto, "mePosxBarrelFullModule");
634  mePosxBarrelFullModule = dbe_->book1D(chisto, chisto, 100, xl, xh);
635  sprintf(chisto, "mePosxBarrelFlippedLadders");
636  mePosxBarrelFlippedLadders = dbe_->book1D(chisto, chisto, 100, xl, xh);
637  sprintf(chisto, "mePosxBarrelNonFlippedLadders");
638  mePosxBarrelNonFlippedLadders = dbe_->book1D(chisto, chisto, 100, xl, xh);
639  sprintf(chisto, "mePosyBarrelHalfModule");
640  mePosyBarrelHalfModule = dbe_->book1D(chisto, chisto, 100, yl, yh);
641  sprintf(chisto, "mePosyBarrelFullModule");
642  mePosyBarrelFullModule = dbe_->book1D(chisto, chisto, 100, yl, yh);
643  sprintf(chisto, "mePosyBarrelFlippedLadders");
644  mePosyBarrelFlippedLadders = dbe_->book1D(chisto, chisto, 100, yl, yh);
645  sprintf(chisto, "mePosyBarrelNonFlippedLadders");
646  mePosyBarrelNonFlippedLadders = dbe_->book1D(chisto, chisto, 100, yl, yh);
647 
648  sprintf(chisto, "meResXvsAlphaBarrelFlippedLadders");
649  meResXvsAlphaBarrelFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh, "");
650  sprintf(chisto, "meResYvsAlphaBarrelFlippedLadders");
651  meResYvsAlphaBarrelFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh, "");
652  sprintf(chisto, "meResXvsBetaBarrelFlippedLadders");
653  meResXvsBetaBarrelFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh, "");
654  sprintf(chisto, "meResYvsBetaBarrelFlippedLadders");
655  meResYvsBetaBarrelFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh, "");
656 
657  sprintf(chisto, "mePullXvsAlphaBarrelFlippedLadders");
658  mePullXvsAlphaBarrelFlippedLadders
659  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh, "");
660  sprintf(chisto, "mePullYvsAlphaBarrelFlippedLadders");
661  mePullYvsAlphaBarrelFlippedLadders
662  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh, "");
663  sprintf(chisto, "mePullXvsBetaBarrelFlippedLadders");
664  mePullXvsBetaBarrelFlippedLadders
665  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh, "");
666  sprintf(chisto, "mePullYvsBetaBarrelFlippedLadders");
667  mePullYvsBetaBarrelFlippedLadders
668  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh, "");
669  sprintf(chisto, "mePullXvsPhiBarrelFlippedLadders");
670  mePullXvsPhiBarrelFlippedLadders
671  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh, "");
672  sprintf(chisto, "mePullYvsPhiBarrelFlippedLadders");
673  mePullYvsPhiBarrelFlippedLadders
674  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh, "");
675  sprintf(chisto, "mePullXvsEtaBarrelFlippedLadders");
676  mePullXvsEtaBarrelFlippedLadders
677  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh, "");
678  sprintf(chisto, "mePullYvsEtaBarrelFlippedLadders");
679  mePullYvsEtaBarrelFlippedLadders
680  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh, "");
681 
682 
683  sprintf(chisto, "meWPullXvsAlphaBarrelFlippedLadders");
684  meWPullXvsAlphaBarrelFlippedLadders
685  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh, "");
686  sprintf(chisto, "meWPullYvsAlphaBarrelFlippedLadders");
687  meWPullYvsAlphaBarrelFlippedLadders
688  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh, "");
689  sprintf(chisto, "meWPullXvsBetaBarrelFlippedLadders");
690  meWPullXvsBetaBarrelFlippedLadders
691  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh, "");
692  sprintf(chisto, "meWPullYvsBetaBarrelFlippedLadders");
693  meWPullYvsBetaBarrelFlippedLadders
694  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh, "");
695 
696  sprintf(chisto, "meResXvsAlphaBarrelNonFlippedLadders");
697  meResXvsAlphaBarrelNonFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh, "");
698  sprintf(chisto, "meResYvsAlphaBarrelNonFlippedLadders");
699  meResYvsAlphaBarrelNonFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh, "");
700  sprintf(chisto, "meResXvsBetaBarrelNonFlippedLadders");
701  meResXvsBetaBarrelNonFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh, "");
702  sprintf(chisto, "meResYvsBetaBarrelNonFlippedLadders");
703  meResYvsBetaBarrelNonFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh, "");
704 
705  sprintf(chisto, "mePullXvsAlphaBarrelNonFlippedLadders");
706  mePullXvsAlphaBarrelNonFlippedLadders
707  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh, "");
708  sprintf(chisto, "mePullYvsAlphaBarrelNonFlippedLadders");
709  mePullYvsAlphaBarrelNonFlippedLadders
710  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh, "");
711  sprintf(chisto, "mePullXvsBetaBarrelNonFlippedLadders");
712  mePullXvsBetaBarrelNonFlippedLadders
713  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh, "");
714  sprintf(chisto, "mePullYvsBetaBarrelNonFlippedLadders");
715  mePullYvsBetaBarrelNonFlippedLadders
716  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh, "");
717  sprintf(chisto, "mePullXvsPhiBarrelNonFlippedLadders");
718  mePullXvsPhiBarrelNonFlippedLadders
719  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh, "");
720  sprintf(chisto, "mePullYvsPhiBarrelNonFlippedLadders");
721  mePullYvsPhiBarrelNonFlippedLadders
722  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh, "");
723  sprintf(chisto, "mePullXvsEtaBarrelNonFlippedLadders");
724  mePullXvsEtaBarrelNonFlippedLadders
725  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh, "");
726  sprintf(chisto, "mePullYvsEtaBarrelNonFlippedLadders");
727  mePullYvsEtaBarrelNonFlippedLadders
728  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh, "");
729 
730 
731  sprintf(chisto, "meWPullXvsAlphaBarrelNonFlippedLadders");
732  meWPullXvsAlphaBarrelNonFlippedLadders
733  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh, "");
734  sprintf(chisto, "meWPullYvsAlphaBarrelNonFlippedLadders");
735  meWPullYvsAlphaBarrelNonFlippedLadders
736  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh, "");
737  sprintf(chisto, "meWPullXvsBetaBarrelNonFlippedLadders");
738  meWPullXvsBetaBarrelNonFlippedLadders
739  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh, "");
740  sprintf(chisto, "meWPullYvsBetaBarrelNonFlippedLadders");
741  meWPullYvsBetaBarrelNonFlippedLadders
742  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh, "");
743 
744 
745  sprintf(chisto, "mePosxZmPanel1");
746  mePosxZmPanel1 = dbe_->book1D(chisto, chisto, 100, xl, xh);
747  sprintf(chisto, "mePosyZmPanel1");
748  mePosyZmPanel1 = dbe_->book1D(chisto, chisto, 100, yl, yh);
749  sprintf(chisto, "meErrxZmPanel1");
750  meErrxZmPanel1 = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
751  sprintf(chisto, "meErryZmPanel1");
752  meErryZmPanel1 = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
753  sprintf(chisto, "meResxZmPanel1");
754  meResxZmPanel1 = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
755  sprintf(chisto, "meResyZmPanel1");
756  meResyZmPanel1 = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
757  sprintf(chisto, "mePullxZmPanel1");
758  mePullxZmPanel1 = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
759  sprintf(chisto, "mePullyZmPanel1");
760  mePullyZmPanel1 = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
761  sprintf(chisto, "meNpixZmPanel1");
762  meNpixZmPanel1 = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
763  sprintf(chisto, "meNxpixZmPanel1");
764  meNxpixZmPanel1 = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
765  sprintf(chisto, "meNypixZmPanel1");
766  meNypixZmPanel1 = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
767  sprintf(chisto, "meChargeZmPanel1");
768  meChargeZmPanel1 = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
769  sprintf(chisto, "meResXvsAlphaZmPanel1");
770  meResXvsAlphaZmPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh, "");
771  sprintf(chisto, "meResYvsAlphaZmPanel1");
772  meResYvsAlphaZmPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh, "");
773  sprintf(chisto, "meResXvsBetaZmPanel1");
774  meResXvsBetaZmPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh, "");
775  sprintf(chisto, "meResYvsBetaZmPanel1");
776  meResYvsBetaZmPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh, "");
777 
778  sprintf(chisto, "mePullXvsAlphaZmPanel1");
779  mePullXvsAlphaZmPanel1
780  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh, "");
781  sprintf(chisto, "mePullYvsAlphaZmPanel1");
782  mePullYvsAlphaZmPanel1
783  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh, "");
784  sprintf(chisto, "mePullXvsBetaZmPanel1");
785  mePullXvsBetaZmPanel1
786  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh, "");
787  sprintf(chisto, "mePullYvsBetaZmPanel1");
788  mePullYvsBetaZmPanel1
789  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh, "");
790  sprintf(chisto, "mePullXvsPhiZmPanel1");
791  mePullXvsPhiZmPanel1
792  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh, "");
793  sprintf(chisto, "mePullYvsPhiZmPanel1");
794  mePullYvsPhiZmPanel1
795  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh, "");
796  sprintf(chisto, "mePullXvsEtaZmPanel1");
797  mePullXvsEtaZmPanel1
798  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh, "");
799  sprintf(chisto, "mePullYvsEtaZmPanel1");
800  mePullYvsEtaZmPanel1
801  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh, "");
802 
803  sprintf(chisto, "meWPullXvsAlphaZmPanel1");
804  meWPullXvsAlphaZmPanel1
805  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh, "");
806  sprintf(chisto, "meWPullYvsAlphaZmPanel1");
807  meWPullYvsAlphaZmPanel1
808  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh, "");
809  sprintf(chisto, "meWPullXvsBetaZmPanel1");
810  meWPullXvsBetaZmPanel1
811  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh, "");
812  sprintf(chisto, "meWPullYvsBetaZmPanel1");
813  meWPullYvsBetaZmPanel1
814  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh, "");
815 
816  sprintf(chisto, "mePosxZpPanel1");
817  mePosxZpPanel1 = dbe_->book1D(chisto, chisto, 100, xl, xh);
818  sprintf(chisto, "mePosyZpPanel1");
819  mePosyZpPanel1 = dbe_->book1D(chisto, chisto, 100, yl, yh);
820  sprintf(chisto, "meErrxZpPanel1");
821  meErrxZpPanel1 = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
822  sprintf(chisto, "meErryZpPanel1");
823  meErryZpPanel1 = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
824  sprintf(chisto, "meResxZpPanel1");
825  meResxZpPanel1 = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
826  sprintf(chisto, "meResyZpPanel1");
827  meResyZpPanel1 = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
828  sprintf(chisto, "mePullxZpPanel1");
829  mePullxZpPanel1 = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
830  sprintf(chisto, "mePullyZpPanel1");
831  mePullyZpPanel1 = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
832  sprintf(chisto, "meNpixZpPanel1");
833  meNpixZpPanel1 = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
834  sprintf(chisto, "meNxpixZpPanel1");
835  meNxpixZpPanel1 = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
836  sprintf(chisto, "meNypixZpPanel1");
837  meNypixZpPanel1 = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
838  sprintf(chisto, "meChargeZpPanel1");
839  meChargeZpPanel1 = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
840  sprintf(chisto, "meResXvsAlphaZpPanel1");
841  meResXvsAlphaZpPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh, "");
842  sprintf(chisto, "meResYvsAlphaZpPanel1");
843  meResYvsAlphaZpPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh, "");
844  sprintf(chisto, "meResXvsBetaZpPanel1");
845  meResXvsBetaZpPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh, "");
846  sprintf(chisto, "meResYvsBetaZpPanel1");
847  meResYvsBetaZpPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh, "");
848 
849  sprintf(chisto, "mePullXvsAlphaZpPanel1");
850  mePullXvsAlphaZpPanel1
851  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh, "");
852  sprintf(chisto, "mePullYvsAlphaZpPanel1");
853  mePullYvsAlphaZpPanel1
854  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh, "");
855  sprintf(chisto, "mePullXvsBetaZpPanel1");
856  mePullXvsBetaZpPanel1
857  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh, "");
858  sprintf(chisto, "mePullYvsBetaZpPanel1");
859  mePullYvsBetaZpPanel1
860  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh, "");
861  sprintf(chisto, "mePullXvsPhiZpPanel1");
862  mePullXvsPhiZpPanel1
863  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh, "");
864  sprintf(chisto, "mePullYvsPhiZpPanel1");
865  mePullYvsPhiZpPanel1
866  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh, "");
867  sprintf(chisto, "mePullXvsEtaZpPanel1");
868  mePullXvsEtaZpPanel1
869  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl, pullxh, "");
870  sprintf(chisto, "mePullYvsEtaZpPanel1");
871  mePullYvsEtaZpPanel1
872  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl, pullyh, "");
873 
874  sprintf(chisto, "meWPullXvsAlphaZpPanel1");
875  meWPullXvsAlphaZpPanel1
876  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh, "");
877  sprintf(chisto, "meWPullYvsAlphaZpPanel1");
878  meWPullYvsAlphaZpPanel1
879  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh, "");
880  sprintf(chisto, "meWPullXvsBetaZpPanel1");
881  meWPullXvsBetaZpPanel1
882  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh, "");
883  sprintf(chisto, "meWPullYvsBetaZpPanel1");
884  meWPullYvsBetaZpPanel1
885  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh, "");
886 
887  sprintf(chisto, "mePosxZmPanel2");
888  mePosxZmPanel2 = dbe_->book1D(chisto, chisto, 100, xl, xh);
889  sprintf(chisto, "mePosyZmPanel2");
890  mePosyZmPanel2 = dbe_->book1D(chisto, chisto, 100, yl, yh);
891  sprintf(chisto, "meErrxZmPanel2");
892  meErrxZmPanel2 = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
893  sprintf(chisto, "meErryZmPanel2");
894  meErryZmPanel2 = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
895  sprintf(chisto, "meResxZmPanel2");
896  meResxZmPanel2 = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
897  sprintf(chisto, "meResyZmPanel2");
898  meResyZmPanel2 = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
899  sprintf(chisto, "mePullxZmPanel2");
900  mePullxZmPanel2 = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
901  sprintf(chisto, "mePullyZmPanel2");
902  mePullyZmPanel2 = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
903  sprintf(chisto, "meNpixZmPanel2");
904  meNpixZmPanel2 = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
905  sprintf(chisto, "meNxpixZmPanel2");
906  meNxpixZmPanel2 = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
907  sprintf(chisto, "meNypixZmPanel2");
908  meNypixZmPanel2 = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
909  sprintf(chisto, "meChargeZmPanel2");
910  meChargeZmPanel2 = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
911  sprintf(chisto, "meResXvsAlphaZmPanel2");
912  meResXvsAlphaZmPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh, "");
913  sprintf(chisto, "meResYvsAlphaZmPanel2");
914  meResYvsAlphaZmPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh, "");
915  sprintf(chisto, "meResXvsBetaZmPanel2");
916  meResXvsBetaZmPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh, "");
917  sprintf(chisto, "meResYvsBetaZmPanel2");
918  meResYvsBetaZmPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh, "");
919 
920  sprintf(chisto, "mePullXvsAlphaZmPanel2");
921  mePullXvsAlphaZmPanel2
922  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh, "");
923  sprintf(chisto, "mePullYvsAlphaZmPanel2");
924  mePullYvsAlphaZmPanel2
925  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh, "");
926  sprintf(chisto, "mePullXvsBetaZmPanel2");
927  mePullXvsBetaZmPanel2
928  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh, "");
929  sprintf(chisto, "mePullYvsBetaZmPanel2");
930  mePullYvsBetaZmPanel2
931  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh, "");
932  sprintf(chisto, "mePullXvsPhiZmPanel2");
933  mePullXvsPhiZmPanel2
934  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh, "");
935  sprintf(chisto, "mePullYvsPhiZmPanel2");
936  mePullYvsPhiZmPanel2
937  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh, "");
938  sprintf(chisto, "mePullXvsEtaZmPanel2");
939  mePullXvsEtaZmPanel2
940  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh, "");
941  sprintf(chisto, "mePullYvsEtaZmPanel2");
942  mePullYvsEtaZmPanel2
943  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh, "");
944 
945  sprintf(chisto, "meWPullXvsAlphaZmPanel2");
946  meWPullXvsAlphaZmPanel2
947  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh, "");
948  sprintf(chisto, "meWPullYvsAlphaZmPanel2");
949  meWPullYvsAlphaZmPanel2
950  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh, "");
951  sprintf(chisto, "meWPullXvsBetaZmPanel2");
952  meWPullXvsBetaZmPanel2
953  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh, "");
954  sprintf(chisto, "meWPullYvsBetaZmPanel2");
955  meWPullYvsBetaZmPanel2
956  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh, "");
957 
958 
959  sprintf(chisto, "mePosxZpPanel2");
960  mePosxZpPanel2 = dbe_->book1D(chisto, chisto, 100, xl, xh);
961  sprintf(chisto, "mePosyZpPanel2");
962  mePosyZpPanel2 = dbe_->book1D(chisto, chisto, 100, yl, yh);
963  sprintf(chisto, "meErrxZpPanel2");
964  meErrxZpPanel2 = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
965  sprintf(chisto, "meErryZpPanel2");
966  meErryZpPanel2 = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
967  sprintf(chisto, "meResxZpPanel2");
968  meResxZpPanel2 = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
969  sprintf(chisto, "meResyZpPanel2");
970  meResyZpPanel2 = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
971  sprintf(chisto, "mePullxZpPanel2");
972  mePullxZpPanel2 = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
973  sprintf(chisto, "mePullyZpPanel2");
974  mePullyZpPanel2 = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
975  sprintf(chisto, "meNpixZpPanel2");
976  meNpixZpPanel2 = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
977  sprintf(chisto, "meNxpixZpPanel2");
978  meNxpixZpPanel2 = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
979  sprintf(chisto, "meNypixZpPanel2");
980  meNypixZpPanel2 = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
981  sprintf(chisto, "meChargeZpPanel2");
982  meChargeZpPanel2 = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
983  sprintf(chisto, "meResXvsAlphaZpPanel2");
984  meResXvsAlphaZpPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh, "");
985  sprintf(chisto, "meResYvsAlphaZpPanel2");
986  meResYvsAlphaZpPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh, "");
987  sprintf(chisto, "meResXvsBetaZpPanel2");
988  meResXvsBetaZpPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh, "");
989  sprintf(chisto, "meResYvsBetaZpPanel2");
990  meResYvsBetaZpPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh, "");
991 
992  sprintf(chisto, "mePullXvsAlphaZpPanel2");
993  mePullXvsAlphaZpPanel2
994  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh, "");
995  sprintf(chisto, "mePullYvsAlphaZpPanel2");
996  mePullYvsAlphaZpPanel2
997  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh, "");
998  sprintf(chisto, "mePullXvsBetaZpPanel2");
999  mePullXvsBetaZpPanel2
1000  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh, "");
1001  sprintf(chisto, "mePullYvsBetaZpPanel2");
1002  mePullYvsBetaZpPanel2
1003  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh, "");
1004  sprintf(chisto, "mePullXvsPhiZpPanel2");
1005  mePullXvsPhiZpPanel2
1006  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh, "");
1007  sprintf(chisto, "mePullYvsPhiZpPanel2");
1008  mePullYvsPhiZpPanel2
1009  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh, "");
1010  sprintf(chisto, "mePullXvsEtaZpPanel2");
1011  mePullXvsEtaZpPanel2
1012  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl, pullxh, "");
1013  sprintf(chisto, "mePullYvsEtaZpPanel2");
1014  mePullYvsEtaZpPanel2
1015  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl, pullyh, "");
1016 
1017  sprintf(chisto, "meWPullXvsAlphaZpPanel2");
1018  meWPullXvsAlphaZpPanel2
1019  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh, "");
1020  sprintf(chisto, "meWPullYvsAlphaZpPanel2");
1021  meWPullYvsAlphaZpPanel2
1022  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh, "");
1023  sprintf(chisto, "meWPullXvsBetaZpPanel2");
1024  meWPullXvsBetaZpPanel2
1025  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh, "");
1026  sprintf(chisto, "meWPullYvsBetaZpPanel2");
1027  meWPullYvsBetaZpPanel2
1028  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh, "");
1029 
1030  // all hits (not only from tracks)
1031  sprintf(chisto, "mePosxBarrel_all_hits");
1032  mePosxBarrel_all_hits = dbe_->book1D(chisto, chisto, 100, xl, xh);
1033  sprintf(chisto, "mePosyBarrel_all_hits");
1034  mePosyBarrel_all_hits = dbe_->book1D(chisto, chisto, 100, yl, yh);
1035 
1036  sprintf(chisto, "mePosxZmPanel1_all_hits");
1037  mePosxZmPanel1_all_hits = dbe_->book1D(chisto, chisto, 100, xl, xh);
1038  sprintf(chisto, "mePosyZmPanel1_all_hits");
1039  mePosyZmPanel1_all_hits = dbe_->book1D(chisto, chisto, 100, yl, yh);
1040  sprintf(chisto, "mePosxZmPanel2_all_hits");
1041  mePosxZmPanel2_all_hits = dbe_->book1D(chisto, chisto, 100, xl, xh);
1042  sprintf(chisto, "mePosyZmPanel2_all_hits");
1043  mePosyZmPanel2_all_hits = dbe_->book1D(chisto, chisto, 100, yl, yh);
1044 
1045  sprintf(chisto, "mePosxZpPanel1_all_hits");
1046  mePosxZpPanel1_all_hits = dbe_->book1D(chisto, chisto, 100, xl, xh);
1047  sprintf(chisto, "mePosyZpPanel1_all_hits");
1048  mePosyZpPanel1_all_hits = dbe_->book1D(chisto, chisto, 100, yl, yh);
1049  sprintf(chisto, "mePosxZpPanel2_all_hits");
1050  mePosxZpPanel2_all_hits = dbe_->book1D(chisto, chisto, 100, xl, xh);
1051  sprintf(chisto, "mePosyZpPanel2_all_hits");
1052  mePosyZpPanel2_all_hits = dbe_->book1D(chisto, chisto, 100, yl, yh);
1053 
1054  // control histograms
1055  meTracksPerEvent = dbe_->book1D("meTracksPerEvent" , "meTracksPerEvent" , 200, 0.0, 200.0);
1056  mePixRecHitsPerTrack = dbe_->book1D("mePixRecHitsPerTrack", "mePixRecHitsPerTrack", 6, 0.0, 6.0);
1057 
1058 }
1059 
1060 // Virtual destructor needed.
1061 SiPixelTrackingRecHitsValid::~SiPixelTrackingRecHitsValid()
1062 {
1063  if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
1064 }
1065 
1066 // Functions that gets called by framework every event
1067 void SiPixelTrackingRecHitsValid::analyze(const edm::Event& e, const edm::EventSetup& es)
1068 {
1069  //Retrieve tracker topology from geometry
1070  edm::ESHandle<TrackerTopology> tTopo;
1071  es.get<IdealGeometryRecord>().get(tTopo);
1072 
1073 
1074  run = e.id().run();
1075  evt = e.id().event();
1076 
1077  // if ( evt%1000 == 0 )
1078  //cout << "evt = " << evt << endl;
1079 
1080  float math_pi = 3.14159265;
1081  float radtodeg = 180.0 / math_pi;
1082 
1083  DetId detId;
1084 
1085  LocalPoint position;
1086  LocalError error;
1087  float mindist = 999999.9;
1088 
1089  std::vector<PSimHit> matched;
1090  TrackerHitAssociator associate(e,conf_);
1091 
1092  edm::ESHandle<TrackerGeometry> pDD;
1093  es.get<TrackerDigiGeometryRecord> ().get (pDD);
1094  const TrackerGeometry* tracker = &(* pDD);
1095 
1096  if ( !MTCCtrack_ )
1097  {
1098  // --------------------------------------- all hits -----------------------------------------------------------
1099  //--- Fetch Pixel RecHits
1100  edm::Handle<SiPixelRecHitCollection> recHitColl;
1101  e.getByToken( siPixelRecHitCollectionToken_, recHitColl );
1102 
1103  //cout <<" ----- Found "
1104  // << const_cast<SiPixelRecHitCollection*>(recHitColl.product())->size()
1105  // << " Pixel RecHits" << std::endl;
1106 
1107  //-----Iterate over detunits
1108  for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++)
1109  {
1110  DetId detId = ((*it)->geographicalId());
1111 
1112  unsigned int subid = detId.subdetId();
1113  if ( !((subid==1) || (subid==2)) )
1114  continue; // end subid if
1115 
1116  SiPixelRecHitCollection::const_iterator match = recHitColl->find(detId);
1117  if (match == recHitColl->end()) continue;
1118 
1119  SiPixelRecHitCollection::DetSet pixelrechitRange = *match;
1120  SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorBegin = pixelrechitRange.begin();
1121  SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorEnd = pixelrechitRange.end();
1122  SiPixelRecHitCollection::DetSet::const_iterator pixeliter = pixelrechitRangeIteratorBegin;
1123  std::vector<PSimHit> matched;
1124 
1125  //----Loop over rechits for this detId
1126  for ( ; pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter)
1127  {
1128  LocalPoint lp = pixeliter->localPosition();
1129  float rechitx = lp.x();
1130  float rechity = lp.y();
1131 
1132  detId = (*it)->geographicalId();
1133  subdetId = (int)detId.subdetId();
1134  if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelBarrel )
1135  {
1136  mePosxBarrel_all_hits->Fill( rechitx );
1137  mePosyBarrel_all_hits->Fill( rechity );
1138  }
1139  else if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelEndcap )
1140  {
1141 
1142  side = tTopo->pxfSide(detId);
1143  disk = tTopo->pxfDisk(detId);
1144  blade = tTopo->pxfBlade(detId);
1145  panel = tTopo->pxfPanel(detId);
1146  plaq = tTopo->pxfModule(detId); // also known as plaquette
1147 
1148  if ( side==1 )
1149  {
1150  if ( panel==1 )
1151  {
1152  mePosxZmPanel1_all_hits->Fill( rechitx );
1153  mePosyZmPanel1_all_hits->Fill( rechity );
1154  }
1155  else if ( panel==2 )
1156  {
1157  mePosxZmPanel2_all_hits->Fill( rechitx );
1158  mePosyZmPanel2_all_hits->Fill( rechity );
1159  }
1160  else edm::LogWarning("SiPixelTrackingRecHitsValid") << "..............................................Wrong panel number !";
1161  } // if ( side==1 )
1162  else if ( side==2 )
1163  {
1164  if ( panel==1 )
1165  {
1166  mePosxZpPanel1_all_hits->Fill( rechitx );
1167  mePosyZpPanel1_all_hits->Fill( rechity );
1168  }
1169  else if ( panel==2 )
1170  {
1171  mePosxZpPanel2_all_hits->Fill( rechitx );
1172  mePosyZpPanel2_all_hits->Fill( rechity );
1173  }
1174  else edm::LogWarning("SiPixelTrackingRecHitsValid")<< "..............................................Wrong panel number !";
1175  } //else if ( side==2 )
1176  else edm::LogWarning("SiPixelTrackingRecHitsValid") << ".......................................................Wrong side !" ;
1177 
1178  } // else if ( detId.subdetId()==PixelSubdetector::PixelEndcap )
1179  else edm::LogWarning("SiPixelTrackingRecHitsValid") << "Pixel rechit collection but we are not in the pixel detector" << (int)detId.subdetId() ;
1180 
1181  }
1182  }
1183  // ------------------------------------------------ all hits ---------------------------------------------------------------
1184 
1185  // Get tracks
1186  edm::Handle<reco::TrackCollection> trackCollection;
1187  e.getByToken( recoTrackCollectionToken_, trackCollection );
1188  const reco::TrackCollection *tracks = trackCollection.product();
1189  reco::TrackCollection::const_iterator tciter;
1190 
1191  int n_tracks = (int)tracks->size(); // number of tracks in this event
1192  meTracksPerEvent->Fill( n_tracks );
1193 
1194  if ( tracks->size() > 0 )
1195  {
1196  // Loop on tracks
1197  for ( tciter=tracks->begin(); tciter!=tracks->end(); tciter++)
1198  {
1199  phi = tciter->momentum().phi() / math_pi*180.0;
1200  eta = tciter->momentum().eta();
1201 
1202  int n_hits = 0;
1203  // First loop on hits: find matched hits
1204  for ( trackingRecHit_iterator it = tciter->recHitsBegin(); it != tciter->recHitsEnd(); it++)
1205  {
1206  const TrackingRecHit &thit = **it;
1207  // Is it a matched hit?
1208  const SiPixelRecHit* matchedhit = dynamic_cast<const SiPixelRecHit*>(&thit);
1209 
1210  if ( matchedhit )
1211  {
1212  ++n_hits;
1213 
1214  layer = -9999;
1215  ladder = -9999;
1216  mod = -9999;
1217  side = -9999;
1218  disk = -9999;
1219  blade = -9999;
1220  panel = -9999;
1221  plaq = -9999;
1222 
1223  rechitx = -9999.9;
1224  rechity = -9999.9;
1225  rechitz = -9999.9;
1226  rechiterrx = -9999.9;
1227  rechiterry = -9999.9;
1228  rechitresx = -9999.9;
1229  rechitresy = -9999.9;
1230  rechitpullx = -9999.9;
1231  rechitpully = -9999.9;
1232 
1233  npix = -9999;
1234  nxpix = -9999;
1235  nypix = -9999;
1236  charge = -9999.9;
1237 
1238  alpha = -9999.9;
1239  beta = -9999.9;
1240 
1241  half = -9999;
1242  flipped = -9999;
1243 
1244  nsimhit = -9999;
1245 
1246  simhitx = -9999.9;
1247  simhity = -9999.9;
1248 
1249  position = (*it)->localPosition();
1250  error = (*it)->localPositionError();
1251 
1252  rechitx = position.x();
1253  rechity = position.y();
1254  rechitz = position.z();
1255  rechiterrx = sqrt(error.xx());
1256  rechiterry = sqrt(error.yy());
1257 
1258  npix = (*matchedhit).cluster()->size();
1259  nxpix = (*matchedhit).cluster()->sizeX();
1260  nypix = (*matchedhit).cluster()->sizeY();
1261  charge = (*matchedhit).cluster()->charge();
1262 
1263  //Association of the rechit to the simhit
1264  matched.clear();
1265  matched = associate.associateHit(*matchedhit);
1266 
1267  nsimhit = (int)matched.size();
1268 
1269  if ( !matched.empty() )
1270  {
1271  mindist = 999999.9;
1272  float distx, disty, dist;
1273  bool found_hit_from_generated_particle = false;
1274 
1275  int n_assoc_muon = 0;
1276 
1277  std::vector<PSimHit>::const_iterator closestit = matched.begin();
1278  for (std::vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); m++)
1279  {
1280  if ( checkType_ )
1281  {
1282  int pid = (*m).particleType();
1283  if ( abs(pid) != genType_ )
1284  continue;
1285 
1286  }
1287 
1288  float simhitx = 0.5 * ( (*m).entryPoint().x() + (*m).exitPoint().x() );
1289  float simhity = 0.5 * ( (*m).entryPoint().y() + (*m).exitPoint().y() );
1290 
1291  distx = fabs(rechitx - simhitx);
1292  disty = fabs(rechity - simhity);
1293  dist = sqrt( distx*distx + disty*disty );
1294 
1295  if ( dist < mindist )
1296  {
1297  n_assoc_muon++;
1298 
1299  mindist = dist;
1300  closestit = m;
1301  found_hit_from_generated_particle = true;
1302  }
1303  } // for (std::vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); m++)
1304 
1305  // This recHit does not have any simHit with the same particleType as the particles generated
1306  // Ignore it as most probably come from delta rays.
1307  if ( checkType_ && !found_hit_from_generated_particle )
1308  continue;
1309 
1310  if ( n_assoc_muon > 1 )
1311  {
1312  edm::LogWarning("SiPixelTrackingRecHitsValid") << " ----- This is not good: n_assoc_muon = " << n_assoc_muon ;
1313  edm::LogWarning("SiPixelTrackingRecHitsValid") << "evt = " << evt ;
1314  }
1315 
1316  pidhit = (*closestit).particleType();
1317 
1318  simhitx = 0.5*( (*closestit).entryPoint().x() + (*closestit).exitPoint().x() );
1319  simhity = 0.5*( (*closestit).entryPoint().y() + (*closestit).exitPoint().y() );
1320 
1321  rechitresx = rechitx - simhitx;
1322  rechitresy = rechity - simhity;
1323  rechitpullx = ( rechitx - simhitx ) / sqrt(error.xx());
1324  rechitpully = ( rechity - simhity ) / sqrt(error.yy());
1325 
1326  float simhitpx = (*closestit).momentumAtEntry().x();
1327  float simhitpy = (*closestit).momentumAtEntry().y();
1328  float simhitpz = (*closestit).momentumAtEntry().z();
1329 
1330  //beta = atan2(simhitpz, simhitpy) * radtodeg;
1331  //alpha = atan2(simhitpz, simhitpx) * radtodeg;
1332 
1333  beta = fabs(atan2(simhitpz, simhitpy)) * radtodeg;
1334  alpha = fabs(atan2(simhitpz, simhitpx)) * radtodeg;
1335 
1336  detId = (*it)->geographicalId();
1337 
1338  subdetId = (int)detId.subdetId();
1339 
1340  if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelBarrel )
1341  {
1342  mePosxBarrel->Fill( rechitx );
1343  mePosyBarrel->Fill( rechity );
1344  meErrxBarrel->Fill( rechiterrx );
1345  meErryBarrel->Fill( rechiterry );
1346  meResxBarrel->Fill( rechitresx );
1347  meResyBarrel->Fill( rechitresy );
1348  mePullxBarrel->Fill( rechitpullx );
1349  mePullyBarrel->Fill( rechitpully );
1350  meNpixBarrel->Fill( npix );
1351  meNxpixBarrel->Fill( nxpix );
1352  meNypixBarrel->Fill( nypix );
1353  meChargeBarrel->Fill( charge );
1354  meResXvsAlphaBarrel->Fill( alpha, fabs(rechitresx) );
1355  meResYvsAlphaBarrel->Fill( alpha, fabs(rechitresy) );
1356  meResXvsBetaBarrel->Fill( beta, fabs(rechitresx) );
1357  meResYvsBetaBarrel->Fill( beta, fabs(rechitresy) );
1358  mePullXvsAlphaBarrel->Fill( alpha, rechitpullx );
1359  mePullYvsAlphaBarrel->Fill( alpha, rechitpully );
1360  mePullXvsBetaBarrel->Fill( beta, rechitpullx );
1361  mePullYvsBetaBarrel->Fill( beta, rechitpully );
1362  mePullXvsPhiBarrel->Fill( phi, rechitpullx );
1363  mePullYvsPhiBarrel->Fill( phi, rechitpully );
1364  mePullXvsEtaBarrel->Fill( eta, rechitpullx );
1365  mePullYvsEtaBarrel->Fill( eta, rechitpully );
1366 
1367  const PixelGeomDetUnit * theGeomDet
1368  = dynamic_cast<const PixelGeomDetUnit*> ( tracker->idToDet(detId) );
1369  //const PixelTopology * topol = (&(theGeomDet->specificTopology()));
1370 
1371  int tmp_nrows = theGeomDet->specificTopology().nrows();
1372 
1373  if ( tmp_nrows == 80 )
1374  {
1375  mePosxBarrelHalfModule->Fill( rechitx );
1376  mePosyBarrelHalfModule->Fill( rechity );
1377  half = 1;
1378  }
1379  else if ( tmp_nrows == 160 )
1380  {
1381  mePosxBarrelFullModule->Fill( rechitx );
1382  mePosyBarrelFullModule->Fill( rechity );
1383  half = 0;
1384  }
1385  else
1386  edm::LogWarning("SiPixelTrackingRecHitsValid") << "-------------------------------------------------- Wrong module size !!!";
1387 
1388  float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,0.)).perp();
1389  float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,1.)).perp();
1390 
1391  if ( tmp2<tmp1 )
1392  { // flipped
1393  mePosxBarrelFlippedLadders->Fill( rechitx );
1394  mePosyBarrelFlippedLadders->Fill( rechity );
1395  flipped = 1;
1396 
1397  meResXvsAlphaBarrelFlippedLadders->Fill( alpha, fabs(rechitresx) );
1398  meResYvsAlphaBarrelFlippedLadders->Fill( alpha, fabs(rechitresy) );
1399  meResXvsBetaBarrelFlippedLadders->Fill( beta, fabs(rechitresx) );
1400  meResYvsBetaBarrelFlippedLadders->Fill( beta, fabs(rechitresy) );
1401  mePullXvsAlphaBarrelFlippedLadders->Fill( alpha, rechitpullx );
1402  mePullYvsAlphaBarrelFlippedLadders->Fill( alpha, rechitpully );
1403  mePullXvsBetaBarrelFlippedLadders->Fill( beta, rechitpullx );
1404  mePullYvsBetaBarrelFlippedLadders->Fill( beta, rechitpully );
1405  mePullXvsPhiBarrelFlippedLadders->Fill( phi, rechitpullx );
1406  mePullYvsPhiBarrelFlippedLadders->Fill( phi, rechitpully );
1407  mePullXvsEtaBarrelFlippedLadders->Fill( eta, rechitpullx );
1408  mePullYvsEtaBarrelFlippedLadders->Fill( eta, rechitpully );
1409 
1410  meWPullXvsAlphaBarrelFlippedLadders->Fill( alpha, fabs(rechitpullx) );
1411  meWPullYvsAlphaBarrelFlippedLadders->Fill( alpha, fabs(rechitpully) );
1412  meWPullXvsBetaBarrelFlippedLadders->Fill( beta, fabs(rechitpullx) );
1413  meWPullYvsBetaBarrelFlippedLadders->Fill( beta, fabs(rechitpully) );
1414  }
1415  else
1416  { // not flipped
1417  mePosxBarrelNonFlippedLadders->Fill( rechitx );
1418  mePosyBarrelNonFlippedLadders->Fill( rechity );
1419  flipped = 0;
1420 
1421  meResXvsAlphaBarrelNonFlippedLadders->Fill( alpha, fabs(rechitresx) );
1422  meResYvsAlphaBarrelNonFlippedLadders->Fill( alpha, fabs(rechitresy) );
1423  meResXvsBetaBarrelNonFlippedLadders->Fill( beta, fabs(rechitresx) );
1424  meResYvsBetaBarrelNonFlippedLadders->Fill( beta, fabs(rechitresy) );
1425  mePullXvsAlphaBarrelNonFlippedLadders->Fill( alpha, rechitpullx );
1426  mePullYvsAlphaBarrelNonFlippedLadders->Fill( alpha, rechitpully );
1427  mePullXvsBetaBarrelNonFlippedLadders->Fill( beta, rechitpullx );
1428  mePullYvsBetaBarrelNonFlippedLadders->Fill( beta, rechitpully );
1429  mePullXvsPhiBarrelNonFlippedLadders->Fill( phi, rechitpullx );
1430  mePullYvsPhiBarrelNonFlippedLadders->Fill( phi, rechitpully );
1431  mePullXvsEtaBarrelNonFlippedLadders->Fill( eta, rechitpullx );
1432  mePullYvsEtaBarrelNonFlippedLadders->Fill( eta, rechitpully );
1433 
1434  meWPullXvsAlphaBarrelNonFlippedLadders->Fill( alpha, fabs(rechitpullx) );
1435  meWPullYvsAlphaBarrelNonFlippedLadders->Fill( alpha, fabs(rechitpully) );
1436  meWPullXvsBetaBarrelNonFlippedLadders->Fill( beta, fabs(rechitpullx) );
1437  meWPullYvsBetaBarrelNonFlippedLadders->Fill( beta, fabs(rechitpully) );
1438  }
1439 
1440 
1441  layer = tTopo->pxbLayer(detId); // Layer: 1,2,3.
1442  ladder = tTopo->pxbLadder(detId); // Ladder: 1-20, 32, 44.
1443  mod = tTopo->pxbModule(detId); // Mod: 1-8.
1444 
1445  mePosxBarrelLayerModule[layer-1][mod-1]->Fill( rechitx );
1446  mePosyBarrelLayerModule[layer-1][mod-1]->Fill( rechity );
1447  meErrxBarrelLayerModule[layer-1][mod-1]->Fill( rechiterrx );
1448  meErryBarrelLayerModule[layer-1][mod-1]->Fill( rechiterry );
1449  meResxBarrelLayerModule[layer-1][mod-1]->Fill( rechitresx );
1450  meResyBarrelLayerModule[layer-1][mod-1]->Fill( rechitresy );
1451  mePullxBarrelLayerModule[layer-1][mod-1]->Fill( rechitpullx );
1452  mePullyBarrelLayerModule[layer-1][mod-1]->Fill( rechitpully );
1453  meNpixBarrelLayerModule[layer-1][mod-1]->Fill( npix );
1454  meNxpixBarrelLayerModule[layer-1][mod-1]->Fill( nxpix );
1455  meNypixBarrelLayerModule[layer-1][mod-1]->Fill( nypix );
1456  meChargeBarrelLayerModule[layer-1][mod-1]->Fill( charge );
1457  meResXvsAlphaBarrelLayerModule[layer-1][mod-1]->Fill( alpha, fabs(rechitresx) );
1458  meResYvsAlphaBarrelLayerModule[layer-1][mod-1]->Fill( alpha, fabs(rechitresy) );
1459  meResXvsBetaBarrelLayerModule[layer-1][mod-1]->Fill( beta, fabs(rechitresx) );
1460  meResYvsBetaBarrelLayerModule[layer-1][mod-1]->Fill( beta, fabs(rechitresy) );
1461  mePullXvsAlphaBarrelLayerModule[layer-1][mod-1]->Fill( alpha, rechitpullx );
1462  mePullYvsAlphaBarrelLayerModule[layer-1][mod-1]->Fill( alpha, rechitpully );
1463  mePullXvsBetaBarrelLayerModule[layer-1][mod-1]->Fill( beta, rechitpullx );
1464  mePullYvsBetaBarrelLayerModule[layer-1][mod-1]->Fill( beta, rechitpully );
1465  mePullXvsPhiBarrelLayerModule[layer-1][mod-1]->Fill( phi, rechitpullx );
1466  mePullYvsPhiBarrelLayerModule[layer-1][mod-1]->Fill( phi, rechitpully );
1467  mePullXvsEtaBarrelLayerModule[layer-1][mod-1]->Fill( eta, rechitpullx );
1468  mePullYvsEtaBarrelLayerModule[layer-1][mod-1]->Fill( eta, rechitpully );
1469 
1470  meResxBarrelLayer[layer-1]->Fill( rechitresx );
1471  meResyBarrelLayer[layer-1]->Fill( rechitresy );
1472  mePullxBarrelLayer[layer-1]->Fill( rechitpullx );
1473  mePullyBarrelLayer[layer-1]->Fill( rechitpully );
1474 
1475  if ( tmp2<tmp1 )
1476  { // flipped
1477  meResXvsAlphaBarrelFlippedLaddersLayer[layer-1]->Fill( alpha, fabs(rechitresx) );
1478  meResYvsAlphaBarrelFlippedLaddersLayer[layer-1]->Fill( alpha, fabs(rechitresy) );
1479  meResXvsBetaBarrelFlippedLaddersLayer[layer-1]->Fill( beta, fabs(rechitresx) );
1480  meResYvsBetaBarrelFlippedLaddersLayer[layer-1]->Fill( beta, fabs(rechitresy) );
1481  }
1482  else
1483  { // not flipped
1484  meResXvsAlphaBarrelNonFlippedLaddersLayer[layer-1]->Fill( alpha, fabs(rechitresx) );
1485  meResYvsAlphaBarrelNonFlippedLaddersLayer[layer-1]->Fill( alpha, fabs(rechitresy) );
1486  meResXvsBetaBarrelNonFlippedLaddersLayer[layer-1]->Fill( beta, fabs(rechitresx) );
1487  meResYvsBetaBarrelNonFlippedLaddersLayer[layer-1]->Fill( beta, fabs(rechitresy) );
1488  }
1489 
1490  }
1491  else if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelEndcap )
1492  {
1493 
1494  side = tTopo->pxfSide(detId);
1495  disk = tTopo->pxfDisk(detId);
1496  blade = tTopo->pxfBlade(detId);
1497  panel = tTopo->pxfPanel(detId);
1498  plaq = tTopo->pxfModule(detId); // also known as plaquette
1499 
1500  if ( side==1 )
1501  {
1502  if ( panel==1 )
1503  {
1504  mePosxZmPanel1->Fill( rechitx );
1505  mePosyZmPanel1->Fill( rechity );
1506  meErrxZmPanel1->Fill( rechiterrx );
1507  meErryZmPanel1->Fill( rechiterry );
1508  meResxZmPanel1->Fill( rechitresx );
1509  meResyZmPanel1->Fill( rechitresy );
1510  mePullxZmPanel1->Fill( rechitpullx );
1511  mePullyZmPanel1->Fill( rechitpully );
1512  meNpixZmPanel1->Fill( npix );
1513  meNxpixZmPanel1->Fill( nxpix );
1514  meNypixZmPanel1->Fill( nypix );
1515  meChargeZmPanel1->Fill( charge );
1516  meResXvsAlphaZmPanel1->Fill( alpha, fabs(rechitresx) );
1517  meResYvsAlphaZmPanel1->Fill( alpha, fabs(rechitresy) );
1518  meResXvsBetaZmPanel1->Fill( beta, fabs(rechitresx) );
1519  meResYvsBetaZmPanel1->Fill( beta, fabs(rechitresy) );
1520  mePullXvsAlphaZmPanel1->Fill( alpha, rechitpullx );
1521  mePullYvsAlphaZmPanel1->Fill( alpha, rechitpully );
1522  mePullXvsBetaZmPanel1->Fill( beta, rechitpullx );
1523  mePullYvsBetaZmPanel1->Fill( beta, rechitpully );
1524  mePullXvsPhiZmPanel1->Fill( phi, rechitpullx );
1525  mePullYvsPhiZmPanel1->Fill( phi, rechitpully );
1526  mePullXvsEtaZmPanel1->Fill( eta, rechitpullx );
1527  mePullYvsEtaZmPanel1->Fill( eta, rechitpully );
1528 
1529  meWPullXvsAlphaZmPanel1->Fill( alpha, fabs(rechitpullx) );
1530  meWPullYvsAlphaZmPanel1->Fill( alpha, fabs(rechitpully) );
1531  meWPullXvsBetaZmPanel1->Fill( beta, fabs(rechitpullx) );
1532  meWPullYvsBetaZmPanel1->Fill( beta, fabs(rechitpully) );
1533 
1534  mePosxZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitx );
1535  mePosyZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechity );
1536  meErrxZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechiterrx );
1537  meErryZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechiterry );
1538  meResxZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitresx );
1539  meResyZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitresy );
1540  mePullxZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitpullx );
1541  mePullyZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitpully );
1542  meNpixZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( npix );
1543  meNxpixZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( nxpix );
1544  meNypixZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( nypix );
1545  meChargeZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( charge );
1546  meResXvsAlphaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresx) );
1547  meResYvsAlphaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresy) );
1548  meResXvsBetaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresx) );
1549  meResYvsBetaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresy) );
1550  mePullXvsAlphaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpullx );
1551  mePullYvsAlphaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpully );
1552  mePullXvsBetaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpullx );
1553  mePullYvsBetaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpully );
1554  mePullXvsPhiZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpullx );
1555  mePullYvsPhiZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpully );
1556  mePullXvsEtaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpullx );
1557  mePullYvsEtaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpully );
1558 
1559  }
1560  else if ( panel==2 )
1561  {
1562  mePosxZmPanel2->Fill( rechitx );
1563  mePosyZmPanel2->Fill( rechity );
1564  meErrxZmPanel2->Fill( rechiterrx );
1565  meErryZmPanel2->Fill( rechiterry );
1566  meResxZmPanel2->Fill( rechitresx );
1567  meResyZmPanel2->Fill( rechitresy );
1568  mePullxZmPanel2->Fill( rechitpullx );
1569  mePullyZmPanel2->Fill( rechitpully );
1570  meNpixZmPanel2->Fill( npix );
1571  meNxpixZmPanel2->Fill( nxpix );
1572  meNypixZmPanel2->Fill( nypix );
1573  meChargeZmPanel2->Fill( charge );
1574  meResXvsAlphaZmPanel2->Fill( alpha, fabs(rechitresx) );
1575  meResYvsAlphaZmPanel2->Fill( alpha, fabs(rechitresy) );
1576  meResXvsBetaZmPanel2->Fill( beta, fabs(rechitresx) );
1577  meResYvsBetaZmPanel2->Fill( beta, fabs(rechitresy) );
1578  mePullXvsAlphaZmPanel2->Fill( alpha, rechitpullx );
1579  mePullYvsAlphaZmPanel2->Fill( alpha, rechitpully );
1580  mePullXvsBetaZmPanel2->Fill( beta, rechitpullx );
1581  mePullYvsBetaZmPanel2->Fill( beta, rechitpully );
1582  mePullXvsPhiZmPanel2->Fill( phi, rechitpullx );
1583  mePullYvsPhiZmPanel2->Fill( phi, rechitpully );
1584  mePullXvsEtaZmPanel2->Fill( eta, rechitpullx );
1585  mePullYvsEtaZmPanel2->Fill( eta, rechitpully );
1586 
1587  meWPullXvsAlphaZmPanel2->Fill( alpha, fabs(rechitpullx) );
1588  meWPullYvsAlphaZmPanel2->Fill( alpha, fabs(rechitpully) );
1589  meWPullXvsBetaZmPanel2->Fill( beta, fabs(rechitpullx) );
1590  meWPullYvsBetaZmPanel2->Fill( beta, fabs(rechitpully) );
1591 
1592  mePosxZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitx );
1593  mePosyZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechity );
1594  meErrxZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechiterrx );
1595  meErryZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechiterry );
1596  meResxZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitresx );
1597  meResyZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitresy );
1598  mePullxZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitpullx );
1599  mePullyZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitpully );
1600  meNpixZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( npix );
1601  meNxpixZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( nxpix );
1602  meNypixZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( nypix );
1603  meChargeZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( charge );
1604  meResXvsAlphaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresx) );
1605  meResYvsAlphaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresy) );
1606  meResXvsBetaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresx) );
1607  meResYvsBetaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresy) );
1608  mePullXvsAlphaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpullx );
1609  mePullYvsAlphaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpully );
1610  mePullXvsBetaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpullx );
1611  mePullYvsBetaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpully );
1612  mePullXvsPhiZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpullx );
1613  mePullYvsPhiZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpully );
1614  mePullXvsEtaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpullx );
1615  mePullYvsEtaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpully );
1616 
1617  }
1618  else edm::LogWarning("SiPixelTrackingRecHitsValid") << "..............................................Wrong panel number !";
1619  } // if ( side==1 )
1620  else if ( side==2 )
1621  {
1622  if ( panel==1 )
1623  {
1624  mePosxZpPanel1->Fill( rechitx );
1625  mePosyZpPanel1->Fill( rechity );
1626  meErrxZpPanel1->Fill( rechiterrx );
1627  meErryZpPanel1->Fill( rechiterry );
1628  meResxZpPanel1->Fill( rechitresx );
1629  meResyZpPanel1->Fill( rechitresy );
1630  mePullxZpPanel1->Fill( rechitpullx );
1631  mePullyZpPanel1->Fill( rechitpully );
1632  meNpixZpPanel1->Fill( npix );
1633  meNxpixZpPanel1->Fill( nxpix );
1634  meNypixZpPanel1->Fill( nypix );
1635  meChargeZpPanel1->Fill( charge );
1636  meResXvsAlphaZpPanel1->Fill( alpha, fabs(rechitresx) );
1637  meResYvsAlphaZpPanel1->Fill( alpha, fabs(rechitresy) );
1638  meResXvsBetaZpPanel1->Fill( beta, fabs(rechitresx) );
1639  meResYvsBetaZpPanel1->Fill( beta, fabs(rechitresy) );
1640  mePullXvsAlphaZpPanel1->Fill( alpha, rechitpullx );
1641  mePullYvsAlphaZpPanel1->Fill( alpha, rechitpully );
1642  mePullXvsBetaZpPanel1->Fill( beta, rechitpullx );
1643  mePullYvsBetaZpPanel1->Fill( beta, rechitpully );
1644  mePullXvsPhiZpPanel1->Fill( phi, rechitpullx );
1645  mePullYvsPhiZpPanel1->Fill( phi, rechitpully );
1646  mePullXvsEtaZpPanel1->Fill( eta, rechitpullx );
1647  mePullYvsEtaZpPanel1->Fill( eta, rechitpully );
1648 
1649  meWPullXvsAlphaZpPanel1->Fill( alpha, fabs(rechitpullx) );
1650  meWPullYvsAlphaZpPanel1->Fill( alpha, fabs(rechitpully) );
1651  meWPullXvsBetaZpPanel1->Fill( beta, fabs(rechitpullx) );
1652  meWPullYvsBetaZpPanel1->Fill( beta, fabs(rechitpully) );
1653 
1654  mePosxZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitx );
1655  mePosyZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechity );
1656  meErrxZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechiterrx );
1657  meErryZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechiterry );
1658  meResxZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitresx );
1659  meResyZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitresy );
1660  mePullxZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitpullx );
1661  mePullyZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitpully );
1662  meNpixZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( npix );
1663  meNxpixZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( nxpix );
1664  meNypixZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( nypix );
1665  meChargeZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( charge );
1666  meResXvsAlphaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresx) );
1667  meResYvsAlphaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresy) );
1668  meResXvsBetaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresx) );
1669  meResYvsBetaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresy) );
1670  mePullXvsAlphaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpullx );
1671  mePullYvsAlphaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpully );
1672  mePullXvsBetaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpullx );
1673  mePullYvsBetaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpully );
1674  mePullXvsPhiZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpullx );
1675  mePullYvsPhiZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpully );
1676  mePullXvsEtaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpullx );
1677  mePullYvsEtaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpully );
1678 
1679  }
1680  else if ( panel==2 )
1681  {
1682  mePosxZpPanel2->Fill( rechitx );
1683  mePosyZpPanel2->Fill( rechity );
1684  meErrxZpPanel2->Fill( rechiterrx );
1685  meErryZpPanel2->Fill( rechiterry );
1686  meResxZpPanel2->Fill( rechitresx );
1687  meResyZpPanel2->Fill( rechitresy );
1688  mePullxZpPanel2->Fill( rechitpullx );
1689  mePullyZpPanel2->Fill( rechitpully );
1690  meNpixZpPanel2->Fill( npix );
1691  meNxpixZpPanel2->Fill( nxpix );
1692  meNypixZpPanel2->Fill( nypix );
1693  meChargeZpPanel2->Fill( charge );
1694  meResXvsAlphaZpPanel2->Fill( alpha, fabs(rechitresx) );
1695  meResYvsAlphaZpPanel2->Fill( alpha, fabs(rechitresy) );
1696  meResXvsBetaZpPanel2->Fill( beta, fabs(rechitresx) );
1697  meResYvsBetaZpPanel2->Fill( beta, fabs(rechitresy) );
1698  mePullXvsAlphaZpPanel2->Fill( alpha, rechitpullx );
1699  mePullYvsAlphaZpPanel2->Fill( alpha, rechitpully );
1700  mePullXvsBetaZpPanel2->Fill( beta, rechitpullx );
1701  mePullYvsBetaZpPanel2->Fill( beta, rechitpully );
1702  mePullXvsPhiZpPanel2->Fill( phi, rechitpullx );
1703  mePullYvsPhiZpPanel2->Fill( phi, rechitpully );
1704  mePullXvsEtaZpPanel2->Fill( eta, rechitpullx );
1705  mePullYvsEtaZpPanel2->Fill( eta, rechitpully );
1706 
1707  meWPullXvsAlphaZpPanel2->Fill( alpha, fabs(rechitpullx) );
1708  meWPullYvsAlphaZpPanel2->Fill( alpha, fabs(rechitpully) );
1709  meWPullXvsBetaZpPanel2->Fill( beta, fabs(rechitpullx) );
1710  meWPullYvsBetaZpPanel2->Fill( beta, fabs(rechitpully) );
1711 
1712  mePosxZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitx );
1713  mePosyZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechity );
1714  meErrxZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechiterrx );
1715  meErryZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechiterry );
1716  meResxZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitresx );
1717  meResyZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitresy );
1718  mePullxZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitpullx );
1719  mePullyZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitpully );
1720  meNpixZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( npix );
1721  meNxpixZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( nxpix );
1722  meNypixZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( nypix );
1723  meChargeZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( charge );
1724  meResXvsAlphaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresx) );
1725  meResYvsAlphaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresy) );
1726  meResXvsBetaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresx) );
1727  meResYvsBetaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresy) );
1728  mePullXvsAlphaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpullx );
1729  mePullYvsAlphaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpully );
1730  mePullXvsBetaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpullx );
1731  mePullYvsBetaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpully );
1732  mePullXvsPhiZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpullx );
1733  mePullYvsPhiZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpully );
1734  mePullXvsEtaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpullx );
1735  mePullYvsEtaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpully );
1736 
1737  }
1738  else edm::LogWarning("SiPixelTrackingRecHitsValid") << "..............................................Wrong panel number !";
1739  } //else if ( side==2 )
1740  else edm::LogWarning("SiPixelTrackingRecHitsValid") << ".......................................................Wrong side !" ;
1741 
1742  } // else if ( detId.subdetId()==PixelSubdetector::PixelEndcap )
1743  else edm::LogWarning("SiPixelTrackingRecHitsValid") << "Pixel rechit but we are not in the pixel detector" << (int)detId.subdetId() ;
1744 
1745  if(debugNtuple_.size()!=0)t_->Fill();
1746 
1747  } // if ( !matched.empty() )
1748  //else
1749  //cout << "---------------- RecHit with no associated SimHit !!! -------------------------- " << endl;
1750 
1751  } // matchedhit.
1752 
1753  } // end of loop on hits
1754 
1755  mePixRecHitsPerTrack->Fill( n_hits );
1756  //cout << "n_hits = " << n_hits << endl;
1757 
1758  } //end of loop on track
1759 
1760  } // tracks > 0.
1761 
1762  } //end of MTCCTrack
1763 
1764 }
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:42
SiPixelTrackingRecHitsValid::meWPullXvsAlphaBarrelNonFlippedLadders
MonitorElement * meWPullXvsAlphaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:191
Surface::toGlobal
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:114
SiPixelTrackingRecHitsValid::siPixelRecHitCollectionToken_
edm::EDGetTokenT< SiPixelRecHitCollection > siPixelRecHitCollectionToken_
Definition: SiPixelTrackingRecHitsValid.h:92
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:44
SiPixelTrackingRecHitsValid::meResXvsAlphaBarrelNonFlippedLaddersLayer
MonitorElement * meResXvsAlphaBarrelNonFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:206
SiPixelTrackingRecHitsValid::meErryBarrelLayerModule
MonitorElement * meErryBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:104
SiPixelTrackingRecHitsValid::meNpixZmPanel1DiskPlaq
MonitorElement * meNpixZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:226
SiPixelTrackingRecHitsValid::meResYvsAlphaBarrelFlippedLaddersLayer
MonitorElement * meResYvsAlphaBarrelFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:203
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
SiPixelTrackingRecHitsValid::mePullYvsEtaBarrelLayerModule
MonitorElement * mePullYvsEtaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:124
i
int i
Definition: DBlmapReader.cc:9
SiPixelTrackingRecHitsValid::meErryZpPanel2DiskPlaq
MonitorElement * meErryZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:300
SiPixelTrackingRecHitsValid::meResyZpPanel1DiskPlaq
MonitorElement * meResyZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:276
LocalError::xx
float xx() const
Definition: LocalError.h:24
SiPixelTrackingRecHitsValid::mePullxZpPanel2DiskPlaq
MonitorElement * mePullxZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:303
SiPixelTrackingRecHitsValid::mePullXvsBetaBarrelLayerModule
MonitorElement * mePullXvsBetaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:119
SiPixelTrackingRecHitsValid::meResyZmPanel2DiskPlaq
MonitorElement * meResyZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:249
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:117
SiPixelTrackingRecHitsValid::mePullXvsAlphaZmPanel1DiskPlaq
MonitorElement * mePullXvsAlphaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:234
SiPixelTrackingRecHitsValid::meResYvsAlphaBarrelNonFlippedLaddersLayer
MonitorElement * meResYvsAlphaBarrelNonFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:207
SiPixelTrackingRecHitsValid::mePullYvsAlphaBarrelNonFlippedLadders
MonitorElement * mePullYvsAlphaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:179
SiPixelTrackingRecHitsValid::mePullYvsAlphaZmPanel1DiskPlaq
MonitorElement * mePullYvsAlphaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:235
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:72
SiPixelTrackingRecHitsValid::meChargeZmPanel2DiskPlaq
MonitorElement * meChargeZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:255
SiPixelTrackingRecHitsValid::mePosxZmPanel1DiskPlaq
MonitorElement * mePosxZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:218
SiPixelTrackingRecHitsValid::mePosyBarrelLayerModule
MonitorElement * mePosyBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:102
SiPixelTrackingRecHitsValid::meResxZpPanel1DiskPlaq
MonitorElement * meResxZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:275
SiPixelTrackingRecHitsValid::meNpixZpPanel2DiskPlaq
MonitorElement * meNpixZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:305
SiPixelTrackingRecHitsValid::meResyZmPanel1DiskPlaq
MonitorElement * meResyZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:223
SiPixelTrackingRecHitsValid::meResXvsBetaZmPanel1DiskPlaq
MonitorElement * meResXvsBetaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:232
SiPixelTrackingRecHitsValid::mePullYvsEtaZpPanel1DiskPlaq
MonitorElement * mePullYvsEtaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:294
SiPixelTrackingRecHitsValid::meWPullXvsBetaBarrelNonFlippedLadders
MonitorElement * meWPullXvsBetaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:193
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434
SiPixelTrackingRecHitsValid::meResYvsBetaZmPanel1DiskPlaq
MonitorElement * meResYvsBetaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:233
SiPixelTrackingRecHitsValid::mePullYvsEtaZmPanel2DiskPlaq
MonitorElement * mePullYvsEtaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:267
SiPixelTrackingRecHitsValid::meErrxZmPanel1DiskPlaq
MonitorElement * meErrxZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:220
SiPixelTrackingRecHitsValid::meResXvsBetaBarrelNonFlippedLaddersLayer
MonitorElement * meResXvsBetaBarrelNonFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:208
SiPixelTrackingRecHitsValid::meErrxZpPanel2DiskPlaq
MonitorElement * meErrxZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:299
SiPixelTrackingRecHitsValid::meResYvsBetaBarrelFlippedLaddersLayer
MonitorElement * meResYvsBetaBarrelFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:205
SiPixelTrackingRecHitsValid::mePullxBarrelLayerModule
MonitorElement * mePullxBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:107
SiPixelTrackingRecHitsValid::mePosyZpPanel1DiskPlaq
MonitorElement * mePosyZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:272
SiPixelTrackingRecHitsValid::meNxpixZpPanel1DiskPlaq
MonitorElement * meNxpixZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:280
SiPixelTrackingRecHitsValid::mePullYvsBetaZpPanel1DiskPlaq
MonitorElement * mePullYvsBetaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:290
SiPixelTrackingRecHitsValid::mePullXvsPhiZmPanel1DiskPlaq
MonitorElement * mePullXvsPhiZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:238
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:10
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
SiPixelTrackingRecHitsValid::meResXvsBetaBarrelFlippedLaddersLayer
MonitorElement * meResXvsBetaBarrelFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:204
SiPixelTrackingRecHitsValid::mePullyZmPanel1DiskPlaq
MonitorElement * mePullyZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:225
SiPixelTrackingRecHitsValid::meErrxBarrelLayerModule
MonitorElement * meErrxBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:103
SiPixelTrackingRecHitsValid::meErryZmPanel2DiskPlaq
MonitorElement * meErryZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:247
PixelTopology::nrows
virtual int nrows() const =0
SiPixelTrackingRecHitsValid::mePosxZmPanel2DiskPlaq
MonitorElement * mePosxZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:244
SiPixelTrackingRecHitsValid::meChargeZpPanel2DiskPlaq
MonitorElement * meChargeZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:308
SiPixelTrackingRecHitsValid::mePullXvsEtaZpPanel2DiskPlaq
MonitorElement * mePullXvsEtaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:319
SiPixelTrackingRecHitsValid::meResYvsAlphaZmPanel1DiskPlaq
MonitorElement * meResYvsAlphaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:231
SiPixelTrackingRecHitsValid::mePullXvsBetaZmPanel1DiskPlaq
MonitorElement * mePullXvsBetaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:236
SiPixelTrackingRecHitsValid::meResXvsAlphaZmPanel2DiskPlaq
MonitorElement * meResXvsAlphaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:256
SiPixelTrackingRecHitsValid::mePullYvsAlphaZmPanel2DiskPlaq
MonitorElement * mePullYvsAlphaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:261
SiPixelTrackingRecHitsValid::mePullxZmPanel1DiskPlaq
MonitorElement * mePullxZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:224
SiPixelTrackingRecHitsValid::mePullyZpPanel2DiskPlaq
MonitorElement * mePullyZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:304
SiPixelTrackingRecHitsValid::mePullXvsAlphaBarrelLayerModule
MonitorElement * mePullXvsAlphaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:117
SiPixelTrackingRecHitsValid::meNpixZpPanel1DiskPlaq
MonitorElement * meNpixZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:279
SiPixelTrackingRecHitsValid::mePosyZmPanel2DiskPlaq
MonitorElement * mePosyZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:245
SiPixelTrackingRecHitsValid::mePullyBarrelLayerModule
MonitorElement * mePullyBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:108
SiPixelTrackingRecHitsValid::mePullXvsPhiZpPanel2DiskPlaq
MonitorElement * mePullXvsPhiZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:317
SiPixelTrackingRecHitsValid::meResyBarrelLayerModule
MonitorElement * meResyBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:106
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:35
SiPixelTrackingRecHitsValid::meWPullYvsBetaBarrelNonFlippedLadders
MonitorElement * meWPullYvsBetaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:194
SiPixelTrackingRecHitsValid::meNypixBarrelLayerModule
MonitorElement * meNypixBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:111
position
static int position[TOTALCHAMBERS][3]
Definition: ReadPGInfo.cc:509
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
SiPixelTrackingRecHitsValid::recoTrackCollectionToken_
edm::EDGetTokenT< reco::TrackCollection > recoTrackCollectionToken_
Definition: SiPixelTrackingRecHitsValid.h:93
SiPixelTrackingRecHitsValid::mePullyZpPanel1DiskPlaq
MonitorElement * mePullyZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:278
SiPixelTrackingRecHitsValid::meNypixZpPanel1DiskPlaq
MonitorElement * meNypixZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:281
SiPixelTrackingRecHitsValid::mePullXvsAlphaZpPanel1DiskPlaq
MonitorElement * mePullXvsAlphaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:287
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:161
SiPixelTrackingRecHitsValid::meErryZmPanel1DiskPlaq
MonitorElement * meErryZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:221
SiPixelTrackingRecHitsValid::mePullxZmPanel2DiskPlaq
MonitorElement * mePullxZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:250
SiPixelTrackingRecHitsValid::mePullXvsEtaZpPanel1DiskPlaq
MonitorElement * mePullXvsEtaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:293
SiPixelTrackingRecHitsValid::mePullXvsPhiZmPanel2DiskPlaq
MonitorElement * mePullXvsPhiZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:264
SiPixelTrackingRecHitsValid::mePullYvsBetaZmPanel2DiskPlaq
MonitorElement * mePullYvsBetaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:263
SiPixelTrackingRecHitsValid::mePullXvsPhiBarrelLayerModule
MonitorElement * mePullXvsPhiBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:121
SiPixelTrackingRecHitsValid::mePullYvsPhiZmPanel2DiskPlaq
MonitorElement * mePullYvsPhiZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:265
SiPixelTrackingRecHitsValid::meWPullYvsAlphaBarrelNonFlippedLadders
MonitorElement * meWPullYvsAlphaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:192
SiPixelTrackingRecHitsValid::meResYvsAlphaZpPanel2DiskPlaq
MonitorElement * meResYvsAlphaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:310
LocalError::yy
float yy() const
Definition: LocalError.h:26
SiPixelTrackingRecHitsValid::meNpixBarrelLayerModule
MonitorElement * meNpixBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:109
SiPixelTrackingRecHitsValid::meResxZpPanel2DiskPlaq
MonitorElement * meResxZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:301
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
SiPixelTrackingRecHitsValid::meChargeZpPanel1DiskPlaq
MonitorElement * meChargeZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:282
SiPixelTrackingRecHitsValid::mePullXvsBetaZpPanel2DiskPlaq
MonitorElement * mePullXvsBetaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:315
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
SiPixelTrackingRecHitsValid::meResYvsAlphaZpPanel1DiskPlaq
MonitorElement * meResYvsAlphaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:284
SiPixelTrackingRecHitsValid::meResxBarrelLayerModule
MonitorElement * meResxBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:105
SiPixelTrackingRecHitsValid::meNxpixZmPanel1DiskPlaq
MonitorElement * meNxpixZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:227
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
SiPixelTrackingRecHitsValid::mePosxZpPanel2DiskPlaq
MonitorElement * mePosxZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:297
j
int j
Definition: DBlmapReader.cc:9
DQMStore::save
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", const uint32_t run=0, SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:2296
patCandidatesForDimuonsSequences_cff.tracker
tuple tracker
Definition: patCandidatesForDimuonsSequences_cff.py:69
SiPixelTrackingRecHitsValid::mePosxZpPanel1DiskPlaq
MonitorElement * mePosxZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:271
SiPixelTrackingRecHitsValid::meNxpixZpPanel2DiskPlaq
MonitorElement * meNxpixZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:306
SiPixelTrackingRecHitsValid::mePullYvsEtaZpPanel2DiskPlaq
MonitorElement * mePullYvsEtaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:320
SiPixelTrackingRecHitsValid::mePullYvsPhiZpPanel1DiskPlaq
MonitorElement * mePullYvsPhiZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:292
SiPixelTrackingRecHitsValid::meResxZmPanel2DiskPlaq
MonitorElement * meResxZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:248
SiPixelTrackingRecHitsValid::meResXvsBetaZmPanel2DiskPlaq
MonitorElement * meResXvsBetaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:258
DetId::subdetId
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:37
SiPixelTrackingRecHitsValid::mePullyZmPanel2DiskPlaq
MonitorElement * mePullyZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:251
SiPixelTrackingRecHitsValid::mePullYvsBetaBarrelLayerModule
MonitorElement * mePullYvsBetaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:120
TrackerGeometry::idToDet
virtual const GeomDet * idToDet(DetId) const
Definition: TrackerGeometry.cc:137
SiPixelTrackingRecHitsValid::meErrxZmPanel2DiskPlaq
MonitorElement * meErrxZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:246
SiPixelTrackingRecHitsValid::meNypixZmPanel1DiskPlaq
MonitorElement * meNypixZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:228
dbe_
DQMStore * dbe_
Definition: PFJetBenchmarkAnalyzer.cc:79
SiPixelTrackingRecHitsValid::mePullXvsEtaZmPanel1DiskPlaq
MonitorElement * mePullXvsEtaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:240
SiPixelTrackingRecHitsValid::meErrxZpPanel1DiskPlaq
MonitorElement * meErrxZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:273
SiPixelTrackingRecHitsValid::meErryZpPanel1DiskPlaq
MonitorElement * meErryZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:274
SiPixelTrackingRecHitsValid::meResXvsBetaZpPanel2DiskPlaq
MonitorElement * meResXvsBetaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:311
SiPixelTrackingRecHitsValid::meResXvsBetaZpPanel1DiskPlaq
MonitorElement * meResXvsBetaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:285
Local3DPoint
Point3DBase< float, LocalTag > Local3DPoint
Definition: LocalPoint.h:9
SiPixelTrackingRecHitsValid::mePullYvsPhiZmPanel1DiskPlaq
MonitorElement * mePullYvsPhiZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:239
DetId
Definition: DetId.h:18
SiPixelTrackingRecHitsValid::meChargeZmPanel1DiskPlaq
MonitorElement * meChargeZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:229
SiPixelTrackingRecHitsValid::mePullYvsAlphaZpPanel2DiskPlaq
MonitorElement * mePullYvsAlphaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:314
SiPixelTrackingRecHitsValid::meNypixZmPanel2DiskPlaq
MonitorElement * meNypixZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:254
SiPixelTrackingRecHitsValid::mePullYvsPhiZpPanel2DiskPlaq
MonitorElement * mePullYvsPhiZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:318
SiPixelTrackingRecHitsValid::mePullYvsBetaZmPanel1DiskPlaq
MonitorElement * mePullYvsBetaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:237
SiPixelTrackingRecHitsValid::mePullXvsEtaBarrelLayerModule
MonitorElement * mePullXvsEtaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:123
SiPixelTrackingRecHitsValid::mePullXvsBetaZmPanel2DiskPlaq
MonitorElement * mePullXvsBetaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:262
testEve_cfg.tracks
tuple tracks
Definition: testEve_cfg.py:39
SiPixelTrackingRecHitsValid::mePullXvsAlphaZpPanel2DiskPlaq
MonitorElement * mePullXvsAlphaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:313
SiPixelTrackingRecHitsValid::meResYvsBetaBarrelLayerModule
MonitorElement * meResYvsBetaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:116
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
SiPixelTrackingRecHitsValid::mePullXvsAlphaBarrelNonFlippedLadders
MonitorElement * mePullXvsAlphaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:178
sysUtil.pid
tuple pid
Definition: sysUtil.py:22
PixelGeomDetUnit::specificTopology
virtual const PixelTopology & specificTopology() const
Returns a reference to the pixel proxy topology.
Definition: PixelGeomDetUnit.cc:19
SiPixelTrackingRecHitsValid::SiPixelTrackingRecHitsValid
SiPixelTrackingRecHitsValid(const edm::ParameterSet &conf)
Definition: SiPixelTrackingRecHitsValid.cc:112
SiPixelTrackingRecHitsValid::mePullYvsAlphaBarrelLayerModule
MonitorElement * mePullYvsAlphaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:118
SiPixelTrackingRecHitsValid::mePullYvsEtaZmPanel1DiskPlaq
MonitorElement * mePullYvsEtaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:241
SiPixelTrackingRecHitsValid::meResXvsAlphaZmPanel1DiskPlaq
MonitorElement * meResXvsAlphaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:230
edm::Handle::product
T const * product() const
Definition: Handle.h:81
SiPixelTrackingRecHitsValid::mePullYvsAlphaZpPanel1DiskPlaq
MonitorElement * mePullYvsAlphaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:288
SiPixelTrackingRecHitsValid::meChargeBarrelLayerModule
MonitorElement * meChargeBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:112
SiPixelTrackingRecHitsValid::meResXvsAlphaZpPanel1DiskPlaq
MonitorElement * meResXvsAlphaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:283
SiPixelTrackingRecHitsValid::mePullYvsBetaZpPanel2DiskPlaq
MonitorElement * mePullYvsBetaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:316
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
SiPixelTrackingRecHitsValid::meResYvsBetaZpPanel1DiskPlaq
MonitorElement * meResYvsBetaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:286
SiPixelTrackingRecHitsValid::mePosxBarrelLayerModule
MonitorElement * mePosxBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:101
edmNew::DetSet::end
iterator end()
Definition: DetSetNew.h:70
SiPixelTrackingRecHitsValid::meResYvsAlphaZmPanel2DiskPlaq
MonitorElement * meResYvsAlphaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:257
SiPixelTrackingRecHitsValid::meResyZpPanel2DiskPlaq
MonitorElement * meResyZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:302
SiPixelTrackingRecHitsValid::mePullXvsBetaZpPanel1DiskPlaq
MonitorElement * mePullXvsBetaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:289
SiPixelTrackingRecHitsValid::meResYvsBetaZpPanel2DiskPlaq
MonitorElement * meResYvsBetaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:312
TrackerHitAssociator::associateHit
std::vector< PSimHit > associateHit(const TrackingRecHit &thit)
Definition: TrackerHitAssociator.cc:132
SiPixelTrackingRecHitsValid::mePullXvsEtaZmPanel2DiskPlaq
MonitorElement * mePullXvsEtaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:266
SiPixelTrackingRecHitsValid::meResYvsAlphaBarrelLayerModule
MonitorElement * meResYvsAlphaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:114
match
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:6
SiPixelTrackingRecHitsValid::mePosyZpPanel2DiskPlaq
MonitorElement * mePosyZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:298
SiPixelTrackingRecHitsValid::mePullXvsAlphaZmPanel2DiskPlaq
MonitorElement * mePullXvsAlphaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:260
SiPixelTrackingRecHitsValid::mePullxZpPanel1DiskPlaq
MonitorElement * mePullxZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:277
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
SiPixelTrackingRecHitsValid::meResXvsAlphaZpPanel2DiskPlaq
MonitorElement * meResXvsAlphaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:309
SiPixelTrackingRecHitsValid::meResXvsAlphaBarrelLayerModule
MonitorElement * meResXvsAlphaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:113
SiPixelTrackingRecHitsValid::meResYvsBetaZmPanel2DiskPlaq
MonitorElement * meResYvsBetaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:259
SiPixelTrackingRecHitsValid::meNpixZmPanel2DiskPlaq
MonitorElement * meNpixZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:252
SiPixelTrackingRecHitsValid::meResXvsAlphaBarrelFlippedLaddersLayer
MonitorElement * meResXvsAlphaBarrelFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:202
SiPixelTrackingRecHitsValid::meResXvsBetaBarrelLayerModule
MonitorElement * meResXvsBetaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:115
SiPixelTrackingRecHitsValid::mePullYvsPhiBarrelLayerModule
MonitorElement * mePullYvsPhiBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:122
SiPixelTrackingRecHitsValid::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: SiPixelTrackingRecHitsValid.cc:1067
SiPixelTrackingRecHitsValid::meNypixZpPanel2DiskPlaq
MonitorElement * meNypixZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:307
SiPixelTrackingRecHitsValid::mePosyZmPanel1DiskPlaq
MonitorElement * mePosyZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:219
SiPixelTrackingRecHitsValid::mePullXvsPhiZpPanel1DiskPlaq
MonitorElement * mePullXvsPhiZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:291
SiPixelTrackingRecHitsValid::meNxpixBarrelLayerModule
MonitorElement * meNxpixBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:110
SiPixelTrackingRecHitsValid::meNxpixZmPanel2DiskPlaq
MonitorElement * meNxpixZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:253
SiPixelRecHit
Pixel Reconstructed Hit.
SiPixelTrackingRecHitsValid::meResxZmPanel1DiskPlaq
MonitorElement * meResxZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:222
edmNew::DetSet::begin
iterator begin()
Definition: DetSetNew.h:67
SiPixelTrackingRecHitsValid::meResYvsBetaBarrelNonFlippedLaddersLayer
MonitorElement * meResYvsBetaBarrelNonFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:209
Generated for CMSSW Reference Manual by doxygen 1.8.5