CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
SiStripRecHitsValid.cc
Go to the documentation of this file.
1 // File: SiStripRecHitsValid.cc
2 // Description: see SiStripRecHitsValid.h
3 // Author: P. Azzi
4 // Creation Date: PA May 2006 Initial version.
5 //
6 //--------------------------------------------
7 #include "Validation/TrackerRecHits/interface/SiStripRecHitsValid.h"
8 #include "SimTracker/TrackerHitAssociation/interface/TrackerHitAssociator.h"
9 
10 //needed for the geometry:
11 #include "DataFormats/DetId/interface/DetId.h"
12 #include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
13 #include "DataFormats/SiStripDetId/interface/TECDetId.h"
14 #include "DataFormats/SiStripDetId/interface/TIBDetId.h"
15 #include "DataFormats/SiStripDetId/interface/TIDDetId.h"
16 #include "DataFormats/SiStripDetId/interface/TOBDetId.h"
17 #include "DataFormats/GeometryVector/interface/LocalPoint.h"
18 #include "DataFormats/GeometryVector/interface/GlobalPoint.h"
19 
20 
21 //--- for RecHit
22 #include "DataFormats/SiStripCluster/interface/SiStripCluster.h"
23 #include "DataFormats/SiStripCluster/interface/SiStripClusterCollection.h"
24 #include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2DCollection.h"
25 #include "DataFormats/TrackerRecHit2D/interface/SiStripMatchedRecHit2DCollection.h"
26 #include "DataFormats/Common/interface/OwnVector.h"
27 #include "DQMServices/Core/interface/DQMStore.h"
28 
29 using namespace std;
30 using namespace edm;
31 
32 namespace helper {
33  struct GetDetId {
34  template<typename X>
35  DetId operator()(const X &x) { return DetId(x.detId()); }
36  };
37 
38  template<typename T>
39  std::pair<typename T::DetSet::const_iterator, typename T::DetSet::const_iterator>
40  getRange(const T &detset, const DetId &id) {
41  typedef std::pair<typename T::DetSet::const_iterator, typename T::DetSet::const_iterator> return_type;
42  typename T::const_iterator match = detset.find(id);
43  if (match == detset.end()) return return_type();
44  typename T::DetSet hits = *match;
45  return return_type(hits.begin(), hits.end());
46  }
47 }
48 
49 
50 
51 //Constructor
52 SiStripRecHitsValid::SiStripRecHitsValid(const ParameterSet& ps) :
53  dbe_(0),
54  conf_(ps),
55  matchedRecHits_( ps.getParameter<edm::InputTag>("matchedRecHits") ),
56  rphiRecHits_( ps.getParameter<edm::InputTag>("rphiRecHits") ),
57  stereoRecHits_( ps.getParameter<edm::InputTag>("stereoRecHits") ) {
58 
59  outputFile_ = ps.getUntrackedParameter<string>("outputFile", "sistriprechitshisto.root");
60  dbe_ = Service<DQMStore>().operator->();
61  //dbe_->showDirStructure();
62  dbe_->setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Strip/SISTRIP");
63 
64  meNumTotRphi = dbe_->book1D("NumTotRphi","Num of RecHits rphi",100, 0, 10000);
65  meNumTotSas = dbe_->book1D("NumTotSas","Num of RecHits sas",100, 0, 10000);
66  meNumTotMatched = dbe_->book1D("NumTotMatched","Num of RecHits rmatched",100, 0, 10000);
67 
68  dbe_->setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Strip/TIB");
69  meNumRphiTIB = dbe_->book1D("NumRphiTIB","Num of RecHits rphi", 100, 0, 1000.);
70  meNumSasTIB = dbe_->book1D("NumSasTIB","Num of RecHits sas", 100, 0, 1000.);
71  meNumMatchedTIB = dbe_->book1D("NumMatchedTIB","Num of RecHits matched", 100, 0, 1000.);
72 
73  //one histo per Layer rphi hits
74  for(int i = 0 ;i<4 ; i++) {
75  Char_t histo[200];
76  sprintf(histo,"Nstp_rphi_layer%dtib",i+1);
77  meNstpRphiTIB[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);
78  sprintf(histo,"Adc_rphi_layer%dtib",i+1);
79  meAdcRphiTIB[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);
80  sprintf(histo,"Posx_rphi_layer%dtib",i+1);
81  mePosxRphiTIB[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);
82  sprintf(histo,"Errx_rphi_layer%dtib",i+1);
83  meErrxRphiTIB[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0,0.01); //<error>~20micron
84  sprintf(histo,"Res_rphi_layer%dtib",i+1);
85  meResRphiTIB[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.02,+0.02);
86  sprintf(histo,"Pull_LF_rphi_layer%dtib",i+1);
87  mePullLFRphiTIB[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
88  sprintf(histo,"Pull_MF_rphi_layer%dtib",i+1);
89  mePullMFRphiTIB[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
90  sprintf(histo,"Chi2_rphi_layer%dtib",i+1);
91  meChi2RphiTIB[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);
92  }
93 
94  //one histo per Layer stereo and matched hits
95  for(int i = 0 ;i<2 ; i++) {
96  Char_t histo[200];
97  sprintf(histo,"Nstp_sas_layer%dtib",i+1);
98  meNstpSasTIB[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);
99  sprintf(histo,"Adc_sas_layer%dtib",i+1);
100  meAdcSasTIB[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);
101  sprintf(histo,"Posx_sas_layer%dtib",i+1);
102  mePosxSasTIB[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);
103  sprintf(histo,"Errx_sas_layer%dtib",i+1);
104  meErrxSasTIB[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0.,0.01);
105  sprintf(histo,"Res_sas_layer%dtib",i+1);
106  meResSasTIB[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.02,+0.02);
107  sprintf(histo,"Pull_LF_sas_layer%dtib",i+1);
108  mePullLFSasTIB[i] = dbe_->book1D(histo,"Pull",100,-4.,4.);
109  sprintf(histo,"Pull_MF_sas_layer%dtib",i+1);
110  mePullMFSasTIB[i] = dbe_->book1D(histo,"Pull",100,-4.,4.);
111  sprintf(histo,"Chi2_sas_layer%dtib",i+1);
112  meChi2SasTIB[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);
113 
114  sprintf(histo,"Posx_matched_layer%dtib",i+1);
115  mePosxMatchedTIB[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0, +6.0);
116  sprintf(histo,"Posy_matched_layer%dtib",i+1);
117  mePosyMatchedTIB[i] = dbe_->book1D(histo,"RecHit y coord.",100,-6.0, +6.0);
118  sprintf(histo,"Errx_matched_layer%dtib",i+1);
119  meErrxMatchedTIB[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0., 0.01);
120  sprintf(histo,"Erry_matched_layer%dtib",i+1);
121  meErryMatchedTIB[i] = dbe_->book1D(histo,"RecHit err(y) coord.",100,0., 0.05);
122  sprintf(histo,"Resx_matched_layer%dtib",i+1);
123  meResxMatchedTIB[i] = dbe_->book1D(histo,"RecHit Res(x) coord.",100,-0.02, +0.02);
124  sprintf(histo,"Resy_matched_layer%dtib",i+1);
125  meResyMatchedTIB[i] = dbe_->book1D(histo,"RecHit Res(y) coord.",100,-1., +1.);
126  sprintf(histo,"Chi2_matched_layer%dtib",i+1);
127  meChi2MatchedTIB[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);
128  }
129 
130  dbe_->setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Strip/TOB");
131  meNumRphiTOB = dbe_->book1D("NumRphiTOB","Num of RecHits rphi", 100, 0, 1000.);
132  meNumSasTOB = dbe_->book1D("NumSasTOB","Num of RecHits sas", 100, 0, 1000.);
133  meNumMatchedTOB = dbe_->book1D("NumMatchedTOB","Num of RecHits matched", 100, 0, 1000.);
134  //one histo per Layer rphi hits
135  for(int i = 0 ;i<6 ; i++) {
136  Char_t histo[200];
137  sprintf(histo,"Nstp_rphi_layer%dtob",i+1);
138  meNstpRphiTOB[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);
139  sprintf(histo,"Adc_rphi_layer%dtob",i+1);
140  meAdcRphiTOB[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);
141  sprintf(histo,"Posx_rphi_layer%dtob",i+1);
142  mePosxRphiTOB[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);
143  sprintf(histo,"Errx_rphi_layer%dtob",i+1);
144  meErrxRphiTOB[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0,0.01);
145  sprintf(histo,"Res_rphi_layer%dtob",i+1);
146  meResRphiTOB[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.02,+0.02);
147  sprintf(histo,"Pull_LF_rphi_layer%dtob",i+1);
148  mePullLFRphiTOB[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
149  sprintf(histo,"Pull_MF_rphi_layer%dtob",i+1);
150  mePullMFRphiTOB[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
151  sprintf(histo,"Chi2_rphi_layer%dtob",i+1);
152  meChi2RphiTOB[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);
153  }
154 
155  //one histo per Layer stereo and matched hits
156  for(int i = 0 ;i<2 ; i++) {
157  Char_t histo[200];
158  sprintf(histo,"Nstp_sas_layer%dtob",i+1);
159  meNstpSasTOB[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);
160  sprintf(histo,"Adc_sas_layer%dtob",i+1);
161  meAdcSasTOB[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);
162  sprintf(histo,"Posx_sas_layer%dtob",i+1);
163  mePosxSasTOB[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);
164  sprintf(histo,"Errx_sas_layer%dtob",i+1);
165  meErrxSasTOB[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0.,0.01);
166  sprintf(histo,"Res_sas_layer%dtob",i+1);
167  meResSasTOB[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.02,+0.02);
168  sprintf(histo,"Pull_LF_sas_layer%dtob",i+1);
169  mePullLFSasTOB[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
170  sprintf(histo,"Pull_MF_sas_layer%dtob",i+1);
171  mePullMFSasTOB[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
172  sprintf(histo,"Chi2_sas_layer%dtob",i+1);
173  meChi2SasTOB[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);
174 
175  sprintf(histo,"Posx_matched_layer%dtob",i+1);
176  mePosxMatchedTOB[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0, +6.0);
177  sprintf(histo,"Posy_matched_layer%dtob",i+1);
178  mePosyMatchedTOB[i] = dbe_->book1D(histo,"RecHit y coord.",100,-10.0, +10.0);
179  sprintf(histo,"Errx_matched_layer%dtob",i+1);
180  meErrxMatchedTOB[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0., 0.01);
181  sprintf(histo,"Erry_matched_layer%dtob",i+1);
182  meErryMatchedTOB[i] = dbe_->book1D(histo,"RecHit err(y) coord.",100,0., 0.1);
183  sprintf(histo,"Resx_matched_layer%dtob",i+1);
184  meResxMatchedTOB[i] = dbe_->book1D(histo,"RecHit Res(x) coord.",100,-0.02, +0.02);
185  sprintf(histo,"Resy_matched_layer%dtob",i+1);
186  meResyMatchedTOB[i] = dbe_->book1D(histo,"RecHit Res(y) coord.",100,-1., +1.);
187  sprintf(histo,"Chi2_matched_layer%dtob",i+1);
188  meChi2MatchedTOB[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);
189  }
190 
191  dbe_->setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Strip/TID");
192  meNumRphiTID = dbe_->book1D("NumRphiTID","Num of RecHits rphi", 100, 0, 1000.);
193  meNumSasTID = dbe_->book1D("NumSasTID","Num of RecHits sas", 100, 0, 1000.);
194  meNumMatchedTID = dbe_->book1D("NumMatchedTID","Num of RecHits matched", 100, 0, 1000.);
195 
196  //one histo per Ring rphi hits: 3 rings, 6 disks, 2 inner rings are glued
197  for(int i = 0 ;i<3 ; i++) {
198  Char_t histo[200];
199  sprintf(histo,"Nstp_rphi_layer%dtid",i+1);
200  meNstpRphiTID[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);
201  sprintf(histo,"Adc_rphi_layer%dtid",i+1);
202  meAdcRphiTID[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);
203  sprintf(histo,"Posx_rphi_layer%dtid",i+1);
204  mePosxRphiTID[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);
205  sprintf(histo,"Errx_rphi_layer%dtid",i+1);
206  meErrxRphiTID[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0,0.5);
207  sprintf(histo,"Res_rphi_layer%dtid",i+1);
208  meResRphiTID[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.5,+0.5);
209  sprintf(histo,"Pull_LF_rphi_layer%dtid",i+1);
210  mePullLFRphiTID[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
211  sprintf(histo,"Pull_MF_rphi_layer%dtid",i+1);
212  mePullMFRphiTID[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
213  sprintf(histo,"Chi2_rphi_layer%dtid",i+1);
214  meChi2RphiTID[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);
215  }
216 
217  //one histo per Ring stereo and matched hits
218  for(int i = 0 ;i<2 ; i++) {
219  Char_t histo[200];
220  sprintf(histo,"Nstp_sas_layer%dtid",i+1);
221  meNstpSasTID[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);
222  sprintf(histo,"Adc_sas_layer%dtid",i+1);
223  meAdcSasTID[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);
224  sprintf(histo,"Posx_sas_layer%dtid",i+1);
225  mePosxSasTID[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);
226  sprintf(histo,"Errx_sas_layer%dtid",i+1);
227  meErrxSasTID[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0.,0.5);
228  sprintf(histo,"Res_sas_layer%dtid",i+1);
229  meResSasTID[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.5,+0.5);
230  sprintf(histo,"Pull_LF_sas_layer%dtid",i+1);
231  mePullLFSasTID[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
232  sprintf(histo,"Pull_MF_sas_layer%dtid",i+1);
233  mePullMFSasTID[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
234  sprintf(histo,"Chi2_sas_layer%dtid",i+1);
235  meChi2SasTID[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);
236 
237  sprintf(histo,"Posx_matched_layer%dtid",i+1);
238  mePosxMatchedTID[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0, +6.0);
239  sprintf(histo,"Posy_matched_layer%dtid",i+1);
240  mePosyMatchedTID[i] = dbe_->book1D(histo,"RecHit y coord.",100,-6.0, +6.0);
241  sprintf(histo,"Errx_matched_layer%dtid",i+1);
242  meErrxMatchedTID[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0., 0.05);
243  sprintf(histo,"Erry_matched_layer%dtid",i+1);
244  meErryMatchedTID[i] = dbe_->book1D(histo,"RecHit err(y) coord.",100,0., 0.1);
245  sprintf(histo,"Resx_matched_layer%dtid",i+1);
246  meResxMatchedTID[i] = dbe_->book1D(histo,"RecHit Res(x) coord.",100,-0.2, +0.2);
247  sprintf(histo,"Resy_matched_layer%dtid",i+1);
248  meResyMatchedTID[i] = dbe_->book1D(histo,"RecHit Res(y) coord.",100,-1., +1.);
249  sprintf(histo,"Chi2_matched_layer%dtid",i+1);
250  meChi2MatchedTID[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);
251  }
252 
253  dbe_->setCurrentFolder("TrackerRecHitsV/TrackerRecHits/Strip/TEC");
254  meNumRphiTEC = dbe_->book1D("NumRphiTEC","Num of RecHits rphi", 100, 0, 1000.);
255  meNumSasTEC = dbe_->book1D("NumSasTEC","Num of RecHits sas", 100, 0, 1000.);
256  meNumMatchedTEC = dbe_->book1D("NumMatchedTEC","Num of RecHits matched", 100, 0, 1000.);
257 
258  //one histo per Ring rphi hits: 7 rings, 18 disks. Innermost 3 rings are same as TID above.
259  for(int i = 0 ;i<7 ; i++) {
260  Char_t histo[200];
261  sprintf(histo,"Nstp_rphi_layer%dtec",i+1);
262  meNstpRphiTEC[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);
263  sprintf(histo,"Adc_rphi_layer%dtec",i+1);
264  meAdcRphiTEC[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);
265  sprintf(histo,"Posx_rphi_layer%dtec",i+1);
266  mePosxRphiTEC[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);
267  sprintf(histo,"Errx_rphi_layer%dtec",i+1);
268  meErrxRphiTEC[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0,0.5);
269  sprintf(histo,"Res_rphi_layer%dtec",i+1);
270  meResRphiTEC[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.5,+0.5);
271  sprintf(histo,"Pull_LF_rphi_layer%dtec",i+1);
272  mePullLFRphiTEC[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
273  sprintf(histo,"Pull_MF_rphi_layer%dtec",i+1);
274  mePullMFRphiTEC[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
275  sprintf(histo,"Chi2_rphi_layer%dtec",i+1);
276  meChi2RphiTEC[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);
277  }
278 
279  //one histo per Layer stereo and matched hits: rings 1,2,5 are double sided
280  for(int i = 0 ;i<5 ; i++) {
281  if(i == 0 || i == 1 || i == 4) {
282  Char_t histo[200];
283  sprintf(histo,"Nstp_sas_layer%dtec",i+1);
284  meNstpSasTEC[i] = dbe_->book1D(histo,"RecHit Cluster Size",10,0.5,10.5);
285  sprintf(histo,"Adc_sas_layer%dtec",i+1);
286  meAdcSasTEC[i] = dbe_->book1D(histo,"RecHit Cluster Charge",100,0.,300.);
287  sprintf(histo,"Posx_sas_layer%dtec",i+1);
288  mePosxSasTEC[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0,+6.0);
289  sprintf(histo,"Errx_sas_layer%dtec",i+1);
290  meErrxSasTEC[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0.,0.5);
291  sprintf(histo,"Res_sas_layer%dtec",i+1);
292  meResSasTEC[i] = dbe_->book1D(histo,"RecHit Residual",100,-0.5,+0.5);
293  sprintf(histo,"Pull_LF_sas_layer%dtec",i+1);
294  mePullLFSasTEC[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
295  sprintf(histo,"Pull_MF_sas_layer%dtec",i+1);
296  mePullMFSasTEC[i] = dbe_->book1D(histo,"Pull",100,-5.,5.);
297  sprintf(histo,"Chi2_sas_layer%dtec",i+1);
298  meChi2SasTEC[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);
299 
300  sprintf(histo,"Posx_matched_layer%dtec",i+1);
301  mePosxMatchedTEC[i] = dbe_->book1D(histo,"RecHit x coord.",100,-6.0, +6.0);
302  sprintf(histo,"Posy_matched_layer%dtec",i+1);
303  mePosyMatchedTEC[i] = dbe_->book1D(histo,"RecHit y coord.",100,-8.0, +8.0);
304  sprintf(histo,"Errx_matched_layer%dtec",i+1);
305  meErrxMatchedTEC[i] = dbe_->book1D(histo,"RecHit err(x) coord.",100,0., 0.05);
306  sprintf(histo,"Erry_matched_layer%dtec",i+1);
307  meErryMatchedTEC[i] = dbe_->book1D(histo,"RecHit err(y) coord.",100,0., 0.1);
308  sprintf(histo,"Resx_matched_layer%dtec",i+1);
309  meResxMatchedTEC[i] = dbe_->book1D(histo,"RecHit Res(x) coord.",100,-0.2, +0.2);
310  sprintf(histo,"Resy_matched_layer%dtec",i+1);
311  meResyMatchedTEC[i] = dbe_->book1D(histo,"RecHit Res(y) coord.",100,-1., +1.);
312  sprintf(histo,"Chi2_matched_layer%dtec",i+1);
313  meChi2MatchedTEC[i] = dbe_->book1D(histo,"RecHit Chi2 test",100,0., 50);
314  }
315  }
316 }
317 
318 
319 SiStripRecHitsValid::~SiStripRecHitsValid(){
320  //if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
321 }
322 
323 void SiStripRecHitsValid::beginJob(){
324 
325 }
326 
327 void SiStripRecHitsValid::endJob() {
328  if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
329 }
330 
331 
332 void SiStripRecHitsValid::analyze(const edm::Event& e, const edm::EventSetup& es) {
333 
334  LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();
335  //cout << " Run = " << e.id().run() << " Event = " << e.id().event() << endl;
336 
337  //--- get RecHits
338 
339  std::string rechitProducer = "SiStripRecHits2D";
340 
341  // Step A: Get Inputs
342  edm::Handle<SiStripMatchedRecHit2DCollection> rechitsmatched;
343  edm::Handle<SiStripRecHit2DCollection> rechitsrphi;
344  edm::Handle<SiStripRecHit2DCollection> rechitsstereo;
345  e.getByLabel(matchedRecHits_, rechitsmatched);
346  e.getByLabel(rphiRecHits_, rechitsrphi);
347  e.getByLabel(stereoRecHits_, rechitsstereo);
348 
349  int numrechitrphi =0;
350  int numrechitsas =0;
351  int numrechitmatched=0;
352 
353  int totTibnumrechitrphi=0;
354  int totTibnumrechitsas=0;
355  int totTibnumrechitmatched=0;
356  int totTobnumrechitrphi=0;
357  int totTobnumrechitsas=0;
358  int totTobnumrechitmatched=0;
359  int totTidnumrechitrphi=0;
360  int totTidnumrechitsas=0;
361  int totTidnumrechitmatched=0;
362  int totTecnumrechitrphi=0;
363  int totTecnumrechitsas=0;
364  int totTecnumrechitmatched=0;
365  int totrechitrphi =0;
366  int totrechitsas =0;
367  int totrechitmatched =0;
368 
369  TrackerHitAssociator associate(e, conf_);
370 
371  edm::ESHandle<TrackerGeometry> pDD;
372  es.get<TrackerDigiGeometryRecord> ().get (pDD);
373  const TrackerGeometry &tracker(*pDD);
374  // FIXME: this using of vector<DetId> is suboptimal, but I don't want to re-write the full class now
375  std::vector<DetId> IDs;
376  IDs.reserve(rechitsrphi->size() + rechitsmatched->size() + rechitsstereo->size());
377  std::transform(rechitsrphi->begin(), rechitsrphi->end(), std::back_inserter(IDs), helper::GetDetId() );
378  std::transform(rechitsstereo->begin(), rechitsstereo->end(), std::back_inserter(IDs), helper::GetDetId() );
379  std::transform(rechitsmatched->begin(), rechitsmatched->end(), std::back_inserter(IDs), helper::GetDetId() );
380  // loop over detunits
381  // for(TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++){
382  for(std::vector<DetId>::const_iterator it = IDs.begin(); it != IDs.end(); ++it ){//loop on rphi detector with hits
383  uint32_t myid=((*it).rawId());
384  DetId detid = ((*it));
385 
386  // initialize here
387  for(int i=0; i<MAXHIT; i++){
388  rechitrphix[i] =0;
389  rechitrphierrx[i] =0;
390  rechitrphiy[i] =0;
391  rechitrphiz[i] =0;
392  rechitsasx[i] =0;
393  rechitsaserrx[i] =0;
394  rechitsasy[i] =0;
395  rechitsasz[i] =0;
396  clusizrphi[i] =0;
397  clusizsas[i] =0;
398  cluchgrphi[i] =0;
399  cluchgsas[i] =0;
400  rechitrphires[i]=-999.;
401  rechitsasres[i]=-999.;
402  rechitrphipullMF[i]=-999.;
403  rechitsaspullMF[i]=-999.;
404  chi2rphi[i] =0;
405  chi2sas[i]=0;
406  rechitmatchedx[i] =0;
407  rechitmatchedy[i] =0;
408  rechitmatchedz[i] =0;
409  rechitmatchederrxx[i] =0;
410  rechitmatchederrxy[i] =0;
411  rechitmatchederryy[i] =0;
412  rechitmatchedresx[i]=-999;
413  rechitmatchedresy[i]=-999;
414  chi2matched[i]=0;
415  }
416 
417  numrechitrphi =0;
418  //loop over rechits-rphi in the same subdetector
419  std::pair<SiStripRecHit2DCollection::DetSet::const_iterator,SiStripRecHit2DCollection::DetSet::const_iterator> rechitrphiRange = helper::getRange(*rechitsrphi, detid);
420  SiStripRecHit2DCollection::DetSet::const_iterator rechitrphiRangeIteratorBegin = rechitrphiRange.first;
421  SiStripRecHit2DCollection::DetSet::const_iterator rechitrphiRangeIteratorEnd = rechitrphiRange.second;
422  SiStripRecHit2DCollection::DetSet::const_iterator iterrphi=rechitrphiRangeIteratorBegin;
423 
424  numrechitrphi = rechitrphiRangeIteratorEnd - rechitrphiRangeIteratorBegin;
425 
426 
427  int i=0;
428  int i2=0;
429 
430  if(numrechitrphi > 0 ){
431  totrechitrphi+=numrechitrphi;
432  for(iterrphi=rechitrphiRangeIteratorBegin; iterrphi!=rechitrphiRangeIteratorEnd;++iterrphi){
433  const GeomDetUnit * det = tracker.idToDetUnit(detid);
434  const StripGeomDetUnit * stripdet=(const StripGeomDetUnit*)(det);
435  const StripTopology &topol=(StripTopology&)stripdet->topology();
436  SiStripRecHit2D const rechit=*iterrphi;
437  LocalPoint position=rechit.localPosition();
438  LocalError error=rechit.localPositionError();
439  MeasurementPoint Mposition;
440  MeasurementError Merror;
441  Mposition = topol.measurementPosition(position);
442  Merror = topol.measurementError(position,error);
443  SiStripRecHit2D::ClusterRef clust=rechit.cluster();
444  int clusiz=0;
445  int totcharge=0;
446  clusiz = clust->amplitudes().size();
447  const std::vector<uint8_t> amplitudes=clust->amplitudes();
448  for(size_t ia=0; ia<amplitudes.size();ia++){
449  totcharge+=amplitudes[ia];
450  }
451  rechitrphix[i] = position.x();
452  rechitrphiy[i] = position.y();
453  rechitrphiz[i] = position.z();
454  rechitrphierrx[i] = error.xx();
455  clusizrphi[i] = clusiz;
456  cluchgrphi[i] = totcharge;
457 
458  matched.clear();
459  matched = associate.associateHit(rechit);
460  float mindist = 999999;
461  float dist = 999999;
462  PSimHit closest;
463  if(!matched.empty()){
464  for(vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); m++){
465  dist = fabs(rechitrphix[i] - (*m).localPosition().x());
466  if(dist<mindist){
467  mindist = dist;
468  closest = (*m);
469  }
470  rechitrphires[i] = rechitrphix[i] - closest.localPosition().x();
471  }
472  rechitrphipullMF[i] = (Mposition.x() - (topol.measurementPosition(closest.localPosition())).x())/sqrt(Merror.uu());
473 
474  //chi2test compare rechit errors with the simhit position ( using null matrix for the simhit).
475  //Can spot problems in the geometry better than a simple residual. (thanks to BorisM)
476  AlgebraicVector rhparameters = rechit.parameters();
477  AlgebraicVector shparameters(2);
478  shparameters[0] = closest.localPosition().x();
479  shparameters[1] = closest.localPosition().y();
480  AlgebraicVector r(rhparameters - shparameters);
481  AlgebraicSymMatrix R = rechit.parametersError();
482  int ierr;
483  R.invert(ierr); // if (ierr != 0) throw exception;
484  double est = R.similarity(r);
485  // std::cout << " ====== Chi2 test rphi hits ====== " << std::endl;
486  // std::cout << "RecHit param. = " << rhparameters << std::endl;
487  // std::cout << "RecHit errors = " << R << std::endl;
488  // std::cout << "SimHit param. = " << shparameters << std::endl;
489  // std::cout << " chi2 = " << est << std::endl;
490  // std::cout << "DEBUG BORIS,filling chi2rphi[i],i: " << i << std::endl;
491  chi2rphi[i2] = est;
492  i2++;
493  }
494  i++;
495  }
496  }
497 
498  //loop over rechits-sas in the same subdetector
499  int j=0;
500  int j2=0;
501  numrechitsas=0;
502  std::pair<SiStripRecHit2DCollection::DetSet::const_iterator,SiStripRecHit2DCollection::DetSet::const_iterator> rechitsasRange = helper::getRange(*rechitsstereo, detid);
503  SiStripRecHit2DCollection::DetSet::const_iterator rechitsasRangeIteratorBegin = rechitsasRange.first;
504  SiStripRecHit2DCollection::DetSet::const_iterator rechitsasRangeIteratorEnd = rechitsasRange.second;
505  SiStripRecHit2DCollection::DetSet::const_iterator itersas=rechitsasRangeIteratorBegin;
506  numrechitsas = rechitsasRangeIteratorEnd - rechitsasRangeIteratorBegin;
507  if(numrechitsas > 0){
508  totrechitsas+=numrechitsas;
509  for(itersas=rechitsasRangeIteratorBegin; itersas!=rechitsasRangeIteratorEnd;++itersas){
510  const GeomDetUnit * det = tracker.idToDetUnit(detid);
511  const StripGeomDetUnit * stripdet=(const StripGeomDetUnit*)(det);
512  const StripTopology &topol=(StripTopology&)stripdet->topology();
513  SiStripRecHit2D const rechit=*itersas;
514  LocalPoint position=rechit.localPosition();
515  LocalError error=rechit.localPositionError();
516  MeasurementPoint Mposition;
517  MeasurementError Merror;
518  Mposition = topol.measurementPosition(position);
519  Merror = topol.measurementError(position,error);
520  SiStripRecHit2D::ClusterRef clust=rechit.cluster(); int clusiz=0;
521  int totcharge=0;
522  clusiz = clust->amplitudes().size();
523  const std::vector<uint8_t> amplitudes=clust->amplitudes();
524  for(size_t ia=0; ia<amplitudes.size();ia++){
525  totcharge+=amplitudes[ia];
526  }
527 
528  rechitsasx[j] = position.x();
529  rechitsasy[j] = position.y();
530  rechitsasz[j] = position.z();
531  rechitsaserrx[j] = error.xx();
532  clusizsas[j] = clusiz;
533  cluchgsas[j] = totcharge;
534 
535  float mindist = 999999;
536  float dist = 999999;
537  PSimHit closest;
538  matched.clear();
539  matched = associate.associateHit(rechit);
540  if(!matched.empty()){
541  for(vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); m++){
542  dist = fabs(rechitsasx[j] - (*m).localPosition().x());
543  if(dist<mindist){
544  mindist = dist;
545  closest = (*m);
546  }
547  rechitsasres[j] = rechitsasx[j] - closest.localPosition().x();
548  }
549  rechitsaspullMF[j] = (Mposition.x() - (topol.measurementPosition(closest.localPosition())).x())/sqrt(Merror.uu());
550  //chi2test compare rechit errors with the simhit position ( using null matrix for the simhit).
551  //Can spot problems in the geometry better than a simple residual. (thanks to BorisM)
552  AlgebraicVector rhparameters = rechit.parameters();
553  AlgebraicVector shparameters(2);
554  shparameters[0] = closest.localPosition().x();
555  shparameters[1] = closest.localPosition().y();
556  AlgebraicVector r(rhparameters - shparameters);
557  AlgebraicSymMatrix R = rechit.parametersError();
558  int ierr;
559  R.invert(ierr); // if (ierr != 0) throw exception;
560  double est = R.similarity(r);
561  // std::cout << " ====== Chi2 test sas hits ====== " << std::endl;
562  // std::cout << "RecHit param. = " << rhparameters << std::endl;
563  // std::cout << "RecHit errors = " << R << std::endl;
564  // std::cout << "SimHit param. = " << shparameters << std::endl;
565  // std::cout << " chi2 = " << est << std::endl;
566  // std::cout << "DEBUG BORIS,filling chi2sas[j],j: " << i << std::endl;
567  chi2sas[j2] = est;
568  j2++;
569  }
570  j++;
571  }
572  }
573 
574  //now matched hits
575 
576  int k=0;
577  int k2=0;
578 
579  //loop over rechits-matched in the same subdetector
580  numrechitmatched=0;
581  std::pair<SiStripMatchedRecHit2DCollection::DetSet::const_iterator,SiStripMatchedRecHit2DCollection::DetSet::const_iterator> rechitmatchedRange = helper::getRange(*rechitsmatched, detid);
582  SiStripMatchedRecHit2DCollection::DetSet::const_iterator rechitmatchedRangeIteratorBegin = rechitmatchedRange.first;
583  SiStripMatchedRecHit2DCollection::DetSet::const_iterator rechitmatchedRangeIteratorEnd = rechitmatchedRange.second;
584  SiStripMatchedRecHit2DCollection::DetSet::const_iterator itermatched=rechitmatchedRangeIteratorBegin;
585  numrechitmatched = rechitmatchedRangeIteratorEnd - rechitmatchedRangeIteratorBegin;
586  if(numrechitmatched > 0){
587  totrechitmatched +=numrechitmatched;
588 
589  for(itermatched=rechitmatchedRangeIteratorBegin; itermatched!=rechitmatchedRangeIteratorEnd;++itermatched){
590  SiStripMatchedRecHit2D const rechit=*itermatched;
591  LocalPoint position=rechit.localPosition();
592  LocalError error=rechit.localPositionError();
593 
594  float mindist = 999999;
595  float distx = 999999;
596  float disty = 999999;
597  float dist = 999999;
598  std::pair<LocalPoint,LocalVector> closestPair;
599  matched.clear();
600  // const SiStripRecHit2D *mono = rechit.monoHit();
601  // const SiStripRecHit2D *st = rechit.stereoHit();
602  // LocalPoint monopos = mono->localPosition();
603  // LocalPoint stpos = st->localPosition();
604 
605  rechitmatchedx[k] = position.x();
606  rechitmatchedy[k] = position.y();
607  rechitmatchedz[k] = position.z();
608  rechitmatchederrxx[k] = error.xx();
609  rechitmatchederrxy[k] = error.xy();
610  rechitmatchederryy[k] = error.yy();
611 
612  // std::cout << " before association " << std::endl;
613  matched = associate.associateHit(rechit);
614  //std::cout << " after association size = " << matched.size() << std::endl;
615 
616  if(!matched.empty()){
617  //project simhit;
618  const GluedGeomDet* gluedDet = (const GluedGeomDet*)tracker.idToDet(rechit.geographicalId());
619  const StripGeomDetUnit* partnerstripdet =(StripGeomDetUnit*) gluedDet->stereoDet();
620  std::pair<LocalPoint,LocalVector> hitPair;
621 
622  //std::cout << " RECHIT position = " << position << std::endl;
623 
624  for(vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); m++){
625  //project simhit;
626  hitPair= projectHit((*m),partnerstripdet,gluedDet->surface());
627  distx = fabs(rechitmatchedx[k] - hitPair.first.x());
628  disty = fabs(rechitmatchedy[k] - hitPair.first.y());
629  dist = sqrt(distx*distx+disty*disty);
630  // std::cout << " Simhit position x = " << hitPair.first.x()
631  // << " y = " << hitPair.first.y() << " dist = " << dist << std::endl;
632  if(dist<mindist){
633  mindist = dist;
634  closestPair = hitPair;
635  }
636  }
637  //std::cout << " Closest position x = " << closestPair.first.x()
638  // << " y = " << closestPair.first.y() << " dist = " << dist << std::endl;
639  rechitmatchedresx[k] = rechitmatchedx[k] - closestPair.first.x();
640  rechitmatchedresy[k] = rechitmatchedy[k] - closestPair.first.y();
641 
642  //chi2test compare rechit errors with the simhit position ( using null matrix for the simhit).
643  //Can spot problems in the geometry better than a simple residual. (thanks to BorisM)
644 
645  AlgebraicVector rhparameters = rechit.parameters();
646  LocalPoint sh = closestPair.first;
647  AlgebraicVector shparameters(2);
648  shparameters[0] = sh.x();
649  shparameters[1] = sh.y();
650  AlgebraicVector r(rhparameters - shparameters);
651  AlgebraicSymMatrix R = rechit.parametersError();
652  int ierr;
653  R.invert(ierr); // if (ierr != 0) throw exception;
654  double est = R.similarity(r);
655  // std::cout << " ====== Chi2 test matched ====== " << std::endl;
656  // std::cout << "RecHit param. = " << rhparameters << std::endl;
657  // std::cout << "RecHit errors = " << R << std::endl;
658  // std::cout << "SimHit param. = " << shparameters << std::endl;
659  // std::cout << " chi2 = " << est << std::endl;
660  chi2matched[k2] = est;
661  k2++;
662  }
663 
664  k++;
665  }
666  }
667 
668  //for each detid
669  if(numrechitrphi>0 || numrechitsas>0 || numrechitmatched){
670  if (detid.subdetId() == int(StripSubdetector::TIB)){
671  TIBDetId tibid(myid);
672  int Tibnumrechitrphi = numrechitrphi;
673  int Tibnumrechitsas = numrechitsas;
674  int Tibnumrechitmatched = numrechitmatched;
675 
676  int Tibnumrechitrphi2 = i2;
677  int Tibnumrechitsas2 = j2;
678  int Tibnumrechitmatched2 = k2;
679 
680  totTibnumrechitrphi +=numrechitrphi;
681  totTibnumrechitsas +=numrechitsas;
682  totTibnumrechitmatched +=numrechitmatched;
683 
684  int ilay = tibid.layer() - 1; //for histogram filling
685 
686  if(tibid.stereo()==0){
687  for(int k = 0; k<Tibnumrechitrphi; k++){
688  meNstpRphiTIB[ilay]->Fill(clusizrphi[k]);
689  meAdcRphiTIB[ilay]->Fill(cluchgrphi[k]);
690  mePosxRphiTIB[ilay]->Fill(rechitrphix[k]);
691  meErrxRphiTIB[ilay]->Fill(sqrt(rechitrphierrx[k]));
692  meResRphiTIB[ilay]->Fill(rechitrphires[k]);
693  mePullLFRphiTIB[ilay]->Fill(rechitrphires[k]/sqrt(rechitrphierrx[k]));
694  mePullMFRphiTIB[ilay]->Fill(rechitrphipullMF[k]);
695 
696  }
697 
698  for(int k = 0; k<Tibnumrechitrphi2; k++){
699  meChi2RphiTIB[ilay]->Fill(chi2rphi[k]);
700  }
701 
702  } else if(tibid.stereo()==1){
703  for(int kk = 0; kk < Tibnumrechitsas; kk++)
704  {
705  meNstpSasTIB[ilay]->Fill(clusizsas[kk]);
706  meAdcSasTIB[ilay]->Fill(cluchgsas[kk]);
707  mePosxSasTIB[ilay]->Fill(rechitsasx[kk]);
708  meErrxSasTIB[ilay]->Fill(sqrt(rechitsaserrx[kk]));
709  meResSasTIB[ilay]->Fill(rechitsasres[kk]);
710  mePullLFSasTIB[ilay]->Fill(rechitsasres[kk]/sqrt(rechitsaserrx[kk]));
711  mePullMFSasTIB[ilay]->Fill(rechitsaspullMF[kk]);
712  }
713  for(int l = 0; l<Tibnumrechitsas2; l++){
714  meChi2SasTIB[ilay]->Fill(chi2sas[l]);
715  }
716 
717  }
718  if(Tibnumrechitmatched>0){
719  for(int kkk = 0; kkk<Tibnumrechitmatched; kkk++)
720  {
721  mePosxMatchedTIB[ilay]->Fill(rechitmatchedx[kkk]);
722  mePosyMatchedTIB[ilay]->Fill(rechitmatchedy[kkk]);
723  meErrxMatchedTIB[ilay]->Fill(sqrt(rechitmatchederrxx[kkk]));
724  meErryMatchedTIB[ilay]->Fill(sqrt(rechitmatchederryy[kkk]));
725  meResxMatchedTIB[ilay]->Fill(rechitmatchedresx[kkk]);
726  meResyMatchedTIB[ilay]->Fill(rechitmatchedresy[kkk]);
727  }
728  for(int l = 0; l<Tibnumrechitmatched2; l++){
729  meChi2MatchedTIB[ilay]->Fill(chi2matched[l]);
730  }
731 
732  }
733  }
734 
735 
736  if (detid.subdetId() == int(StripSubdetector::TOB)){
737  TOBDetId tobid(myid);
738  int Tobnumrechitrphi = numrechitrphi;
739  int Tobnumrechitsas = numrechitsas;
740  int Tobnumrechitmatched = numrechitmatched;
741  totTobnumrechitrphi +=numrechitrphi;
742  totTobnumrechitsas +=numrechitsas;
743  totTobnumrechitmatched +=numrechitmatched;
744 
745  int Tobnumrechitrphi2 = i2;
746  int Tobnumrechitsas2 = j2;
747  int Tobnumrechitmatched2 = k2;
748 
749  int ilay = tobid.layer() - 1; //for histogram filling
750 
751  if(tobid.stereo()==0){
752  for(int k = 0; k<Tobnumrechitrphi; k++){
753  meNstpRphiTOB[ilay]->Fill(clusizrphi[k]);
754  meAdcRphiTOB[ilay]->Fill(cluchgrphi[k]);
755  mePosxRphiTOB[ilay]->Fill(rechitrphix[k]);
756  meErrxRphiTOB[ilay]->Fill(sqrt(rechitrphierrx[k]));
757  meResRphiTOB[ilay]->Fill(rechitrphires[k]);
758  mePullLFRphiTOB[ilay]->Fill(rechitrphires[k]/sqrt(rechitrphierrx[k]));
759  mePullMFRphiTOB[ilay]->Fill(rechitrphipullMF[k]);
760  }
761  for(int l = 0; l<Tobnumrechitrphi2; l++){
762  meChi2RphiTOB[ilay]->Fill(chi2rphi[l]);
763  }
764  } else if(tobid.stereo()==1){
765  for(int kk = 0; kk < Tobnumrechitsas; kk++)
766  {
767  meNstpSasTOB[ilay]->Fill(clusizsas[kk]);
768  meAdcSasTOB[ilay]->Fill(cluchgsas[kk]);
769  mePosxSasTOB[ilay]->Fill(rechitsasx[kk]);
770  meErrxSasTOB[ilay]->Fill(sqrt(rechitsaserrx[kk]));
771  meResSasTOB[ilay]->Fill(rechitsasres[kk]);
772  mePullLFSasTOB[ilay]->Fill(rechitsasres[kk]/sqrt(rechitsaserrx[kk]));
773  mePullMFSasTOB[ilay]->Fill(rechitsaspullMF[kk]);
774  }
775  for(int l = 0; l<Tobnumrechitsas2; l++){
776  meChi2SasTOB[ilay]->Fill(chi2sas[l]);
777  }
778 
779  }
780  if(Tobnumrechitmatched>0){
781  for(int kkk = 0; kkk<Tobnumrechitmatched; kkk++)
782  {
783  mePosxMatchedTOB[ilay]->Fill(rechitmatchedx[kkk]);
784  mePosyMatchedTOB[ilay]->Fill(rechitmatchedy[kkk]);
785  meErrxMatchedTOB[ilay]->Fill(sqrt(rechitmatchederrxx[kkk]));
786  meErryMatchedTOB[ilay]->Fill(sqrt(rechitmatchederryy[kkk]));
787  meResxMatchedTOB[ilay]->Fill(rechitmatchedresx[kkk]);
788  meResyMatchedTOB[ilay]->Fill(rechitmatchedresy[kkk]);
789  }
790  for(int l = 0; l<Tobnumrechitmatched2; l++){
791  meChi2MatchedTOB[ilay]->Fill(chi2matched[l]);
792  }
793 
794  }
795  }
796  if (detid.subdetId() == int(StripSubdetector::TID)){
797  TIDDetId tidid(myid);
798  int Tidnumrechitrphi = numrechitrphi;
799  int Tidnumrechitsas = numrechitsas;
800  int Tidnumrechitmatched = numrechitmatched;
801  totTidnumrechitrphi +=numrechitrphi;
802  totTidnumrechitsas +=numrechitsas;
803  totTidnumrechitmatched +=numrechitmatched;
804 
805  int Tidnumrechitrphi2 = i2;
806  int Tidnumrechitsas2 = j2;
807  int Tidnumrechitmatched2 = k2;
808 
809  int ilay = tidid.ring() - 1; //for histogram filling
810 
811  if(tidid.stereo()==0){
812  for(int k = 0; k<Tidnumrechitrphi; k++){
813  meNstpRphiTID[ilay]->Fill(clusizrphi[k]);
814  meAdcRphiTID[ilay]->Fill(cluchgrphi[k]);
815  mePosxRphiTID[ilay]->Fill(rechitrphix[k]);
816  meErrxRphiTID[ilay]->Fill(sqrt(rechitrphierrx[k]));
817  meResRphiTID[ilay]->Fill(rechitrphires[k]);
818  mePullLFRphiTID[ilay]->Fill(rechitrphires[k]/sqrt(rechitrphierrx[k]));
819  mePullMFRphiTID[ilay]->Fill(rechitrphipullMF[k]);
820  }
821  for(int l = 0; l<Tidnumrechitrphi2; l++){
822  meChi2RphiTID[ilay]->Fill(chi2rphi[l]);
823  }
824 
825  } else if(tidid.stereo()==1){
826  for(int kk = 0; kk < Tidnumrechitsas; kk++)
827  {
828  meNstpSasTID[ilay]->Fill(clusizsas[kk]);
829  meAdcSasTID[ilay]->Fill(cluchgsas[kk]);
830  mePosxSasTID[ilay]->Fill(rechitsasx[kk]);
831  meErrxSasTID[ilay]->Fill(sqrt(rechitsaserrx[kk]));
832  meResSasTID[ilay]->Fill(rechitsasres[kk]);
833  mePullLFSasTID[ilay]->Fill(rechitsasres[kk]/sqrt(rechitsaserrx[kk]));
834  mePullMFSasTID[ilay]->Fill(rechitsaspullMF[kk]);
835  }
836  for(int l = 0; l<Tidnumrechitsas2; l++){
837  meChi2SasTID[ilay]->Fill(chi2sas[l]);
838  }
839 
840  }
841  if(Tidnumrechitmatched>0){
842  for(int kkk = 0; kkk<Tidnumrechitmatched; kkk++)
843  {
844  mePosxMatchedTID[ilay]->Fill(rechitmatchedx[kkk]);
845  mePosyMatchedTID[ilay]->Fill(rechitmatchedy[kkk]);
846  meErrxMatchedTID[ilay]->Fill(sqrt(rechitmatchederrxx[kkk]));
847  meErryMatchedTID[ilay]->Fill(sqrt(rechitmatchederryy[kkk]));
848  meResxMatchedTID[ilay]->Fill(rechitmatchedresx[kkk]);
849  meResyMatchedTID[ilay]->Fill(rechitmatchedresy[kkk]);
850  }
851  for(int l = 0; l<Tidnumrechitmatched2; l++){
852  meChi2MatchedTID[ilay]->Fill(chi2matched[l]);
853  }
854 
855  }
856  }
857  if (detid.subdetId() == int(StripSubdetector::TEC)){
858  TECDetId tecid(myid);
859  int Tecnumrechitrphi = numrechitrphi;
860  int Tecnumrechitsas = numrechitsas;
861  int Tecnumrechitmatched = numrechitmatched;
862  totTecnumrechitrphi +=numrechitrphi;
863  totTecnumrechitsas +=numrechitsas;
864  totTecnumrechitmatched +=numrechitmatched;
865 
866  int Tecnumrechitrphi2 = i2;
867  int Tecnumrechitsas2 = j2;
868  int Tecnumrechitmatched2 = k2;
869 
870  int ilay = tecid.ring() - 1; //for histogram filling
871 
872  if(tecid.stereo()==0){
873  for(int k = 0; k<Tecnumrechitrphi; k++){
874  meNstpRphiTEC[ilay]->Fill(clusizrphi[k]);
875  meAdcRphiTEC[ilay]->Fill(cluchgrphi[k]);
876  mePosxRphiTEC[ilay]->Fill(rechitrphix[k]);
877  meErrxRphiTEC[ilay]->Fill(sqrt(rechitrphierrx[k]));
878  meResRphiTEC[ilay]->Fill(rechitrphires[k]);
879  mePullLFRphiTEC[ilay]->Fill(rechitrphires[k]/sqrt(rechitrphierrx[k]));
880  mePullMFRphiTEC[ilay]->Fill(rechitrphipullMF[k]);
881  }
882  for(int l = 0; l<Tecnumrechitrphi2; l++){
883  meChi2RphiTEC[ilay]->Fill(chi2rphi[l]);
884  }
885 
886  } else if(tecid.stereo()==1){
887  for(int kk = 0; kk < Tecnumrechitsas; kk++)
888  {
889  meNstpSasTEC[ilay]->Fill(clusizsas[kk]);
890  meAdcSasTEC[ilay]->Fill(cluchgsas[kk]);
891  mePosxSasTEC[ilay]->Fill(rechitsasx[kk]);
892  meErrxSasTEC[ilay]->Fill(sqrt(rechitsaserrx[kk]));
893  meResSasTEC[ilay]->Fill(rechitsasres[kk]);
894  mePullLFSasTEC[ilay]->Fill(rechitsasres[kk]/sqrt(rechitsaserrx[kk]));
895  mePullMFSasTEC[ilay]->Fill(rechitsaspullMF[kk]);
896  }
897  for(int l = 0; l<Tecnumrechitsas2; l++){
898  meChi2SasTEC[ilay]->Fill(chi2sas[l]);
899  }
900 
901  }
902  if(Tecnumrechitmatched>0){
903  for(int kkk = 0; kkk<Tecnumrechitmatched; kkk++)
904  {
905  mePosxMatchedTEC[ilay]->Fill(rechitmatchedx[kkk]);
906  mePosyMatchedTEC[ilay]->Fill(rechitmatchedy[kkk]);
907  meErrxMatchedTEC[ilay]->Fill(sqrt(rechitmatchederrxx[kkk]));
908  meErryMatchedTEC[ilay]->Fill(sqrt(rechitmatchederryy[kkk]));
909  meResxMatchedTEC[ilay]->Fill(rechitmatchedresx[kkk]);
910  meResyMatchedTEC[ilay]->Fill(rechitmatchedresy[kkk]);
911  }
912  for(int l = 0; l<Tecnumrechitmatched2; l++){
913  meChi2MatchedTEC[ilay]->Fill(chi2matched[l]);
914  }
915 
916  }
917  }
918 
919  }
920  }
921 
922  //now fill the cumulative histograms of the hits
923 
924  meNumRphiTIB->Fill(totTibnumrechitrphi);
925  meNumSasTIB->Fill(totTibnumrechitsas);
926  meNumMatchedTIB->Fill(totTibnumrechitmatched);
927  meNumRphiTOB->Fill(totTobnumrechitrphi);
928  meNumSasTOB->Fill(totTobnumrechitsas);
929  meNumMatchedTOB->Fill(totTobnumrechitmatched);
930  meNumRphiTID->Fill(totTidnumrechitrphi);
931  meNumSasTID->Fill(totTidnumrechitsas);
932  meNumMatchedTID->Fill(totTidnumrechitmatched);
933  meNumRphiTEC->Fill(totTecnumrechitrphi);
934  meNumSasTEC->Fill(totTecnumrechitsas);
935  meNumMatchedTEC->Fill(totTecnumrechitmatched);
936 
937  meNumTotRphi->Fill(totrechitrphi);
938  meNumTotSas->Fill(totrechitsas);
939  meNumTotMatched->Fill(totrechitmatched);
940 
941 
942 }
943 
944 
945 //needed by to do the residual for matched hits
946 std::pair<LocalPoint,LocalVector> SiStripRecHitsValid::projectHit( const PSimHit& hit, const StripGeomDetUnit* stripDet,
947  const BoundPlane& plane)
948 {
949  // const StripGeomDetUnit* stripDet = dynamic_cast<const StripGeomDetUnit*>(hit.det());
950  //if (stripDet == 0) throw MeasurementDetException("HitMatcher hit is not on StripGeomDetUnit");
951 
952  const StripTopology& topol = stripDet->specificTopology();
953  GlobalPoint globalpos= stripDet->surface().toGlobal(hit.localPosition());
954  LocalPoint localHit = plane.toLocal(globalpos);
955  //track direction
956  LocalVector locdir=hit.localDirection();
957  //rotate track in new frame
958 
959  GlobalVector globaldir= stripDet->surface().toGlobal(locdir);
960  LocalVector dir=plane.toLocal(globaldir);
961  float scale = -localHit.z() / dir.z();
962 
963  LocalPoint projectedPos = localHit + scale*dir;
964 
965  // std::cout << "projectedPos " << projectedPos << std::endl;
966 
967  float selfAngle = topol.stripAngle( topol.strip( hit.localPosition()));
968 
969  LocalVector stripDir( sin(selfAngle), cos(selfAngle), 0); // vector along strip in hit frame
970 
971  LocalVector localStripDir( plane.toLocal(stripDet->surface().toGlobal( stripDir)));
972 
973  return std::pair<LocalPoint,LocalVector>( projectedPos, localStripDir);
974 }
975 
976 
SiStripRecHitsValid::meNstpRphiTEC
MonitorElement * meNstpRphiTEC[7]
Definition: SiStripRecHitsValid.h:159
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:42
SiStripRecHitsValid::meNumTotRphi
MonitorElement * meNumTotRphi
Definition: SiStripRecHitsValid.h:69
RecHit2DLocalPos::parametersError
virtual AlgebraicSymMatrix parametersError() const
Definition: RecHit2DLocalPos.cc:23
SiStripRecHitsValid::meNstpSasTEC
MonitorElement * meNstpSasTEC[5]
Definition: SiStripRecHitsValid.h:167
Surface::toGlobal
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:78
SiStripRecHitsValid::meResRphiTEC
MonitorElement * meResRphiTEC[7]
Definition: SiStripRecHitsValid.h:163
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:44
SiStripRecHitsValid::meChi2RphiTID
MonitorElement * meChi2RphiTID[3]
Definition: SiStripRecHitsValid.h:141
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
LocalError::xx
float xx() const
Definition: LocalError.h:19
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:105
SiStripRecHitsValid::mePullLFSasTIB
MonitorElement * mePullLFSasTIB[4]
Definition: SiStripRecHitsValid.h:99
SiStripRecHitsValid::meErrxSasTEC
MonitorElement * meErrxSasTEC[5]
Definition: SiStripRecHitsValid.h:170
BaseSiTrackerRecHit2DLocalPos::localPositionError
virtual LocalError localPositionError() const
Definition: BaseSiTrackerRecHit2DLocalPos.cc:15
SiStripRecHitsValid::meAdcRphiTIB
MonitorElement * meAdcRphiTIB[4]
Definition: SiStripRecHitsValid.h:87
SiStripRecHitsValid::meErrxMatchedTIB
MonitorElement * meErrxMatchedTIB[2]
Definition: SiStripRecHitsValid.h:104
TOBDetId::layer
unsigned int layer() const
layer id
Definition: TOBDetId.h:39
SiStripRecHitsValid::mePullLFRphiTID
MonitorElement * mePullLFRphiTID[3]
Definition: SiStripRecHitsValid.h:139
StripTopology::stripAngle
virtual float stripAngle(float strip) const =0
SiStripRecHitsValid::meChi2RphiTEC
MonitorElement * meChi2RphiTEC[7]
Definition: SiStripRecHitsValid.h:166
SiStripDetId::stereo
virtual uint32_t stereo() const
Definition: SiStripDetId.h:150
SiStripRecHitsValid::rechitmatchedresy
float rechitmatchedresy[MAXHIT]
Definition: SiStripRecHitsValid.h:219
SiStripRecHitsValid::mePosyMatchedTEC
MonitorElement * mePosyMatchedTEC[5]
Definition: SiStripRecHitsValid.h:176
SiStripRecHitsValid::meNumRphiTIB
MonitorElement * meNumRphiTIB
Definition: SiStripRecHitsValid.h:72
SiStripRecHitsValid::meErrxMatchedTOB
MonitorElement * meErrxMatchedTOB[2]
Definition: SiStripRecHitsValid.h:128
SiStripRecHitsValid::mePosxRphiTID
MonitorElement * mePosxRphiTID[3]
Definition: SiStripRecHitsValid.h:136
SiStripRecHitsValid::meChi2MatchedTID
MonitorElement * meChi2MatchedTID[2]
Definition: SiStripRecHitsValid.h:157
SiStripRecHitsValid::meErryMatchedTID
MonitorElement * meErryMatchedTID[2]
Definition: SiStripRecHitsValid.h:154
SiStripRecHitsValid::matched
std::vector< PSimHit > matched
Definition: SiStripRecHitsValid.h:183
SiStripRecHitsValid::meNstpRphiTIB
MonitorElement * meNstpRphiTIB[4]
Definition: SiStripRecHitsValid.h:86
SiStripRecHitsValid::meResyMatchedTID
MonitorElement * meResyMatchedTID[2]
Definition: SiStripRecHitsValid.h:156
SiStripRecHitsValid::meNumMatchedTOB
MonitorElement * meNumMatchedTOB
Definition: SiStripRecHitsValid.h:77
SiStripRecHitsValid::rechitmatchedresx
float rechitmatchedresx[MAXHIT]
Definition: SiStripRecHitsValid.h:218
DQMStore::save
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:1898
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
SiStripRecHitsValid::meNumSasTOB
MonitorElement * meNumSasTOB
Definition: SiStripRecHitsValid.h:76
SiStripRecHitsValid::meErrxRphiTIB
MonitorElement * meErrxRphiTIB[4]
Definition: SiStripRecHitsValid.h:89
SiStripRecHitsValid::rechitmatchederrxx
float rechitmatchederrxx[MAXHIT]
Definition: SiStripRecHitsValid.h:214
PV3DBase::y
T y() const
Definition: PV3DBase.h:57
SiStripRecHitsValid::meResxMatchedTEC
MonitorElement * meResxMatchedTEC[5]
Definition: SiStripRecHitsValid.h:179
SiStripRecHitsValid::meNstpSasTOB
MonitorElement * meNstpSasTOB[2]
Definition: SiStripRecHitsValid.h:118
SiStripRecHitsValid::mePullLFSasTEC
MonitorElement * mePullLFSasTEC[5]
Definition: SiStripRecHitsValid.h:172
SiStripRecHitsValid::meResSasTID
MonitorElement * meResSasTID[2]
Definition: SiStripRecHitsValid.h:146
X
#define X(str)
Definition: MuonsGrabber.cc:49
SiStripRecHitsValid::meNumMatchedTID
MonitorElement * meNumMatchedTID
Definition: SiStripRecHitsValid.h:80
SiStripRecHitsValid::mePullMFRphiTOB
MonitorElement * mePullMFRphiTOB[6]
Definition: SiStripRecHitsValid.h:116
SiStripRecHitsValid::meAdcRphiTEC
MonitorElement * meAdcRphiTEC[7]
Definition: SiStripRecHitsValid.h:160
SiStripRecHitsValid::meChi2MatchedTIB
MonitorElement * meChi2MatchedTIB[2]
Definition: SiStripRecHitsValid.h:108
SiStripRecHitsValid::meNumSasTIB
MonitorElement * meNumSasTIB
Definition: SiStripRecHitsValid.h:73
trackerHits.histo
tuple histo
Definition: trackerHits.py:12
SiStripRecHitsValid::meResxMatchedTID
MonitorElement * meResxMatchedTID[2]
Definition: SiStripRecHitsValid.h:155
SiStripRecHitsValid::meNstpRphiTID
MonitorElement * meNstpRphiTID[3]
Definition: SiStripRecHitsValid.h:134
StripGeomDetUnit::specificTopology
virtual const StripTopology & specificTopology() const
Returns a reference to the strip proxy topology.
Definition: StripGeomDetUnit.cc:19
SiStripRecHitsValid::mePosxRphiTOB
MonitorElement * mePosxRphiTOB[6]
Definition: SiStripRecHitsValid.h:112
position
static int position[TOTALCHAMBERS][3]
Definition: ReadPGInfo.cc:509
BaseSiTrackerRecHit2DLocalPos::localPosition
virtual LocalPoint localPosition() const
Definition: BaseSiTrackerRecHit2DLocalPos.cc:10
SiStripRecHitsValid::mePullLFRphiTIB
MonitorElement * mePullLFRphiTIB[4]
Definition: SiStripRecHitsValid.h:91
SiStripRecHitsValid::meNumRphiTID
MonitorElement * meNumRphiTID
Definition: SiStripRecHitsValid.h:78
StripTopology::strip
virtual float strip(const LocalPoint &) const =0
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:147
SiStripRecHitsValid::rechitmatchederryy
float rechitmatchederryy[MAXHIT]
Definition: SiStripRecHitsValid.h:216
SiStripRecHitsValid::meErrxRphiTEC
MonitorElement * meErrxRphiTEC[7]
Definition: SiStripRecHitsValid.h:162
SiStripRecHitsValid::rechitmatchederrxy
float rechitmatchederrxy[MAXHIT]
Definition: SiStripRecHitsValid.h:215
SiStripRecHitsValid::meChi2RphiTOB
MonitorElement * meChi2RphiTOB[6]
Definition: SiStripRecHitsValid.h:117
SiStripRecHitsValid::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: SiStripRecHitsValid.cc:332
SiStripRecHitsValid::matchedRecHits_
edm::InputTag matchedRecHits_
Definition: SiStripRecHitsValid.h:223
LocalError::xy
float xy() const
Definition: LocalError.h:20
SiStripRecHitsValid::meResSasTOB
MonitorElement * meResSasTOB[2]
Definition: SiStripRecHitsValid.h:122
SiStripRecHitsValid::mePullMFRphiTID
MonitorElement * mePullMFRphiTID[3]
Definition: SiStripRecHitsValid.h:140
SiStripRecHitsValid::mePosxRphiTEC
MonitorElement * mePosxRphiTEC[7]
Definition: SiStripRecHitsValid.h:161
SiStripRecHitsValid::meResRphiTOB
MonitorElement * meResRphiTOB[6]
Definition: SiStripRecHitsValid.h:114
SiStripRecHitsValid::mePosxSasTID
MonitorElement * mePosxSasTID[2]
Definition: SiStripRecHitsValid.h:144
SiStripRecHitsValid::meNumMatchedTIB
MonitorElement * meNumMatchedTIB
Definition: SiStripRecHitsValid.h:74
SiStripRecHitsValid::rechitsaspullMF
float rechitsaspullMF[MAXHIT]
Definition: SiStripRecHitsValid.h:205
Topology::measurementError
virtual MeasurementError measurementError(const LocalPoint &, const LocalError &) const =0
LocalError::yy
float yy() const
Definition: LocalError.h:21
PSimHit::localPosition
Local3DPoint localPosition() const
Definition: PSimHit.h:44
SiStripRecHitsValid::meNumTotMatched
MonitorElement * meNumTotMatched
Definition: SiStripRecHitsValid.h:71
SiStripRecHitsValid::meErryMatchedTEC
MonitorElement * meErryMatchedTEC[5]
Definition: SiStripRecHitsValid.h:178
SiStripRecHitsValid::mePullMFSasTID
MonitorElement * mePullMFSasTID[2]
Definition: SiStripRecHitsValid.h:148
SiStripRecHitsValid::meResyMatchedTOB
MonitorElement * meResyMatchedTOB[2]
Definition: SiStripRecHitsValid.h:131
SiStripRecHitsValid::mePullMFRphiTEC
MonitorElement * mePullMFRphiTEC[7]
Definition: SiStripRecHitsValid.h:165
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:28
GloballyPositioned::toLocal
LocalPoint toLocal(const GlobalPoint &gp) const
Definition: GloballyPositioned.h:118
PV3DBase::z
T z() const
Definition: PV3DBase.h:58
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
SiStripRecHitsValid::meErrxSasTIB
MonitorElement * meErrxSasTIB[4]
Definition: SiStripRecHitsValid.h:97
StripGeomDetUnit::topology
virtual const Topology & topology() const
Returns a reference to the strip proxy topology.
Definition: StripGeomDetUnit.cc:17
SiStripRecHitsValid::meChi2SasTOB
MonitorElement * meChi2SasTOB[2]
Definition: SiStripRecHitsValid.h:125
TIDDetId::ring
unsigned int ring() const
ring id
Definition: TIDDetId.h:55
SiStripRecHitsValid::meAdcRphiTOB
MonitorElement * meAdcRphiTOB[6]
Definition: SiStripRecHitsValid.h:111
MeasurementError::uu
float uu() const
Definition: MeasurementError.h:16
SiStripRecHitsValid::meChi2SasTEC
MonitorElement * meChi2SasTEC[5]
Definition: SiStripRecHitsValid.h:174
j
int j
Definition: DBlmapReader.cc:9
SiStripRecHitsValid::meErrxSasTOB
MonitorElement * meErrxSasTOB[2]
Definition: SiStripRecHitsValid.h:121
SiStripRecHitsValid::mePosxMatchedTOB
MonitorElement * mePosxMatchedTOB[2]
Definition: SiStripRecHitsValid.h:126
SiStripRecHitsValid::rechitmatchedz
float rechitmatchedz[MAXHIT]
Definition: SiStripRecHitsValid.h:213
SiStripRecHitsValid::mePosyMatchedTIB
MonitorElement * mePosyMatchedTIB[2]
Definition: SiStripRecHitsValid.h:103
patCandidatesForDimuonsSequences_cff.tracker
tuple tracker
Definition: patCandidatesForDimuonsSequences_cff.py:69
Topology::measurementPosition
virtual MeasurementPoint measurementPosition(const LocalPoint &) const =0
SiStripRecHitsValid::meResRphiTIB
MonitorElement * meResRphiTIB[4]
Definition: SiStripRecHitsValid.h:90
SiStripRecHitsValid::rechitmatchedx
float rechitmatchedx[MAXHIT]
Definition: SiStripRecHitsValid.h:211
SiStripRecHitsValid::meChi2RphiTIB
MonitorElement * meChi2RphiTIB[4]
Definition: SiStripRecHitsValid.h:93
SiStripRecHitsValid::meErrxMatchedTEC
MonitorElement * meErrxMatchedTEC[5]
Definition: SiStripRecHitsValid.h:177
SiStripRecHitsValid::meChi2SasTIB
MonitorElement * meChi2SasTIB[4]
Definition: SiStripRecHitsValid.h:101
DetId::subdetId
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:39
SiStripRecHitsValid::mePosyMatchedTID
MonitorElement * mePosyMatchedTID[2]
Definition: SiStripRecHitsValid.h:152
SiStripRecHitsValid::mePullLFRphiTEC
MonitorElement * mePullLFRphiTEC[7]
Definition: SiStripRecHitsValid.h:164
helper::getRange
std::pair< typename T::DetSet::const_iterator, typename T::DetSet::const_iterator > getRange(const T &detset, const DetId &id)
Definition: SiStripRecHitsValid.cc:40
SiStripRecHitsValid::meAdcSasTIB
MonitorElement * meAdcSasTIB[4]
Definition: SiStripRecHitsValid.h:95
TrackerGeometry::idToDet
virtual const GeomDet * idToDet(DetId) const
Definition: TrackerGeometry.cc:135
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:355
SiStripRecHitsValid::mePosxSasTIB
MonitorElement * mePosxSasTIB[4]
Definition: SiStripRecHitsValid.h:96
dbe_
DQMStore * dbe_
Definition: PFJetBenchmarkAnalyzer.cc:78
SiStripRecHitsValid::meChi2MatchedTOB
MonitorElement * meChi2MatchedTOB[2]
Definition: SiStripRecHitsValid.h:132
gen::k
int k[5][pyjets_maxn]
Definition: Pythia6Hadronizer.cc:86
SiStripRecHitsValid::meResSasTEC
MonitorElement * meResSasTEC[5]
Definition: SiStripRecHitsValid.h:171
SiStripRecHitsValid::mePosxRphiTIB
MonitorElement * mePosxRphiTIB[4]
Definition: SiStripRecHitsValid.h:88
SiStripRecHitsValid::MAXHIT
static const int MAXHIT
Definition: SiStripRecHitsValid.h:189
PSimHit::localDirection
LocalVector localDirection() const
Obsolete. Same as momentumAtEntry().unit(), for backward compatibility.
Definition: PSimHit.h:52
SiStripRecHitsValid::meNstpSasTIB
MonitorElement * meNstpSasTIB[4]
Definition: SiStripRecHitsValid.h:94
SiStripRecHitsValid::SiStripRecHitsValid
SiStripRecHitsValid(const edm::ParameterSet &conf)
Definition: SiStripRecHitsValid.cc:52
SiStripRecHitsValid::mePullLFSasTOB
MonitorElement * mePullLFSasTOB[2]
Definition: SiStripRecHitsValid.h:123
DetId
Definition: DetId.h:20
AlgebraicVector
CLHEP::HepVector AlgebraicVector
Definition: AlgebraicObjects.h:12
SiStripRecHitsValid::meErrxRphiTOB
MonitorElement * meErrxRphiTOB[6]
Definition: SiStripRecHitsValid.h:113
SiStripRecHitsValid::meResyMatchedTEC
MonitorElement * meResyMatchedTEC[5]
Definition: SiStripRecHitsValid.h:180
SiStripRecHitsValid::meErrxRphiTID
MonitorElement * meErrxRphiTID[3]
Definition: SiStripRecHitsValid.h:137
SiStripRecHitsValid::meAdcSasTID
MonitorElement * meAdcSasTID[2]
Definition: SiStripRecHitsValid.h:143
SiStripRecHitsValid::meResxMatchedTIB
MonitorElement * meResxMatchedTIB[2]
Definition: SiStripRecHitsValid.h:106
SiStripRecHitsValid::meErryMatchedTIB
MonitorElement * meErryMatchedTIB[2]
Definition: SiStripRecHitsValid.h:105
SiStripRecHitsValid::meNumRphiTEC
MonitorElement * meNumRphiTEC
Definition: SiStripRecHitsValid.h:81
SiStripRecHitsValid::mePosxSasTOB
MonitorElement * mePosxSasTOB[2]
Definition: SiStripRecHitsValid.h:120
SiStripRecHitsValid::meNstpSasTID
MonitorElement * meNstpSasTID[2]
Definition: SiStripRecHitsValid.h:142
RecHit2DLocalPos::parameters
virtual AlgebraicVector parameters() const
Definition: RecHit2DLocalPos.cc:12
SiStripRecHitsValid::mePullLFSasTID
MonitorElement * mePullLFSasTID[2]
Definition: SiStripRecHitsValid.h:147
SiStripRecHitsValid::meAdcSasTEC
MonitorElement * meAdcSasTEC[5]
Definition: SiStripRecHitsValid.h:168
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
SiStripRecHitsValid::mePullLFRphiTOB
MonitorElement * mePullLFRphiTOB[6]
Definition: SiStripRecHitsValid.h:115
SiStripRecHit2D::cluster
ClusterRef const & cluster() const
Definition: SiStripRecHit2D.h:32
SiStripRecHitsValid::mePullMFSasTEC
MonitorElement * mePullMFSasTEC[5]
Definition: SiStripRecHitsValid.h:173
SiStripRecHitsValid::meNumTotSas
MonitorElement * meNumTotSas
Definition: SiStripRecHitsValid.h:70
SiStripRecHitsValid::meNumSasTEC
MonitorElement * meNumSasTEC
Definition: SiStripRecHitsValid.h:82
SiStripRecHitsValid::rechitrphipullMF
float rechitrphipullMF[MAXHIT]
Definition: SiStripRecHitsValid.h:196
SiStripRecHitsValid::projectHit
std::pair< LocalPoint, LocalVector > projectHit(const PSimHit &hit, const StripGeomDetUnit *stripDet, const BoundPlane &plane)
Definition: SiStripRecHitsValid.cc:946
SiStripRecHitsValid::rechitrphierrx
float rechitrphierrx[MAXHIT]
Definition: SiStripRecHitsValid.h:191
SiStripRecHitsValid::meErryMatchedTOB
MonitorElement * meErryMatchedTOB[2]
Definition: SiStripRecHitsValid.h:129
SiStripRecHitsValid::meChi2MatchedTEC
MonitorElement * meChi2MatchedTEC[5]
Definition: SiStripRecHitsValid.h:181
TIBDetId::layer
unsigned int layer() const
layer id
Definition: TIBDetId.h:41
SiStripRecHitsValid::mePosxMatchedTEC
MonitorElement * mePosxMatchedTEC[5]
Definition: SiStripRecHitsValid.h:175
SiStripRecHitsValid::mePosxSasTEC
MonitorElement * mePosxSasTEC[5]
Definition: SiStripRecHitsValid.h:169
TrackerGeometry::idToDetUnit
virtual const GeomDetUnit * idToDetUnit(DetId) const
Return the pointer to the GeomDetUnit corresponding to a given DetId.
Definition: TrackerGeometry.cc:124
SiStripRecHitsValid::rechitmatchedy
float rechitmatchedy[MAXHIT]
Definition: SiStripRecHitsValid.h:212
SiStripRecHitsValid::conf_
edm::ParameterSet conf_
Definition: SiStripRecHitsValid.h:186
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
SiStripRecHitsValid::stereoRecHits_
edm::InputTag stereoRecHits_
Definition: SiStripRecHitsValid.h:223
SiStripRecHitsValid::mePullMFSasTOB
MonitorElement * mePullMFSasTOB[2]
Definition: SiStripRecHitsValid.h:124
SiStripRecHitsValid::mePullMFRphiTIB
MonitorElement * mePullMFRphiTIB[4]
Definition: SiStripRecHitsValid.h:92
SiStripRecHitsValid::mePosxMatchedTID
MonitorElement * mePosxMatchedTID[2]
Definition: SiStripRecHitsValid.h:151
TECDetId::ring
unsigned int ring() const
ring id
Definition: TECDetId.h:71
AlgebraicSymMatrix
CLHEP::HepSymMatrix AlgebraicSymMatrix
Definition: AlgebraicObjects.h:14
SiStripRecHitsValid::meResxMatchedTOB
MonitorElement * meResxMatchedTOB[2]
Definition: SiStripRecHitsValid.h:130
SiStripRecHitsValid::meResSasTIB
MonitorElement * meResSasTIB[4]
Definition: SiStripRecHitsValid.h:98
TrackerHitAssociator::associateHit
std::vector< PSimHit > associateHit(const TrackingRecHit &thit)
Definition: TrackerHitAssociator.cc:120
SiStripRecHitsValid::mePullMFSasTIB
MonitorElement * mePullMFSasTIB[4]
Definition: SiStripRecHitsValid.h:100
SiStripRecHitsValid::meErrxMatchedTID
MonitorElement * meErrxMatchedTID[2]
Definition: SiStripRecHitsValid.h:153
match
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:6
TrackingRecHit::geographicalId
DetId geographicalId() const
Definition: TrackingRecHit.h:53
dir
dbl *** dir
Definition: mlp_gen.cc:35
x
Definition: DDAxes.h:10
SiStripRecHitsValid::meNumSasTID
MonitorElement * meNumSasTID
Definition: SiStripRecHitsValid.h:79
SiStripRecHitsValid::meResyMatchedTIB
MonitorElement * meResyMatchedTIB[2]
Definition: SiStripRecHitsValid.h:107
SiStripRecHitsValid::meAdcRphiTID
MonitorElement * meAdcRphiTID[3]
Definition: SiStripRecHitsValid.h:135
SiStripRecHitsValid::meNumRphiTOB
MonitorElement * meNumRphiTOB
Definition: SiStripRecHitsValid.h:75
PV2DBase::x
T x() const
Definition: PV2DBase.h:39
T
long double T
Definition: Basic3DVectorLD.h:57
SiStripRecHitsValid::meChi2SasTID
MonitorElement * meChi2SasTID[2]
Definition: SiStripRecHitsValid.h:149
PV3DBase::x
T x() const
Definition: PV3DBase.h:56
SiStripRecHitsValid::meNumMatchedTEC
MonitorElement * meNumMatchedTEC
Definition: SiStripRecHitsValid.h:83
helper::GetDetId::operator()
DetId operator()(const X &x)
Definition: SiStripRecHitsValid.cc:35
GeomDet::surface
virtual const BoundPlane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
SiStripRecHitsValid::meErrxSasTID
MonitorElement * meErrxSasTID[2]
Definition: SiStripRecHitsValid.h:145
SiStripRecHitsValid::mePosxMatchedTIB
MonitorElement * mePosxMatchedTIB[2]
Definition: SiStripRecHitsValid.h:102
GluedGeomDet::stereoDet
const GeomDetUnit * stereoDet() const
Definition: GluedGeomDet.h:21
SiStripRecHitsValid::meResRphiTID
MonitorElement * meResRphiTID[3]
Definition: SiStripRecHitsValid.h:138
SiStripRecHitsValid::meNstpRphiTOB
MonitorElement * meNstpRphiTOB[6]
Definition: SiStripRecHitsValid.h:110
SiStripRecHitsValid::mePosyMatchedTOB
MonitorElement * mePosyMatchedTOB[2]
Definition: SiStripRecHitsValid.h:127
SiStripRecHitsValid::meAdcSasTOB
MonitorElement * meAdcSasTOB[2]
Definition: SiStripRecHitsValid.h:119
Generated for CMSSW Reference Manual by doxygen 1.8.5