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