CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
SiPixelHitEfficiencyModule.cc
Go to the documentation of this file.
1 // Package: SiPixelMonitorTrack
2 // Class: SiPixelHitEfficiencyModule
3 //
4 // class SiPixelHitEfficiencyModule SiPixelHitEfficiencyModule.cc
5 // DQM/SiPixelMonitorTrack/src/SiPixelHitEfficiencyModule.cc
6 //
7 // Description: SiPixel hit efficiency data quality monitoring modules
8 // Implementation: prototype -> improved -> never final - end of the 1st step
9 //
10 // Original Authors: Romain Rougny & Luca Mucibello
11 // Created: Mar Nov 10 13:29:00 CET nbinangle9
12 
13 
14 
15 #include <string>
16 #include <iostream>
17 
18 #include "DQMServices/Core/interface/DQMStore.h"
19 #include "FWCore/Framework/interface/ESHandle.h"
20 #include "FWCore/ParameterSet/interface/ParameterSet.h"
21 #include "FWCore/ServiceRegistry/interface/Service.h"
22 #include "FWCore/Utilities/interface/InputTag.h"
23 #include "DQM/SiPixelCommon/interface/SiPixelHistogramId.h"
24 #include "DQM/SiPixelMonitorTrack/interface/SiPixelHitEfficiencyModule.h"
25 
26 // Data Formats
27 #include "DataFormats/SiPixelDetId/interface/PixelBarrelName.h"
28 #include "DataFormats/SiPixelDetId/interface/PixelBarrelNameUpgrade.h"
29 #include "DataFormats/SiPixelDetId/interface/PixelEndcapName.h"
30 #include "DataFormats/SiPixelDetId/interface/PixelEndcapNameUpgrade.h"
31 #include "DataFormats/DetId/interface/DetId.h"
32 #include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
33 #include "DataFormats/VertexReco/interface/Vertex.h"
34 #include "DataFormats/VertexReco/interface/VertexFwd.h"
35 #include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
36 #include "Geometry/Records/interface/IdealGeometryRecord.h"
37 
38 
39 using namespace std;
40 
41 
42 SiPixelHitEfficiencyModule::SiPixelHitEfficiencyModule() : id_(0) {
43  bBookTracks = true;
44 }
45 
46 
47 SiPixelHitEfficiencyModule::SiPixelHitEfficiencyModule(uint32_t id) : id_(id) {
48  bBookTracks = true;
49 }
50 
51 
52 SiPixelHitEfficiencyModule::~SiPixelHitEfficiencyModule() {
53 
54 }
55 
56 
57 void SiPixelHitEfficiencyModule::book(const edm::ParameterSet& iConfig, edm::EventSetup const & iSetup, DQMStore::IBooker & iBooker,int type, bool isUpgrade) {
58 
59  edm::ESHandle<TrackerTopology> tTopoHandle;
60  iSetup.get<IdealGeometryRecord>().get(tTopoHandle);
61  const TrackerTopology *pTT = tTopoHandle.product();
62 
63  bool barrel = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
64  bool endcap = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
65  bool isHalfModule = false;
66  if(barrel){
67  isHalfModule = PixelBarrelName(DetId(id_),pTT,isUpgrade).isHalfModule();
68  }
69 
70  edm::InputTag src = iConfig.getParameter<edm::InputTag>("src");
71  debug_ = iConfig.getUntrackedParameter<bool>("debug",false);
72  updateEfficiencies = iConfig.getUntrackedParameter<bool>("updateEfficiencies",false);
73  std::string hisID;
74 
75  int nbinangle = 28;
76  int nbinX = 45;
77  int nbinY = 40;
78 
79  if(type==0){
80  SiPixelHistogramId* theHistogramId = new SiPixelHistogramId(src.label());
81 
82  if(updateEfficiencies){
83  //EFFICIENCY
84  hisID = theHistogramId->setHistoId("efficiency",id_);
85  meEfficiency_ = iBooker.book1D(hisID,"Hit efficiency",1,0,1.);
86  meEfficiency_->setAxisTitle("Hit efficiency",1);
87 
88  hisID = theHistogramId->setHistoId("efficiencyX",id_);
89  meEfficiencyX_ = iBooker.book1D(hisID,"Hit efficiency in X",nbinX,-1.5,1.5);
90  meEfficiencyX_->setAxisTitle("Hit efficiency in X",1);
91 
92  hisID = theHistogramId->setHistoId("efficiencyY",id_);
93  meEfficiencyY_ = iBooker.book1D(hisID,"Hit efficiency in Y",nbinY,-4.,4.);
94  meEfficiencyY_->setAxisTitle("Hit efficiency in Y",1);
95 
96  hisID = theHistogramId->setHistoId("efficiencyAlpha",id_);
97  meEfficiencyAlpha_ = iBooker.book1D(hisID,"Hit efficiency in Alpha",nbinangle,-3.5,3.5);
98  meEfficiencyAlpha_->setAxisTitle("Hit efficiency in Alpha",1);
99 
100  hisID = theHistogramId->setHistoId("efficiencyBeta",id_);
101  meEfficiencyBeta_ = iBooker.book1D(hisID,"Hit efficiency in Beta",nbinangle,-3.5,3.5);
102  meEfficiencyBeta_->setAxisTitle("Hit efficiency in Beta",1);
103  }
104 
105  //VALID
106  hisID = theHistogramId->setHistoId("valid",id_);
107  meValid_ = iBooker.book1D(hisID,"# Valid hits",1,0,1.);
108  meValid_->setAxisTitle("# Valid hits",1);
109 
110  //MISSING
111  hisID = theHistogramId->setHistoId("missing",id_);
112  meMissing_ = iBooker.book1D(hisID,"# Missing hits",1,0,1.);
113  meMissing_->setAxisTitle("# Missing hits",1);
114 
115  delete theHistogramId;
116  }
117 
118  if(type==1 && barrel){
119  uint32_t DBladder;
120  DBladder = PixelBarrelName(DetId(id_),pTT,isUpgrade).ladderName();
121  char sladder[80]; sprintf(sladder,"Ladder_%02i",DBladder);
122  hisID = src.label() + "_" + sladder;
123  if(isHalfModule) hisID += "H";
124  else hisID += "F";
125 
126  if(updateEfficiencies){
127  //EFFICIENCY
128  meEfficiencyLad_ = iBooker.book1D("efficiency_"+hisID,"Hit efficiency",1,0,1.);
129  meEfficiencyLad_->setAxisTitle("Hit efficiency",1);
130 
131  meEfficiencyXLad_ = iBooker.book1D("efficiencyX_"+hisID,"Hit efficiency in X",nbinX,-1.5,1.5);
132  meEfficiencyXLad_->setAxisTitle("Hit efficiency in X",1);
133 
134  meEfficiencyYLad_ = iBooker.book1D("efficiencyY_"+hisID,"Hit efficiency in Y",nbinY,-4.,4.);
135  meEfficiencyYLad_->setAxisTitle("Hit efficiency in Y",1);
136 
137  meEfficiencyAlphaLad_ = iBooker.book1D("efficiencyAlpha_"+hisID,"Hit efficiency in Alpha",nbinangle,-3.5,3.5);
138  meEfficiencyAlphaLad_->setAxisTitle("Hit efficiency in Alpha",1);
139 
140  meEfficiencyBetaLad_ = iBooker.book1D("efficiencyBeta_"+hisID,"Hit efficiency in Beta",nbinangle,-3.5,3.5);
141  meEfficiencyBetaLad_->setAxisTitle("Hit efficiency in Beta",1);
142  }
143 
144  //VALID
145  meValidLad_ = iBooker.book1D("valid_"+hisID,"# Valid hits",1,0,1.);
146  meValidLad_->setAxisTitle("# Valid hits",1);
147 
148  meValidXLad_ = iBooker.book1D("validX_"+hisID,"# Valid hits in X",nbinX,-1.5,1.5);
149  meValidXLad_->setAxisTitle("# Valid hits in X",1);
150 
151  meValidYLad_ = iBooker.book1D("validY_"+hisID,"# Valid hits in Y",nbinY,-4.,4.);
152  meValidYLad_->setAxisTitle("# Valid hits in Y",1);
153 
154  meValidModLad_ = iBooker.book1D("validMod_"+hisID,"# Valid hits on Module",20,1,21.);
155  meValidModLad_->setAxisTitle("# Valid hits on Module",1);
156 
157  meValidAlphaLad_ = iBooker.book1D("validAlpha_"+hisID,"# Valid hits in Alpha",nbinangle,-3.5,3.5);
158  meValidAlphaLad_->setAxisTitle("# Valid hits in Alpha",1);
159 
160  meValidBetaLad_ = iBooker.book1D("validBeta_"+hisID,"# Valid hits in Beta",nbinangle,-3.5,3.5);
161  meValidBetaLad_->setAxisTitle("# Valid hits in Beta",1);
162 
163  //MISSING
164  meMissingLad_ = iBooker.book1D("missing_"+hisID,"# Missing hits",1,0,1.);
165  meMissingLad_->setAxisTitle("# Missing hits",1);
166 
167  meMissingXLad_ = iBooker.book1D("missingX_"+hisID,"# Missing hits in X",nbinX,-1.5,1.5);
168  meMissingXLad_->setAxisTitle("# Missing hits in X",1);
169 
170  meMissingYLad_ = iBooker.book1D("missingY_"+hisID,"# Missing hits in Y",nbinY,-4.,4.);
171  meMissingYLad_->setAxisTitle("# Missing hits in Y",1);
172 
173  meMissingModLad_ = iBooker.book1D("missingMod_"+hisID,"# Missing hits on Module",20,1,21.);
174  meMissingModLad_->setAxisTitle("# Missing hits on Module",1);
175 
176  meMissingAlphaLad_ = iBooker.book1D("missingAlpha_"+hisID,"# Missing hits in Alpha",nbinangle,-3.5,3.5);
177  meMissingAlphaLad_->setAxisTitle("# Missing hits in Alpha",1);
178 
179  meMissingBetaLad_ = iBooker.book1D("missingBeta_"+hisID,"# Missing hits in Beta",nbinangle,-3.5,3.5);
180  meMissingBetaLad_->setAxisTitle("# Missing hits in Beta",1);
181  }
182 
183  if(type==2 && barrel){
184  uint32_t DBlayer;
185  DBlayer = PixelBarrelName(DetId(id_),pTT,isUpgrade).layerName();
186  char slayer[80]; sprintf(slayer,"Layer_%i",DBlayer);
187  hisID = src.label() + "_" + slayer;
188 
189 
190  if(updateEfficiencies){
191  //EFFICIENCY
192  meEfficiencyLay_ = iBooker.book1D("efficiency_"+hisID,"Hit efficiency",1,0,1.);
193  meEfficiencyLay_->setAxisTitle("Hit efficiency",1);
194 
195  meEfficiencyXLay_ = iBooker.book1D("efficiencyX_"+hisID,"Hit efficiency in X",nbinX,-1.5,1.5);
196  meEfficiencyXLay_->setAxisTitle("Hit efficiency in X",1);
197 
198  meEfficiencyYLay_ = iBooker.book1D("efficiencyY_"+hisID,"Hit efficiency in Y",nbinY,-4.,4.);
199  meEfficiencyYLay_->setAxisTitle("Hit efficiency in Y",1);
200 
201  meEfficiencyAlphaLay_ = iBooker.book1D("efficiencyAlpha_"+hisID,"Hit efficiency in Alpha",nbinangle,-3.5,3.5);
202  meEfficiencyAlphaLay_->setAxisTitle("Hit efficiency in Alpha",1);
203 
204  meEfficiencyBetaLay_ = iBooker.book1D("efficiencyBeta_"+hisID,"Hit efficiency in Beta",nbinangle,-3.5,3.5);
205  meEfficiencyBetaLay_->setAxisTitle("Hit efficiency in Beta",1);
206  }
207 
208  //VALID
209  meValidLay_ = iBooker.book1D("valid_"+hisID,"# Valid hits",1,0,1.);
210  meValidLay_->setAxisTitle("# Valid hits",1);
211 
212  meValidXLay_ = iBooker.book1D("validX_"+hisID,"# Valid hits in X",nbinX,-1.5,1.5);
213  meValidXLay_->setAxisTitle("# Valid hits in X",1);
214 
215  meValidYLay_ = iBooker.book1D("validY_"+hisID,"# Valid hits in Y",nbinY,-4.,4.);
216  meValidYLay_->setAxisTitle("# Valid hits in Y",1);
217 
218  meValidAlphaLay_ = iBooker.book1D("validAlpha_"+hisID,"# Valid hits in Alpha",nbinangle,-3.5,3.5);
219  meValidAlphaLay_->setAxisTitle("# Valid hits in Alpha",1);
220 
221  meValidBetaLay_ = iBooker.book1D("validBeta_"+hisID,"# Valid hits in Beta",nbinangle,-3.5,3.5);
222  meValidBetaLay_->setAxisTitle("# Valid hits in Beta",1);
223 
224  //MISSING
225  meMissingLay_ = iBooker.book1D("missing_"+hisID,"# Missing hits",1,0,1.);
226  meMissingLay_->setAxisTitle("# Missing hits",1);
227 
228  meMissingXLay_ = iBooker.book1D("missingX_"+hisID,"# Missing hits in X",nbinX,-1.5,1.5);
229  meMissingXLay_->setAxisTitle("# Missing hits in X",1);
230 
231  meMissingYLay_ = iBooker.book1D("missingY_"+hisID,"# Missing hits in Y",nbinY,-4.,4.);
232  meMissingYLay_->setAxisTitle("# Missing hits in Y",1);
233 
234  meMissingAlphaLay_ = iBooker.book1D("missingAlpha_"+hisID,"# Missing hits in Alpha",nbinangle,-3.5,3.5);
235  meMissingAlphaLay_->setAxisTitle("# Missing hits in Alpha",1);
236 
237  meMissingBetaLay_ = iBooker.book1D("missingBeta_"+hisID,"# Missing hits in Beta",nbinangle,-3.5,3.5);
238  meMissingBetaLay_->setAxisTitle("# Missing hits in Beta",1);
239  }
240 
241  if(type==3 && barrel){
242  uint32_t DBmodule;
243  DBmodule = PixelBarrelName(DetId(id_),pTT,isUpgrade).moduleName();
244  char smodule[80]; sprintf(smodule,"Ring_%i",DBmodule);
245  hisID = src.label() + "_" + smodule;
246 
247  if(updateEfficiencies){
248  //EFFICIENCY
249  meEfficiencyPhi_ = iBooker.book1D("efficiency_"+hisID,"Hit efficiency",1,0,1.);
250  meEfficiencyPhi_->setAxisTitle("Hit efficiency",1);
251 
252  meEfficiencyXPhi_ = iBooker.book1D("efficiencyX_"+hisID,"Hit efficiency in X",nbinX,-1.5,1.5);
253  meEfficiencyXPhi_->setAxisTitle("Hit efficiency in X",1);
254 
255  meEfficiencyYPhi_ = iBooker.book1D("efficiencyY_"+hisID,"Hit efficiency in Y",nbinY,-4.,4.);
256  meEfficiencyYPhi_->setAxisTitle("Hit efficiency in Y",1);
257 
258  meEfficiencyAlphaPhi_ = iBooker.book1D("efficiencyAlpha_"+hisID,"Hit efficiency in Alpha",nbinangle,-3.5,3.5);
259  meEfficiencyAlphaPhi_->setAxisTitle("Hit efficiency in Alpha",1);
260 
261  meEfficiencyBetaPhi_ = iBooker.book1D("efficiencyBeta_"+hisID,"Hit efficiency in Beta",nbinangle,-3.5,3.5);
262  meEfficiencyBetaPhi_->setAxisTitle("Hit efficiency in Beta",1);
263  }
264 
265  //VALID
266  meValidPhi_ = iBooker.book1D("valid_"+hisID,"# Valid hits",1,0,1.);
267  meValidPhi_->setAxisTitle("# Valid hits",1);
268 
269  meValidXPhi_ = iBooker.book1D("validX_"+hisID,"# Valid hits in X",nbinX,-1.5,1.5);
270  meValidXPhi_->setAxisTitle("# Valid hits in X",1);
271 
272  meValidYPhi_ = iBooker.book1D("validY_"+hisID,"# Valid hits in Y",nbinY,-4.,4.);
273  meValidYPhi_->setAxisTitle("# Valid hits in Y",1);
274 
275  meValidAlphaPhi_ = iBooker.book1D("validAlpha_"+hisID,"# Valid hits in Alpha",nbinangle,-3.5,3.5);
276  meValidAlphaPhi_->setAxisTitle("# Valid hits in Alpha",1);
277 
278  meValidBetaPhi_ = iBooker.book1D("validBeta_"+hisID,"# Valid hits in Beta",nbinangle,-3.5,3.5);
279  meValidBetaPhi_->setAxisTitle("# Valid hits in Beta",1);
280 
281  //MISSING
282  meMissingPhi_ = iBooker.book1D("missing_"+hisID,"# Missing hits",1,0,1.);
283  meMissingPhi_->setAxisTitle("# Missing hits",1);
284 
285  meMissingXPhi_ = iBooker.book1D("missingX_"+hisID,"# Missing hits in X",nbinX,-1.5,1.5);
286  meMissingXPhi_->setAxisTitle("# Missing hits in X",1);
287 
288  meMissingYPhi_ = iBooker.book1D("missingY_"+hisID,"# Missing hits in Y",nbinY,-4.,4.);
289  meMissingYPhi_->setAxisTitle("# Missing hits in Y",1);
290 
291  meMissingAlphaPhi_ = iBooker.book1D("missingAlpha_"+hisID,"# Missing hits in Alpha",nbinangle,-3.5,3.5);
292  meMissingAlphaPhi_->setAxisTitle("# Missing hits in Alpha",1);
293 
294  meMissingBetaPhi_ = iBooker.book1D("missingBeta_"+hisID,"# Missing hits in Beta",nbinangle,-3.5,3.5);
295  meMissingBetaPhi_->setAxisTitle("# Missing hits in Beta",1);
296  }
297 
298  if(type==4 && endcap){
299  uint32_t blade;
300  blade= PixelEndcapName(DetId(id_),pTT,isUpgrade).bladeName();
301 
302  char sblade[80]; sprintf(sblade, "Blade_%02i",blade);
303  hisID = src.label() + "_" + sblade;
304 
305  if(updateEfficiencies){
306  //EFFICIENCY
307  meEfficiencyBlade_ = iBooker.book1D("efficiency_"+hisID,"Hit efficiency",1,0,1.);
308  meEfficiencyBlade_->setAxisTitle("Hit efficiency",1);
309 
310  meEfficiencyXBlade_ = iBooker.book1D("efficiencyX_"+hisID,"Hit efficiency in X",nbinX,-1.5,1.5);
311  meEfficiencyXBlade_->setAxisTitle("Hit efficiency in X",1);
312 
313  meEfficiencyYBlade_ = iBooker.book1D("efficiencyY_"+hisID,"Hit efficiency in Y",nbinY,-4.,4.);
314  meEfficiencyYBlade_->setAxisTitle("Hit efficiency in Y",1);
315 
316  meEfficiencyAlphaBlade_ = iBooker.book1D("efficiencyAlpha_"+hisID,"Hit efficiency in Alpha",nbinangle,-3.5,3.5);
317  meEfficiencyAlphaBlade_->setAxisTitle("Hit efficiency in Alpha",1);
318 
319  meEfficiencyBetaBlade_ = iBooker.book1D("efficiencyBeta_"+hisID,"Hit efficiency in Beta",nbinangle,-3.5,3.5);
320  meEfficiencyBetaBlade_->setAxisTitle("Hit efficiency in Beta",1);
321  }
322 
323  //VALID
324  meValidBlade_ = iBooker.book1D("valid_"+hisID,"# Valid hits",1,0,1.);
325  meValidBlade_->setAxisTitle("# Valid hits",1);
326 
327  meValidXBlade_ = iBooker.book1D("validX_"+hisID,"# Valid hits in X",nbinX,-1.5,1.5);
328  meValidXBlade_->setAxisTitle("# Valid hits in X",1);
329 
330  meValidYBlade_ = iBooker.book1D("validY_"+hisID,"# Valid hits in Y",nbinY,-4.,4.);
331  meValidYBlade_->setAxisTitle("# Valid hits in Y",1);
332 
333  meValidAlphaBlade_ = iBooker.book1D("validAlpha_"+hisID,"# Valid hits in Alpha",nbinangle,-3.5,3.5);
334  meValidAlphaBlade_->setAxisTitle("# Valid hits in Alpha",1);
335 
336  meValidBetaBlade_ = iBooker.book1D("validBeta_"+hisID,"# Valid hits in Beta",nbinangle,-3.5,3.5);
337  meValidBetaBlade_->setAxisTitle("# Valid hits in Beta",1);
338 
339  //MISSING
340  meMissingBlade_ = iBooker.book1D("missing_"+hisID,"# Missing hits",1,0,1.);
341  meMissingBlade_->setAxisTitle("# Missing hits",1);
342 
343  meMissingXBlade_ = iBooker.book1D("missingX_"+hisID,"# Missing hits in X",nbinX,-1.5,1.5);
344  meMissingXBlade_->setAxisTitle("# Missing hits in X",1);
345 
346  meMissingYBlade_ = iBooker.book1D("missingY_"+hisID,"# Missing hits in Y",nbinY,-4.,4.);
347  meMissingYBlade_->setAxisTitle("# Missing hits in Y",1);
348 
349  meMissingAlphaBlade_ = iBooker.book1D("missingAlpha_"+hisID,"# Missing hits in Alpha",nbinangle,-3.5,3.5);
350  meMissingAlphaBlade_->setAxisTitle("# Missing hits in Alpha",1);
351 
352  meMissingBetaBlade_ = iBooker.book1D("missingBeta_"+hisID,"# Missing hits in Beta",nbinangle,-3.5,3.5);
353  meMissingBetaBlade_->setAxisTitle("# Missing hits in Beta",1);
354  }
355 
356  if(type==5 && endcap){
357  uint32_t disk;
358  disk = PixelEndcapName(DetId(id_),pTT,isUpgrade).diskName();
359 
360  char sdisk[80]; sprintf(sdisk, "Disk_%i",disk);
361  hisID = src.label() + "_" + sdisk;
362 
363  if(updateEfficiencies){
364  //EFFICIENCY
365  meEfficiencyDisk_ = iBooker.book1D("efficiency_"+hisID,"Hit efficiency",1,0,1.);
366  meEfficiencyDisk_->setAxisTitle("Hit efficiency",1);
367 
368  meEfficiencyXDisk_ = iBooker.book1D("efficiencyX_"+hisID,"Hit efficiency in X",nbinX,-1.5,1.5);
369  meEfficiencyXDisk_->setAxisTitle("Hit efficiency in X",1);
370 
371  meEfficiencyYDisk_ = iBooker.book1D("efficiencyY_"+hisID,"Hit efficiency in Y",nbinY,-4.,4.);
372  meEfficiencyYDisk_->setAxisTitle("Hit efficiency in Y",1);
373 
374  meEfficiencyAlphaDisk_ = iBooker.book1D("efficiencyAlpha_"+hisID,"Hit efficiency in Alpha",nbinangle,-3.5,3.5);
375  meEfficiencyAlphaDisk_->setAxisTitle("Hit efficiency in Alpha",1);
376 
377  meEfficiencyBetaDisk_ = iBooker.book1D("efficiencyBeta_"+hisID,"Hit efficiency in Beta",nbinangle,-3.5,3.5);
378  meEfficiencyBetaDisk_->setAxisTitle("Hit efficiency in Beta",1);
379  }
380 
381  //VALID
382  meValidDisk_ = iBooker.book1D("valid_"+hisID,"# Valid hits",1,0,1.);
383  meValidDisk_->setAxisTitle("# Valid hits",1);
384 
385  meValidXDisk_ = iBooker.book1D("validX_"+hisID,"# Valid hits in X",nbinX,-1.5,1.5);
386  meValidXDisk_->setAxisTitle("# Valid hits in X",1);
387 
388  meValidYDisk_ = iBooker.book1D("validY_"+hisID,"# Valid hits in Y",nbinY,-4.,4.);
389  meValidYDisk_->setAxisTitle("# Valid hits in Y",1);
390 
391  meValidAlphaDisk_ = iBooker.book1D("validAlpha_"+hisID,"# Valid hits in Alpha",nbinangle,-3.5,3.5);
392  meValidAlphaDisk_->setAxisTitle("# Valid hits in Alpha",1);
393 
394  meValidBetaDisk_ = iBooker.book1D("validBeta_"+hisID,"# Valid hits in Beta",nbinangle,-3.5,3.5);
395  meValidBetaDisk_->setAxisTitle("# Valid hits in Beta",1);
396 
397  //MISSING
398  meMissingDisk_ = iBooker.book1D("missing_"+hisID,"# Missing hits",1,0,1.);
399  meMissingDisk_->setAxisTitle("# Missing hits",1);
400 
401  meMissingXDisk_ = iBooker.book1D("missingX_"+hisID,"# Missing hits in X",nbinX,-1.5,1.5);
402  meMissingXDisk_->setAxisTitle("# Missing hits in X",1);
403 
404  meMissingYDisk_ = iBooker.book1D("missingY_"+hisID,"# Missing hits in Y",nbinY,-4.,4.);
405  meMissingYDisk_->setAxisTitle("# Missing hits in Y",1);
406 
407  meMissingAlphaDisk_ = iBooker.book1D("missingAlpha_"+hisID,"# Missing hits in Alpha",nbinangle,-3.5,3.5);
408  meMissingAlphaDisk_->setAxisTitle("# Missing hits in Alpha",1);
409 
410  meMissingBetaDisk_ = iBooker.book1D("missingBeta_"+hisID,"# Missing hits in Beta",nbinangle,-3.5,3.5);
411  meMissingBetaDisk_->setAxisTitle("# Missing hits in Beta",1);
412  }
413 
414 
415  if(type==6 && endcap){
416  uint32_t panel;
417  uint32_t module;
418  panel= PixelEndcapName(DetId(id_),pTT,isUpgrade).pannelName();
419  module= PixelEndcapName(DetId(id_),pTT,isUpgrade).plaquetteName();
420 
421  char slab[80]; sprintf(slab, "Panel_%i_Ring_%i",panel, module);
422  hisID = src.label() + "_" + slab;
423 
424  if(updateEfficiencies){
425  //EFFICIENCY
426  meEfficiencyRing_ = iBooker.book1D("efficiency_"+hisID,"Hit efficiency",1,0,1.);
427  meEfficiencyRing_->setAxisTitle("Hit efficiency",1);
428 
429  meEfficiencyXRing_ = iBooker.book1D("efficiencyX_"+hisID,"Hit efficiency in X",nbinX,-1.5,1.5);
430  meEfficiencyXRing_->setAxisTitle("Hit efficiency in X",1);
431 
432  meEfficiencyYRing_ = iBooker.book1D("efficiencyY_"+hisID,"Hit efficiency in Y",nbinY,-4.,4.);
433  meEfficiencyYRing_->setAxisTitle("Hit efficiency in Y",1);
434 
435  meEfficiencyAlphaRing_ = iBooker.book1D("efficiencyAlpha_"+hisID,"Hit efficiency in Alpha",nbinangle,-3.5,3.5);
436  meEfficiencyAlphaRing_->setAxisTitle("Hit efficiency in Alpha",1);
437 
438  meEfficiencyBetaRing_ = iBooker.book1D("efficiencyBeta_"+hisID,"Hit efficiency in Beta",nbinangle,-3.5,3.5);
439  meEfficiencyBetaRing_->setAxisTitle("Hit efficiency in Beta",1);
440  }
441 
442  //VALID
443  meValidRing_ = iBooker.book1D("valid_"+hisID,"# Valid hits",1,0,1.);
444  meValidRing_->setAxisTitle("# Valid hits",1);
445 
446  meValidXRing_ = iBooker.book1D("validX_"+hisID,"# Valid hits in X",nbinX,-1.5,1.5);
447  meValidXRing_->setAxisTitle("# Valid hits in X",1);
448 
449  meValidYRing_ = iBooker.book1D("validY_"+hisID,"# Valid hits in Y",nbinY,-4.,4.);
450  meValidYRing_->setAxisTitle("# Valid hits in Y",1);
451 
452  meValidAlphaRing_ = iBooker.book1D("validAlpha_"+hisID,"# Valid hits in Alpha",nbinangle,-3.5,3.5);
453  meValidAlphaRing_->setAxisTitle("# Valid hits in Alpha",1);
454 
455  meValidBetaRing_ = iBooker.book1D("validBeta_"+hisID,"# Valid hits in Beta",nbinangle,-3.5,3.5);
456  meValidBetaRing_->setAxisTitle("# Valid hits in Beta",1);
457 
458  //MISSING
459  meMissingRing_ = iBooker.book1D("missing_"+hisID,"# Missing hits",1,0,1.);
460  meMissingRing_->setAxisTitle("# Missing hits",1);
461 
462  meMissingXRing_ = iBooker.book1D("missingX_"+hisID,"# Missing hits in X",nbinX,-1.5,1.5);
463  meMissingXRing_->setAxisTitle("# Missing hits in X",1);
464 
465  meMissingYRing_ = iBooker.book1D("missingY_"+hisID,"# Missing hits in Y",nbinY,-4.,4.);
466  meMissingYRing_->setAxisTitle("# Missing hits in Y",1);
467 
468  meMissingAlphaRing_ = iBooker.book1D("missingAlpha_"+hisID,"# Missing hits in Alpha",nbinangle,-3.5,3.5);
469  meMissingAlphaRing_->setAxisTitle("# Missing hits in Alpha",1);
470 
471  meMissingBetaRing_ = iBooker.book1D("missingBeta_"+hisID,"# Missing hits in Beta",nbinangle,-3.5,3.5);
472  meMissingBetaRing_->setAxisTitle("# Missing hits in Beta",1);
473  }
474 }
475 
476 
477 void SiPixelHitEfficiencyModule::fill(const LocalTrajectoryParameters& ltp, bool isHitValid, bool modon, bool ladon, bool layon, bool phion, bool bladeon, bool diskon, bool ringon) {
478 
479  bool barrel = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
480  bool endcap = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
481 
482  LocalVector localDir = ltp.momentum()/ltp.momentum().mag();
483  float prediction_alpha = atan2(localDir.z(), localDir.x());
484  float prediction_beta = atan2(localDir.z(), localDir.y());
485  float prediction_x = ltp.position().x();
486  float prediction_y = ltp.position().y();
487  //CS - this will probably break with isUpgrade
488  int imod = PXBDetId(DetId(id_)).module();
489 
490  if(isHitValid){
491  if(modon){
492  meValid_->Fill(0.5);
493  }
494  if(barrel && ladon){
495  meValidLad_->Fill(0.5);
496  meValidXLad_->Fill(prediction_x);
497  meValidYLad_->Fill(prediction_y);
498  meValidModLad_->Fill(imod);
499  meValidAlphaLad_->Fill(prediction_alpha);
500  meValidBetaLad_->Fill(prediction_beta);
501  }
502  if(barrel && layon){
503  meValidLay_->Fill(0.5);
504  meValidXLay_->Fill(prediction_x);
505  meValidYLay_->Fill(prediction_y);
506  meValidAlphaLay_->Fill(prediction_alpha);
507  meValidBetaLay_->Fill(prediction_beta);
508  }
509  if(barrel && phion){
510  meValidPhi_->Fill(0.5);
511  meValidXPhi_->Fill(prediction_x);
512  meValidYPhi_->Fill(prediction_y);
513  meValidAlphaPhi_->Fill(prediction_alpha);
514  meValidBetaPhi_->Fill(prediction_beta);
515  }
516  if(endcap && bladeon){
517  meValidBlade_->Fill(0.5);
518  meValidXBlade_->Fill(prediction_x);
519  meValidYBlade_->Fill(prediction_y);
520  meValidAlphaBlade_->Fill(prediction_alpha);
521  meValidBetaBlade_->Fill(prediction_beta);
522  }
523  if(endcap && diskon){
524  meValidDisk_->Fill(0.5);
525  meValidXDisk_->Fill(prediction_x);
526  meValidYDisk_->Fill(prediction_y);
527  meValidAlphaDisk_->Fill(prediction_alpha);
528  meValidBetaDisk_->Fill(prediction_beta);
529  }
530  if(endcap && ringon){
531  meValidRing_->Fill(0.5);
532  meValidXRing_->Fill(prediction_x);
533  meValidYRing_->Fill(prediction_y);
534  meValidAlphaRing_->Fill(prediction_alpha);
535  meValidBetaRing_->Fill(prediction_beta);
536  }
537  }
538  else {
539  if(modon){
540  meMissing_->Fill(0.5);
541  }
542  if(barrel && ladon){
543  meMissingLad_->Fill(0.5);
544  meMissingXLad_->Fill(prediction_x);
545  meMissingYLad_->Fill(prediction_y);
546  meMissingModLad_->Fill(imod);
547  meMissingAlphaLad_->Fill(prediction_alpha);
548  meMissingBetaLad_->Fill(prediction_beta);
549  }
550  if(barrel && layon){
551  meMissingLay_->Fill(0.5);
552  meMissingXLay_->Fill(prediction_x);
553  meMissingYLay_->Fill(prediction_y);
554  meMissingAlphaLay_->Fill(prediction_alpha);
555  meMissingBetaLay_->Fill(prediction_beta);
556  }
557  if(barrel && phion){
558  meMissingPhi_->Fill(0.5);
559  meMissingXPhi_->Fill(prediction_x);
560  meMissingYPhi_->Fill(prediction_y);
561  meMissingAlphaPhi_->Fill(prediction_alpha);
562  meMissingBetaPhi_->Fill(prediction_beta);
563  }
564  if(endcap && bladeon){
565  meMissingBlade_->Fill(0.5);
566  meMissingXBlade_->Fill(prediction_x);
567  meMissingYBlade_->Fill(prediction_y);
568  meMissingAlphaBlade_->Fill(prediction_alpha);
569  meMissingBetaBlade_->Fill(prediction_beta);
570  }
571  if(endcap && diskon){
572  meMissingDisk_->Fill(0.5);
573  meMissingXDisk_->Fill(prediction_x);
574  meMissingYDisk_->Fill(prediction_y);
575  meMissingAlphaDisk_->Fill(prediction_alpha);
576  meMissingBetaDisk_->Fill(prediction_beta);
577  }
578  if(endcap && ringon){
579  meMissingRing_->Fill(0.5);
580  meMissingXRing_->Fill(prediction_x);
581  meMissingYRing_->Fill(prediction_y);
582  meMissingAlphaRing_->Fill(prediction_alpha);
583  meMissingBetaRing_->Fill(prediction_beta);
584  }
585  }
586 
587  if(updateEfficiencies)
588  computeEfficiencies(modon, ladon, layon, phion, bladeon, diskon, ringon);
589 }
590 
591 void SiPixelHitEfficiencyModule::computeEfficiencies(bool modon, bool ladon, bool layon, bool phion, bool bladeon, bool diskon, bool ringon) {
592 
593  if(debug_)
594  std::cout<<"Now Filling histos for detid "<<id_<<std::endl;
595 
596  bool barrel = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
597  bool endcap = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
598 
599  if(modon){
600  meEfficiency_->setBinContent(1,(eff(meValid_->getBinContent(1),meMissing_->getBinContent(1))).first);
601  meEfficiency_->setBinError(1,(eff(meValid_->getBinContent(1),meMissing_->getBinContent(1))).second);
602  }
603  if(ladon && barrel){
604  meEfficiencyLad_->setBinContent(1,(eff(meValidLad_->getBinContent(1),meMissingLad_->getBinContent(1))).first);
605  meEfficiencyLad_->setBinError(1,(eff(meValidLad_->getBinContent(1),meMissingLad_->getBinContent(1))).second);
606  for(int i=1;i<=meValidXLad_->getNbinsX();++i){
607  meEfficiencyXLad_->setBinContent(i,(eff(meValidXLad_->getBinContent(i),meMissingXLad_->getBinContent(i))).first);
608  meEfficiencyXLad_->setBinError(i,(eff(meValidXLad_->getBinContent(i),meMissingXLad_->getBinContent(i))).second);
609  }
610  for(int i=1;i<=meValidYLad_->getNbinsX();++i){
611  meEfficiencyYLad_->setBinContent(i,(eff(meValidYLad_->getBinContent(i),meMissingYLad_->getBinContent(i))).first);
612  meEfficiencyYLad_->setBinError(i,(eff(meValidYLad_->getBinContent(i),meMissingYLad_->getBinContent(i))).second);
613  }
614  for(int i=1;i<=meValidAlphaLad_->getNbinsX();++i){
615  meEfficiencyAlphaLad_->setBinContent(i,(eff(meValidAlphaLad_->getBinContent(i),meMissingAlphaLad_->getBinContent(i))).first);
616  meEfficiencyAlphaLad_->setBinError(i,(eff(meValidAlphaLad_->getBinContent(i),meMissingAlphaLad_->getBinContent(i))).second);
617  }
618  for(int i=1;i<=meValidBetaLad_->getNbinsX();++i){
619  meEfficiencyBetaLad_->setBinContent(i,(eff(meValidBetaLad_->getBinContent(i),meMissingBetaLad_->getBinContent(i))).first);
620  meEfficiencyBetaLad_->setBinError(i,(eff(meValidBetaLad_->getBinContent(i),meMissingBetaLad_->getBinContent(i))).second);
621  }
622  }
623 
624  if(layon && barrel){
625  meEfficiencyLay_->setBinContent(1,(eff(meValidLay_->getBinContent(1),meMissingLay_->getBinContent(1))).first);
626  meEfficiencyLay_->setBinError(1,(eff(meValidLay_->getBinContent(1),meMissingLay_->getBinContent(1))).second);
627  for(int i=1;i<=meValidXLay_->getNbinsX();++i){
628  meEfficiencyXLay_->setBinContent(i,(eff(meValidXLay_->getBinContent(i),meMissingXLay_->getBinContent(i))).first);
629  meEfficiencyXLay_->setBinError(i,(eff(meValidXLay_->getBinContent(i),meMissingXLay_->getBinContent(i))).second);
630  }
631  for(int i=1;i<=meValidYLay_->getNbinsX();++i){
632  meEfficiencyYLay_->setBinContent(i,(eff(meValidYLay_->getBinContent(i),meMissingYLay_->getBinContent(i))).first);
633  meEfficiencyYLay_->setBinError(i,(eff(meValidYLay_->getBinContent(i),meMissingYLay_->getBinContent(i))).second);
634  }
635  for(int i=1;i<=meValidAlphaLay_->getNbinsX();++i){
636  meEfficiencyAlphaLay_->setBinContent(i,(eff(meValidAlphaLay_->getBinContent(i),meMissingAlphaLay_->getBinContent(i))).first);
637  meEfficiencyAlphaLay_->setBinError(i,(eff(meValidAlphaLay_->getBinContent(i),meMissingAlphaLay_->getBinContent(i))).second);
638  }
639  for(int i=1;i<=meValidBetaLay_->getNbinsX();++i){
640  meEfficiencyBetaLay_->setBinContent(i,(eff(meValidBetaLay_->getBinContent(i),meMissingBetaLay_->getBinContent(i))).first);
641  meEfficiencyBetaLay_->setBinError(i,(eff(meValidBetaLay_->getBinContent(i),meMissingBetaLay_->getBinContent(i))).second);
642  }
643  }
644 
645  if(phion && barrel){
646  meEfficiencyPhi_->setBinContent(1,(eff(meValidPhi_->getBinContent(1),meMissingPhi_->getBinContent(1))).first);
647  meEfficiencyPhi_->setBinError(1,(eff(meValidPhi_->getBinContent(1),meMissingPhi_->getBinContent(1))).second);
648  for(int i=1;i<=meValidXPhi_->getNbinsX();++i){
649  meEfficiencyXPhi_->setBinContent(i,(eff(meValidXPhi_->getBinContent(i),meMissingXPhi_->getBinContent(i))).first);
650  meEfficiencyXPhi_->setBinError(i,(eff(meValidXPhi_->getBinContent(i),meMissingXPhi_->getBinContent(i))).second);
651  }
652  for(int i=1;i<=meValidYPhi_->getNbinsX();++i){
653  meEfficiencyYPhi_->setBinContent(i,(eff(meValidYPhi_->getBinContent(i),meMissingYPhi_->getBinContent(i))).first);
654  meEfficiencyYPhi_->setBinError(i,(eff(meValidYPhi_->getBinContent(i),meMissingYPhi_->getBinContent(i))).second);
655  }
656  for(int i=1;i<=meValidAlphaPhi_->getNbinsX();++i){
657  meEfficiencyAlphaPhi_->setBinContent(i,(eff(meValidAlphaPhi_->getBinContent(i),meMissingAlphaPhi_->getBinContent(i))).first);
658  meEfficiencyAlphaPhi_->setBinError(i,(eff(meValidAlphaPhi_->getBinContent(i),meMissingAlphaPhi_->getBinContent(i))).second);
659  }
660  for(int i=1;i<=meValidBetaPhi_->getNbinsX();++i){
661  meEfficiencyBetaPhi_->setBinContent(i,(eff(meValidBetaPhi_->getBinContent(i),meMissingBetaPhi_->getBinContent(i))).first);
662  meEfficiencyBetaPhi_->setBinError(i,(eff(meValidBetaPhi_->getBinContent(i),meMissingBetaPhi_->getBinContent(i))).second);
663  }
664  }
665  if(bladeon && endcap){
666  meEfficiencyBlade_->setBinContent(1,(eff(meValidBlade_->getBinContent(1),meMissingBlade_->getBinContent(1))).first);
667  meEfficiencyBlade_->setBinError(1,(eff(meValidBlade_->getBinContent(1),meMissingBlade_->getBinContent(1))).second);
668  for(int i=1;i<=meValidXBlade_->getNbinsX();++i){
669  meEfficiencyXBlade_->setBinContent(i,(eff(meValidXBlade_->getBinContent(i),meMissingXBlade_->getBinContent(i))).first);
670  meEfficiencyXBlade_->setBinError(i,(eff(meValidXBlade_->getBinContent(i),meMissingXBlade_->getBinContent(i))).second);
671  }
672  for(int i=1;i<=meValidYBlade_->getNbinsX();++i){
673  meEfficiencyYBlade_->setBinContent(i,(eff(meValidYBlade_->getBinContent(i),meMissingYBlade_->getBinContent(i))).first);
674  meEfficiencyYBlade_->setBinError(i,(eff(meValidYBlade_->getBinContent(i),meMissingYBlade_->getBinContent(i))).second);
675  }
676  for(int i=1;i<=meValidAlphaBlade_->getNbinsX();++i){
677  meEfficiencyAlphaBlade_->setBinContent(i,(eff(meValidAlphaBlade_->getBinContent(i),meMissingAlphaBlade_->getBinContent(i))).first);
678  meEfficiencyAlphaBlade_->setBinError(i,(eff(meValidAlphaBlade_->getBinContent(i),meMissingAlphaBlade_->getBinContent(i))).second);
679  }
680  for(int i=1;i<=meValidBetaBlade_->getNbinsX();++i){
681  meEfficiencyBetaBlade_->setBinContent(i,(eff(meValidBetaBlade_->getBinContent(i),meMissingBetaBlade_->getBinContent(i))).first);
682  meEfficiencyBetaBlade_->setBinError(i,(eff(meValidBetaBlade_->getBinContent(i),meMissingBetaBlade_->getBinContent(i))).second);
683  }
684  }
685  if(diskon && endcap){
686  meEfficiencyDisk_->setBinContent(1,(eff(meValidDisk_->getBinContent(1),meMissingDisk_->getBinContent(1))).first);
687  meEfficiencyDisk_->setBinError(1,(eff(meValidDisk_->getBinContent(1),meMissingDisk_->getBinContent(1))).second);
688  for(int i=1;i<=meValidXDisk_->getNbinsX();++i){
689  meEfficiencyXDisk_->setBinContent(i,(eff(meValidXDisk_->getBinContent(i),meMissingXDisk_->getBinContent(i))).first);
690  meEfficiencyXDisk_->setBinError(i,(eff(meValidXDisk_->getBinContent(i),meMissingXDisk_->getBinContent(i))).second);
691  }
692  for(int i=1;i<=meValidYDisk_->getNbinsX();++i){
693  meEfficiencyYDisk_->setBinContent(i,(eff(meValidYDisk_->getBinContent(i),meMissingYDisk_->getBinContent(i))).first);
694  meEfficiencyYDisk_->setBinError(i,(eff(meValidYDisk_->getBinContent(i),meMissingYDisk_->getBinContent(i))).second);
695  }
696  for(int i=1;i<=meValidAlphaDisk_->getNbinsX();++i){
697  meEfficiencyAlphaDisk_->setBinContent(i,(eff(meValidAlphaDisk_->getBinContent(i),meMissingAlphaDisk_->getBinContent(i))).first);
698  meEfficiencyAlphaDisk_->setBinError(i,(eff(meValidAlphaDisk_->getBinContent(i),meMissingAlphaDisk_->getBinContent(i))).second);
699  }
700  for(int i=1;i<=meValidBetaDisk_->getNbinsX();++i){
701  meEfficiencyBetaDisk_->setBinContent(i,(eff(meValidBetaDisk_->getBinContent(i),meMissingBetaDisk_->getBinContent(i))).first);
702  meEfficiencyBetaDisk_->setBinError(i,(eff(meValidBetaDisk_->getBinContent(i),meMissingBetaDisk_->getBinContent(i))).second);
703  }
704  }
705  if(ringon && endcap){
706  meEfficiencyRing_->setBinContent(1,(eff(meValidRing_->getBinContent(1),meMissingRing_->getBinContent(1))).first);
707  meEfficiencyRing_->setBinError(1,(eff(meValidRing_->getBinContent(1),meMissingRing_->getBinContent(1))).second);
708  for(int i=1;i<=meValidXRing_->getNbinsX();++i){
709  meEfficiencyXRing_->setBinContent(i,(eff(meValidXRing_->getBinContent(i),meMissingXRing_->getBinContent(i))).first);
710  meEfficiencyXRing_->setBinError(i,(eff(meValidXRing_->getBinContent(i),meMissingXRing_->getBinContent(i))).second);
711  }
712  for(int i=1;i<=meValidYRing_->getNbinsX();++i){
713  meEfficiencyYRing_->setBinContent(i,(eff(meValidYRing_->getBinContent(i),meMissingYRing_->getBinContent(i))).first);
714  meEfficiencyYRing_->setBinError(i,(eff(meValidYRing_->getBinContent(i),meMissingYRing_->getBinContent(i))).second);
715  }
716  for(int i=1;i<=meValidAlphaRing_->getNbinsX();++i){
717  meEfficiencyAlphaRing_->setBinContent(i,(eff(meValidAlphaRing_->getBinContent(i),meMissingAlphaRing_->getBinContent(i))).first);
718  meEfficiencyAlphaRing_->setBinError(i,(eff(meValidAlphaRing_->getBinContent(i),meMissingAlphaRing_->getBinContent(i))).second);
719  }
720  for(int i=1;i<=meValidBetaRing_->getNbinsX();++i){
721  meEfficiencyBetaRing_->setBinContent(i,(eff(meValidBetaRing_->getBinContent(i),meMissingBetaRing_->getBinContent(i))).first);
722  meEfficiencyBetaRing_->setBinError(i,(eff(meValidBetaRing_->getBinContent(i),meMissingBetaRing_->getBinContent(i))).second);
723  }
724  }
725 
726 }
727 
728 std::pair<double,double> SiPixelHitEfficiencyModule::eff(double nValid, double nMissing){
729  double efficiency = 0, error = 0 ;
730  if(nValid+nMissing!=0){
731  efficiency=nValid/(nValid+nMissing);
732  error=sqrt(efficiency*(1.-efficiency)/(nValid+nMissing));
733  }
734  return make_pair(efficiency,error);
735 }
PixelEndcapName::plaquetteName
int plaquetteName() const
plaquetteId (in pannel)
Definition: PixelEndcapName.h:53
SiPixelHitEfficiencyModule::book
void book(const edm::ParameterSet &, edm::EventSetup const &, DQMStore::IBooker &, int type=0, bool isUpgrade=false)
Definition: SiPixelHitEfficiencyModule.cc:57
type
type
Definition: HCALResponse.h:21
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:844
Reference_intrackfit_cff.barrel
barrel
Definition: Reference_intrackfit_cff.py:37
PixelBarrelName::moduleName
int moduleName() const
module id (index in z)
Definition: PixelBarrelName.h:48
Reference_intrackfit_cff.endcap
endcap
Definition: Reference_intrackfit_cff.py:82
LocalTrajectoryParameters::position
LocalPoint position() const
Local x and y position coordinates.
Definition: LocalTrajectoryParameters.h:84
SiPixelHitEfficiencyModule::eff
std::pair< double, double > eff(double nValid, double nMissing)
Definition: SiPixelHitEfficiencyModule.cc:728
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
PixelEndcapName::bladeName
int bladeName() const
blade id
Definition: PixelEndcapName.h:47
PV3DBase::mag
T mag() const
Definition: PV3DBase.h:67
PixelBarrelName::isHalfModule
bool isHalfModule() const
full or half module
Definition: PixelBarrelName.cc:399
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
PXBDetId::module
unsigned int module() const
det id
Definition: PXBDetId.h:43
MonitorElement::setBinError
void setBinError(int binx, double error)
set uncertainty on content of bin (1-D)
Definition: MonitorElement.cc:869
alcazmumu_cfi.src
tuple src
Definition: alcazmumu_cfi.py:30
LocalTrajectoryParameters::momentum
LocalVector momentum() const
Momentum vector in the local frame.
Definition: LocalTrajectoryParameters.h:89
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
SiPixelHistogramId::setHistoId
std::string setHistoId(std::string variable, uint32_t &rawId)
Set Histogram Id.
Definition: SiPixelHistogramId.cc:39
DetId
Definition: DetId.h:18
PixelBarrelName::ladderName
int ladderName() const
ladder id (index in phi)
Definition: PixelBarrelName.h:51
SiPixelHitEfficiencyModule::computeEfficiencies
void computeEfficiencies(bool modon=true, bool ladon=true, bool layon=true, bool phion=true, bool bladeon=true, bool diskon=true, bool ringon=true)
Definition: SiPixelHitEfficiencyModule.cc:591
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
PixelBarrelName::layerName
int layerName() const
layer id
Definition: PixelBarrelName.h:45
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:42
MonitorElement::getBinContent
double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:743
SiPixelHitEfficiencyModule::fill
void fill(const LocalTrajectoryParameters &ltp, bool isHitValid, bool modon=true, bool ladon=true, bool layon=true, bool phion=true, bool bladeon=true, bool diskon=true, bool ringon=true)
Definition: SiPixelHitEfficiencyModule.cc:477
PixelEndcapName::pannelName
int pannelName() const
pannel id
Definition: PixelEndcapName.h:50
MonitorElement::getNbinsX
int getNbinsX(void) const
get # of bins in X-axis
Definition: MonitorElement.cc:725
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
PixelEndcapName::diskName
int diskName() const
disk id
Definition: PixelEndcapName.h:44
postValidation_cfi.efficiency
tuple efficiency
Definition: postValidation_cfi.py:73
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:947
Generated for CMSSW Reference Manual by doxygen 1.8.5