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/TrackFwd.h"
25 #include "DataFormats/TrackReco/interface/TrackExtra.h"
26 #include "DataFormats/DetId/interface/DetId.h"
27 #include "DataFormats/GeometryVector/interface/LocalPoint.h"
28 #include "DataFormats/GeometryVector/interface/GlobalPoint.h"
29 
30 #include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHit.h"
31 #include "DataFormats/SiPixelDetId/interface/PXBDetId.h"
32 #include "DataFormats/SiPixelDetId/interface/PXFDetId.h"
33 #include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
34 
35 #include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
36 
37 #include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHitCollection.h"
38 
39 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
40 
41 #include "Geometry/TrackerGeometryBuilder/interface/PixelTopologyBuilder.h"
42 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
43 #include "Geometry/CommonDetUnit/interface/GeomDetType.h"
44 #include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"
45 
46 
47 #include <TTree.h>
48 #include <TFile.h>
49 
50 using namespace std;
51 using namespace edm;
52 
53 // End job: write and close the ntuple file
54 void SiPixelTrackingRecHitsValid::endJob()
55 {
56  if(debugNtuple_.size()!=0){
57  tfile_->Write();
58  tfile_->Close();
59  }
60 }
61 
62 void SiPixelTrackingRecHitsValid::beginJob()
63 {
64 
65  if(debugNtuple_.size()!=0){
66  tfile_ = new TFile (debugNtuple_.c_str() , "RECREATE");
67 
68  t_ = new TTree("Ntuple", "Ntuple");
69  int bufsize = 64000;
70 
71  t_->Branch("subdetId", &subdetId, "subdetId/I", bufsize);
72 
73  t_->Branch("layer" , &layer , "layer/I" , bufsize);
74  t_->Branch("ladder", &ladder, "ladder/I", bufsize);
75  t_->Branch("mod" , &mod , "mod/I" , bufsize);
76  t_->Branch("side" , &side , "side/I" , bufsize);
77  t_->Branch("disk" , &disk , "disk/I" , bufsize);
78  t_->Branch("blade" , &blade , "blade/I" , bufsize);
79  t_->Branch("panel" , &panel , "panel/I" , bufsize);
80  t_->Branch("plaq" , &plaq , "plaq/I" , bufsize);
81 
82  t_->Branch("rechitx" , &rechitx , "rechitx/F" , bufsize);
83  t_->Branch("rechity" , &rechity , "rechity/F" , bufsize);
84  t_->Branch("rechitz" , &rechitz , "rechitz/F" , bufsize);
85  t_->Branch("rechiterrx" , &rechiterrx , "rechiterrx/F" , bufsize);
86  t_->Branch("rechiterry" , &rechiterry , "rechiterry/F" , bufsize);
87  t_->Branch("rechitresx" , &rechitresx , "rechitresx/F" , bufsize);
88  t_->Branch("rechitresy" , &rechitresy , "rechitresy/F" , bufsize);
89  t_->Branch("rechitpullx", &rechitpullx, "rechitpullx/F", bufsize);
90  t_->Branch("rechitpully", &rechitpully, "rechitpully/F", bufsize);
91 
92  t_->Branch("npix" , &npix , "npix/I" , bufsize);
93  t_->Branch("nxpix" , &nxpix , "nxpix/I" , bufsize);
94  t_->Branch("nypix" , &nypix , "nypix/I" , bufsize);
95  t_->Branch("charge", &charge, "charge/F", bufsize);
96 
97  t_->Branch("alpha", &alpha, "alpha/F", bufsize);
98  t_->Branch("beta" , &beta , "beta/F" , bufsize);
99 
100  t_->Branch("phi", &phi, "phi/F", bufsize);
101  t_->Branch("eta", &eta, "eta/F", bufsize);
102 
103  t_->Branch("half" , &half , "half/I" , bufsize);
104  t_->Branch("flipped", &flipped, "flipped/I", bufsize);
105 
106  t_->Branch("simhitx", &simhitx, "simhitx/F", bufsize);
107  t_->Branch("simhity", &simhity, "simhity/F", bufsize);
108 
109  t_->Branch("nsimhit", &nsimhit, "nsimhit/I", bufsize);
110  t_->Branch("pidhit" , &pidhit , "pidhit/I" , bufsize);
111 
112  t_->Branch("evt", &evt, "evt/I", bufsize);
113  t_->Branch("run", &run, "run/I", bufsize);
114  }
115 
116 }
117 
118 SiPixelTrackingRecHitsValid::SiPixelTrackingRecHitsValid(const ParameterSet& ps):conf_(ps), dbe_(0), tfile_(0), t_(0)
119 {
120  //Read config file
121  MTCCtrack_ = ps.getParameter<bool>("MTCCtrack");
122  outputFile_ = ps.getUntrackedParameter<string>("outputFile", "pixeltrackingrechitshisto.root");
123  src_ = ps.getUntrackedParameter<std::string>( "src" );
124  builderName_ = ps.getParameter<std::string>("TTRHBuilder");
125  checkType_ = ps.getParameter<bool>("checkType");
126  genType_ = ps.getParameter<int>("genType");
127  debugNtuple_=ps.getUntrackedParameter<string>("debugNtuple", "SiPixelTrackingRecHitsValid_Ntuple.root");
128 
129  // Book histograms
130  dbe_ = Service<DQMStore>().operator->();
131  //dbe_->showDirStructure();
132 
133  //float math_pi = 3.14159265;
134  //float radtodeg = 180.0 / math_pi;
135 
136  // Histogram ranges (low and high)
137 
138  float xl = -1.0;
139  float xh = 1.0;
140  float errxl = 0.0;
141  float errxh = 0.003;
142  float resxl = -0.02;
143  float resxh = 0.02;
144  float pullxl = -10.0;
145  float pullxh = 10.0;
146 
147  float yl = -4.0;
148  float yh = 4.0;
149  float erryl = 0.0;
150  float erryh = 0.010;
151  float resyl = -0.04;
152  float resyh = 0.04;
153  float pullyl = -10.0;
154  float pullyh = 10.0;
155 
156  float barrel_alphal = 80.0;
157  float barrel_alphah = 100.0;
158  float barrel_betal = 10.0;
159  float barrel_betah = 170.0;
160  //float barrel_phil = -180.0;
161  //float barrel_phih = 180.0;
162  //float barrel_etal = -2.5;
163  //float barrel_etah = 2.5;
164 
165  float forward_p1_alphal = 100.0;
166  float forward_p1_alphah = 115.0;
167  float forward_p2_alphal = 65.0;
168  float forward_p2_alphah = 80.0;
169  float forward_neg_betal = 67.0;
170  float forward_neg_betah = 73.0;
171  float forward_pos_betal = 107.0;
172  float forward_pos_betah = 113.0;
173  //float forward_phil = -180.0;
174  //float forward_phih = 180.0;
175  //float forward_neg_etal = -2.5;
176  //float forward_neg_etah = -1.5;
177  //float forward_pos_etal = 1.5;
178  //float forward_pos_etah = 2.5;
179 
180  // special ranges for pulls
181  float pull_barrel_alphal = 80.0;
182  float pull_barrel_alphah = 100.0;
183  float pull_barrel_betal = 10.0;
184  float pull_barrel_betah = 170.0;
185  float pull_barrel_phil = -180.0;
186  float pull_barrel_phih = 180.0;
187  float pull_barrel_etal = -2.4;
188  float pull_barrel_etah = 2.4;
189 
190  float pull_forward_p1_alphal = 100.0;
191  float pull_forward_p1_alphah = 112.0;
192  float pull_forward_p2_alphal = 68.0;
193  float pull_forward_p2_alphah = 80.0;
194  float pull_forward_neg_betal = 68.0;
195  float pull_forward_neg_betah = 72.0;
196  float pull_forward_pos_betal = 108.0;
197  float pull_forward_pos_betah = 112.0;
198  float pull_forward_phil = -180.0;
199  float pull_forward_phih = 180.0;
200  float pull_forward_neg_etal = -2.4;
201  float pull_forward_neg_etah = -1.4;
202  float pull_forward_pos_etal = 1.5;
203  float pull_forward_pos_etah = 2.5;
204 
205  int npixl = 0;
206  int npixh = 20;
207  int nxpixl = 0;
208  int nxpixh = 10;
209  int nypixl = 0;
210  int nypixh = 20;
211 
212  float barrel_chargel = 0.0;
213  float barrel_chargeh = 250000.0;
214  float forward_chargel = 0.0;
215  float forward_chargeh = 100000.0;
216 
217  dbe_->setCurrentFolder("Tracking/TrackingRecHits/Pixel/Histograms_per_ring-layer_or_disk-plaquette");
218 
219  // Pixel barrel has 3 layers and 8 rings; book a histogram for each module given by the (layer, ring) pair
220  for (int i=0; i<3 ; i++) // loop over layers
221  {
222  Char_t chisto[100];
223 
224  sprintf(chisto, "meResxBarrelLayer_%d", i+1);
225  meResxBarrelLayer[i] = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
226  sprintf(chisto, "meResyBarrelLayer_%d", i+1);
227  meResyBarrelLayer[i] = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
228  sprintf(chisto, "mePullxBarrelLayer_%d", i+1);
229  mePullxBarrelLayer[i] = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
230  sprintf(chisto, "mePullyBarrelLayer_%d", i+1);
231  mePullyBarrelLayer[i] = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
232 
233  sprintf(chisto, "meResXvsAlphaBarrelFlippedLaddersLayer_%d", i+1);
234  meResXvsAlphaBarrelFlippedLaddersLayer[i] = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh, "");
235  sprintf(chisto, "meResYvsAlphaBarrelFlippedLaddersLayer_%d", i+1);
236  meResYvsAlphaBarrelFlippedLaddersLayer[i] = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh, "");
237  sprintf(chisto, "meResXvsBetaBarrelFlippedLaddersLayer_%d", i+1);
238  meResXvsBetaBarrelFlippedLaddersLayer[i] = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh, "");
239  sprintf(chisto, "meResYvsBetaBarrelFlippedLaddersLayer_%d", i+1);
240  meResYvsBetaBarrelFlippedLaddersLayer[i] = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh, "");
241 
242  sprintf(chisto, "meResXvsAlphaBarrelNonFlippedLaddersLayer_%d", i+1);
243  meResXvsAlphaBarrelNonFlippedLaddersLayer[i]
244  = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh, "");
245  sprintf(chisto, "meResYvsAlphaBarrelNonFlippedLaddersLayer_%d", i+1);
246  meResYvsAlphaBarrelNonFlippedLaddersLayer[i]
247  = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh, "");
248  sprintf(chisto, "meResXvsBetaBarrelNonFlippedLaddersLayer_%d", i+1);
249  meResXvsBetaBarrelNonFlippedLaddersLayer[i]
250  = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh, "");
251  sprintf(chisto, "meResYvsBetaBarrelNonFlippedLaddersLayer_%d", i+1);
252  meResYvsBetaBarrelNonFlippedLaddersLayer[i]
253  = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh, "");
254 
255  for (int j=0; j<8; j++) // loop over rings
256  {
257  sprintf(chisto, "mePosxBarrelLayerModule_%d_%d", i+1, j+1);
258  mePosxBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, xl, xh);
259  sprintf(chisto, "mePosyBarrelLayerModule_%d_%d", i+1, j+1);
260  mePosyBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, yl, yh);
261  sprintf(chisto, "meErrxBarrelLayerModule_%d_%d", i+1, j+1);
262  meErrxBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
263  sprintf(chisto, "meErryBarrelLayerModule_%d_%d", i+1, j+1);
264  meErryBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
265  sprintf(chisto, "meResxBarrelLayerModule_%d_%d", i+1, j+1);
266  meResxBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
267  sprintf(chisto, "meResyBarrelLayerModule_%d_%d", i+1, j+1);
268  meResyBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
269  sprintf(chisto, "mePullxBarrelLayerModule_%d_%d", i+1, j+1);
270  mePullxBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
271  sprintf(chisto, "mePullyBarrelLayerModule_%d_%d", i+1, j+1);
272  mePullyBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
273  sprintf(chisto, "meNpixBarrelLayerModule_%d_%d", i+1, j+1);
274  meNpixBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
275  sprintf(chisto, "meNxpixBarrelLayerModule_%d_%d", i+1, j+1);
276  meNxpixBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
277  sprintf(chisto, "meNypixBarrelLayerModule_%d_%d", i+1, j+1);
278  meNypixBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
279  sprintf(chisto, "meChargeBarrelLayerModule_%d_%d", i+1, j+1);
280  meChargeBarrelLayerModule[i][j] = dbe_->book1D(chisto, chisto, 100, barrel_chargel, barrel_chargeh);
281 
282  sprintf(chisto, "meResXvsAlphaBarrelLayerModule_%d_%d", i+1, j+1);
283  meResXvsAlphaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh, "");
284  sprintf(chisto, "meResYvsAlphaBarrelLayerModule_%d_%d", i+1, j+1);
285  meResYvsAlphaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh, "");
286  sprintf(chisto, "meResXvsBetaBarrelLayerModule_%d_%d", i+1, j+1);
287  meResXvsBetaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh, "");
288  sprintf(chisto, "meResYvsBetaBarrelLayerModule_%d_%d", i+1, j+1);
289  meResYvsBetaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh, "");
290 
291  sprintf(chisto, "mePullXvsAlphaBarrelLayerModule_%d_%d", i+1, j+1);
292  mePullXvsAlphaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah,
293  100, pullxl, pullxh, "");
294  sprintf(chisto, "mePullYvsAlphaBarrelLayerModule_%d_%d", i+1, j+1);
295  mePullYvsAlphaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah,
296  100, pullyl, pullyh, "");
297  sprintf(chisto, "mePullXvsBetaBarrelLayerModule_%d_%d", i+1, j+1);
298  mePullXvsBetaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah,
299  100, pullxl, pullxh, "");
300  sprintf(chisto, "mePullYvsBetaBarrelLayerModule_%d_%d", i+1, j+1);
301  mePullYvsBetaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah,
302  100, pullyl, pullyh, "");
303  sprintf(chisto, "mePullXvsPhiBarrelLayerModule_%d_%d", i+1, j+1);
304  mePullXvsPhiBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100,
305  pullxl, pullxh, "");
306  sprintf(chisto, "mePullYvsPhiBarrelLayerModule_%d_%d", i+1, j+1);
307  mePullYvsPhiBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih,
308  100, pullyl, pullyh, "");
309  sprintf(chisto, "mePullXvsEtaBarrelLayerModule_%d_%d", i+1, j+1);
310  mePullXvsEtaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah,
311  100, pullxl, pullxh, "");
312  sprintf(chisto, "mePullYvsEtaBarrelLayerModule_%d_%d", i+1, j+1);
313  mePullYvsEtaBarrelLayerModule[i][j] = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah,
314  100, pullyl, pullyh, "");
315  } // for (int j=0; j<8; j++) // loop over rings
316 
317  } // for (int i=0; i<3 ; i++) // loop over layers
318 
319  // Pixel forward detector has 2 disks, 2 panels and either 3 or 4 plaquettes
320  // Panel 1 has 4 plaquettes
321  // Panel 2 has 3 plaquettes
322 
323  // Panel 1: 2 disks, 4 plaquets
324  // Panel 2: 2 disks, 3 plaquets
325  for (int i=0; i<2 ; i++) // loop over disks
326  for (int j=0; j<4; j++) // loop over plaquettes
327  {
328  Char_t chisto[100];
329 
330  sprintf(chisto, "mePosxZmPanel1DiskPlaq_%d_%d", i+1, j+1);
331  mePosxZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, xl, xh);
332  sprintf(chisto, "mePosyZmPanel1DiskPlaq_%d_%d", i+1, j+1);
333  mePosyZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, yl, yh);
334  sprintf(chisto, "meErrxZmPanel1DiskPlaq_%d_%d", i+1, j+1);
335  meErrxZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
336  sprintf(chisto, "meErryZmPanel1DiskPlaq_%d_%d", i+1, j+1);
337  meErryZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
338  sprintf(chisto, "meResxZmPanel1DiskPlaq_%d_%d", i+1, j+1);
339  meResxZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
340  sprintf(chisto, "meResyZmPanel1DiskPlaq_%d_%d", i+1, j+1);
341  meResyZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
342  sprintf(chisto, "mePullxZmPanel1DiskPlaq_%d_%d", i+1, j+1);
343  mePullxZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
344  sprintf(chisto, "mePullyZmPanel1DiskPlaq_%d_%d", i+1, j+1);
345  mePullyZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
346  sprintf(chisto, "meNpixZmPanel1DiskPlaq_%d_%d", i+1, j+1);
347  meNpixZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
348  sprintf(chisto, "meNxpixZmPanel1DiskPlaq_%d_%d", i+1, j+1);
349  meNxpixZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
350  sprintf(chisto, "meNypixZmPanel1DiskPlaq_%d_%d", i+1, j+1);
351  meNypixZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
352  sprintf(chisto, "meChargeZmPanel1DiskPlaq_%d_%d", i+1, j+1);
353  meChargeZmPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
354 
355  sprintf(chisto, "meResXvsAlphaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
356  meResXvsAlphaZmPanel1DiskPlaq[i][j]
357  = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh, "");
358  sprintf(chisto, "meResYvsAlphaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
359  meResYvsAlphaZmPanel1DiskPlaq[i][j]
360  = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh, "");
361  sprintf(chisto, "meResXvsBetaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
362  meResXvsBetaZmPanel1DiskPlaq[i][j]
363  = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh, "");
364  sprintf(chisto, "meResYvsBetaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
365  meResYvsBetaZmPanel1DiskPlaq[i][j]
366  = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh, "");
367 
368  sprintf(chisto, "mePullXvsAlphaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
369  mePullXvsAlphaZmPanel1DiskPlaq[i][j]
370  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl,pullxh, "");
371  sprintf(chisto, "mePullYvsAlphaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
372  mePullYvsAlphaZmPanel1DiskPlaq[i][j]
373  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl,pullyh, "");
374  sprintf(chisto, "mePullXvsBetaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
375  mePullXvsBetaZmPanel1DiskPlaq[i][j]
376  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl,pullxh, "");
377  sprintf(chisto, "mePullYvsBetaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
378  mePullYvsBetaZmPanel1DiskPlaq[i][j]
379  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl,pullyh, "");
380  sprintf(chisto, "mePullXvsPhiZmPanel1DiskPlaq_%d_%d", i+1, j+1);
381  mePullXvsPhiZmPanel1DiskPlaq[i][j]
382  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl,pullxh, "");
383  sprintf(chisto, "mePullYvsPhiZmPanel1DiskPlaq_%d_%d", i+1, j+1);
384  mePullYvsPhiZmPanel1DiskPlaq[i][j]
385  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl,pullyh, "");
386  sprintf(chisto, "mePullXvsEtaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
387  mePullXvsEtaZmPanel1DiskPlaq[i][j]
388  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl,pullxh, "");
389  sprintf(chisto, "mePullYvsEtaZmPanel1DiskPlaq_%d_%d", i+1, j+1);
390  mePullYvsEtaZmPanel1DiskPlaq[i][j]
391  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl,pullyh, "");
392 
393  sprintf(chisto, "mePosxZpPanel1DiskPlaq_%d_%d", i+1, j+1);
394  mePosxZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, xl, xh);
395  sprintf(chisto, "mePosyZpPanel1DiskPlaq_%d_%d", i+1, j+1);
396  mePosyZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, yl, yh);
397  sprintf(chisto, "meErrxZpPanel1DiskPlaq_%d_%d", i+1, j+1);
398  meErrxZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
399  sprintf(chisto, "meErryZpPanel1DiskPlaq_%d_%d", i+1, j+1);
400  meErryZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
401  sprintf(chisto, "meResxZpPanel1DiskPlaq_%d_%d", i+1, j+1);
402  meResxZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
403  sprintf(chisto, "meResyZpPanel1DiskPlaq_%d_%d", i+1, j+1);
404  meResyZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
405  sprintf(chisto, "mePullxZpPanel1DiskPlaq_%d_%d", i+1, j+1);
406  mePullxZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
407  sprintf(chisto, "mePullyZpPanel1DiskPlaq_%d_%d", i+1, j+1);
408  mePullyZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
409  sprintf(chisto, "meNpixZpPanel1DiskPlaq_%d_%d", i+1, j+1);
410  meNpixZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
411  sprintf(chisto, "meNxpixZpPanel1DiskPlaq_%d_%d", i+1, j+1);
412  meNxpixZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
413  sprintf(chisto, "meNypixZpPanel1DiskPlaq_%d_%d", i+1, j+1);
414  meNypixZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
415  sprintf(chisto, "meChargeZpPanel1DiskPlaq_%d_%d", i+1, j+1);
416  meChargeZpPanel1DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
417  sprintf(chisto, "meResXvsAlphaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
418  meResXvsAlphaZpPanel1DiskPlaq[i][j]
419  = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh, "");
420  sprintf(chisto, "meResYvsAlphaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
421  meResYvsAlphaZpPanel1DiskPlaq[i][j]
422  = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh, "");
423  sprintf(chisto, "meResXvsBetaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
424  meResXvsBetaZpPanel1DiskPlaq[i][j]
425  = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh, "");
426  sprintf(chisto, "meResYvsBetaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
427  meResYvsBetaZpPanel1DiskPlaq[i][j]
428  = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh, "");
429 
430  sprintf(chisto, "mePullXvsAlphaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
431  mePullXvsAlphaZpPanel1DiskPlaq[i][j]
432  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl,pullxh, "");
433  sprintf(chisto, "mePullYvsAlphaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
434  mePullYvsAlphaZpPanel1DiskPlaq[i][j]
435  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl,pullyh, "");
436  sprintf(chisto, "mePullXvsBetaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
437  mePullXvsBetaZpPanel1DiskPlaq[i][j]
438  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl,pullxh, "");
439  sprintf(chisto, "mePullYvsBetaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
440  mePullYvsBetaZpPanel1DiskPlaq[i][j]
441  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl,pullyh, "");
442  sprintf(chisto, "mePullXvsPhiZpPanel1DiskPlaq_%d_%d", i+1, j+1);
443  mePullXvsPhiZpPanel1DiskPlaq[i][j]
444  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl,pullxh, "");
445  sprintf(chisto, "mePullYvsPhiZpPanel1DiskPlaq_%d_%d", i+1, j+1);
446  mePullYvsPhiZpPanel1DiskPlaq[i][j]
447  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl,pullyh, "");
448  sprintf(chisto, "mePullXvsEtaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
449  mePullXvsEtaZpPanel1DiskPlaq[i][j]
450  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl,pullxh, "");
451  sprintf(chisto, "mePullYvsEtaZpPanel1DiskPlaq_%d_%d", i+1, j+1);
452  mePullYvsEtaZpPanel1DiskPlaq[i][j]
453  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl,pullyh, "");
454 
455  if ( j>2 ) continue; // panel 2 has only 3 plaquettes
456 
457  sprintf(chisto, "mePosxZmPanel2DiskPlaq_%d_%d", i+1, j+1);
458  mePosxZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, xl, xh);
459  sprintf(chisto, "mePosyZmPanel2DiskPlaq_%d_%d", i+1, j+1);
460  mePosyZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, yl, yh);
461  sprintf(chisto, "meErrxZmPanel2DiskPlaq_%d_%d", i+1, j+1);
462  meErrxZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
463  sprintf(chisto, "meErryZmPanel2DiskPlaq_%d_%d", i+1, j+1);
464  meErryZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
465  sprintf(chisto, "meResxZmPanel2DiskPlaq_%d_%d", i+1, j+1);
466  meResxZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
467  sprintf(chisto, "meResyZmPanel2DiskPlaq_%d_%d", i+1, j+1);
468  meResyZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
469  sprintf(chisto, "mePullxZmPanel2DiskPlaq_%d_%d", i+1, j+1);
470  mePullxZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
471  sprintf(chisto, "mePullyZmPanel2DiskPlaq_%d_%d", i+1, j+1);
472  mePullyZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
473  sprintf(chisto, "meNpixZmPanel2DiskPlaq_%d_%d", i+1, j+1);
474  meNpixZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
475  sprintf(chisto, "meNxpixZmPanel2DiskPlaq_%d_%d", i+1, j+1);
476  meNxpixZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
477  sprintf(chisto, "meNypixZmPanel2DiskPlaq_%d_%d", i+1, j+1);
478  meNypixZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
479  sprintf(chisto, "meChargeZmPanel2DiskPlaq_%d_%d", i+1, j+1);
480  meChargeZmPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
481  sprintf(chisto, "meResXvsAlphaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
482  meResXvsAlphaZmPanel2DiskPlaq[i][j]
483  = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh, "");
484  sprintf(chisto, "meResYvsAlphaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
485  meResYvsAlphaZmPanel2DiskPlaq[i][j]
486  = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh, "");
487  sprintf(chisto, "meResXvsBetaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
488  meResXvsBetaZmPanel2DiskPlaq[i][j]
489  = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh, "");
490  sprintf(chisto, "meResYvsBetaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
491  meResYvsBetaZmPanel2DiskPlaq[i][j]
492  = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh, "");
493 
494  sprintf(chisto, "mePullXvsAlphaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
495  mePullXvsAlphaZmPanel2DiskPlaq[i][j]
496  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl,pullxh, "");
497  sprintf(chisto, "mePullYvsAlphaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
498  mePullYvsAlphaZmPanel2DiskPlaq[i][j]
499  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl,pullyh, "");
500  sprintf(chisto, "mePullXvsBetaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
501  mePullXvsBetaZmPanel2DiskPlaq[i][j]
502  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl,pullxh, "");
503  sprintf(chisto, "mePullYvsBetaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
504  mePullYvsBetaZmPanel2DiskPlaq[i][j]
505  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl,pullyh, "");
506  sprintf(chisto, "mePullXvsPhiZmPanel2DiskPlaq_%d_%d", i+1, j+1);
507  mePullXvsPhiZmPanel2DiskPlaq[i][j]
508  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl,pullxh, "");
509  sprintf(chisto, "mePullYvsPhiZmPanel2DiskPlaq_%d_%d", i+1, j+1);
510  mePullYvsPhiZmPanel2DiskPlaq[i][j]
511  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl,pullyh, "");
512  sprintf(chisto, "mePullXvsEtaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
513  mePullXvsEtaZmPanel2DiskPlaq[i][j]
514  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh, "");
515  sprintf(chisto, "mePullYvsEtaZmPanel2DiskPlaq_%d_%d", i+1, j+1);
516  mePullYvsEtaZmPanel2DiskPlaq[i][j]
517  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh, "");
518 
519  sprintf(chisto, "mePosxZpPanel2DiskPlaq_%d_%d", i+1, j+1);
520  mePosxZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, xl, xh);
521  sprintf(chisto, "mePosyZpPanel2DiskPlaq_%d_%d", i+1, j+1);
522  mePosyZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, yl, yh);
523  sprintf(chisto, "meErrxZpPanel2DiskPlaq_%d_%d", i+1, j+1);
524  meErrxZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
525  sprintf(chisto, "meErryZpPanel2DiskPlaq_%d_%d", i+1, j+1);
526  meErryZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
527  sprintf(chisto, "meResxZpPanel2DiskPlaq_%d_%d", i+1, j+1);
528  meResxZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
529  sprintf(chisto, "meResyZpPanel2DiskPlaq_%d_%d", i+1, j+1);
530  meResyZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
531  sprintf(chisto, "mePullxZpPanel2DiskPlaq_%d_%d", i+1, j+1);
532  mePullxZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
533  sprintf(chisto, "mePullyZpPanel2DiskPlaq_%d_%d", i+1, j+1);
534  mePullyZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
535  sprintf(chisto, "meNpixZpPanel2DiskPlaq_%d_%d", i+1, j+1);
536  meNpixZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
537  sprintf(chisto, "meNxpixZpPanel2DiskPlaq_%d_%d", i+1, j+1);
538  meNxpixZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
539  sprintf(chisto, "meNypixZpPanel2DiskPlaq_%d_%d", i+1, j+1);
540  meNypixZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
541  sprintf(chisto, "meChargeZpPanel2DiskPlaq_%d_%d", i+1, j+1);
542  meChargeZpPanel2DiskPlaq[i][j] = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
543  sprintf(chisto, "meResXvsAlphaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
544  meResXvsAlphaZpPanel2DiskPlaq[i][j]
545  = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh, "");
546  sprintf(chisto, "meResYvsAlphaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
547  meResYvsAlphaZpPanel2DiskPlaq[i][j]
548  = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh, "");
549  sprintf(chisto, "meResXvsBetaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
550  meResXvsBetaZpPanel2DiskPlaq[i][j]
551  = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh, "");
552  sprintf(chisto, "meResYvsBetaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
553  meResYvsBetaZpPanel2DiskPlaq[i][j]
554  = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh, "");
555 
556  sprintf(chisto, "mePullXvsAlphaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
557  mePullXvsAlphaZpPanel2DiskPlaq[i][j]
558  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl,pullxh, "");
559  sprintf(chisto, "mePullYvsAlphaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
560  mePullYvsAlphaZpPanel2DiskPlaq[i][j]
561  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl,pullyh, "");
562  sprintf(chisto, "mePullXvsBetaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
563  mePullXvsBetaZpPanel2DiskPlaq[i][j]
564  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl,pullxh, "");
565  sprintf(chisto, "mePullYvsBetaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
566  mePullYvsBetaZpPanel2DiskPlaq[i][j]
567  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl,pullyh, "");
568  sprintf(chisto, "mePullXvsPhiZpPanel2DiskPlaq_%d_%d", i+1, j+1);
569  mePullXvsPhiZpPanel2DiskPlaq[i][j]
570  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl,pullxh, "");
571  sprintf(chisto, "mePullYvsPhiZpPanel2DiskPlaq_%d_%d", i+1, j+1);
572  mePullYvsPhiZpPanel2DiskPlaq[i][j]
573  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl,pullyh, "");
574  sprintf(chisto, "mePullXvsEtaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
575  mePullXvsEtaZpPanel2DiskPlaq[i][j]
576  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl,pullxh, "");
577  sprintf(chisto, "mePullYvsEtaZpPanel2DiskPlaq_%d_%d", i+1, j+1);
578  mePullYvsEtaZpPanel2DiskPlaq[i][j]
579  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl,pullyh, "");
580 
581  } // for (int j=0; j<4; j++) // loop over plaquettes
582 
583  dbe_->setCurrentFolder("Tracking/TrackingRecHits/Pixel/Histograms_all");
584 
585  Char_t chisto[100];
586  sprintf(chisto, "mePosxBarrel");
587  mePosxBarrel = dbe_->book1D(chisto, chisto, 100, xl, xh);
588  sprintf(chisto, "mePosyBarrel");
589  mePosyBarrel = dbe_->book1D(chisto, chisto, 100, yl, yh);
590  sprintf(chisto, "meErrxBarrel");
591  meErrxBarrel = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
592  sprintf(chisto, "meErryBarrel");
593  meErryBarrel = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
594  sprintf(chisto, "meResxBarrel");
595  meResxBarrel = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
596  sprintf(chisto, "meResyBarrel");
597  meResyBarrel = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
598  sprintf(chisto, "mePullxBarrel");
599  mePullxBarrel = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
600  sprintf(chisto, "mePullyBarrel");
601  mePullyBarrel = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
602  sprintf(chisto, "meNpixBarrel");
603  meNpixBarrel = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
604  sprintf(chisto, "meNxpixBarrel");
605  meNxpixBarrel = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
606  sprintf(chisto, "meNypixBarrel");
607  meNypixBarrel = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
608  sprintf(chisto, "meChargeBarrel");
609  meChargeBarrel = dbe_->book1D(chisto, chisto, 100, barrel_chargel, barrel_chargeh);
610  sprintf(chisto, "meResXvsAlphaBarrel");
611  meResXvsAlphaBarrel = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh, "");
612  sprintf(chisto, "meResYvsAlphaBarrel");
613  meResYvsAlphaBarrel = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh, "");
614  sprintf(chisto, "meResXvsBetaBarrel");
615  meResXvsBetaBarrel = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh, "");
616  sprintf(chisto, "meResYvsBetaBarrel");
617  meResYvsBetaBarrel = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh, "");
618 
619  sprintf(chisto, "mePullXvsAlphaBarrel");
620  mePullXvsAlphaBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh, "");
621  sprintf(chisto, "mePullYvsAlphaBarrel");
622  mePullYvsAlphaBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh, "");
623  sprintf(chisto, "mePullXvsBetaBarrel");
624  mePullXvsBetaBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh, "");
625  sprintf(chisto, "mePullYvsBetaBarrel");
626  mePullYvsBetaBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh, "");
627  sprintf(chisto, "mePullXvsPhiBarrel");
628  mePullXvsPhiBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh, "");
629  sprintf(chisto, "mePullYvsPhiBarrel");
630  mePullYvsPhiBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh, "");
631  sprintf(chisto, "mePullXvsEtaBarrel");
632  mePullXvsEtaBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh, "");
633  sprintf(chisto, "mePullYvsEtaBarrel");
634  mePullYvsEtaBarrel = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh, "");
635 
636  sprintf(chisto, "mePosxBarrelHalfModule");
637  mePosxBarrelHalfModule = dbe_->book1D(chisto, chisto, 100, xl, xh);
638  sprintf(chisto, "mePosxBarrelFullModule");
639  mePosxBarrelFullModule = dbe_->book1D(chisto, chisto, 100, xl, xh);
640  sprintf(chisto, "mePosxBarrelFlippedLadders");
641  mePosxBarrelFlippedLadders = dbe_->book1D(chisto, chisto, 100, xl, xh);
642  sprintf(chisto, "mePosxBarrelNonFlippedLadders");
643  mePosxBarrelNonFlippedLadders = dbe_->book1D(chisto, chisto, 100, xl, xh);
644  sprintf(chisto, "mePosyBarrelHalfModule");
645  mePosyBarrelHalfModule = dbe_->book1D(chisto, chisto, 100, yl, yh);
646  sprintf(chisto, "mePosyBarrelFullModule");
647  mePosyBarrelFullModule = dbe_->book1D(chisto, chisto, 100, yl, yh);
648  sprintf(chisto, "mePosyBarrelFlippedLadders");
649  mePosyBarrelFlippedLadders = dbe_->book1D(chisto, chisto, 100, yl, yh);
650  sprintf(chisto, "mePosyBarrelNonFlippedLadders");
651  mePosyBarrelNonFlippedLadders = dbe_->book1D(chisto, chisto, 100, yl, yh);
652 
653  sprintf(chisto, "meResXvsAlphaBarrelFlippedLadders");
654  meResXvsAlphaBarrelFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh, "");
655  sprintf(chisto, "meResYvsAlphaBarrelFlippedLadders");
656  meResYvsAlphaBarrelFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh, "");
657  sprintf(chisto, "meResXvsBetaBarrelFlippedLadders");
658  meResXvsBetaBarrelFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh, "");
659  sprintf(chisto, "meResYvsBetaBarrelFlippedLadders");
660  meResYvsBetaBarrelFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh, "");
661 
662  sprintf(chisto, "mePullXvsAlphaBarrelFlippedLadders");
663  mePullXvsAlphaBarrelFlippedLadders
664  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh, "");
665  sprintf(chisto, "mePullYvsAlphaBarrelFlippedLadders");
666  mePullYvsAlphaBarrelFlippedLadders
667  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh, "");
668  sprintf(chisto, "mePullXvsBetaBarrelFlippedLadders");
669  mePullXvsBetaBarrelFlippedLadders
670  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh, "");
671  sprintf(chisto, "mePullYvsBetaBarrelFlippedLadders");
672  mePullYvsBetaBarrelFlippedLadders
673  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh, "");
674  sprintf(chisto, "mePullXvsPhiBarrelFlippedLadders");
675  mePullXvsPhiBarrelFlippedLadders
676  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh, "");
677  sprintf(chisto, "mePullYvsPhiBarrelFlippedLadders");
678  mePullYvsPhiBarrelFlippedLadders
679  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh, "");
680  sprintf(chisto, "mePullXvsEtaBarrelFlippedLadders");
681  mePullXvsEtaBarrelFlippedLadders
682  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh, "");
683  sprintf(chisto, "mePullYvsEtaBarrelFlippedLadders");
684  mePullYvsEtaBarrelFlippedLadders
685  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh, "");
686 
687 
688  sprintf(chisto, "meWPullXvsAlphaBarrelFlippedLadders");
689  meWPullXvsAlphaBarrelFlippedLadders
690  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh, "");
691  sprintf(chisto, "meWPullYvsAlphaBarrelFlippedLadders");
692  meWPullYvsAlphaBarrelFlippedLadders
693  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh, "");
694  sprintf(chisto, "meWPullXvsBetaBarrelFlippedLadders");
695  meWPullXvsBetaBarrelFlippedLadders
696  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh, "");
697  sprintf(chisto, "meWPullYvsBetaBarrelFlippedLadders");
698  meWPullYvsBetaBarrelFlippedLadders
699  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh, "");
700 
701  sprintf(chisto, "meResXvsAlphaBarrelNonFlippedLadders");
702  meResXvsAlphaBarrelNonFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resxh, "");
703  sprintf(chisto, "meResYvsAlphaBarrelNonFlippedLadders");
704  meResYvsAlphaBarrelNonFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_alphal, barrel_alphah, 100, 0.0, resyh, "");
705  sprintf(chisto, "meResXvsBetaBarrelNonFlippedLadders");
706  meResXvsBetaBarrelNonFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resxh, "");
707  sprintf(chisto, "meResYvsBetaBarrelNonFlippedLadders");
708  meResYvsBetaBarrelNonFlippedLadders = dbe_->bookProfile(chisto, chisto, 20, barrel_betal, barrel_betah, 100, 0.0, resyh, "");
709 
710  sprintf(chisto, "mePullXvsAlphaBarrelNonFlippedLadders");
711  mePullXvsAlphaBarrelNonFlippedLadders
712  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh, "");
713  sprintf(chisto, "mePullYvsAlphaBarrelNonFlippedLadders");
714  mePullYvsAlphaBarrelNonFlippedLadders
715  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh, "");
716  sprintf(chisto, "mePullXvsBetaBarrelNonFlippedLadders");
717  mePullXvsBetaBarrelNonFlippedLadders
718  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh, "");
719  sprintf(chisto, "mePullYvsBetaBarrelNonFlippedLadders");
720  mePullYvsBetaBarrelNonFlippedLadders
721  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh, "");
722  sprintf(chisto, "mePullXvsPhiBarrelNonFlippedLadders");
723  mePullXvsPhiBarrelNonFlippedLadders
724  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullxl, pullxh, "");
725  sprintf(chisto, "mePullYvsPhiBarrelNonFlippedLadders");
726  mePullYvsPhiBarrelNonFlippedLadders
727  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_phil, pull_barrel_phih, 100, pullyl, pullyh, "");
728  sprintf(chisto, "mePullXvsEtaBarrelNonFlippedLadders");
729  mePullXvsEtaBarrelNonFlippedLadders
730  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullxl, pullxh, "");
731  sprintf(chisto, "mePullYvsEtaBarrelNonFlippedLadders");
732  mePullYvsEtaBarrelNonFlippedLadders
733  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_etal, pull_barrel_etah, 100, pullyl, pullyh, "");
734 
735 
736  sprintf(chisto, "meWPullXvsAlphaBarrelNonFlippedLadders");
737  meWPullXvsAlphaBarrelNonFlippedLadders
738  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullxl, pullxh, "");
739  sprintf(chisto, "meWPullYvsAlphaBarrelNonFlippedLadders");
740  meWPullYvsAlphaBarrelNonFlippedLadders
741  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_alphal, pull_barrel_alphah, 100, pullyl, pullyh, "");
742  sprintf(chisto, "meWPullXvsBetaBarrelNonFlippedLadders");
743  meWPullXvsBetaBarrelNonFlippedLadders
744  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullxl, pullxh, "");
745  sprintf(chisto, "meWPullYvsBetaBarrelNonFlippedLadders");
746  meWPullYvsBetaBarrelNonFlippedLadders
747  = dbe_->bookProfile(chisto, chisto, 20, pull_barrel_betal, pull_barrel_betah, 100, pullyl, pullyh, "");
748 
749 
750  sprintf(chisto, "mePosxZmPanel1");
751  mePosxZmPanel1 = dbe_->book1D(chisto, chisto, 100, xl, xh);
752  sprintf(chisto, "mePosyZmPanel1");
753  mePosyZmPanel1 = dbe_->book1D(chisto, chisto, 100, yl, yh);
754  sprintf(chisto, "meErrxZmPanel1");
755  meErrxZmPanel1 = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
756  sprintf(chisto, "meErryZmPanel1");
757  meErryZmPanel1 = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
758  sprintf(chisto, "meResxZmPanel1");
759  meResxZmPanel1 = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
760  sprintf(chisto, "meResyZmPanel1");
761  meResyZmPanel1 = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
762  sprintf(chisto, "mePullxZmPanel1");
763  mePullxZmPanel1 = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
764  sprintf(chisto, "mePullyZmPanel1");
765  mePullyZmPanel1 = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
766  sprintf(chisto, "meNpixZmPanel1");
767  meNpixZmPanel1 = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
768  sprintf(chisto, "meNxpixZmPanel1");
769  meNxpixZmPanel1 = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
770  sprintf(chisto, "meNypixZmPanel1");
771  meNypixZmPanel1 = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
772  sprintf(chisto, "meChargeZmPanel1");
773  meChargeZmPanel1 = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
774  sprintf(chisto, "meResXvsAlphaZmPanel1");
775  meResXvsAlphaZmPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh, "");
776  sprintf(chisto, "meResYvsAlphaZmPanel1");
777  meResYvsAlphaZmPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh, "");
778  sprintf(chisto, "meResXvsBetaZmPanel1");
779  meResXvsBetaZmPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh, "");
780  sprintf(chisto, "meResYvsBetaZmPanel1");
781  meResYvsBetaZmPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh, "");
782 
783  sprintf(chisto, "mePullXvsAlphaZmPanel1");
784  mePullXvsAlphaZmPanel1
785  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh, "");
786  sprintf(chisto, "mePullYvsAlphaZmPanel1");
787  mePullYvsAlphaZmPanel1
788  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh, "");
789  sprintf(chisto, "mePullXvsBetaZmPanel1");
790  mePullXvsBetaZmPanel1
791  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh, "");
792  sprintf(chisto, "mePullYvsBetaZmPanel1");
793  mePullYvsBetaZmPanel1
794  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh, "");
795  sprintf(chisto, "mePullXvsPhiZmPanel1");
796  mePullXvsPhiZmPanel1
797  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh, "");
798  sprintf(chisto, "mePullYvsPhiZmPanel1");
799  mePullYvsPhiZmPanel1
800  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh, "");
801  sprintf(chisto, "mePullXvsEtaZmPanel1");
802  mePullXvsEtaZmPanel1
803  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh, "");
804  sprintf(chisto, "mePullYvsEtaZmPanel1");
805  mePullYvsEtaZmPanel1
806  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh, "");
807 
808  sprintf(chisto, "meWPullXvsAlphaZmPanel1");
809  meWPullXvsAlphaZmPanel1
810  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh, "");
811  sprintf(chisto, "meWPullYvsAlphaZmPanel1");
812  meWPullYvsAlphaZmPanel1
813  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh, "");
814  sprintf(chisto, "meWPullXvsBetaZmPanel1");
815  meWPullXvsBetaZmPanel1
816  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh, "");
817  sprintf(chisto, "meWPullYvsBetaZmPanel1");
818  meWPullYvsBetaZmPanel1
819  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh, "");
820 
821  sprintf(chisto, "mePosxZpPanel1");
822  mePosxZpPanel1 = dbe_->book1D(chisto, chisto, 100, xl, xh);
823  sprintf(chisto, "mePosyZpPanel1");
824  mePosyZpPanel1 = dbe_->book1D(chisto, chisto, 100, yl, yh);
825  sprintf(chisto, "meErrxZpPanel1");
826  meErrxZpPanel1 = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
827  sprintf(chisto, "meErryZpPanel1");
828  meErryZpPanel1 = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
829  sprintf(chisto, "meResxZpPanel1");
830  meResxZpPanel1 = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
831  sprintf(chisto, "meResyZpPanel1");
832  meResyZpPanel1 = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
833  sprintf(chisto, "mePullxZpPanel1");
834  mePullxZpPanel1 = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
835  sprintf(chisto, "mePullyZpPanel1");
836  mePullyZpPanel1 = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
837  sprintf(chisto, "meNpixZpPanel1");
838  meNpixZpPanel1 = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
839  sprintf(chisto, "meNxpixZpPanel1");
840  meNxpixZpPanel1 = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
841  sprintf(chisto, "meNypixZpPanel1");
842  meNypixZpPanel1 = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
843  sprintf(chisto, "meChargeZpPanel1");
844  meChargeZpPanel1 = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
845  sprintf(chisto, "meResXvsAlphaZpPanel1");
846  meResXvsAlphaZpPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resxh, "");
847  sprintf(chisto, "meResYvsAlphaZpPanel1");
848  meResYvsAlphaZpPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_p1_alphal, forward_p1_alphah, 100, 0.0, resyh, "");
849  sprintf(chisto, "meResXvsBetaZpPanel1");
850  meResXvsBetaZpPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh, "");
851  sprintf(chisto, "meResYvsBetaZpPanel1");
852  meResYvsBetaZpPanel1 = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh, "");
853 
854  sprintf(chisto, "mePullXvsAlphaZpPanel1");
855  mePullXvsAlphaZpPanel1
856  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh, "");
857  sprintf(chisto, "mePullYvsAlphaZpPanel1");
858  mePullYvsAlphaZpPanel1
859  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh, "");
860  sprintf(chisto, "mePullXvsBetaZpPanel1");
861  mePullXvsBetaZpPanel1
862  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh, "");
863  sprintf(chisto, "mePullYvsBetaZpPanel1");
864  mePullYvsBetaZpPanel1
865  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh, "");
866  sprintf(chisto, "mePullXvsPhiZpPanel1");
867  mePullXvsPhiZpPanel1
868  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh, "");
869  sprintf(chisto, "mePullYvsPhiZpPanel1");
870  mePullYvsPhiZpPanel1
871  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh, "");
872  sprintf(chisto, "mePullXvsEtaZpPanel1");
873  mePullXvsEtaZpPanel1
874  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl, pullxh, "");
875  sprintf(chisto, "mePullYvsEtaZpPanel1");
876  mePullYvsEtaZpPanel1
877  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl, pullyh, "");
878 
879  sprintf(chisto, "meWPullXvsAlphaZpPanel1");
880  meWPullXvsAlphaZpPanel1
881  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullxl, pullxh, "");
882  sprintf(chisto, "meWPullYvsAlphaZpPanel1");
883  meWPullYvsAlphaZpPanel1
884  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p1_alphal, pull_forward_p1_alphah, 100, pullyl, pullyh, "");
885  sprintf(chisto, "meWPullXvsBetaZpPanel1");
886  meWPullXvsBetaZpPanel1
887  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh, "");
888  sprintf(chisto, "meWPullYvsBetaZpPanel1");
889  meWPullYvsBetaZpPanel1
890  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh, "");
891 
892  sprintf(chisto, "mePosxZmPanel2");
893  mePosxZmPanel2 = dbe_->book1D(chisto, chisto, 100, xl, xh);
894  sprintf(chisto, "mePosyZmPanel2");
895  mePosyZmPanel2 = dbe_->book1D(chisto, chisto, 100, yl, yh);
896  sprintf(chisto, "meErrxZmPanel2");
897  meErrxZmPanel2 = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
898  sprintf(chisto, "meErryZmPanel2");
899  meErryZmPanel2 = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
900  sprintf(chisto, "meResxZmPanel2");
901  meResxZmPanel2 = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
902  sprintf(chisto, "meResyZmPanel2");
903  meResyZmPanel2 = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
904  sprintf(chisto, "mePullxZmPanel2");
905  mePullxZmPanel2 = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
906  sprintf(chisto, "mePullyZmPanel2");
907  mePullyZmPanel2 = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
908  sprintf(chisto, "meNpixZmPanel2");
909  meNpixZmPanel2 = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
910  sprintf(chisto, "meNxpixZmPanel2");
911  meNxpixZmPanel2 = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
912  sprintf(chisto, "meNypixZmPanel2");
913  meNypixZmPanel2 = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
914  sprintf(chisto, "meChargeZmPanel2");
915  meChargeZmPanel2 = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
916  sprintf(chisto, "meResXvsAlphaZmPanel2");
917  meResXvsAlphaZmPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh, "");
918  sprintf(chisto, "meResYvsAlphaZmPanel2");
919  meResYvsAlphaZmPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh, "");
920  sprintf(chisto, "meResXvsBetaZmPanel2");
921  meResXvsBetaZmPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resxh, "");
922  sprintf(chisto, "meResYvsBetaZmPanel2");
923  meResYvsBetaZmPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_neg_betal, forward_neg_betah, 100, 0.0, resyh, "");
924 
925  sprintf(chisto, "mePullXvsAlphaZmPanel2");
926  mePullXvsAlphaZmPanel2
927  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh, "");
928  sprintf(chisto, "mePullYvsAlphaZmPanel2");
929  mePullYvsAlphaZmPanel2
930  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh, "");
931  sprintf(chisto, "mePullXvsBetaZmPanel2");
932  mePullXvsBetaZmPanel2
933  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh, "");
934  sprintf(chisto, "mePullYvsBetaZmPanel2");
935  mePullYvsBetaZmPanel2
936  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh, "");
937  sprintf(chisto, "mePullXvsPhiZmPanel2");
938  mePullXvsPhiZmPanel2
939  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh, "");
940  sprintf(chisto, "mePullYvsPhiZmPanel2");
941  mePullYvsPhiZmPanel2
942  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh, "");
943  sprintf(chisto, "mePullXvsEtaZmPanel2");
944  mePullXvsEtaZmPanel2
945  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullxl, pullxh, "");
946  sprintf(chisto, "mePullYvsEtaZmPanel2");
947  mePullYvsEtaZmPanel2
948  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_etal, pull_forward_neg_etah, 100, pullyl, pullyh, "");
949 
950  sprintf(chisto, "meWPullXvsAlphaZmPanel2");
951  meWPullXvsAlphaZmPanel2
952  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh, "");
953  sprintf(chisto, "meWPullYvsAlphaZmPanel2");
954  meWPullYvsAlphaZmPanel2
955  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh, "");
956  sprintf(chisto, "meWPullXvsBetaZmPanel2");
957  meWPullXvsBetaZmPanel2
958  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullxl, pullxh, "");
959  sprintf(chisto, "meWPullYvsBetaZmPanel2");
960  meWPullYvsBetaZmPanel2
961  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_neg_betal, pull_forward_neg_betah, 100, pullyl, pullyh, "");
962 
963 
964  sprintf(chisto, "mePosxZpPanel2");
965  mePosxZpPanel2 = dbe_->book1D(chisto, chisto, 100, xl, xh);
966  sprintf(chisto, "mePosyZpPanel2");
967  mePosyZpPanel2 = dbe_->book1D(chisto, chisto, 100, yl, yh);
968  sprintf(chisto, "meErrxZpPanel2");
969  meErrxZpPanel2 = dbe_->book1D(chisto, chisto, 100, errxl, errxh);
970  sprintf(chisto, "meErryZpPanel2");
971  meErryZpPanel2 = dbe_->book1D(chisto, chisto, 100, erryl, erryh);
972  sprintf(chisto, "meResxZpPanel2");
973  meResxZpPanel2 = dbe_->book1D(chisto, chisto, 100, resxl, resxh);
974  sprintf(chisto, "meResyZpPanel2");
975  meResyZpPanel2 = dbe_->book1D(chisto, chisto, 100, resyl, resyh);
976  sprintf(chisto, "mePullxZpPanel2");
977  mePullxZpPanel2 = dbe_->book1D(chisto, chisto, 100, pullxl, pullxh);
978  sprintf(chisto, "mePullyZpPanel2");
979  mePullyZpPanel2 = dbe_->book1D(chisto, chisto, 100, pullyl, pullyh);
980  sprintf(chisto, "meNpixZpPanel2");
981  meNpixZpPanel2 = dbe_->book1D(chisto, chisto, 100, npixl, npixh);
982  sprintf(chisto, "meNxpixZpPanel2");
983  meNxpixZpPanel2 = dbe_->book1D(chisto, chisto, 100, nxpixl, nxpixh);
984  sprintf(chisto, "meNypixZpPanel2");
985  meNypixZpPanel2 = dbe_->book1D(chisto, chisto, 100, nypixl, nypixh);
986  sprintf(chisto, "meChargeZpPanel2");
987  meChargeZpPanel2 = dbe_->book1D(chisto, chisto, 100, forward_chargel, forward_chargeh);
988  sprintf(chisto, "meResXvsAlphaZpPanel2");
989  meResXvsAlphaZpPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resxh, "");
990  sprintf(chisto, "meResYvsAlphaZpPanel2");
991  meResYvsAlphaZpPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_p2_alphal, forward_p2_alphah, 100, 0.0, resyh, "");
992  sprintf(chisto, "meResXvsBetaZpPanel2");
993  meResXvsBetaZpPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resxh, "");
994  sprintf(chisto, "meResYvsBetaZpPanel2");
995  meResYvsBetaZpPanel2 = dbe_->bookProfile(chisto, chisto, 20, forward_pos_betal, forward_pos_betah, 100, 0.0, resyh, "");
996 
997  sprintf(chisto, "mePullXvsAlphaZpPanel2");
998  mePullXvsAlphaZpPanel2
999  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh, "");
1000  sprintf(chisto, "mePullYvsAlphaZpPanel2");
1001  mePullYvsAlphaZpPanel2
1002  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh, "");
1003  sprintf(chisto, "mePullXvsBetaZpPanel2");
1004  mePullXvsBetaZpPanel2
1005  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh, "");
1006  sprintf(chisto, "mePullYvsBetaZpPanel2");
1007  mePullYvsBetaZpPanel2
1008  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh, "");
1009  sprintf(chisto, "mePullXvsPhiZpPanel2");
1010  mePullXvsPhiZpPanel2
1011  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullxl, pullxh, "");
1012  sprintf(chisto, "mePullYvsPhiZpPanel2");
1013  mePullYvsPhiZpPanel2
1014  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_phil, pull_forward_phih, 100, pullyl, pullyh, "");
1015  sprintf(chisto, "mePullXvsEtaZpPanel2");
1016  mePullXvsEtaZpPanel2
1017  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullxl, pullxh, "");
1018  sprintf(chisto, "mePullYvsEtaZpPanel2");
1019  mePullYvsEtaZpPanel2
1020  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_etal, pull_forward_pos_etah, 100, pullyl, pullyh, "");
1021 
1022  sprintf(chisto, "meWPullXvsAlphaZpPanel2");
1023  meWPullXvsAlphaZpPanel2
1024  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullxl, pullxh, "");
1025  sprintf(chisto, "meWPullYvsAlphaZpPanel2");
1026  meWPullYvsAlphaZpPanel2
1027  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_p2_alphal, pull_forward_p2_alphah, 100, pullyl, pullyh, "");
1028  sprintf(chisto, "meWPullXvsBetaZpPanel2");
1029  meWPullXvsBetaZpPanel2
1030  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullxl, pullxh, "");
1031  sprintf(chisto, "meWPullYvsBetaZpPanel2");
1032  meWPullYvsBetaZpPanel2
1033  = dbe_->bookProfile(chisto, chisto, 20, pull_forward_pos_betal, pull_forward_pos_betah, 100, pullyl, pullyh, "");
1034 
1035  // all hits (not only from tracks)
1036  sprintf(chisto, "mePosxBarrel_all_hits");
1037  mePosxBarrel_all_hits = dbe_->book1D(chisto, chisto, 100, xl, xh);
1038  sprintf(chisto, "mePosyBarrel_all_hits");
1039  mePosyBarrel_all_hits = dbe_->book1D(chisto, chisto, 100, yl, yh);
1040 
1041  sprintf(chisto, "mePosxZmPanel1_all_hits");
1042  mePosxZmPanel1_all_hits = dbe_->book1D(chisto, chisto, 100, xl, xh);
1043  sprintf(chisto, "mePosyZmPanel1_all_hits");
1044  mePosyZmPanel1_all_hits = dbe_->book1D(chisto, chisto, 100, yl, yh);
1045  sprintf(chisto, "mePosxZmPanel2_all_hits");
1046  mePosxZmPanel2_all_hits = dbe_->book1D(chisto, chisto, 100, xl, xh);
1047  sprintf(chisto, "mePosyZmPanel2_all_hits");
1048  mePosyZmPanel2_all_hits = dbe_->book1D(chisto, chisto, 100, yl, yh);
1049 
1050  sprintf(chisto, "mePosxZpPanel1_all_hits");
1051  mePosxZpPanel1_all_hits = dbe_->book1D(chisto, chisto, 100, xl, xh);
1052  sprintf(chisto, "mePosyZpPanel1_all_hits");
1053  mePosyZpPanel1_all_hits = dbe_->book1D(chisto, chisto, 100, yl, yh);
1054  sprintf(chisto, "mePosxZpPanel2_all_hits");
1055  mePosxZpPanel2_all_hits = dbe_->book1D(chisto, chisto, 100, xl, xh);
1056  sprintf(chisto, "mePosyZpPanel2_all_hits");
1057  mePosyZpPanel2_all_hits = dbe_->book1D(chisto, chisto, 100, yl, yh);
1058 
1059  // control histograms
1060  meTracksPerEvent = dbe_->book1D("meTracksPerEvent" , "meTracksPerEvent" , 200, 0.0, 200.0);
1061  mePixRecHitsPerTrack = dbe_->book1D("mePixRecHitsPerTrack", "mePixRecHitsPerTrack", 6, 0.0, 6.0);
1062 
1063 }
1064 
1065 // Virtual destructor needed.
1066 SiPixelTrackingRecHitsValid::~SiPixelTrackingRecHitsValid()
1067 {
1068  if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
1069 }
1070 
1071 // Functions that gets called by framework every event
1072 void SiPixelTrackingRecHitsValid::analyze(const edm::Event& e, const edm::EventSetup& es)
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.getByLabel( "siPixelRecHits", 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  PXFDetId fdetid(detId);
1142  side = fdetid.side();
1143  disk = fdetid.disk();
1144  blade = fdetid.blade();
1145  panel = fdetid.panel();
1146  plaq = fdetid.module(); // 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 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 LogWarning("SiPixelTrackingRecHitsValid")<< "..............................................Wrong panel number !";
1175  } //else if ( side==2 )
1176  else LogWarning("SiPixelTrackingRecHitsValid") << ".......................................................Wrong side !" ;
1177 
1178  } // else if ( detId.subdetId()==PixelSubdetector::PixelEndcap )
1179  else 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.getByLabel(src_, 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  vector<PSimHit>::const_iterator closestit = matched.begin();
1278  for (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 (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  LogWarning("SiPixelTrackingRecHitsValid") << " ----- This is not good: n_assoc_muon = " << n_assoc_muon ;
1313  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  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  PXBDetId bdetid(detId);
1441  layer = bdetid.layer(); // Layer: 1,2,3.
1442  ladder = bdetid.ladder(); // Ladder: 1-20, 32, 44.
1443  mod = bdetid.module(); // 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  PXFDetId fdetid(detId);
1494  side = fdetid.side();
1495  disk = fdetid.disk();
1496  blade = fdetid.blade();
1497  panel = fdetid.panel();
1498  plaq = fdetid.module(); // 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 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 LogWarning("SiPixelTrackingRecHitsValid") << "..............................................Wrong panel number !";
1739  } //else if ( side==2 )
1740  else LogWarning("SiPixelTrackingRecHitsValid") << ".......................................................Wrong side !" ;
1741 
1742  } // else if ( detId.subdetId()==PixelSubdetector::PixelEndcap )
1743  else 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:189
beta
const double beta
Definition: CosmicMuonParameters.h:69
Surface::toGlobal
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:78
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:204
SiPixelTrackingRecHitsValid::meErryBarrelLayerModule
MonitorElement * meErryBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:102
SiPixelTrackingRecHitsValid::meNpixZmPanel1DiskPlaq
MonitorElement * meNpixZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:224
SiPixelTrackingRecHitsValid::meResYvsAlphaBarrelFlippedLaddersLayer
MonitorElement * meResYvsAlphaBarrelFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:201
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
SiPixelTrackingRecHitsValid::mePullYvsEtaBarrelLayerModule
MonitorElement * mePullYvsEtaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:122
i
int i
Definition: DBlmapReader.cc:9
SiPixelTrackingRecHitsValid::meErryZpPanel2DiskPlaq
MonitorElement * meErryZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:298
SiPixelTrackingRecHitsValid::meResyZpPanel1DiskPlaq
MonitorElement * meResyZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:274
LocalError::xx
float xx() const
Definition: LocalError.h:19
SiPixelTrackingRecHitsValid::mePullxZpPanel2DiskPlaq
MonitorElement * mePullxZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:301
alpha
float alpha
Definition: AMPTWrapper.h:95
SiPixelTrackingRecHitsValid::mePullXvsBetaBarrelLayerModule
MonitorElement * mePullXvsBetaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:117
SiPixelTrackingRecHitsValid::meResyZmPanel2DiskPlaq
MonitorElement * meResyZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:247
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:105
SiPixelTrackingRecHitsValid::mePullXvsAlphaZmPanel1DiskPlaq
MonitorElement * mePullXvsAlphaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:232
PXFDetId::panel
unsigned int panel() const
panel id
Definition: PXFDetId.h:52
SiPixelTrackingRecHitsValid::meResYvsAlphaBarrelNonFlippedLaddersLayer
MonitorElement * meResYvsAlphaBarrelNonFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:205
SiPixelTrackingRecHitsValid::mePullYvsAlphaBarrelNonFlippedLadders
MonitorElement * mePullYvsAlphaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:177
SiPixelTrackingRecHitsValid::mePullYvsAlphaZmPanel1DiskPlaq
MonitorElement * mePullYvsAlphaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:233
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:66
SiPixelTrackingRecHitsValid::meChargeZmPanel2DiskPlaq
MonitorElement * meChargeZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:253
SiPixelTrackingRecHitsValid::mePosxZmPanel1DiskPlaq
MonitorElement * mePosxZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:216
SiPixelTrackingRecHitsValid::mePosyBarrelLayerModule
MonitorElement * mePosyBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:100
SiPixelTrackingRecHitsValid::meResxZpPanel1DiskPlaq
MonitorElement * meResxZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:273
SiPixelTrackingRecHitsValid::meNpixZpPanel2DiskPlaq
MonitorElement * meNpixZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:303
SiPixelTrackingRecHitsValid::meResyZmPanel1DiskPlaq
MonitorElement * meResyZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:221
SiPixelTrackingRecHitsValid::meResXvsBetaZmPanel1DiskPlaq
MonitorElement * meResXvsBetaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:230
SiPixelTrackingRecHitsValid::mePullYvsEtaZpPanel1DiskPlaq
MonitorElement * mePullYvsEtaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:292
SiPixelTrackingRecHitsValid::meWPullXvsBetaBarrelNonFlippedLadders
MonitorElement * meWPullXvsBetaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:191
SiPixelTrackingRecHitsValid::meResYvsBetaZmPanel1DiskPlaq
MonitorElement * meResYvsBetaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:231
SiPixelTrackingRecHitsValid::mePullYvsEtaZmPanel2DiskPlaq
MonitorElement * mePullYvsEtaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:265
SiPixelTrackingRecHitsValid::meErrxZmPanel1DiskPlaq
MonitorElement * meErrxZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:218
SiPixelTrackingRecHitsValid::meResXvsBetaBarrelNonFlippedLaddersLayer
MonitorElement * meResXvsBetaBarrelNonFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:206
SiPixelTrackingRecHitsValid::meErrxZpPanel2DiskPlaq
MonitorElement * meErrxZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:297
DQMStore::save
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:1898
SiPixelTrackingRecHitsValid::meResYvsBetaBarrelFlippedLaddersLayer
MonitorElement * meResYvsBetaBarrelFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:203
SiPixelTrackingRecHitsValid::mePullxBarrelLayerModule
MonitorElement * mePullxBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:105
SiPixelTrackingRecHitsValid::mePosyZpPanel1DiskPlaq
MonitorElement * mePosyZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:270
SiPixelTrackingRecHitsValid::meNxpixZpPanel1DiskPlaq
MonitorElement * meNxpixZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:278
SiPixelTrackingRecHitsValid::mePullYvsBetaZpPanel1DiskPlaq
MonitorElement * mePullYvsBetaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:288
SiPixelTrackingRecHitsValid::mePullXvsPhiZmPanel1DiskPlaq
MonitorElement * mePullXvsPhiZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:236
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:10
PV3DBase::y
T y() const
Definition: PV3DBase.h:57
SiPixelTrackingRecHitsValid::meResXvsBetaBarrelFlippedLaddersLayer
MonitorElement * meResXvsBetaBarrelFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:202
SiPixelTrackingRecHitsValid::mePullyZmPanel1DiskPlaq
MonitorElement * mePullyZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:223
SiPixelTrackingRecHitsValid::meErrxBarrelLayerModule
MonitorElement * meErrxBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:101
SiPixelTrackingRecHitsValid::meErryZmPanel2DiskPlaq
MonitorElement * meErryZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:245
PixelTopology::nrows
virtual int nrows() const =0
SiPixelTrackingRecHitsValid::mePosxZmPanel2DiskPlaq
MonitorElement * mePosxZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:242
abs
#define abs(x)
Definition: mlp_lapack.h:159
SiPixelTrackingRecHitsValid::meChargeZpPanel2DiskPlaq
MonitorElement * meChargeZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:306
PXBDetId::ladder
unsigned int ladder() const
ladder id
Definition: PXBDetId.h:39
SiPixelTrackingRecHitsValid::mePullXvsEtaZpPanel2DiskPlaq
MonitorElement * mePullXvsEtaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:317
SiPixelTrackingRecHitsValid::meResYvsAlphaZmPanel1DiskPlaq
MonitorElement * meResYvsAlphaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:229
SiPixelTrackingRecHitsValid::mePullXvsBetaZmPanel1DiskPlaq
MonitorElement * mePullXvsBetaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:234
SiPixelTrackingRecHitsValid::meResXvsAlphaZmPanel2DiskPlaq
MonitorElement * meResXvsAlphaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:254
SiPixelTrackingRecHitsValid::mePullYvsAlphaZmPanel2DiskPlaq
MonitorElement * mePullYvsAlphaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:259
SiPixelTrackingRecHitsValid::mePullxZmPanel1DiskPlaq
MonitorElement * mePullxZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:222
SiPixelTrackingRecHitsValid::mePullyZpPanel2DiskPlaq
MonitorElement * mePullyZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:302
eta
T eta() const
Definition: Basic3DVectorLD.h:168
SiPixelTrackingRecHitsValid::mePullXvsAlphaBarrelLayerModule
MonitorElement * mePullXvsAlphaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:115
SiPixelTrackingRecHitsValid::meNpixZpPanel1DiskPlaq
MonitorElement * meNpixZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:277
SiPixelTrackingRecHitsValid::mePosyZmPanel2DiskPlaq
MonitorElement * mePosyZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:243
SiPixelTrackingRecHitsValid::mePullyBarrelLayerModule
MonitorElement * mePullyBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:106
SiPixelTrackingRecHitsValid::mePullXvsPhiZpPanel2DiskPlaq
MonitorElement * mePullXvsPhiZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:315
SiPixelTrackingRecHitsValid::meResyBarrelLayerModule
MonitorElement * meResyBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:104
SiPixelTrackingRecHitsValid::meWPullYvsBetaBarrelNonFlippedLadders
MonitorElement * meWPullYvsBetaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:192
DeDxDiscriminatorTools::charge
double charge(const std::vector< uint8_t > &Ampls)
Definition: DeDxDiscriminatorTools.cc:43
SiPixelTrackingRecHitsValid::meNypixBarrelLayerModule
MonitorElement * meNypixBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:109
position
static int position[TOTALCHAMBERS][3]
Definition: ReadPGInfo.cc:509
PXBDetId::layer
unsigned int layer() const
layer id
Definition: PXBDetId.h:35
SiPixelTrackingRecHitsValid::mePullyZpPanel1DiskPlaq
MonitorElement * mePullyZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:276
SiPixelTrackingRecHitsValid::meNypixZpPanel1DiskPlaq
MonitorElement * meNypixZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:279
SiPixelTrackingRecHitsValid::mePullXvsAlphaZpPanel1DiskPlaq
MonitorElement * mePullXvsAlphaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:285
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:147
SiPixelTrackingRecHitsValid::meErryZmPanel1DiskPlaq
MonitorElement * meErryZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:219
SiPixelTrackingRecHitsValid::mePullxZmPanel2DiskPlaq
MonitorElement * mePullxZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:248
PXFDetId::blade
unsigned int blade() const
blade id
Definition: PXFDetId.h:48
SiPixelTrackingRecHitsValid::mePullXvsEtaZpPanel1DiskPlaq
MonitorElement * mePullXvsEtaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:291
SiPixelTrackingRecHitsValid::mePullXvsPhiZmPanel2DiskPlaq
MonitorElement * mePullXvsPhiZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:262
SiPixelTrackingRecHitsValid::mePullYvsBetaZmPanel2DiskPlaq
MonitorElement * mePullYvsBetaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:261
SiPixelTrackingRecHitsValid::mePullXvsPhiBarrelLayerModule
MonitorElement * mePullXvsPhiBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:119
SiPixelTrackingRecHitsValid::mePullYvsPhiZmPanel2DiskPlaq
MonitorElement * mePullYvsPhiZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:263
SiPixelTrackingRecHitsValid::meWPullYvsAlphaBarrelNonFlippedLadders
MonitorElement * meWPullYvsAlphaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:190
SiPixelTrackingRecHitsValid::meResYvsAlphaZpPanel2DiskPlaq
MonitorElement * meResYvsAlphaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:308
LocalError::yy
float yy() const
Definition: LocalError.h:21
SiPixelTrackingRecHitsValid::meNpixBarrelLayerModule
MonitorElement * meNpixBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:107
SiPixelTrackingRecHitsValid::meResxZpPanel2DiskPlaq
MonitorElement * meResxZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:299
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:28
SiPixelTrackingRecHitsValid::meChargeZpPanel1DiskPlaq
MonitorElement * meChargeZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:280
SiPixelTrackingRecHitsValid::mePullXvsBetaZpPanel2DiskPlaq
MonitorElement * mePullXvsBetaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:313
PV3DBase::z
T z() const
Definition: PV3DBase.h:58
SiPixelTrackingRecHitsValid::meResYvsAlphaZpPanel1DiskPlaq
MonitorElement * meResYvsAlphaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:282
SiPixelTrackingRecHitsValid::meResxBarrelLayerModule
MonitorElement * meResxBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:103
SiPixelTrackingRecHitsValid::meNxpixZmPanel1DiskPlaq
MonitorElement * meNxpixZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:225
PXBDetId::module
unsigned int module() const
det id
Definition: PXBDetId.h:43
SiPixelTrackingRecHitsValid::mePosxZpPanel2DiskPlaq
MonitorElement * mePosxZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:295
j
int j
Definition: DBlmapReader.cc:9
PXFDetId::module
unsigned int module() const
det id
Definition: PXFDetId.h:56
patCandidatesForDimuonsSequences_cff.tracker
tuple tracker
Definition: patCandidatesForDimuonsSequences_cff.py:69
SiPixelTrackingRecHitsValid::mePosxZpPanel1DiskPlaq
MonitorElement * mePosxZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:269
SiPixelTrackingRecHitsValid::meNxpixZpPanel2DiskPlaq
MonitorElement * meNxpixZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:304
SiPixelTrackingRecHitsValid::mePullYvsEtaZpPanel2DiskPlaq
MonitorElement * mePullYvsEtaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:318
SiPixelTrackingRecHitsValid::mePullYvsPhiZpPanel1DiskPlaq
MonitorElement * mePullYvsPhiZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:290
SiPixelTrackingRecHitsValid::meResxZmPanel2DiskPlaq
MonitorElement * meResxZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:246
SiPixelTrackingRecHitsValid::meResXvsBetaZmPanel2DiskPlaq
MonitorElement * meResXvsBetaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:256
DetId::subdetId
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:39
SiPixelTrackingRecHitsValid::mePullyZmPanel2DiskPlaq
MonitorElement * mePullyZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:249
SiPixelTrackingRecHitsValid::mePullYvsBetaBarrelLayerModule
MonitorElement * mePullYvsBetaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:118
TrackerGeometry::idToDet
virtual const GeomDet * idToDet(DetId) const
Definition: TrackerGeometry.cc:135
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:355
SiPixelTrackingRecHitsValid::meErrxZmPanel2DiskPlaq
MonitorElement * meErrxZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:244
SiPixelTrackingRecHitsValid::meNypixZmPanel1DiskPlaq
MonitorElement * meNypixZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:226
PXFDetId::disk
unsigned int disk() const
disk id
Definition: PXFDetId.h:43
dbe_
DQMStore * dbe_
Definition: PFJetBenchmarkAnalyzer.cc:78
SiPixelTrackingRecHitsValid::mePullXvsEtaZmPanel1DiskPlaq
MonitorElement * mePullXvsEtaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:238
SiPixelTrackingRecHitsValid::meErrxZpPanel1DiskPlaq
MonitorElement * meErrxZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:271
SiPixelTrackingRecHitsValid::meErryZpPanel1DiskPlaq
MonitorElement * meErryZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:272
SiPixelTrackingRecHitsValid::meResXvsBetaZpPanel2DiskPlaq
MonitorElement * meResXvsBetaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:309
SiPixelTrackingRecHitsValid::meResXvsBetaZpPanel1DiskPlaq
MonitorElement * meResXvsBetaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:283
Local3DPoint
Point3DBase< float, LocalTag > Local3DPoint
Definition: LocalPoint.h:9
SiPixelTrackingRecHitsValid::mePullYvsPhiZmPanel1DiskPlaq
MonitorElement * mePullYvsPhiZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:237
DetId
Definition: DetId.h:20
SiPixelTrackingRecHitsValid::meChargeZmPanel1DiskPlaq
MonitorElement * meChargeZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:227
SiPixelTrackingRecHitsValid::mePullYvsAlphaZpPanel2DiskPlaq
MonitorElement * mePullYvsAlphaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:312
SiPixelTrackingRecHitsValid::meNypixZmPanel2DiskPlaq
MonitorElement * meNypixZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:252
SiPixelTrackingRecHitsValid::mePullYvsPhiZpPanel2DiskPlaq
MonitorElement * mePullYvsPhiZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:316
SiPixelTrackingRecHitsValid::mePullYvsBetaZmPanel1DiskPlaq
MonitorElement * mePullYvsBetaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:235
SiPixelTrackingRecHitsValid::mePullXvsEtaBarrelLayerModule
MonitorElement * mePullXvsEtaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:121
SiPixelTrackingRecHitsValid::mePullXvsBetaZmPanel2DiskPlaq
MonitorElement * mePullXvsBetaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:260
testEve_cfg.tracks
tuple tracks
Definition: testEve_cfg.py:39
SiPixelTrackingRecHitsValid::mePullXvsAlphaZpPanel2DiskPlaq
MonitorElement * mePullXvsAlphaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:311
SiPixelTrackingRecHitsValid::meResYvsBetaBarrelLayerModule
MonitorElement * meResYvsBetaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:114
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
SiPixelTrackingRecHitsValid::mePullXvsAlphaBarrelNonFlippedLadders
MonitorElement * mePullXvsAlphaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:176
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:118
SiPixelTrackingRecHitsValid::mePullYvsAlphaBarrelLayerModule
MonitorElement * mePullYvsAlphaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:116
SiPixelTrackingRecHitsValid::mePullYvsEtaZmPanel1DiskPlaq
MonitorElement * mePullYvsEtaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:239
SiPixelTrackingRecHitsValid::meResXvsAlphaZmPanel1DiskPlaq
MonitorElement * meResXvsAlphaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:228
edm::Handle::product
T const * product() const
Definition: Handle.h:74
SiPixelTrackingRecHitsValid::mePullYvsAlphaZpPanel1DiskPlaq
MonitorElement * mePullYvsAlphaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:286
SiPixelTrackingRecHitsValid::meChargeBarrelLayerModule
MonitorElement * meChargeBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:110
SiPixelTrackingRecHitsValid::meResXvsAlphaZpPanel1DiskPlaq
MonitorElement * meResXvsAlphaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:281
SiPixelTrackingRecHitsValid::mePullYvsBetaZpPanel2DiskPlaq
MonitorElement * mePullYvsBetaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:314
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
SiPixelTrackingRecHitsValid::meResYvsBetaZpPanel1DiskPlaq
MonitorElement * meResYvsBetaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:284
SiPixelTrackingRecHitsValid::mePosxBarrelLayerModule
MonitorElement * mePosxBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:99
edmNew::DetSet::end
iterator end()
Definition: DetSetNew.h:59
SiPixelTrackingRecHitsValid::meResYvsAlphaZmPanel2DiskPlaq
MonitorElement * meResYvsAlphaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:255
SiPixelTrackingRecHitsValid::meResyZpPanel2DiskPlaq
MonitorElement * meResyZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:300
SiPixelTrackingRecHitsValid::mePullXvsBetaZpPanel1DiskPlaq
MonitorElement * mePullXvsBetaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:287
SiPixelTrackingRecHitsValid::meResYvsBetaZpPanel2DiskPlaq
MonitorElement * meResYvsBetaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:310
TrackerHitAssociator::associateHit
std::vector< PSimHit > associateHit(const TrackingRecHit &thit)
Definition: TrackerHitAssociator.cc:120
SiPixelTrackingRecHitsValid::mePullXvsEtaZmPanel2DiskPlaq
MonitorElement * mePullXvsEtaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:264
SiPixelTrackingRecHitsValid::meResYvsAlphaBarrelLayerModule
MonitorElement * meResYvsAlphaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:112
PXFDetId::side
unsigned int side() const
positive or negative id
Definition: PXFDetId.h:38
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:296
SiPixelTrackingRecHitsValid::mePullXvsAlphaZmPanel2DiskPlaq
MonitorElement * mePullXvsAlphaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:258
SiPixelTrackingRecHitsValid::mePullxZpPanel1DiskPlaq
MonitorElement * mePullxZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:275
PV3DBase::x
T x() const
Definition: PV3DBase.h:56
SiPixelTrackingRecHitsValid::meResXvsAlphaZpPanel2DiskPlaq
MonitorElement * meResXvsAlphaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:307
mod
T mod(const T &a, const T &b)
Definition: ecalDccMap.h:4
SiPixelTrackingRecHitsValid::meResXvsAlphaBarrelLayerModule
MonitorElement * meResXvsAlphaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:111
SiPixelTrackingRecHitsValid::meResYvsBetaZmPanel2DiskPlaq
MonitorElement * meResYvsBetaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:257
SiPixelTrackingRecHitsValid::meNpixZmPanel2DiskPlaq
MonitorElement * meNpixZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:250
GeomDet::surface
virtual const BoundPlane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
SiPixelTrackingRecHitsValid::meResXvsAlphaBarrelFlippedLaddersLayer
MonitorElement * meResXvsAlphaBarrelFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:200
SiPixelTrackingRecHitsValid::meResXvsBetaBarrelLayerModule
MonitorElement * meResXvsBetaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:113
SiPixelTrackingRecHitsValid::mePullYvsPhiBarrelLayerModule
MonitorElement * mePullYvsPhiBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:120
SiPixelTrackingRecHitsValid::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: SiPixelTrackingRecHitsValid.cc:1072
SiPixelTrackingRecHitsValid::meNypixZpPanel2DiskPlaq
MonitorElement * meNypixZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:305
SiPixelTrackingRecHitsValid::mePosyZmPanel1DiskPlaq
MonitorElement * mePosyZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:217
SiPixelTrackingRecHitsValid::mePullXvsPhiZpPanel1DiskPlaq
MonitorElement * mePullXvsPhiZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:289
SiPixelTrackingRecHitsValid::meNxpixBarrelLayerModule
MonitorElement * meNxpixBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:108
SiPixelTrackingRecHitsValid::meNxpixZmPanel2DiskPlaq
MonitorElement * meNxpixZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:251
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:27
SiPixelTrackingRecHitsValid::meResxZmPanel1DiskPlaq
MonitorElement * meResxZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:220
phi
Definition: DDAxes.h:10
edmNew::DetSet::begin
iterator begin()
Definition: DetSetNew.h:56
SiPixelTrackingRecHitsValid::meResYvsBetaBarrelNonFlippedLaddersLayer
MonitorElement * meResYvsBetaBarrelNonFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:207
Generated for CMSSW Reference Manual by doxygen 1.8.5