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/TrackerCommon/interface/TrackerTopology.h"
32 #include "Geometry/Records/interface/IdealGeometryRecord.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  //Retrieve tracker topology from geometry
1075  edm::ESHandle<TrackerTopology> tTopo;
1076  es.get<IdealGeometryRecord>().get(tTopo);
1077 
1078 
1079  run = e.id().run();
1080  evt = e.id().event();
1081 
1082  // if ( evt%1000 == 0 )
1083  //cout << "evt = " << evt << endl;
1084 
1085  float math_pi = 3.14159265;
1086  float radtodeg = 180.0 / math_pi;
1087 
1088  DetId detId;
1089 
1090  LocalPoint position;
1091  LocalError error;
1092  float mindist = 999999.9;
1093 
1094  std::vector<PSimHit> matched;
1095  TrackerHitAssociator associate(e,conf_);
1096 
1097  edm::ESHandle<TrackerGeometry> pDD;
1098  es.get<TrackerDigiGeometryRecord> ().get (pDD);
1099  const TrackerGeometry* tracker = &(* pDD);
1100 
1101  if ( !MTCCtrack_ )
1102  {
1103  // --------------------------------------- all hits -----------------------------------------------------------
1104  //--- Fetch Pixel RecHits
1105  edm::Handle<SiPixelRecHitCollection> recHitColl;
1106  e.getByLabel( "siPixelRecHits", recHitColl);
1107 
1108  //cout <<" ----- Found "
1109  // << const_cast<SiPixelRecHitCollection*>(recHitColl.product())->size()
1110  // << " Pixel RecHits" << std::endl;
1111 
1112  //-----Iterate over detunits
1113  for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++)
1114  {
1115  DetId detId = ((*it)->geographicalId());
1116 
1117  unsigned int subid = detId.subdetId();
1118  if ( !((subid==1) || (subid==2)) )
1119  continue; // end subid if
1120 
1121  SiPixelRecHitCollection::const_iterator match = recHitColl->find(detId);
1122  if (match == recHitColl->end()) continue;
1123 
1124  SiPixelRecHitCollection::DetSet pixelrechitRange = *match;
1125  SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorBegin = pixelrechitRange.begin();
1126  SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorEnd = pixelrechitRange.end();
1127  SiPixelRecHitCollection::DetSet::const_iterator pixeliter = pixelrechitRangeIteratorBegin;
1128  std::vector<PSimHit> matched;
1129 
1130  //----Loop over rechits for this detId
1131  for ( ; pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter)
1132  {
1133  LocalPoint lp = pixeliter->localPosition();
1134  float rechitx = lp.x();
1135  float rechity = lp.y();
1136 
1137  detId = (*it)->geographicalId();
1138  subdetId = (int)detId.subdetId();
1139  if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelBarrel )
1140  {
1141  mePosxBarrel_all_hits->Fill( rechitx );
1142  mePosyBarrel_all_hits->Fill( rechity );
1143  }
1144  else if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelEndcap )
1145  {
1146 
1147  side = tTopo->pxfSide(detId);
1148  disk = tTopo->pxfDisk(detId);
1149  blade = tTopo->pxfBlade(detId);
1150  panel = tTopo->pxfPanel(detId);
1151  plaq = tTopo->pxfModule(detId); // also known as plaquette
1152 
1153  if ( side==1 )
1154  {
1155  if ( panel==1 )
1156  {
1157  mePosxZmPanel1_all_hits->Fill( rechitx );
1158  mePosyZmPanel1_all_hits->Fill( rechity );
1159  }
1160  else if ( panel==2 )
1161  {
1162  mePosxZmPanel2_all_hits->Fill( rechitx );
1163  mePosyZmPanel2_all_hits->Fill( rechity );
1164  }
1165  else LogWarning("SiPixelTrackingRecHitsValid") << "..............................................Wrong panel number !";
1166  } // if ( side==1 )
1167  else if ( side==2 )
1168  {
1169  if ( panel==1 )
1170  {
1171  mePosxZpPanel1_all_hits->Fill( rechitx );
1172  mePosyZpPanel1_all_hits->Fill( rechity );
1173  }
1174  else if ( panel==2 )
1175  {
1176  mePosxZpPanel2_all_hits->Fill( rechitx );
1177  mePosyZpPanel2_all_hits->Fill( rechity );
1178  }
1179  else LogWarning("SiPixelTrackingRecHitsValid")<< "..............................................Wrong panel number !";
1180  } //else if ( side==2 )
1181  else LogWarning("SiPixelTrackingRecHitsValid") << ".......................................................Wrong side !" ;
1182 
1183  } // else if ( detId.subdetId()==PixelSubdetector::PixelEndcap )
1184  else LogWarning("SiPixelTrackingRecHitsValid") << "Pixel rechit collection but we are not in the pixel detector" << (int)detId.subdetId() ;
1185 
1186  }
1187  }
1188  // ------------------------------------------------ all hits ---------------------------------------------------------------
1189 
1190  // Get tracks
1191  edm::Handle<reco::TrackCollection> trackCollection;
1192  e.getByLabel(src_, trackCollection);
1193  const reco::TrackCollection *tracks = trackCollection.product();
1194  reco::TrackCollection::const_iterator tciter;
1195 
1196  int n_tracks = (int)tracks->size(); // number of tracks in this event
1197  meTracksPerEvent->Fill( n_tracks );
1198 
1199  if ( tracks->size() > 0 )
1200  {
1201  // Loop on tracks
1202  for ( tciter=tracks->begin(); tciter!=tracks->end(); tciter++)
1203  {
1204  phi = tciter->momentum().phi() / math_pi*180.0;
1205  eta = tciter->momentum().eta();
1206 
1207  int n_hits = 0;
1208  // First loop on hits: find matched hits
1209  for ( trackingRecHit_iterator it = tciter->recHitsBegin(); it != tciter->recHitsEnd(); it++)
1210  {
1211  const TrackingRecHit &thit = **it;
1212  // Is it a matched hit?
1213  const SiPixelRecHit* matchedhit = dynamic_cast<const SiPixelRecHit*>(&thit);
1214 
1215  if ( matchedhit )
1216  {
1217  ++n_hits;
1218 
1219  layer = -9999;
1220  ladder = -9999;
1221  mod = -9999;
1222  side = -9999;
1223  disk = -9999;
1224  blade = -9999;
1225  panel = -9999;
1226  plaq = -9999;
1227 
1228  rechitx = -9999.9;
1229  rechity = -9999.9;
1230  rechitz = -9999.9;
1231  rechiterrx = -9999.9;
1232  rechiterry = -9999.9;
1233  rechitresx = -9999.9;
1234  rechitresy = -9999.9;
1235  rechitpullx = -9999.9;
1236  rechitpully = -9999.9;
1237 
1238  npix = -9999;
1239  nxpix = -9999;
1240  nypix = -9999;
1241  charge = -9999.9;
1242 
1243  alpha = -9999.9;
1244  beta = -9999.9;
1245 
1246  half = -9999;
1247  flipped = -9999;
1248 
1249  nsimhit = -9999;
1250 
1251  simhitx = -9999.9;
1252  simhity = -9999.9;
1253 
1254  position = (*it)->localPosition();
1255  error = (*it)->localPositionError();
1256 
1257  rechitx = position.x();
1258  rechity = position.y();
1259  rechitz = position.z();
1260  rechiterrx = sqrt(error.xx());
1261  rechiterry = sqrt(error.yy());
1262 
1263  npix = (*matchedhit).cluster()->size();
1264  nxpix = (*matchedhit).cluster()->sizeX();
1265  nypix = (*matchedhit).cluster()->sizeY();
1266  charge = (*matchedhit).cluster()->charge();
1267 
1268  //Association of the rechit to the simhit
1269  matched.clear();
1270  matched = associate.associateHit(*matchedhit);
1271 
1272  nsimhit = (int)matched.size();
1273 
1274  if ( !matched.empty() )
1275  {
1276  mindist = 999999.9;
1277  float distx, disty, dist;
1278  bool found_hit_from_generated_particle = false;
1279 
1280  int n_assoc_muon = 0;
1281 
1282  vector<PSimHit>::const_iterator closestit = matched.begin();
1283  for (vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); m++)
1284  {
1285  if ( checkType_ )
1286  {
1287  int pid = (*m).particleType();
1288  if ( abs(pid) != genType_ )
1289  continue;
1290 
1291  }
1292 
1293  float simhitx = 0.5 * ( (*m).entryPoint().x() + (*m).exitPoint().x() );
1294  float simhity = 0.5 * ( (*m).entryPoint().y() + (*m).exitPoint().y() );
1295 
1296  distx = fabs(rechitx - simhitx);
1297  disty = fabs(rechity - simhity);
1298  dist = sqrt( distx*distx + disty*disty );
1299 
1300  if ( dist < mindist )
1301  {
1302  n_assoc_muon++;
1303 
1304  mindist = dist;
1305  closestit = m;
1306  found_hit_from_generated_particle = true;
1307  }
1308  } // for (vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); m++)
1309 
1310  // This recHit does not have any simHit with the same particleType as the particles generated
1311  // Ignore it as most probably come from delta rays.
1312  if ( checkType_ && !found_hit_from_generated_particle )
1313  continue;
1314 
1315  if ( n_assoc_muon > 1 )
1316  {
1317  LogWarning("SiPixelTrackingRecHitsValid") << " ----- This is not good: n_assoc_muon = " << n_assoc_muon ;
1318  LogWarning("SiPixelTrackingRecHitsValid") << "evt = " << evt ;
1319  }
1320 
1321  pidhit = (*closestit).particleType();
1322 
1323  simhitx = 0.5*( (*closestit).entryPoint().x() + (*closestit).exitPoint().x() );
1324  simhity = 0.5*( (*closestit).entryPoint().y() + (*closestit).exitPoint().y() );
1325 
1326  rechitresx = rechitx - simhitx;
1327  rechitresy = rechity - simhity;
1328  rechitpullx = ( rechitx - simhitx ) / sqrt(error.xx());
1329  rechitpully = ( rechity - simhity ) / sqrt(error.yy());
1330 
1331  float simhitpx = (*closestit).momentumAtEntry().x();
1332  float simhitpy = (*closestit).momentumAtEntry().y();
1333  float simhitpz = (*closestit).momentumAtEntry().z();
1334 
1335  //beta = atan2(simhitpz, simhitpy) * radtodeg;
1336  //alpha = atan2(simhitpz, simhitpx) * radtodeg;
1337 
1338  beta = fabs(atan2(simhitpz, simhitpy)) * radtodeg;
1339  alpha = fabs(atan2(simhitpz, simhitpx)) * radtodeg;
1340 
1341  detId = (*it)->geographicalId();
1342 
1343  subdetId = (int)detId.subdetId();
1344 
1345  if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelBarrel )
1346  {
1347  mePosxBarrel->Fill( rechitx );
1348  mePosyBarrel->Fill( rechity );
1349  meErrxBarrel->Fill( rechiterrx );
1350  meErryBarrel->Fill( rechiterry );
1351  meResxBarrel->Fill( rechitresx );
1352  meResyBarrel->Fill( rechitresy );
1353  mePullxBarrel->Fill( rechitpullx );
1354  mePullyBarrel->Fill( rechitpully );
1355  meNpixBarrel->Fill( npix );
1356  meNxpixBarrel->Fill( nxpix );
1357  meNypixBarrel->Fill( nypix );
1358  meChargeBarrel->Fill( charge );
1359  meResXvsAlphaBarrel->Fill( alpha, fabs(rechitresx) );
1360  meResYvsAlphaBarrel->Fill( alpha, fabs(rechitresy) );
1361  meResXvsBetaBarrel->Fill( beta, fabs(rechitresx) );
1362  meResYvsBetaBarrel->Fill( beta, fabs(rechitresy) );
1363  mePullXvsAlphaBarrel->Fill( alpha, rechitpullx );
1364  mePullYvsAlphaBarrel->Fill( alpha, rechitpully );
1365  mePullXvsBetaBarrel->Fill( beta, rechitpullx );
1366  mePullYvsBetaBarrel->Fill( beta, rechitpully );
1367  mePullXvsPhiBarrel->Fill( phi, rechitpullx );
1368  mePullYvsPhiBarrel->Fill( phi, rechitpully );
1369  mePullXvsEtaBarrel->Fill( eta, rechitpullx );
1370  mePullYvsEtaBarrel->Fill( eta, rechitpully );
1371 
1372  const PixelGeomDetUnit * theGeomDet
1373  = dynamic_cast<const PixelGeomDetUnit*> ( tracker->idToDet(detId) );
1374  //const PixelTopology * topol = (&(theGeomDet->specificTopology()));
1375 
1376  int tmp_nrows = theGeomDet->specificTopology().nrows();
1377 
1378  if ( tmp_nrows == 80 )
1379  {
1380  mePosxBarrelHalfModule->Fill( rechitx );
1381  mePosyBarrelHalfModule->Fill( rechity );
1382  half = 1;
1383  }
1384  else if ( tmp_nrows == 160 )
1385  {
1386  mePosxBarrelFullModule->Fill( rechitx );
1387  mePosyBarrelFullModule->Fill( rechity );
1388  half = 0;
1389  }
1390  else
1391  LogWarning("SiPixelTrackingRecHitsValid") << "-------------------------------------------------- Wrong module size !!!";
1392 
1393  float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,0.)).perp();
1394  float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,1.)).perp();
1395 
1396  if ( tmp2<tmp1 )
1397  { // flipped
1398  mePosxBarrelFlippedLadders->Fill( rechitx );
1399  mePosyBarrelFlippedLadders->Fill( rechity );
1400  flipped = 1;
1401 
1402  meResXvsAlphaBarrelFlippedLadders->Fill( alpha, fabs(rechitresx) );
1403  meResYvsAlphaBarrelFlippedLadders->Fill( alpha, fabs(rechitresy) );
1404  meResXvsBetaBarrelFlippedLadders->Fill( beta, fabs(rechitresx) );
1405  meResYvsBetaBarrelFlippedLadders->Fill( beta, fabs(rechitresy) );
1406  mePullXvsAlphaBarrelFlippedLadders->Fill( alpha, rechitpullx );
1407  mePullYvsAlphaBarrelFlippedLadders->Fill( alpha, rechitpully );
1408  mePullXvsBetaBarrelFlippedLadders->Fill( beta, rechitpullx );
1409  mePullYvsBetaBarrelFlippedLadders->Fill( beta, rechitpully );
1410  mePullXvsPhiBarrelFlippedLadders->Fill( phi, rechitpullx );
1411  mePullYvsPhiBarrelFlippedLadders->Fill( phi, rechitpully );
1412  mePullXvsEtaBarrelFlippedLadders->Fill( eta, rechitpullx );
1413  mePullYvsEtaBarrelFlippedLadders->Fill( eta, rechitpully );
1414 
1415  meWPullXvsAlphaBarrelFlippedLadders->Fill( alpha, fabs(rechitpullx) );
1416  meWPullYvsAlphaBarrelFlippedLadders->Fill( alpha, fabs(rechitpully) );
1417  meWPullXvsBetaBarrelFlippedLadders->Fill( beta, fabs(rechitpullx) );
1418  meWPullYvsBetaBarrelFlippedLadders->Fill( beta, fabs(rechitpully) );
1419  }
1420  else
1421  { // not flipped
1422  mePosxBarrelNonFlippedLadders->Fill( rechitx );
1423  mePosyBarrelNonFlippedLadders->Fill( rechity );
1424  flipped = 0;
1425 
1426  meResXvsAlphaBarrelNonFlippedLadders->Fill( alpha, fabs(rechitresx) );
1427  meResYvsAlphaBarrelNonFlippedLadders->Fill( alpha, fabs(rechitresy) );
1428  meResXvsBetaBarrelNonFlippedLadders->Fill( beta, fabs(rechitresx) );
1429  meResYvsBetaBarrelNonFlippedLadders->Fill( beta, fabs(rechitresy) );
1430  mePullXvsAlphaBarrelNonFlippedLadders->Fill( alpha, rechitpullx );
1431  mePullYvsAlphaBarrelNonFlippedLadders->Fill( alpha, rechitpully );
1432  mePullXvsBetaBarrelNonFlippedLadders->Fill( beta, rechitpullx );
1433  mePullYvsBetaBarrelNonFlippedLadders->Fill( beta, rechitpully );
1434  mePullXvsPhiBarrelNonFlippedLadders->Fill( phi, rechitpullx );
1435  mePullYvsPhiBarrelNonFlippedLadders->Fill( phi, rechitpully );
1436  mePullXvsEtaBarrelNonFlippedLadders->Fill( eta, rechitpullx );
1437  mePullYvsEtaBarrelNonFlippedLadders->Fill( eta, rechitpully );
1438 
1439  meWPullXvsAlphaBarrelNonFlippedLadders->Fill( alpha, fabs(rechitpullx) );
1440  meWPullYvsAlphaBarrelNonFlippedLadders->Fill( alpha, fabs(rechitpully) );
1441  meWPullXvsBetaBarrelNonFlippedLadders->Fill( beta, fabs(rechitpullx) );
1442  meWPullYvsBetaBarrelNonFlippedLadders->Fill( beta, fabs(rechitpully) );
1443  }
1444 
1445 
1446  layer = tTopo->pxbLayer(detId); // Layer: 1,2,3.
1447  ladder = tTopo->pxbLadder(detId); // Ladder: 1-20, 32, 44.
1448  mod = tTopo->pxbModule(detId); // Mod: 1-8.
1449 
1450  mePosxBarrelLayerModule[layer-1][mod-1]->Fill( rechitx );
1451  mePosyBarrelLayerModule[layer-1][mod-1]->Fill( rechity );
1452  meErrxBarrelLayerModule[layer-1][mod-1]->Fill( rechiterrx );
1453  meErryBarrelLayerModule[layer-1][mod-1]->Fill( rechiterry );
1454  meResxBarrelLayerModule[layer-1][mod-1]->Fill( rechitresx );
1455  meResyBarrelLayerModule[layer-1][mod-1]->Fill( rechitresy );
1456  mePullxBarrelLayerModule[layer-1][mod-1]->Fill( rechitpullx );
1457  mePullyBarrelLayerModule[layer-1][mod-1]->Fill( rechitpully );
1458  meNpixBarrelLayerModule[layer-1][mod-1]->Fill( npix );
1459  meNxpixBarrelLayerModule[layer-1][mod-1]->Fill( nxpix );
1460  meNypixBarrelLayerModule[layer-1][mod-1]->Fill( nypix );
1461  meChargeBarrelLayerModule[layer-1][mod-1]->Fill( charge );
1462  meResXvsAlphaBarrelLayerModule[layer-1][mod-1]->Fill( alpha, fabs(rechitresx) );
1463  meResYvsAlphaBarrelLayerModule[layer-1][mod-1]->Fill( alpha, fabs(rechitresy) );
1464  meResXvsBetaBarrelLayerModule[layer-1][mod-1]->Fill( beta, fabs(rechitresx) );
1465  meResYvsBetaBarrelLayerModule[layer-1][mod-1]->Fill( beta, fabs(rechitresy) );
1466  mePullXvsAlphaBarrelLayerModule[layer-1][mod-1]->Fill( alpha, rechitpullx );
1467  mePullYvsAlphaBarrelLayerModule[layer-1][mod-1]->Fill( alpha, rechitpully );
1468  mePullXvsBetaBarrelLayerModule[layer-1][mod-1]->Fill( beta, rechitpullx );
1469  mePullYvsBetaBarrelLayerModule[layer-1][mod-1]->Fill( beta, rechitpully );
1470  mePullXvsPhiBarrelLayerModule[layer-1][mod-1]->Fill( phi, rechitpullx );
1471  mePullYvsPhiBarrelLayerModule[layer-1][mod-1]->Fill( phi, rechitpully );
1472  mePullXvsEtaBarrelLayerModule[layer-1][mod-1]->Fill( eta, rechitpullx );
1473  mePullYvsEtaBarrelLayerModule[layer-1][mod-1]->Fill( eta, rechitpully );
1474 
1475  meResxBarrelLayer[layer-1]->Fill( rechitresx );
1476  meResyBarrelLayer[layer-1]->Fill( rechitresy );
1477  mePullxBarrelLayer[layer-1]->Fill( rechitpullx );
1478  mePullyBarrelLayer[layer-1]->Fill( rechitpully );
1479 
1480  if ( tmp2<tmp1 )
1481  { // flipped
1482  meResXvsAlphaBarrelFlippedLaddersLayer[layer-1]->Fill( alpha, fabs(rechitresx) );
1483  meResYvsAlphaBarrelFlippedLaddersLayer[layer-1]->Fill( alpha, fabs(rechitresy) );
1484  meResXvsBetaBarrelFlippedLaddersLayer[layer-1]->Fill( beta, fabs(rechitresx) );
1485  meResYvsBetaBarrelFlippedLaddersLayer[layer-1]->Fill( beta, fabs(rechitresy) );
1486  }
1487  else
1488  { // not flipped
1489  meResXvsAlphaBarrelNonFlippedLaddersLayer[layer-1]->Fill( alpha, fabs(rechitresx) );
1490  meResYvsAlphaBarrelNonFlippedLaddersLayer[layer-1]->Fill( alpha, fabs(rechitresy) );
1491  meResXvsBetaBarrelNonFlippedLaddersLayer[layer-1]->Fill( beta, fabs(rechitresx) );
1492  meResYvsBetaBarrelNonFlippedLaddersLayer[layer-1]->Fill( beta, fabs(rechitresy) );
1493  }
1494 
1495  }
1496  else if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelEndcap )
1497  {
1498 
1499  side = tTopo->pxfSide(detId);
1500  disk = tTopo->pxfDisk(detId);
1501  blade = tTopo->pxfBlade(detId);
1502  panel = tTopo->pxfPanel(detId);
1503  plaq = tTopo->pxfModule(detId); // also known as plaquette
1504 
1505  if ( side==1 )
1506  {
1507  if ( panel==1 )
1508  {
1509  mePosxZmPanel1->Fill( rechitx );
1510  mePosyZmPanel1->Fill( rechity );
1511  meErrxZmPanel1->Fill( rechiterrx );
1512  meErryZmPanel1->Fill( rechiterry );
1513  meResxZmPanel1->Fill( rechitresx );
1514  meResyZmPanel1->Fill( rechitresy );
1515  mePullxZmPanel1->Fill( rechitpullx );
1516  mePullyZmPanel1->Fill( rechitpully );
1517  meNpixZmPanel1->Fill( npix );
1518  meNxpixZmPanel1->Fill( nxpix );
1519  meNypixZmPanel1->Fill( nypix );
1520  meChargeZmPanel1->Fill( charge );
1521  meResXvsAlphaZmPanel1->Fill( alpha, fabs(rechitresx) );
1522  meResYvsAlphaZmPanel1->Fill( alpha, fabs(rechitresy) );
1523  meResXvsBetaZmPanel1->Fill( beta, fabs(rechitresx) );
1524  meResYvsBetaZmPanel1->Fill( beta, fabs(rechitresy) );
1525  mePullXvsAlphaZmPanel1->Fill( alpha, rechitpullx );
1526  mePullYvsAlphaZmPanel1->Fill( alpha, rechitpully );
1527  mePullXvsBetaZmPanel1->Fill( beta, rechitpullx );
1528  mePullYvsBetaZmPanel1->Fill( beta, rechitpully );
1529  mePullXvsPhiZmPanel1->Fill( phi, rechitpullx );
1530  mePullYvsPhiZmPanel1->Fill( phi, rechitpully );
1531  mePullXvsEtaZmPanel1->Fill( eta, rechitpullx );
1532  mePullYvsEtaZmPanel1->Fill( eta, rechitpully );
1533 
1534  meWPullXvsAlphaZmPanel1->Fill( alpha, fabs(rechitpullx) );
1535  meWPullYvsAlphaZmPanel1->Fill( alpha, fabs(rechitpully) );
1536  meWPullXvsBetaZmPanel1->Fill( beta, fabs(rechitpullx) );
1537  meWPullYvsBetaZmPanel1->Fill( beta, fabs(rechitpully) );
1538 
1539  mePosxZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitx );
1540  mePosyZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechity );
1541  meErrxZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechiterrx );
1542  meErryZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechiterry );
1543  meResxZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitresx );
1544  meResyZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitresy );
1545  mePullxZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitpullx );
1546  mePullyZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitpully );
1547  meNpixZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( npix );
1548  meNxpixZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( nxpix );
1549  meNypixZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( nypix );
1550  meChargeZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( charge );
1551  meResXvsAlphaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresx) );
1552  meResYvsAlphaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresy) );
1553  meResXvsBetaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresx) );
1554  meResYvsBetaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresy) );
1555  mePullXvsAlphaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpullx );
1556  mePullYvsAlphaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpully );
1557  mePullXvsBetaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpullx );
1558  mePullYvsBetaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpully );
1559  mePullXvsPhiZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpullx );
1560  mePullYvsPhiZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpully );
1561  mePullXvsEtaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpullx );
1562  mePullYvsEtaZmPanel1DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpully );
1563 
1564  }
1565  else if ( panel==2 )
1566  {
1567  mePosxZmPanel2->Fill( rechitx );
1568  mePosyZmPanel2->Fill( rechity );
1569  meErrxZmPanel2->Fill( rechiterrx );
1570  meErryZmPanel2->Fill( rechiterry );
1571  meResxZmPanel2->Fill( rechitresx );
1572  meResyZmPanel2->Fill( rechitresy );
1573  mePullxZmPanel2->Fill( rechitpullx );
1574  mePullyZmPanel2->Fill( rechitpully );
1575  meNpixZmPanel2->Fill( npix );
1576  meNxpixZmPanel2->Fill( nxpix );
1577  meNypixZmPanel2->Fill( nypix );
1578  meChargeZmPanel2->Fill( charge );
1579  meResXvsAlphaZmPanel2->Fill( alpha, fabs(rechitresx) );
1580  meResYvsAlphaZmPanel2->Fill( alpha, fabs(rechitresy) );
1581  meResXvsBetaZmPanel2->Fill( beta, fabs(rechitresx) );
1582  meResYvsBetaZmPanel2->Fill( beta, fabs(rechitresy) );
1583  mePullXvsAlphaZmPanel2->Fill( alpha, rechitpullx );
1584  mePullYvsAlphaZmPanel2->Fill( alpha, rechitpully );
1585  mePullXvsBetaZmPanel2->Fill( beta, rechitpullx );
1586  mePullYvsBetaZmPanel2->Fill( beta, rechitpully );
1587  mePullXvsPhiZmPanel2->Fill( phi, rechitpullx );
1588  mePullYvsPhiZmPanel2->Fill( phi, rechitpully );
1589  mePullXvsEtaZmPanel2->Fill( eta, rechitpullx );
1590  mePullYvsEtaZmPanel2->Fill( eta, rechitpully );
1591 
1592  meWPullXvsAlphaZmPanel2->Fill( alpha, fabs(rechitpullx) );
1593  meWPullYvsAlphaZmPanel2->Fill( alpha, fabs(rechitpully) );
1594  meWPullXvsBetaZmPanel2->Fill( beta, fabs(rechitpullx) );
1595  meWPullYvsBetaZmPanel2->Fill( beta, fabs(rechitpully) );
1596 
1597  mePosxZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitx );
1598  mePosyZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechity );
1599  meErrxZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechiterrx );
1600  meErryZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechiterry );
1601  meResxZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitresx );
1602  meResyZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitresy );
1603  mePullxZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitpullx );
1604  mePullyZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitpully );
1605  meNpixZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( npix );
1606  meNxpixZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( nxpix );
1607  meNypixZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( nypix );
1608  meChargeZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( charge );
1609  meResXvsAlphaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresx) );
1610  meResYvsAlphaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresy) );
1611  meResXvsBetaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresx) );
1612  meResYvsBetaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresy) );
1613  mePullXvsAlphaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpullx );
1614  mePullYvsAlphaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpully );
1615  mePullXvsBetaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpullx );
1616  mePullYvsBetaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpully );
1617  mePullXvsPhiZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpullx );
1618  mePullYvsPhiZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpully );
1619  mePullXvsEtaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpullx );
1620  mePullYvsEtaZmPanel2DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpully );
1621 
1622  }
1623  else LogWarning("SiPixelTrackingRecHitsValid") << "..............................................Wrong panel number !";
1624  } // if ( side==1 )
1625  else if ( side==2 )
1626  {
1627  if ( panel==1 )
1628  {
1629  mePosxZpPanel1->Fill( rechitx );
1630  mePosyZpPanel1->Fill( rechity );
1631  meErrxZpPanel1->Fill( rechiterrx );
1632  meErryZpPanel1->Fill( rechiterry );
1633  meResxZpPanel1->Fill( rechitresx );
1634  meResyZpPanel1->Fill( rechitresy );
1635  mePullxZpPanel1->Fill( rechitpullx );
1636  mePullyZpPanel1->Fill( rechitpully );
1637  meNpixZpPanel1->Fill( npix );
1638  meNxpixZpPanel1->Fill( nxpix );
1639  meNypixZpPanel1->Fill( nypix );
1640  meChargeZpPanel1->Fill( charge );
1641  meResXvsAlphaZpPanel1->Fill( alpha, fabs(rechitresx) );
1642  meResYvsAlphaZpPanel1->Fill( alpha, fabs(rechitresy) );
1643  meResXvsBetaZpPanel1->Fill( beta, fabs(rechitresx) );
1644  meResYvsBetaZpPanel1->Fill( beta, fabs(rechitresy) );
1645  mePullXvsAlphaZpPanel1->Fill( alpha, rechitpullx );
1646  mePullYvsAlphaZpPanel1->Fill( alpha, rechitpully );
1647  mePullXvsBetaZpPanel1->Fill( beta, rechitpullx );
1648  mePullYvsBetaZpPanel1->Fill( beta, rechitpully );
1649  mePullXvsPhiZpPanel1->Fill( phi, rechitpullx );
1650  mePullYvsPhiZpPanel1->Fill( phi, rechitpully );
1651  mePullXvsEtaZpPanel1->Fill( eta, rechitpullx );
1652  mePullYvsEtaZpPanel1->Fill( eta, rechitpully );
1653 
1654  meWPullXvsAlphaZpPanel1->Fill( alpha, fabs(rechitpullx) );
1655  meWPullYvsAlphaZpPanel1->Fill( alpha, fabs(rechitpully) );
1656  meWPullXvsBetaZpPanel1->Fill( beta, fabs(rechitpullx) );
1657  meWPullYvsBetaZpPanel1->Fill( beta, fabs(rechitpully) );
1658 
1659  mePosxZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitx );
1660  mePosyZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechity );
1661  meErrxZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechiterrx );
1662  meErryZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechiterry );
1663  meResxZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitresx );
1664  meResyZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitresy );
1665  mePullxZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitpullx );
1666  mePullyZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( rechitpully );
1667  meNpixZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( npix );
1668  meNxpixZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( nxpix );
1669  meNypixZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( nypix );
1670  meChargeZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( charge );
1671  meResXvsAlphaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresx) );
1672  meResYvsAlphaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresy) );
1673  meResXvsBetaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresx) );
1674  meResYvsBetaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresy) );
1675  mePullXvsAlphaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpullx );
1676  mePullYvsAlphaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpully );
1677  mePullXvsBetaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpullx );
1678  mePullYvsBetaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpully );
1679  mePullXvsPhiZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpullx );
1680  mePullYvsPhiZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpully );
1681  mePullXvsEtaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpullx );
1682  mePullYvsEtaZpPanel1DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpully );
1683 
1684  }
1685  else if ( panel==2 )
1686  {
1687  mePosxZpPanel2->Fill( rechitx );
1688  mePosyZpPanel2->Fill( rechity );
1689  meErrxZpPanel2->Fill( rechiterrx );
1690  meErryZpPanel2->Fill( rechiterry );
1691  meResxZpPanel2->Fill( rechitresx );
1692  meResyZpPanel2->Fill( rechitresy );
1693  mePullxZpPanel2->Fill( rechitpullx );
1694  mePullyZpPanel2->Fill( rechitpully );
1695  meNpixZpPanel2->Fill( npix );
1696  meNxpixZpPanel2->Fill( nxpix );
1697  meNypixZpPanel2->Fill( nypix );
1698  meChargeZpPanel2->Fill( charge );
1699  meResXvsAlphaZpPanel2->Fill( alpha, fabs(rechitresx) );
1700  meResYvsAlphaZpPanel2->Fill( alpha, fabs(rechitresy) );
1701  meResXvsBetaZpPanel2->Fill( beta, fabs(rechitresx) );
1702  meResYvsBetaZpPanel2->Fill( beta, fabs(rechitresy) );
1703  mePullXvsAlphaZpPanel2->Fill( alpha, rechitpullx );
1704  mePullYvsAlphaZpPanel2->Fill( alpha, rechitpully );
1705  mePullXvsBetaZpPanel2->Fill( beta, rechitpullx );
1706  mePullYvsBetaZpPanel2->Fill( beta, rechitpully );
1707  mePullXvsPhiZpPanel2->Fill( phi, rechitpullx );
1708  mePullYvsPhiZpPanel2->Fill( phi, rechitpully );
1709  mePullXvsEtaZpPanel2->Fill( eta, rechitpullx );
1710  mePullYvsEtaZpPanel2->Fill( eta, rechitpully );
1711 
1712  meWPullXvsAlphaZpPanel2->Fill( alpha, fabs(rechitpullx) );
1713  meWPullYvsAlphaZpPanel2->Fill( alpha, fabs(rechitpully) );
1714  meWPullXvsBetaZpPanel2->Fill( beta, fabs(rechitpullx) );
1715  meWPullYvsBetaZpPanel2->Fill( beta, fabs(rechitpully) );
1716 
1717  mePosxZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitx );
1718  mePosyZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechity );
1719  meErrxZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechiterrx );
1720  meErryZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechiterry );
1721  meResxZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitresx );
1722  meResyZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitresy );
1723  mePullxZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitpullx );
1724  mePullyZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( rechitpully );
1725  meNpixZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( npix );
1726  meNxpixZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( nxpix );
1727  meNypixZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( nypix );
1728  meChargeZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( charge );
1729  meResXvsAlphaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresx) );
1730  meResYvsAlphaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, fabs(rechitresy) );
1731  meResXvsBetaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresx) );
1732  meResYvsBetaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, fabs(rechitresy) );
1733  mePullXvsAlphaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpullx );
1734  mePullYvsAlphaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( alpha, rechitpully );
1735  mePullXvsBetaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpullx );
1736  mePullYvsBetaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( beta, rechitpully );
1737  mePullXvsPhiZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpullx );
1738  mePullYvsPhiZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( phi, rechitpully );
1739  mePullXvsEtaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpullx );
1740  mePullYvsEtaZpPanel2DiskPlaq[disk-1][plaq-1]->Fill( eta, rechitpully );
1741 
1742  }
1743  else LogWarning("SiPixelTrackingRecHitsValid") << "..............................................Wrong panel number !";
1744  } //else if ( side==2 )
1745  else LogWarning("SiPixelTrackingRecHitsValid") << ".......................................................Wrong side !" ;
1746 
1747  } // else if ( detId.subdetId()==PixelSubdetector::PixelEndcap )
1748  else LogWarning("SiPixelTrackingRecHitsValid") << "Pixel rechit but we are not in the pixel detector" << (int)detId.subdetId() ;
1749 
1750  if(debugNtuple_.size()!=0)t_->Fill();
1751 
1752  } // if ( !matched.empty() )
1753  //else
1754  //cout << "---------------- RecHit with no associated SimHit !!! -------------------------- " << endl;
1755 
1756  } // matchedhit.
1757 
1758  } // end of loop on hits
1759 
1760  mePixRecHitsPerTrack->Fill( n_hits );
1761  //cout << "n_hits = " << n_hits << endl;
1762 
1763  } //end of loop on track
1764 
1765  } // tracks > 0.
1766 
1767  } //end of MTCCTrack
1768 
1769 }
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:42
SiPixelTrackingRecHitsValid::meWPullXvsAlphaBarrelNonFlippedLadders
MonitorElement * meWPullXvsAlphaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:188
beta
const double beta
Definition: CosmicMuonParameters.h:69
Surface::toGlobal
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:114
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:203
SiPixelTrackingRecHitsValid::meErryBarrelLayerModule
MonitorElement * meErryBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:101
SiPixelTrackingRecHitsValid::meNpixZmPanel1DiskPlaq
MonitorElement * meNpixZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:223
SiPixelTrackingRecHitsValid::meResYvsAlphaBarrelFlippedLaddersLayer
MonitorElement * meResYvsAlphaBarrelFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:200
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
SiPixelTrackingRecHitsValid::mePullYvsEtaBarrelLayerModule
MonitorElement * mePullYvsEtaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:121
i
int i
Definition: DBlmapReader.cc:9
SiPixelTrackingRecHitsValid::meErryZpPanel2DiskPlaq
MonitorElement * meErryZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:297
SiPixelTrackingRecHitsValid::meResyZpPanel1DiskPlaq
MonitorElement * meResyZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:273
LocalError::xx
float xx() const
Definition: LocalError.h:24
SiPixelTrackingRecHitsValid::mePullxZpPanel2DiskPlaq
MonitorElement * mePullxZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:300
alpha
float alpha
Definition: AMPTWrapper.h:95
SiPixelTrackingRecHitsValid::mePullXvsBetaBarrelLayerModule
MonitorElement * mePullXvsBetaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:116
SiPixelTrackingRecHitsValid::meResyZmPanel2DiskPlaq
MonitorElement * meResyZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:246
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:107
SiPixelTrackingRecHitsValid::mePullXvsAlphaZmPanel1DiskPlaq
MonitorElement * mePullXvsAlphaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:231
SiPixelTrackingRecHitsValid::meResYvsAlphaBarrelNonFlippedLaddersLayer
MonitorElement * meResYvsAlphaBarrelNonFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:204
SiPixelTrackingRecHitsValid::mePullYvsAlphaBarrelNonFlippedLadders
MonitorElement * mePullYvsAlphaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:176
SiPixelTrackingRecHitsValid::mePullYvsAlphaZmPanel1DiskPlaq
MonitorElement * mePullYvsAlphaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:232
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:72
SiPixelTrackingRecHitsValid::meChargeZmPanel2DiskPlaq
MonitorElement * meChargeZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:252
SiPixelTrackingRecHitsValid::mePosxZmPanel1DiskPlaq
MonitorElement * mePosxZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:215
SiPixelTrackingRecHitsValid::mePosyBarrelLayerModule
MonitorElement * mePosyBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:99
SiPixelTrackingRecHitsValid::meResxZpPanel1DiskPlaq
MonitorElement * meResxZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:272
SiPixelTrackingRecHitsValid::meNpixZpPanel2DiskPlaq
MonitorElement * meNpixZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:302
SiPixelTrackingRecHitsValid::meResyZmPanel1DiskPlaq
MonitorElement * meResyZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:220
SiPixelTrackingRecHitsValid::meResXvsBetaZmPanel1DiskPlaq
MonitorElement * meResXvsBetaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:229
SiPixelTrackingRecHitsValid::mePullYvsEtaZpPanel1DiskPlaq
MonitorElement * mePullYvsEtaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:291
SiPixelTrackingRecHitsValid::meWPullXvsBetaBarrelNonFlippedLadders
MonitorElement * meWPullXvsBetaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:190
SiPixelTrackingRecHitsValid::meResYvsBetaZmPanel1DiskPlaq
MonitorElement * meResYvsBetaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:230
SiPixelTrackingRecHitsValid::mePullYvsEtaZmPanel2DiskPlaq
MonitorElement * mePullYvsEtaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:264
SiPixelTrackingRecHitsValid::meErrxZmPanel1DiskPlaq
MonitorElement * meErrxZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:217
SiPixelTrackingRecHitsValid::meResXvsBetaBarrelNonFlippedLaddersLayer
MonitorElement * meResXvsBetaBarrelNonFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:205
SiPixelTrackingRecHitsValid::meErrxZpPanel2DiskPlaq
MonitorElement * meErrxZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:296
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:2118
SiPixelTrackingRecHitsValid::meResYvsBetaBarrelFlippedLaddersLayer
MonitorElement * meResYvsBetaBarrelFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:202
SiPixelTrackingRecHitsValid::mePullxBarrelLayerModule
MonitorElement * mePullxBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:104
SiPixelTrackingRecHitsValid::mePosyZpPanel1DiskPlaq
MonitorElement * mePosyZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:269
SiPixelTrackingRecHitsValid::meNxpixZpPanel1DiskPlaq
MonitorElement * meNxpixZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:277
SiPixelTrackingRecHitsValid::mePullYvsBetaZpPanel1DiskPlaq
MonitorElement * mePullYvsBetaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:287
SiPixelTrackingRecHitsValid::mePullXvsPhiZmPanel1DiskPlaq
MonitorElement * mePullXvsPhiZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:235
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:10
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
SiPixelTrackingRecHitsValid::meResXvsBetaBarrelFlippedLaddersLayer
MonitorElement * meResXvsBetaBarrelFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:201
SiPixelTrackingRecHitsValid::mePullyZmPanel1DiskPlaq
MonitorElement * mePullyZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:222
SiPixelTrackingRecHitsValid::meErrxBarrelLayerModule
MonitorElement * meErrxBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:100
SiPixelTrackingRecHitsValid::meErryZmPanel2DiskPlaq
MonitorElement * meErryZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:244
PixelTopology::nrows
virtual int nrows() const =0
SiPixelTrackingRecHitsValid::mePosxZmPanel2DiskPlaq
MonitorElement * mePosxZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:241
abs
#define abs(x)
Definition: mlp_lapack.h:159
SiPixelTrackingRecHitsValid::meChargeZpPanel2DiskPlaq
MonitorElement * meChargeZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:305
SiPixelTrackingRecHitsValid::mePullXvsEtaZpPanel2DiskPlaq
MonitorElement * mePullXvsEtaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:316
SiPixelTrackingRecHitsValid::meResYvsAlphaZmPanel1DiskPlaq
MonitorElement * meResYvsAlphaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:228
SiPixelTrackingRecHitsValid::mePullXvsBetaZmPanel1DiskPlaq
MonitorElement * mePullXvsBetaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:233
SiPixelTrackingRecHitsValid::meResXvsAlphaZmPanel2DiskPlaq
MonitorElement * meResXvsAlphaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:253
SiPixelTrackingRecHitsValid::mePullYvsAlphaZmPanel2DiskPlaq
MonitorElement * mePullYvsAlphaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:258
SiPixelTrackingRecHitsValid::mePullxZmPanel1DiskPlaq
MonitorElement * mePullxZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:221
SiPixelTrackingRecHitsValid::mePullyZpPanel2DiskPlaq
MonitorElement * mePullyZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:301
eta
T eta() const
Definition: Basic3DVectorLD.h:179
SiPixelTrackingRecHitsValid::mePullXvsAlphaBarrelLayerModule
MonitorElement * mePullXvsAlphaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:114
SiPixelTrackingRecHitsValid::meNpixZpPanel1DiskPlaq
MonitorElement * meNpixZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:276
SiPixelTrackingRecHitsValid::mePosyZmPanel2DiskPlaq
MonitorElement * mePosyZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:242
SiPixelTrackingRecHitsValid::mePullyBarrelLayerModule
MonitorElement * mePullyBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:105
SiPixelTrackingRecHitsValid::mePullXvsPhiZpPanel2DiskPlaq
MonitorElement * mePullXvsPhiZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:314
SiPixelTrackingRecHitsValid::meResyBarrelLayerModule
MonitorElement * meResyBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:103
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:35
SiPixelTrackingRecHitsValid::meWPullYvsBetaBarrelNonFlippedLadders
MonitorElement * meWPullYvsBetaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:191
DeDxDiscriminatorTools::charge
double charge(const std::vector< uint8_t > &Ampls)
Definition: DeDxDiscriminatorTools.cc:43
SiPixelTrackingRecHitsValid::meNypixBarrelLayerModule
MonitorElement * meNypixBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:108
position
static int position[TOTALCHAMBERS][3]
Definition: ReadPGInfo.cc:509
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
SiPixelTrackingRecHitsValid::mePullyZpPanel1DiskPlaq
MonitorElement * mePullyZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:275
SiPixelTrackingRecHitsValid::meNypixZpPanel1DiskPlaq
MonitorElement * meNypixZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:278
SiPixelTrackingRecHitsValid::mePullXvsAlphaZpPanel1DiskPlaq
MonitorElement * mePullXvsAlphaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:284
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:152
SiPixelTrackingRecHitsValid::meErryZmPanel1DiskPlaq
MonitorElement * meErryZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:218
SiPixelTrackingRecHitsValid::mePullxZmPanel2DiskPlaq
MonitorElement * mePullxZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:247
SiPixelTrackingRecHitsValid::mePullXvsEtaZpPanel1DiskPlaq
MonitorElement * mePullXvsEtaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:290
SiPixelTrackingRecHitsValid::mePullXvsPhiZmPanel2DiskPlaq
MonitorElement * mePullXvsPhiZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:261
SiPixelTrackingRecHitsValid::mePullYvsBetaZmPanel2DiskPlaq
MonitorElement * mePullYvsBetaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:260
SiPixelTrackingRecHitsValid::mePullXvsPhiBarrelLayerModule
MonitorElement * mePullXvsPhiBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:118
SiPixelTrackingRecHitsValid::mePullYvsPhiZmPanel2DiskPlaq
MonitorElement * mePullYvsPhiZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:262
SiPixelTrackingRecHitsValid::meWPullYvsAlphaBarrelNonFlippedLadders
MonitorElement * meWPullYvsAlphaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:189
SiPixelTrackingRecHitsValid::meResYvsAlphaZpPanel2DiskPlaq
MonitorElement * meResYvsAlphaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:307
LocalError::yy
float yy() const
Definition: LocalError.h:26
SiPixelTrackingRecHitsValid::meNpixBarrelLayerModule
MonitorElement * meNpixBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:106
SiPixelTrackingRecHitsValid::meResxZpPanel2DiskPlaq
MonitorElement * meResxZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:298
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
SiPixelTrackingRecHitsValid::meChargeZpPanel1DiskPlaq
MonitorElement * meChargeZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:279
SiPixelTrackingRecHitsValid::mePullXvsBetaZpPanel2DiskPlaq
MonitorElement * mePullXvsBetaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:312
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
SiPixelTrackingRecHitsValid::meResYvsAlphaZpPanel1DiskPlaq
MonitorElement * meResYvsAlphaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:281
SiPixelTrackingRecHitsValid::meResxBarrelLayerModule
MonitorElement * meResxBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:102
SiPixelTrackingRecHitsValid::meNxpixZmPanel1DiskPlaq
MonitorElement * meNxpixZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:224
SiPixelTrackingRecHitsValid::mePosxZpPanel2DiskPlaq
MonitorElement * mePosxZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:294
j
int j
Definition: DBlmapReader.cc:9
patCandidatesForDimuonsSequences_cff.tracker
tuple tracker
Definition: patCandidatesForDimuonsSequences_cff.py:69
SiPixelTrackingRecHitsValid::mePosxZpPanel1DiskPlaq
MonitorElement * mePosxZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:268
SiPixelTrackingRecHitsValid::meNxpixZpPanel2DiskPlaq
MonitorElement * meNxpixZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:303
SiPixelTrackingRecHitsValid::mePullYvsEtaZpPanel2DiskPlaq
MonitorElement * mePullYvsEtaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:317
SiPixelTrackingRecHitsValid::mePullYvsPhiZpPanel1DiskPlaq
MonitorElement * mePullYvsPhiZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:289
SiPixelTrackingRecHitsValid::meResxZmPanel2DiskPlaq
MonitorElement * meResxZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:245
SiPixelTrackingRecHitsValid::meResXvsBetaZmPanel2DiskPlaq
MonitorElement * meResXvsBetaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:255
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:248
SiPixelTrackingRecHitsValid::mePullYvsBetaBarrelLayerModule
MonitorElement * mePullYvsBetaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:117
TrackerGeometry::idToDet
virtual const GeomDet * idToDet(DetId) const
Definition: TrackerGeometry.cc:137
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:361
SiPixelTrackingRecHitsValid::meErrxZmPanel2DiskPlaq
MonitorElement * meErrxZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:243
SiPixelTrackingRecHitsValid::meNypixZmPanel1DiskPlaq
MonitorElement * meNypixZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:225
dbe_
DQMStore * dbe_
Definition: PFJetBenchmarkAnalyzer.cc:78
SiPixelTrackingRecHitsValid::mePullXvsEtaZmPanel1DiskPlaq
MonitorElement * mePullXvsEtaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:237
SiPixelTrackingRecHitsValid::meErrxZpPanel1DiskPlaq
MonitorElement * meErrxZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:270
SiPixelTrackingRecHitsValid::meErryZpPanel1DiskPlaq
MonitorElement * meErryZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:271
SiPixelTrackingRecHitsValid::meResXvsBetaZpPanel2DiskPlaq
MonitorElement * meResXvsBetaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:308
SiPixelTrackingRecHitsValid::meResXvsBetaZpPanel1DiskPlaq
MonitorElement * meResXvsBetaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:282
Local3DPoint
Point3DBase< float, LocalTag > Local3DPoint
Definition: LocalPoint.h:9
SiPixelTrackingRecHitsValid::mePullYvsPhiZmPanel1DiskPlaq
MonitorElement * mePullYvsPhiZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:236
DetId
Definition: DetId.h:20
SiPixelTrackingRecHitsValid::meChargeZmPanel1DiskPlaq
MonitorElement * meChargeZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:226
SiPixelTrackingRecHitsValid::mePullYvsAlphaZpPanel2DiskPlaq
MonitorElement * mePullYvsAlphaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:311
SiPixelTrackingRecHitsValid::meNypixZmPanel2DiskPlaq
MonitorElement * meNypixZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:251
SiPixelTrackingRecHitsValid::mePullYvsPhiZpPanel2DiskPlaq
MonitorElement * mePullYvsPhiZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:315
SiPixelTrackingRecHitsValid::mePullYvsBetaZmPanel1DiskPlaq
MonitorElement * mePullYvsBetaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:234
SiPixelTrackingRecHitsValid::mePullXvsEtaBarrelLayerModule
MonitorElement * mePullXvsEtaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:120
SiPixelTrackingRecHitsValid::mePullXvsBetaZmPanel2DiskPlaq
MonitorElement * mePullXvsBetaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:259
testEve_cfg.tracks
tuple tracks
Definition: testEve_cfg.py:39
SiPixelTrackingRecHitsValid::mePullXvsAlphaZpPanel2DiskPlaq
MonitorElement * mePullXvsAlphaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:310
SiPixelTrackingRecHitsValid::meResYvsBetaBarrelLayerModule
MonitorElement * meResYvsBetaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:113
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
SiPixelTrackingRecHitsValid::mePullXvsAlphaBarrelNonFlippedLadders
MonitorElement * mePullXvsAlphaBarrelNonFlippedLadders
Definition: SiPixelTrackingRecHitsValid.h:175
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:115
SiPixelTrackingRecHitsValid::mePullYvsEtaZmPanel1DiskPlaq
MonitorElement * mePullYvsEtaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:238
SiPixelTrackingRecHitsValid::meResXvsAlphaZmPanel1DiskPlaq
MonitorElement * meResXvsAlphaZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:227
edm::Handle::product
T const * product() const
Definition: Handle.h:74
SiPixelTrackingRecHitsValid::mePullYvsAlphaZpPanel1DiskPlaq
MonitorElement * mePullYvsAlphaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:285
SiPixelTrackingRecHitsValid::meChargeBarrelLayerModule
MonitorElement * meChargeBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:109
SiPixelTrackingRecHitsValid::meResXvsAlphaZpPanel1DiskPlaq
MonitorElement * meResXvsAlphaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:280
SiPixelTrackingRecHitsValid::mePullYvsBetaZpPanel2DiskPlaq
MonitorElement * mePullYvsBetaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:313
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
SiPixelTrackingRecHitsValid::meResYvsBetaZpPanel1DiskPlaq
MonitorElement * meResYvsBetaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:283
SiPixelTrackingRecHitsValid::mePosxBarrelLayerModule
MonitorElement * mePosxBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:98
edmNew::DetSet::end
iterator end()
Definition: DetSetNew.h:59
SiPixelTrackingRecHitsValid::meResYvsAlphaZmPanel2DiskPlaq
MonitorElement * meResYvsAlphaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:254
SiPixelTrackingRecHitsValid::meResyZpPanel2DiskPlaq
MonitorElement * meResyZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:299
SiPixelTrackingRecHitsValid::mePullXvsBetaZpPanel1DiskPlaq
MonitorElement * mePullXvsBetaZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:286
SiPixelTrackingRecHitsValid::meResYvsBetaZpPanel2DiskPlaq
MonitorElement * meResYvsBetaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:309
TrackerHitAssociator::associateHit
std::vector< PSimHit > associateHit(const TrackingRecHit &thit)
Definition: TrackerHitAssociator.cc:132
SiPixelTrackingRecHitsValid::mePullXvsEtaZmPanel2DiskPlaq
MonitorElement * mePullXvsEtaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:263
SiPixelTrackingRecHitsValid::meResYvsAlphaBarrelLayerModule
MonitorElement * meResYvsAlphaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:111
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:295
SiPixelTrackingRecHitsValid::mePullXvsAlphaZmPanel2DiskPlaq
MonitorElement * mePullXvsAlphaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:257
SiPixelTrackingRecHitsValid::mePullxZpPanel1DiskPlaq
MonitorElement * mePullxZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:274
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
SiPixelTrackingRecHitsValid::meResXvsAlphaZpPanel2DiskPlaq
MonitorElement * meResXvsAlphaZpPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:306
mod
T mod(const T &a, const T &b)
Definition: ecalDccMap.h:4
SiPixelTrackingRecHitsValid::meResXvsAlphaBarrelLayerModule
MonitorElement * meResXvsAlphaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:110
SiPixelTrackingRecHitsValid::meResYvsBetaZmPanel2DiskPlaq
MonitorElement * meResYvsBetaZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:256
SiPixelTrackingRecHitsValid::meNpixZmPanel2DiskPlaq
MonitorElement * meNpixZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:249
SiPixelTrackingRecHitsValid::meResXvsAlphaBarrelFlippedLaddersLayer
MonitorElement * meResXvsAlphaBarrelFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:199
SiPixelTrackingRecHitsValid::meResXvsBetaBarrelLayerModule
MonitorElement * meResXvsBetaBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:112
SiPixelTrackingRecHitsValid::mePullYvsPhiBarrelLayerModule
MonitorElement * mePullYvsPhiBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:119
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:304
SiPixelTrackingRecHitsValid::mePosyZmPanel1DiskPlaq
MonitorElement * mePosyZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:216
SiPixelTrackingRecHitsValid::mePullXvsPhiZpPanel1DiskPlaq
MonitorElement * mePullXvsPhiZpPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:288
SiPixelTrackingRecHitsValid::meNxpixBarrelLayerModule
MonitorElement * meNxpixBarrelLayerModule[3][8]
Definition: SiPixelTrackingRecHitsValid.h:107
SiPixelTrackingRecHitsValid::meNxpixZmPanel2DiskPlaq
MonitorElement * meNxpixZmPanel2DiskPlaq[2][3]
Definition: SiPixelTrackingRecHitsValid.h:250
SiPixelRecHit
Pixel Reconstructed Hit.
SiPixelTrackingRecHitsValid::meResxZmPanel1DiskPlaq
MonitorElement * meResxZmPanel1DiskPlaq[2][4]
Definition: SiPixelTrackingRecHitsValid.h:219
phi
Definition: DDAxes.h:10
edmNew::DetSet::begin
iterator begin()
Definition: DetSetNew.h:56
SiPixelTrackingRecHitsValid::meResYvsBetaBarrelNonFlippedLaddersLayer
MonitorElement * meResYvsBetaBarrelNonFlippedLaddersLayer[3]
Definition: SiPixelTrackingRecHitsValid.h:206
Generated for CMSSW Reference Manual by doxygen 1.8.5