65 dbe_->setCurrentFolder(
"Tracking/TrackingRecHits/Strip/ALL");
69 sprintf(histo,
"Errx_LF");
70 meErrxLF = dbe_->book1D(histo,
"RecHit err(x) Local Frame coord.",100,0,0.005);
73 sprintf(histo,
"Errx_MF");
74 meErrxMF = dbe_->book1D(histo,
"RecHit err(x) Meas. Frame coord.",100,0,0.5);
76 sprintf(histo,
"Res_LF");
77 meResLF = dbe_->book1D(histo,
"RecHit Residual",100,-0.02,+0.02);
79 sprintf(histo,
"Res_MF");
80 meResMF = dbe_->book1D(histo,
"RecHit Residual",100,-2,+2);
82 sprintf(histo,
"Pull_LF");
83 mePullLF = dbe_->book1D(histo,
"Pull",100,-5.,5.);
85 sprintf(histo,
"Pull_MF");
86 mePullMF = dbe_->book1D(histo,
"Pull",100,-5.,5.);
88 sprintf(histo,
"Category");
89 meCategory = dbe_->book1D(histo,
"Cluster Category",10,0.,10.);
91 sprintf(histo,
"Trackwidth");
92 meTrackwidth = dbe_->book1D(histo,
"Track width",100,0.,4.);
94 sprintf(histo,
"Expectedwidth");
97 sprintf(histo,
"Clusterwidth");
100 sprintf(histo,
"Trackanglealpha");
103 sprintf(histo,
"Trackanglebeta");
106 sprintf(histo,
"ErrxMFTrackwidthProfile_WClus1");
109 sprintf(histo,
"ErrxMFTrackwidthProfile_WClus2");
112 sprintf(histo,
"ErrxMFTrackwidthProfile_WClus3");
115 sprintf(histo,
"ErrxMFTrackwidthProfile_WClus4");
118 sprintf(histo,
"ResMFTrackwidthProfile_WClus1");
121 sprintf(histo,
"ResMFTrackwidthProfile_WClus2");
123 sprintf(histo,
"ResMFTrackwidthProfile_WClus21");
125 sprintf(histo,
"ResMFTrackwidthProfile_WClus22");
127 sprintf(histo,
"ResMFTrackwidthProfile_WClus23");
132 sprintf(histo,
"ResMFTrackwidthProfile_WClus3");
135 sprintf(histo,
"ResMFTrackwidthProfile_WClus4");
138 sprintf(histo,
"ErrxMFTrackwidthProfile");
141 sprintf(histo,
"ErrxMFTrackwidthProfile_Category1");
143 sprintf(histo,
"ErrxMFTrackwidthProfile_Category2");
146 sprintf(histo,
"ErrxMFTrackwidthProfile_Category3");
149 sprintf(histo,
"ErrxMFTrackwidthProfile_Category4");
152 sprintf(histo,
"ErrxMFClusterwidthProfile_Category1");
155 sprintf(histo,
"ErrxMFAngleProfile");
156 meErrxMFAngleProfile = dbe_->bookProfile(histo,
"Resolution Angle Profile", 100, 0., 60.,100, 0.,2.,
"s");
158 dbe_->setCurrentFolder(
"Tracking/TrackingRecHits/Strip/TIB");
160 for(
int i = 0 ;
i<4 ;
i++) {
162 sprintf(histo,
"Nstp_rphi_layer%dtib",
i+1);
163 meNstpRphiTIB[
i] = dbe_->book1D(histo,
"RecHit Cluster Size",20,0.5,20.5);
164 sprintf(histo,
"Adc_rphi_layer%dtib",
i+1);
165 meAdcRphiTIB[
i] = dbe_->book1D(histo,
"RecHit Cluster Charge",100,0.,300.);
166 sprintf(histo,
"Posx_rphi_layer%dtib",
i+1);
167 mePosxRphiTIB[
i] = dbe_->book1D(histo,
"RecHit x coord.",100,-6.0,+6.0);
169 sprintf(histo,
"Errx_LF_rphi_layer%dtib",
i+1);
170 meErrxLFRphiTIB[
i] = dbe_->book1D(histo,
"RecHit err(x) Local Frame coord.",100,0,0.005);
171 sprintf(histo,
"Errx_MF_rphi_layer%dtib",
i+1);
172 meErrxMFRphiTIB[
i] = dbe_->book1D(histo,
"RecHit err(x) Meas. Frame coord.",100,0,0.5);
174 sprintf(histo,
"Res_LF_rphi_layer%dtib",
i+1);
175 meResLFRphiTIB[
i] = dbe_->book1D(histo,
"RecHit Residual",100,-0.02,+0.02);
176 sprintf(histo,
"Res_MF_rphi_layer%dtib",
i+1);
179 sprintf(histo,
"Pull_LF_rphi_layer%dtib",
i+1);
181 sprintf(histo,
"Pull_MF_rphi_layer%dtib",
i+1);
184 sprintf(histo,
"Trackangle_rphi_layer%dtib",
i+1);
187 sprintf(histo,
"Trackanglebeta_rphi_layer%dtib",
i+1);
190 sprintf(histo,
"Trackangle2_rphi_layer%dtib",
i+1);
193 sprintf(histo,
"PullTrackangleProfile_rphi_layer%dtib",
i+1);
196 sprintf(histo,
"Trackwidth_rphi_layer%dtib",
i+1);
199 sprintf(histo,
"Expectedwidth_rphi_layer%dtib",
i+1);
202 sprintf(histo,
"Clusterwidth_rphi_layer%dtib",
i+1);
205 sprintf(histo,
"Category_rphi_layer%dtib",
i+1);
208 sprintf(histo,
"PullTrackwidthProfile_rphi_layer%dtib",
i+1);
211 sprintf(histo,
"PullTrackwidthProfile_Category1_rphi_layer%dtib",
i+1);
214 sprintf(histo,
"PullTrackwidthProfile_Category2_rphi_layer%dtib",
i+1);
217 sprintf(histo,
"PullTrackwidthProfile_Category3_rphi_layer%dtib",
i+1);
220 sprintf(histo,
"PullTrackwidthProfile_Category4_rphi_layer%dtib",
i+1);
223 sprintf(histo,
"ErrxMFTrackwidthProfile_rphi_layer%dtib",
i+1);
227 sprintf(histo,
"ErrxMFTrackwidthProfile_WClus1_rphi_layer%dtib",
i+1);
230 sprintf(histo,
"ErrxMFTrackwidthProfile_WClus2_rphi_layer%dtib",
i+1);
233 sprintf(histo,
"ErrxMFTrackwidthProfile_WClus3_rphi_layer%dtib",
i+1);
236 sprintf(histo,
"ErrxMFTrackwidthProfile_WClus4_rphi_layer%dtib",
i+1);
239 sprintf(histo,
"ResMFTrackwidthProfile_WClus1_rphi_layer%dtib",
i+1);
241 sprintf(histo,
"ResMFTrackwidthProfile_WClus2_rphi_layer%dtib",
i+1);
244 sprintf(histo,
"ResMFTrackwidthProfile_WClus3_rphi_layer%dtib",
i+1);
247 sprintf(histo,
"ResMFTrackwidthProfile_WClus4_rphi_layer%dtib",
i+1);
252 sprintf(histo,
"ErrxMFTrackwidthProfile_Category1_rphi_layer%dtib",
i+1);
255 sprintf(histo,
"ErrxMFTrackwidthProfile_Category2_rphi_layer%dtib",
i+1);
258 sprintf(histo,
"ErrxMFTrackwidthProfile_Category3_rphi_layer%dtib",
i+1);
261 sprintf(histo,
"ErrxMFTrackwidthProfile_Category4_rphi_layer%dtib",
i+1);
264 sprintf(histo,
"ErrxMFAngleProfile_rphi_layer%dtib",
i+1);
267 sprintf(histo,
"ErrxMFClusterwidthProfile_Category1_rphi_layer%dtib",
i+1);
273 for(
int i = 0 ;
i<2 ;
i++) {
275 sprintf(histo,
"Nstp_sas_layer%dtib",
i+1);
276 meNstpSasTIB[
i] = dbe_->book1D(histo,
"RecHit Cluster Size",20,0.5,20.5);
277 sprintf(histo,
"Adc_sas_layer%dtib",
i+1);
278 meAdcSasTIB[
i] = dbe_->book1D(histo,
"RecHit Cluster Charge",100,0.,300.);
279 sprintf(histo,
"Posx_sas_layer%dtib",
i+1);
280 mePosxSasTIB[
i] = dbe_->book1D(histo,
"RecHit x coord.",100,-6.0,+6.0);
282 sprintf(histo,
"Errx_LF_sas_layer%dtib",
i+1);
283 meErrxLFSasTIB[
i] = dbe_->book1D(histo,
"RecHit err(x) coord.",100,0.,0.005);
284 sprintf(histo,
"Errx_MF_sas_layer%dtib",
i+1);
285 meErrxMFSasTIB[
i] = dbe_->book1D(histo,
"RecHit err(x) coord.",100,0.,0.5);
287 sprintf(histo,
"Res_LF_sas_layer%dtib",
i+1);
288 meResLFSasTIB[
i] = dbe_->book1D(histo,
"RecHit Residual",100,-0.02,+0.02);
289 sprintf(histo,
"Res_MF_sas_layer%dtib",
i+1);
290 meResMFSasTIB[
i] = dbe_->book1D(histo,
"RecHit Residual",100,-2,+2);
292 sprintf(histo,
"Pull_LF_sas_layer%dtib",
i+1);
294 sprintf(histo,
"Pull_MF_sas_layer%dtib",
i+1);
297 sprintf(histo,
"Trackangle_sas_layer%dtib",
i+1);
300 sprintf(histo,
"Trackanglebeta_sas_layer%dtib",
i+1);
303 sprintf(histo,
"PullTrackangleProfile_sas_layer%dtib",
i+1);
306 sprintf(histo,
"Trackwidth_sas_layer%dtib",
i+1);
309 sprintf(histo,
"Expectedwidth_sas_layer%dtib",
i+1);
312 sprintf(histo,
"Clusterwidth_sas_layer%dtib",
i+1);
315 sprintf(histo,
"Category_sas_layer%dtib",
i+1);
318 sprintf(histo,
"PullTrackwidthProfile_sas_layer%dtib",
i+1);
321 sprintf(histo,
"PullTrackwidthProfile_Category1_sas_layer%dtib",
i+1);
324 sprintf(histo,
"PullTrackwidthProfile_Category2_sas_layer%dtib",
i+1);
327 sprintf(histo,
"PullTrackwidthProfile_Category3_sas_layer%dtib",
i+1);
330 sprintf(histo,
"PullTrackwidthProfile_Category4_sas_layer%dtib",
i+1);
333 sprintf(histo,
"ErrxMFTrackwidthProfile_sas_layer%dtib",
i+1);
336 sprintf(histo,
"ErrxMFTrackwidthProfile_Category1_sas_layer%dtib",
i+1);
338 sprintf(histo,
"ErrxMFTrackwidthProfile_Category2_sas_layer%dtib",
i+1);
340 sprintf(histo,
"ErrxMFTrackwidthProfile_Category3_sas_layer%dtib",
i+1);
342 sprintf(histo,
"ErrxMFTrackwidthProfile_Category4_sas_layer%dtib",
i+1);
345 sprintf(histo,
"ErrxMFAngleProfile_sas_layer%dtib",
i+1);
348 sprintf(histo,
"ErrxMFClusterwidthProfile_Category1_sas_layer%dtib",
i+1);
353 sprintf(histo,
"Posx_matched_layer%dtib",
i+1);
355 sprintf(histo,
"Posy_matched_layer%dtib",
i+1);
357 sprintf(histo,
"Errx_matched_layer%dtib",
i+1);
359 sprintf(histo,
"Erry_matched_layer%dtib",
i+1);
361 sprintf(histo,
"Resx_matched_layer%dtib",
i+1);
362 meResxMatchedTIB[
i] = dbe_->book1D(histo,
"RecHit Res(x) coord.",100,-0.02, +0.02);
363 sprintf(histo,
"Resy_matched_layer%dtib",
i+1);
365 sprintf(histo,
"Pullx_matched_layer%dtib",
i+1);
367 sprintf(histo,
"Pully_matched_layer%dtib",
i+1);
371 dbe_->setCurrentFolder(
"Tracking/TrackingRecHits/Strip/TOB");
373 for(
int i = 0 ;
i<6 ;
i++) {
375 sprintf(histo,
"Nstp_rphi_layer%dtob",
i+1);
376 meNstpRphiTOB[
i] = dbe_->book1D(histo,
"RecHit Cluster Size",20,0.5,20.5);
377 sprintf(histo,
"Adc_rphi_layer%dtob",
i+1);
378 meAdcRphiTOB[
i] = dbe_->book1D(histo,
"RecHit Cluster Charge",100,0.,300.);
379 sprintf(histo,
"Posx_rphi_layer%dtob",
i+1);
380 mePosxRphiTOB[
i] = dbe_->book1D(histo,
"RecHit x coord.",100,-6.0,+6.0);
382 sprintf(histo,
"Errx_LF_rphi_layer%dtob",
i+1);
384 sprintf(histo,
"Errx_MF_rphi_layer%dtob",
i+1);
387 sprintf(histo,
"Res_LF_rphi_layer%dtob",
i+1);
388 meResLFRphiTOB[
i] = dbe_->book1D(histo,
"RecHit Residual",100,-0.02,+0.02);
389 sprintf(histo,
"Res_MF_rphi_layer%dtob",
i+1);
392 sprintf(histo,
"Pull_LF_rphi_layer%dtob",
i+1);
394 sprintf(histo,
"Pull_MF_rphi_layer%dtob",
i+1);
397 sprintf(histo,
"Trackangle_rphi_layer%dtob",
i+1);
400 sprintf(histo,
"Trackanglebeta_rphi_layer%dtob",
i+1);
403 sprintf(histo,
"PullTrackangleProfile_rphi_layer%dtob",
i+1);
406 sprintf(histo,
"Trackwidth_rphi_layer%dtob",
i+1);
409 sprintf(histo,
"Expectedwidth_rphi_layer%dtob",
i+1);
412 sprintf(histo,
"Clusterwidth_rphi_layer%dtob",
i+1);
415 sprintf(histo,
"Category_rphi_layer%dtob",
i+1);
418 sprintf(histo,
"PullTrackwidthProfile_rphi_layer%dtob",
i+1);
421 sprintf(histo,
"PullTrackwidthProfile_Category1_rphi_layer%dtob",
i+1);
424 sprintf(histo,
"PullTrackwidthProfile_Category2_rphi_layer%dtob",
i+1);
427 sprintf(histo,
"PullTrackwidthProfile_Category3_rphi_layer%dtob",
i+1);
430 sprintf(histo,
"PullTrackwidthProfile_Category4_rphi_layer%dtob",
i+1);
433 sprintf(histo,
"ErrxMFTrackwidthProfile_rphi_layer%dtob",
i+1);
436 sprintf(histo,
"ErrxMFTrackwidthProfile_WClus1_rphi_layer%dtob",
i+1);
439 sprintf(histo,
"ErrxMFTrackwidthProfile_WClus2_rphi_layer%dtob",
i+1);
442 sprintf(histo,
"ErrxMFTrackwidthProfile_WClus3_rphi_layer%dtob",
i+1);
445 sprintf(histo,
"ErrxMFTrackwidthProfile_WClus4_rphi_layer%dtob",
i+1);
448 sprintf(histo,
"ResMFTrackwidthProfile_WClus1_rphi_layer%dtob",
i+1);
450 sprintf(histo,
"ResMFTrackwidthProfile_WClus2_rphi_layer%dtob",
i+1);
453 sprintf(histo,
"ResMFTrackwidthProfile_WClus3_rphi_layer%dtob",
i+1);
456 sprintf(histo,
"ResMFTrackwidthProfile_WClus4_rphi_layer%dtob",
i+1);
461 sprintf(histo,
"ErrxMFTrackwidthProfile_Category1_rphi_layer%dtob",
i+1);
464 sprintf(histo,
"ErrxMFTrackwidthProfile_Category2_rphi_layer%dtob",
i+1);
467 sprintf(histo,
"ErrxMFTrackwidthProfile_Category3_rphi_layer%dtob",
i+1);
470 sprintf(histo,
"ErrxMFTrackwidthProfile_Category4_rphi_layer%dtob",
i+1);
473 sprintf(histo,
"ErrxMFAngleProfile_rphi_layer%dtob",
i+1);
476 sprintf(histo,
"ErrxMFClusterwidthProfile_Category1_rphi_layer%dtob",
i+1);
482 for(
int i = 0 ;
i<2 ;
i++) {
484 sprintf(histo,
"Nstp_sas_layer%dtob",
i+1);
485 meNstpSasTOB[
i] = dbe_->book1D(histo,
"RecHit Cluster Size",20,0.5,20.5);
486 sprintf(histo,
"Adc_sas_layer%dtob",
i+1);
487 meAdcSasTOB[
i] = dbe_->book1D(histo,
"RecHit Cluster Charge",100,0.,300.);
488 sprintf(histo,
"Posx_sas_layer%dtob",
i+1);
489 mePosxSasTOB[
i] = dbe_->book1D(histo,
"RecHit x coord.",100,-6.0,+6.0);
491 sprintf(histo,
"Errx_LF_sas_layer%dtob",
i+1);
492 meErrxLFSasTOB[
i] = dbe_->book1D(histo,
"RecHit err(x) coord.",100,0.,0.01);
493 sprintf(histo,
"Errx_MF_sas_layer%dtob",
i+1);
494 meErrxMFSasTOB[
i] = dbe_->book1D(histo,
"RecHit err(x) coord.",100,0.,0.5);
496 sprintf(histo,
"Res_LF_sas_layer%dtob",
i+1);
497 meResLFSasTOB[
i] = dbe_->book1D(histo,
"RecHit Residual",100,-0.02,+0.02);
498 sprintf(histo,
"Res_MF_sas_layer%dtob",
i+1);
499 meResMFSasTOB[
i] = dbe_->book1D(histo,
"RecHit Residual",100,-2,2);
501 sprintf(histo,
"Pull_LF_sas_layer%dtob",
i+1);
503 sprintf(histo,
"Pull_MF_sas_layer%dtob",
i+1);
506 sprintf(histo,
"Trackangle_sas_layer%dtob",
i+1);
509 sprintf(histo,
"Trackanglebeta_sas_layer%dtob",
i+1);
512 sprintf(histo,
"PullTrackangleProfile_sas_layer%dtob",
i+1);
515 sprintf(histo,
"Trackwidth_sas_layer%dtob",
i+1);
518 sprintf(histo,
"Expectedwidth_sas_layer%dtob",
i+1);
521 sprintf(histo,
"Clusterwidth_sas_layer%dtob",
i+1);
524 sprintf(histo,
"Category_sas_layer%dtob",
i+1);
527 sprintf(histo,
"PullTrackwidthProfile_sas_layer%dtob",
i+1);
530 sprintf(histo,
"PullTrackwidthProfile_Category1_sas_layer%dtob",
i+1);
533 sprintf(histo,
"PullTrackwidthProfile_Category2_sas_layer%dtob",
i+1);
536 sprintf(histo,
"PullTrackwidthProfile_Category3_sas_layer%dtob",
i+1);
539 sprintf(histo,
"PullTrackwidthProfile_Category4_sas_layer%dtob",
i+1);
542 sprintf(histo,
"ErrxMFTrackwidthProfile_sas_layer%dtob",
i+1);
545 sprintf(histo,
"ErrxMFTrackwidthProfile_Category1_sas_layer%dtob",
i+1);
547 sprintf(histo,
"ErrxMFTrackwidthProfile_Category2_sas_layer%dtob",
i+1);
549 sprintf(histo,
"ErrxMFTrackwidthProfile_Category3_sas_layer%dtob",
i+1);
551 sprintf(histo,
"ErrxMFTrackwidthProfile_Category4_sas_layer%dtob",
i+1);
554 sprintf(histo,
"ErrxMFAngleProfile_sas_layer%dtob",
i+1);
557 sprintf(histo,
"ErrxMFClusterwidthProfile_Category1_sas_layer%dtob",
i+1);
560 sprintf(histo,
"Posx_matched_layer%dtob",
i+1);
562 sprintf(histo,
"Posy_matched_layer%dtob",
i+1);
564 sprintf(histo,
"Errx_matched_layer%dtob",
i+1);
566 sprintf(histo,
"Erry_matched_layer%dtob",
i+1);
568 sprintf(histo,
"Resx_matched_layer%dtob",
i+1);
569 meResxMatchedTOB[
i] = dbe_->book1D(histo,
"RecHit Res(x) coord.",100,-0.02, +0.02);
570 sprintf(histo,
"Resy_matched_layer%dtob",
i+1);
572 sprintf(histo,
"Pullx_matched_layer%dtob",
i+1);
574 sprintf(histo,
"Pully_matched_layer%dtob",
i+1);
578 dbe_->setCurrentFolder(
"Tracking/TrackingRecHits/Strip/TID");
580 for(
int i = 0 ;
i<3 ;
i++) {
582 sprintf(histo,
"Nstp_rphi_layer%dtid",
i+1);
583 meNstpRphiTID[
i] = dbe_->book1D(histo,
"RecHit Cluster Size",20,0.5,20.5);
584 sprintf(histo,
"Adc_rphi_layer%dtid",
i+1);
585 meAdcRphiTID[
i] = dbe_->book1D(histo,
"RecHit Cluster Charge",100,0.,300.);
586 sprintf(histo,
"Posx_rphi_layer%dtid",
i+1);
587 mePosxRphiTID[
i] = dbe_->book1D(histo,
"RecHit x coord.",100,-6.0,+6.0);
588 sprintf(histo,
"Errx_LF_rphi_layer%dtid",
i+1);
590 sprintf(histo,
"Errx_MF_rphi_layer%dtid",
i+1);
592 sprintf(histo,
"Res_LF_rphi_layer%dtid",
i+1);
593 meResLFRphiTID[
i] = dbe_->book1D(histo,
"RecHit Residual",100,-0.5,+0.5);
594 sprintf(histo,
"Res_MF_rphi_layer%dtid",
i+1);
596 sprintf(histo,
"Pull_LF_rphi_layer%dtid",
i+1);
598 sprintf(histo,
"Pull_MF_rphi_layer%dtid",
i+1);
600 sprintf(histo,
"Trackangle_rphi_layer%dtid",
i+1);
603 sprintf(histo,
"Trackanglebeta_rphi_layer%dtid",
i+1);
606 sprintf(histo,
"PullTrackangleProfile_rphi_layer%dtid",
i+1);
609 sprintf(histo,
"Trackwidth_rphi_layer%dtid",
i+1);
612 sprintf(histo,
"Expectedwidth_rphi_layer%dtid",
i+1);
615 sprintf(histo,
"Clusterwidth_rphi_layer%dtid",
i+1);
618 sprintf(histo,
"Category_rphi_layer%dtid",
i+1);
621 sprintf(histo,
"PullTrackwidthProfile_rphi_layer%dtid",
i+1);
624 sprintf(histo,
"PullTrackwidthProfile_Category1_rphi_layer%dtid",
i+1);
627 sprintf(histo,
"PullTrackwidthProfile_Category2_rphi_layer%dtid",
i+1);
630 sprintf(histo,
"PullTrackwidthProfile_Category3_rphi_layer%dtid",
i+1);
633 sprintf(histo,
"PullTrackwidthProfile_Category4_rphi_layer%dtid",
i+1);
636 sprintf(histo,
"ErrxMFTrackwidthProfile_rphi_layer%dtid",
i+1);
639 sprintf(histo,
"ErrxMFTrackwidthProfile_Category1_rphi_layer%dtid",
i+1);
641 sprintf(histo,
"ErrxMFTrackwidthProfile_Category2_rphi_layer%dtid",
i+1);
643 sprintf(histo,
"ErrxMFTrackwidthProfile_Category3_rphi_layer%dtid",
i+1);
645 sprintf(histo,
"ErrxMFTrackwidthProfile_Category4_rphi_layer%dtid",
i+1);
648 sprintf(histo,
"ErrxMFAngleProfile_rphi_layer%dtid",
i+1);
651 sprintf(histo,
"ErrxMFClusterwidthProfile_Category1_rphi_layer%dtid",
i+1);
657 for(
int i = 0 ;
i<2 ;
i++) {
659 sprintf(histo,
"Nstp_sas_layer%dtid",
i+1);
660 meNstpSasTID[
i] = dbe_->book1D(histo,
"RecHit Cluster Size",20,0.5,20.5);
661 sprintf(histo,
"Adc_sas_layer%dtid",
i+1);
662 meAdcSasTID[
i] = dbe_->book1D(histo,
"RecHit Cluster Charge",100,0.,300.);
663 sprintf(histo,
"Posx_sas_layer%dtid",
i+1);
664 mePosxSasTID[
i] = dbe_->book1D(histo,
"RecHit x coord.",100,-6.0,+6.0);
665 sprintf(histo,
"Errx_LF_sas_layer%dtid",
i+1);
666 meErrxLFSasTID[
i] = dbe_->book1D(histo,
"RecHit err(x) coord.",100,0.,0.5);
667 sprintf(histo,
"Errx_MF_sas_layer%dtid",
i+1);
668 meErrxMFSasTID[
i] = dbe_->book1D(histo,
"RecHit err(x) coord.",100,0.,0.5);
669 sprintf(histo,
"Res_LF_sas_layer%dtid",
i+1);
670 meResLFSasTID[
i] = dbe_->book1D(histo,
"RecHit Residual",100,-0.5,+0.5);
671 sprintf(histo,
"Res_MF_sas_layer%dtid",
i+1);
672 meResMFSasTID[
i] = dbe_->book1D(histo,
"RecHit Residual",100,-2,2);
673 sprintf(histo,
"Pull_LF_sas_layer%dtid",
i+1);
675 sprintf(histo,
"Pull_MF_sas_layer%dtid",
i+1);
677 sprintf(histo,
"Trackangle_sas_layer%dtid",
i+1);
679 sprintf(histo,
"Trackanglebeta_sas_layer%dtid",
i+1);
682 sprintf(histo,
"PullTrackangleProfile_sas_layer%dtid",
i+1);
685 sprintf(histo,
"Trackwidth_sas_layer%dtid",
i+1);
688 sprintf(histo,
"Expectedwidth_sas_layer%dtid",
i+1);
691 sprintf(histo,
"Clusterwidth_sas_layer%dtid",
i+1);
694 sprintf(histo,
"Category_sas_layer%dtid",
i+1);
697 sprintf(histo,
"PullTrackwidthProfile_sas_layer%dtid",
i+1);
700 sprintf(histo,
"PullTrackwidthProfile_Category1_sas_layer%dtid",
i+1);
703 sprintf(histo,
"PullTrackwidthProfile_Category2_sas_layer%dtid",
i+1);
706 sprintf(histo,
"PullTrackwidthProfile_Category3_sas_layer%dtid",
i+1);
709 sprintf(histo,
"PullTrackwidthProfile_Category4_sas_layer%dtid",
i+1);
712 sprintf(histo,
"ErrxMFTrackwidthProfile_sas_layer%dtid",
i+1);
715 sprintf(histo,
"ErrxMFTrackwidthProfile_Category1_sas_layer%dtid",
i+1);
717 sprintf(histo,
"ErrxMFTrackwidthProfile_Category2_sas_layer%dtid",
i+1);
719 sprintf(histo,
"ErrxMFTrackwidthProfile_Category3_sas_layer%dtid",
i+1);
721 sprintf(histo,
"ErrxMFTrackwidthProfile_Category4_sas_layer%dtid",
i+1);
724 sprintf(histo,
"ErrxMFAngleProfile_sas_layer%dtid",
i+1);
727 sprintf(histo,
"ErrxMFClusterwidthProfile_Category1_sas_layer%dtid",
i+1);
730 sprintf(histo,
"Posx_matched_layer%dtid",
i+1);
732 sprintf(histo,
"Posy_matched_layer%dtid",
i+1);
734 sprintf(histo,
"Errx_matched_layer%dtid",
i+1);
736 sprintf(histo,
"Erry_matched_layer%dtid",
i+1);
738 sprintf(histo,
"Resx_matched_layer%dtid",
i+1);
740 sprintf(histo,
"Resy_matched_layer%dtid",
i+1);
742 sprintf(histo,
"Pullx_matched_layer%dtid",
i+1);
744 sprintf(histo,
"Pully_matched_layer%dtid",
i+1);
748 dbe_->setCurrentFolder(
"Tracking/TrackingRecHits/Strip/TEC");
750 for(
int i = 0 ;
i<7 ;
i++) {
752 sprintf(histo,
"Nstp_rphi_layer%dtec",
i+1);
753 meNstpRphiTEC[
i] = dbe_->book1D(histo,
"RecHit Cluster Size",20,0.5,20.5);
754 sprintf(histo,
"Adc_rphi_layer%dtec",
i+1);
755 meAdcRphiTEC[
i] = dbe_->book1D(histo,
"RecHit Cluster Charge",100,0.,300.);
756 sprintf(histo,
"Posx_rphi_layer%dtec",
i+1);
757 mePosxRphiTEC[
i] = dbe_->book1D(histo,
"RecHit x coord.",100,-6.0,+6.0);
759 sprintf(histo,
"Errx_LF_rphi_layer%dtec",
i+1);
761 sprintf(histo,
"Errx_MF_rphi_layer%dtec",
i+1);
764 sprintf(histo,
"Res_LF_rphi_layer%dtec",
i+1);
765 meResLFRphiTEC[
i] = dbe_->book1D(histo,
"RecHit Residual",100,-0.5,+0.5);
766 sprintf(histo,
"Res_MF_rphi_layer%dtec",
i+1);
769 sprintf(histo,
"Pull_LF_rphi_layer%dtec",
i+1);
771 sprintf(histo,
"Pull_MF_rphi_layer%dtec",
i+1);
774 sprintf(histo,
"Trackangle_rphi_layer%dtec",
i+1);
777 sprintf(histo,
"Trackanglebeta_rphi_layer%dtec",
i+1);
780 sprintf(histo,
"PullTrackangleProfile_rphi_layer%dtec",
i+1);
783 sprintf(histo,
"Trackwidth_rphi_layer%dtec",
i+1);
786 sprintf(histo,
"Expectedwidth_rphi_layer%dtec",
i+1);
789 sprintf(histo,
"Clusterwidth_rphi_layer%dtec",
i+1);
792 sprintf(histo,
"Category_rphi_layer%dtec",
i+1);
795 sprintf(histo,
"PullTrackwidthProfile_rphi_layer%dtec",
i+1);
798 sprintf(histo,
"PullTrackwidthProfile_Category1_rphi_layer%dtec",
i+1);
801 sprintf(histo,
"PullTrackwidthProfile_Category2_rphi_layer%dtec",
i+1);
804 sprintf(histo,
"PullTrackwidthProfile_Category3_rphi_layer%dtec",
i+1);
807 sprintf(histo,
"PullTrackwidthProfile_Category4_rphi_layer%dtec",
i+1);
810 sprintf(histo,
"ErrxMFTrackwidthProfile_rphi_layer%dtec",
i+1);
813 sprintf(histo,
"ErrxMFTrackwidthProfile_Category1_rphi_layer%dtec",
i+1);
815 sprintf(histo,
"ErrxMFTrackwidthProfile_Category2_rphi_layer%dtec",
i+1);
817 sprintf(histo,
"ErrxMFTrackwidthProfile_Category3_rphi_layer%dtec",
i+1);
819 sprintf(histo,
"ErrxMFTrackwidthProfile_Category4_rphi_layer%dtec",
i+1);
822 sprintf(histo,
"ErrxMFAngleProfile_rphi_layer%dtec",
i+1);
825 sprintf(histo,
"ErrxMFClusterwidthProfile_Category1_rphi_layer%dtec",
i+1);
831 for(
int i = 0 ;
i<5 ;
i++) {
832 if(
i == 0 ||
i == 1 ||
i == 4) {
834 sprintf(histo,
"Nstp_sas_layer%dtec",
i+1);
835 meNstpSasTEC[
i] = dbe_->book1D(histo,
"RecHit Cluster Size",20,0.5,20.5);
836 sprintf(histo,
"Adc_sas_layer%dtec",
i+1);
837 meAdcSasTEC[
i] = dbe_->book1D(histo,
"RecHit Cluster Charge",100,0.,300.);
838 sprintf(histo,
"Posx_sas_layer%dtec",
i+1);
839 mePosxSasTEC[
i] = dbe_->book1D(histo,
"RecHit x coord.",100,-6.0,+6.0);
840 sprintf(histo,
"Errx_LF_sas_layer%dtec",
i+1);
841 meErrxLFSasTEC[
i] = dbe_->book1D(histo,
"RecHit err(x) coord.",100,0.,0.5);
842 sprintf(histo,
"Errx_MF_sas_layer%dtec",
i+1);
843 meErrxMFSasTEC[
i] = dbe_->book1D(histo,
"RecHit err(x) coord.",100,0.,0.5);
844 sprintf(histo,
"Res_LF_sas_layer%dtec",
i+1);
845 meResLFSasTEC[
i] = dbe_->book1D(histo,
"RecHit Residual",100,-0.5,+0.5);
846 sprintf(histo,
"Res_MF_sas_layer%dtec",
i+1);
847 meResMFSasTEC[
i] = dbe_->book1D(histo,
"RecHit Residual",100,-2,+2);
848 sprintf(histo,
"Pull_LF_sas_layer%dtec",
i+1);
850 sprintf(histo,
"Pull_MF_sas_layer%dtec",
i+1);
852 sprintf(histo,
"Trackangle_sas_layer%dtec",
i+1);
854 sprintf(histo,
"Trackanglebeta_sas_layer%dtec",
i+1);
857 sprintf(histo,
"PullTrackangleProfile_sas_layer%dtec",
i+1);
860 sprintf(histo,
"Trackwidth_sas_layer%dtec",
i+1);
863 sprintf(histo,
"Expectedwidth_sas_layer%dtec",
i+1);
866 sprintf(histo,
"Clusterwidth_sas_layer%dtec",
i+1);
869 sprintf(histo,
"Category_sas_layer%dtec",
i+1);
872 sprintf(histo,
"PullTrackwidthProfile_sas_layer%dtec",
i+1);
875 sprintf(histo,
"PullTrackwidthProfile_Category1_sas_layer%dtec",
i+1);
878 sprintf(histo,
"PullTrackwidthProfile_Category2_sas_layer%dtec",
i+1);
881 sprintf(histo,
"PullTrackwidthProfile_Category3_sas_layer%dtec",
i+1);
884 sprintf(histo,
"PullTrackwidthProfile_Category4_sas_layer%dtec",
i+1);
887 sprintf(histo,
"ErrxMFTrackwidthProfile_sas_layer%dtec",
i+1);
890 sprintf(histo,
"ErrxMFTrackwidthProfile_Category1_sas_layer%dtec",
i+1);
892 sprintf(histo,
"ErrxMFTrackwidthProfile_Category2_sas_layer%dtec",
i+1);
894 sprintf(histo,
"ErrxMFTrackwidthProfile_Category3_sas_layer%dtec",
i+1);
896 sprintf(histo,
"ErrxMFTrackwidthProfile_Category4_sas_layer%dtec",
i+1);
899 sprintf(histo,
"ErrxMFAngleProfile_sas_layer%dtec",
i+1);
902 sprintf(histo,
"ErrxMFClusterwidthProfile_Category1_sas_layer%dtec",
i+1);
905 sprintf(histo,
"Posx_matched_layer%dtec",
i+1);
907 sprintf(histo,
"Posy_matched_layer%dtec",
i+1);
909 sprintf(histo,
"Errx_matched_layer%dtec",
i+1);
911 sprintf(histo,
"Erry_matched_layer%dtec",
i+1);
913 sprintf(histo,
"Resx_matched_layer%dtec",
i+1);
915 sprintf(histo,
"Resy_matched_layer%dtec",
i+1);
917 sprintf(histo,
"Pullx_matched_layer%dtec",
i+1);
919 sprintf(histo,
"Pully_matched_layer%dtec",
i+1);
1205 int isrechitrphi = 0;
1206 int isrechitsas = 0;
1207 int isrechitmatched = 0;
1226 float mindist = 999999;
1228 std::vector<PSimHit> matched;
1242 tracker2=&(* estracker);
1264 edm::LogVerbatim(
"TrajectoryAnalyzer") <<
"trajColl->size(): " << trajCollectionHandle->size() ;
1268 for(vector<Trajectory>::const_iterator it = trajCollectionHandle->begin(); it!=trajCollectionHandle->end();it++){
1270 edm::LogVerbatim(
"TrajectoryAnalyzer") <<
"this traj has " << it->foundHits() <<
" valid hits" <<
" , "
1271 <<
"isValid: " << it->isValid() ;
1273 vector<TrajectoryMeasurement> tmColl = it->measurements();
1274 for(vector<TrajectoryMeasurement>::const_iterator itTraj = tmColl.begin(); itTraj!=tmColl.end(); itTraj++){
1275 if(! itTraj->updatedState().isValid())
continue;
1283 if ( itTraj->updatedState().globalMomentum().perp() < 1.)
continue;
1288 DetId detid2 = itTraj->recHit()->geographicalId();
1291 const SiStripMatchedRecHit2D* matchedhit=
dynamic_cast<const SiStripMatchedRecHit2D*
>((*thit2).hit());
1292 const SiStripRecHit2D* hit2d=
dynamic_cast<const SiStripRecHit2D*
>((*thit2).hit());
1293 const SiStripRecHit1D* hit1d=
dynamic_cast<const SiStripRecHit1D*
>((*thit2).hit());
1299 detid = (thit)->geographicalId();
1300 myid=((thit)->geographicalId()).rawId();
1304 isrechitmatched = 0;
1308 isrechitmatched = 1;
1310 position = (thit)->localPosition();
1312 error = (thit)->localPositionError();
1324 std::pair<LocalPoint,LocalVector> closestPair;
1327 if(!matched.empty()){
1331 std::pair<LocalPoint,LocalVector> hitPair;
1332 for(vector<PSimHit>::const_iterator
m=matched.begin();
m<matched.end();
m++){
1337 dist = distx*distx+disty*disty;
1338 if(
sqrt(dist)<mindist){
1340 closestPair = hitPair;
1352 if(isrechitmatched){
1356 int Tibisrechitmatched = isrechitmatched;
1357 int ilay = tTopo->tibLayer(myid) - 1;
1358 if(Tibisrechitmatched>0){
1372 int Tobisrechitmatched = isrechitmatched;
1373 int ilay = tTopo->tobLayer(myid) - 1;
1374 if(Tobisrechitmatched>0){
1387 int Tidisrechitmatched = isrechitmatched;
1388 int ilay = tTopo->tidRing(myid) - 1;
1389 if(Tidisrechitmatched>0){
1402 int Tecisrechitmatched = isrechitmatched;
1403 int ilay = tTopo->tecRing(myid) - 1;
1404 if(Tecisrechitmatched>0){
1461 const SiStripRecHit2D *monohit;
1466 auto hm =matchedhit->monoHit();
1478 if(monotkdir.
z()!=0){
1479 anglealpha = atan(monotkdir.
x()/monotkdir.
z())*180/
TMath::Pi();
1480 anglebeta = atan(monotkdir.
y()/monotkdir.
z())*180/
TMath::Pi();
1487 position = monohit->localPosition();
1488 error = monohit->localPositionError();
1499 float tanalpha =
tan(anglealpha/57.3);
1501 float tanalphaL = drift.
x()/drift.
z();
1510 Wtrack = fabs((thickness/pitch)*tanalpha - (thickness/pitch)*tanalphaL);
1512 float SLorentz = 0.5*(thickness/pitch)*tanalphaL;
1515 int Sp = int(position.
x()/pitch+SLorentz+0.5*Wtrack);
1516 int Sm = int(position.
x()/pitch+SLorentz-0.5*Wtrack);
1523 SiStripRecHit2D::ClusterRef cluster=monohit->cluster();
1527 clusiz = cluster->amplitudes().size();
1529 const std::vector<uint8_t> amplitudes=cluster->amplitudes();
1530 for(
size_t ia=0; ia<amplitudes.size();ia++){
1531 totcharge+=amplitudes[ia];
1547 if(!matched.empty()){
1550 for(vector<PSimHit>::const_iterator
m=matched.begin();
m<matched.end();
m++){
1551 dist =
abs((monohit)->localPosition().
x() - (*m).localPosition().x());
1573 clusterWidth = clusiz;
1575 if (clusterWidth > Wexp + 2) {
1577 }
else if (Wexp == 1) {
1579 }
else if (clusterWidth <= Wexp) {
1587 auto s =matchedhit->stereoHit();
1588 const SiStripRecHit2D *stereohit=&
s;
1591 gtrkdir=gdet->
toGlobal(trackdirection);
1593 if(stereotkdir.
z()!=0){
1594 anglealpha = atan(stereotkdir.
x()/stereotkdir.
z())*180/
TMath::Pi();
1595 anglebeta = atan(stereotkdir.
y()/stereotkdir.
z())*180/
TMath::Pi();
1602 SiStripRecHit2D::ClusterRef cluster=stereohit->cluster();
1608 position = stereohit->localPosition();
1610 error = stereohit->localPositionError();
1621 float tanalpha =
tan(anglealpha/57.3);
1623 float tanalphaL = drift.
x()/drift.
z();
1627 Wtrack = fabs((thickness/pitch)*tanalpha - (thickness/pitch)*tanalphaL);
1629 float SLorentz = 0.5*(thickness/pitch)*tanalphaL;
1631 int Sp = int(position.
x()/pitch+SLorentz+0.5*Wtrack);
1632 int Sm = int(position.
x()/pitch+SLorentz-0.5*Wtrack);
1638 clusiz = cluster->amplitudes().size();
1639 const std::vector<uint8_t> amplitudes=cluster->amplitudes();
1640 for(
size_t ia=0; ia<amplitudes.size();ia++){
1641 totcharge+=amplitudes[ia];
1657 if(!matched.empty()){
1659 for(vector<PSimHit>::const_iterator
m=matched.begin();
m<matched.end();
m++){
1660 dist =
abs((stereohit)->localPosition().
x() - (*m).localPosition().x());
1683 clusterWidth = clusiz;
1685 if (clusterWidth > Wexp + 2) {
1687 }
else if (Wexp == 1) {
1689 }
else if (clusterWidth <= Wexp) {
1708 if (StripSubdet.
stereo() == 0){
1716 float anglealpha = atan(trackdirection.
x()/trackdirection.
z())*180/
TMath::Pi();
1717 float anglebeta = atan(trackdirection.
y()/trackdirection.
z())*180/
TMath::Pi();
1720 SiStripRecHit1D::ClusterRef cluster=hit1d->cluster();
1733 float tanalpha =
tan(anglealpha/57.3);
1735 float tanalphaL = drift.
x()/drift.
z();
1742 Wtrack = fabs((thickness/pitch)*tanalpha - (thickness/pitch)*tanalphaL);
1745 float SLorentz = 0.5*(thickness/pitch)*tanalphaL;
1747 int Sp = int(position.
x()/pitch+SLorentz+0.5*Wtrack);
1748 int Sm = int(position.
x()/pitch+SLorentz-0.5*Wtrack);
1753 clusiz = cluster->amplitudes().size();
1755 const std::vector<uint8_t> amplitudes=cluster->amplitudes();
1756 for(
size_t ia=0; ia<amplitudes.size();ia++){
1757 totcharge+=amplitudes[ia];
1774 if(!matched.empty()){
1776 for(vector<PSimHit>::const_iterator
m=matched.begin();
m<matched.end();
m++){
1777 dist =
abs((hit1d)->localPosition().
x() - (*m).localPosition().x());
1793 clusterWidth = clusiz;
1795 if (clusterWidth > Wexp + 2) {
1797 }
else if (Wexp == 1) {
1799 }
else if (clusterWidth <= Wexp) {
1881 if (StripSubdet.
stereo() == 1){
1890 float anglealpha = atan(trackdirection.
x()/trackdirection.
z())*180/
TMath::Pi();
1891 float anglebeta = atan(trackdirection.
y()/trackdirection.
z())*180/
TMath::Pi();
1894 SiStripRecHit1D::ClusterRef cluster=hit1d->cluster();
1909 float tanalpha =
tan(anglealpha/57.3);
1911 float tanalphaL = drift.
x()/drift.
z();
1915 Wtrack = fabs((thickness/pitch)*tanalpha - (thickness/pitch)*tanalphaL);
1917 float SLorentz = 0.5*(thickness/pitch)*tanalphaL;
1919 int Sp = int(position.
x()/pitch+SLorentz+0.5*Wtrack);
1920 int Sm = int(position.
x()/pitch+SLorentz-0.5*Wtrack);
1925 clusiz = cluster->amplitudes().size();
1926 const std::vector<uint8_t> amplitudes=cluster->amplitudes();
1927 for(
size_t ia=0; ia<amplitudes.size();ia++){
1928 totcharge+=amplitudes[ia];
1944 if(!matched.empty()){
1946 for(vector<PSimHit>::const_iterator
m=matched.begin();
m<matched.end();
m++){
1947 dist =
abs((hit1d)->localPosition().
x() - (*m).localPosition().x());
1965 clusterWidth = clusiz;
1967 if (clusterWidth > Wexp + 2) {
1969 }
else if (Wexp == 1) {
1971 }
else if (clusterWidth <= Wexp) {
1985 if (StripSubdet.
stereo() == 0){
1993 float anglealpha = atan(trackdirection.
x()/trackdirection.
z())*180/
TMath::Pi();
1994 float anglebeta = atan(trackdirection.
y()/trackdirection.
z())*180/
TMath::Pi();
1996 SiStripRecHit2D::ClusterRef cluster=hit2d->cluster();
2009 float tanalpha =
tan(anglealpha/57.3);
2011 float tanalphaL = drift.
x()/drift.
z();
2018 Wtrack = fabs((thickness/pitch)*tanalpha - (thickness/pitch)*tanalphaL);
2021 float SLorentz = 0.5*(thickness/pitch)*tanalphaL;
2023 int Sp = int(position.
x()/pitch+SLorentz+0.5*Wtrack);
2024 int Sm = int(position.
x()/pitch+SLorentz-0.5*Wtrack);
2029 clusiz = cluster->amplitudes().size();
2031 const std::vector<uint8_t> amplitudes=cluster->amplitudes();
2032 for(
size_t ia=0; ia<amplitudes.size();ia++){
2033 totcharge+=amplitudes[ia];
2050 if(!matched.empty()){
2052 for(vector<PSimHit>::const_iterator
m=matched.begin();
m<matched.end();
m++){
2053 dist =
abs((hit2d)->localPosition().
x() - (*m).localPosition().x());
2069 clusterWidth = clusiz;
2071 if (clusterWidth > Wexp + 2) {
2073 }
else if (Wexp == 1) {
2075 }
else if (clusterWidth <= Wexp) {
2157 if (StripSubdet.
stereo() == 1){
2166 float anglealpha = atan(trackdirection.
x()/trackdirection.
z())*180/
TMath::Pi();
2167 float anglebeta = atan(trackdirection.
y()/trackdirection.
z())*180/
TMath::Pi();
2169 SiStripRecHit2D::ClusterRef cluster=hit2d->cluster();
2184 float tanalpha =
tan(anglealpha/57.3);
2186 float tanalphaL = drift.
x()/drift.
z();
2190 Wtrack = fabs((thickness/pitch)*tanalpha - (thickness/pitch)*tanalphaL);
2192 float SLorentz = 0.5*(thickness/pitch)*tanalphaL;
2194 int Sp = int(position.
x()/pitch+SLorentz+0.5*Wtrack);
2195 int Sm = int(position.
x()/pitch+SLorentz-0.5*Wtrack);
2200 clusiz = cluster->amplitudes().size();
2201 const std::vector<uint8_t> amplitudes=cluster->amplitudes();
2202 for(
size_t ia=0; ia<amplitudes.size();ia++){
2203 totcharge+=amplitudes[ia];
2219 if(!matched.empty()){
2221 for(vector<PSimHit>::const_iterator
m=matched.begin();
m<matched.end();
m++){
2222 dist =
abs((hit2d)->localPosition().
x() - (*m).localPosition().x());
2240 clusterWidth = clusiz;
2242 if (clusterWidth > Wexp + 2) {
2244 }
else if (Wexp == 1) {
2246 }
else if (clusterWidth <= Wexp) {
2260 float CutThickness=0.04;
2263 if(isrechitrphi>0 || isrechitsas>0){
2266 if (isrechitrphi>0) {
2418 if (isrechitsas>0) {
2553 int Tibisrechitrphi = isrechitrphi;
2554 int Tibisrechitsas = isrechitsas;
2556 int ilay = tTopo->tibLayer(myid) - 1;
2558 if(Tibisrechitrphi!=0){
2624 if(Tibisrechitsas!=0){
2673 int Tobisrechitrphi = isrechitrphi;
2674 int Tobisrechitsas = isrechitsas;
2675 int ilay = tTopo->tobLayer(myid) - 1;
2676 if(Tobisrechitrphi!=0){
2741 if(Tobisrechitsas!=0){
2791 int Tidisrechitrphi = isrechitrphi;
2792 int Tidisrechitsas = isrechitsas;
2793 int ilay = tTopo->tidRing(myid) - 1;
2794 if(Tidisrechitrphi!=0){
2840 if(Tidisrechitsas!=0){
2889 int Tecisrechitrphi = isrechitrphi;
2890 int Tecisrechitsas = isrechitsas;
2891 int ilay = tTopo->tecRing(myid) - 1;
2892 if(Tecisrechitrphi!=0){
2938 if(Tecisrechitsas!=0){
3007 LocalPoint localHit = plane.toLocal(globalpos);
3014 float scale = -localHit.
z() / dir.
z();
3026 return std::pair<LocalPoint,LocalVector>( projectedPos, localStripDir);
MonitorElement * mePullTrackwidthProfileCategory1SasTEC[5]
MonitorElement * meErrxMFTrackwidthProfileRphiTEC[7]
MonitorElement * mePullTrackwidthProfileSasTOB[2]
MonitorElement * meResyMatchedTEC[5]
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
MonitorElement * mePullTrackwidthProfileCategory2RphiTID[3]
MonitorElement * mePullTrackwidthProfileRphiTID[3]
T getParameter(std::string const &) const
SiStripTrackingRecHitsValid(const edm::ParameterSet &conf)
T getUntrackedParameter(std::string const &, T const &) const
virtual ~SiStripTrackingRecHitsValid()
MonitorElement * meErrxMFTrackwidthProfileCategory2RphiTIB[4]
MonitorElement * mePullMFSasTOB[2]
MonitorElement * mePullMFRphiTEC[7]
MonitorElement * meTrackangleRphiTEC[7]
MonitorElement * meExpectedwidthRphiTEC[7]
MonitorElement * meResyMatchedTIB[2]
MonitorElement * meCategoryRphiTIB[4]
MonitorElement * meErrxMatchedTID[2]
MonitorElement * meTrackwidthRphiTID[3]
MonitorElement * meErrxMatchedTEC[5]
MonitorElement * mePullTrackwidthProfileCategory3SasTIB[4]
MonitorElement * meErrxMFClusterwidthProfileCategory1RphiTIB[4]
MonitorElement * meExpectedwidthSasTOB[2]
MonitorElement * meErrxMFTrackwidthProfileWClus2
MonitorElement * mePullTrackwidthProfileCategory4RphiTEC[7]
MonitorElement * mePullTrackwidthProfileCategory2SasTID[2]
MonitorElement * meResMFTrackwidthProfileWclus4RphiTOB[6]
MonitorElement * meErrxMFTrackwidthProfileCategory3SasTOB[2]
MonitorElement * mePullyMatchedTIB[2]
MonitorElement * meAdcRphiTID[3]
MonitorElement * meErrxMFClusterwidthProfileCategory1SasTID[2]
MonitorElement * meResyMatchedTOB[2]
MonitorElement * mePullTrackwidthProfileRphiTEC[7]
virtual float stripAngle(float strip) const =0
const GeomDetUnit * monoDet() const
MonitorElement * mePosyMatchedTEC[5]
MonitorElement * meResMFSasTOB[2]
MonitorElement * meErrxMFTrackwidthProfileCategory4RphiTIB[4]
MonitorElement * meResLFSasTEC[5]
MonitorElement * meErrxMFTrackwidthProfileCategory1SasTID[2]
MonitorElement * mePosxMatchedTIB[2]
MonitorElement * meTrackwidthSasTIB[4]
MonitorElement * meTrackwidthRphiTEC[7]
MonitorElement * meResLFSasTID[2]
MonitorElement * meNstpSasTOB[2]
MonitorElement * mePullTrackwidthProfileCategory1RphiTOB[6]
MonitorElement * meExpectedwidthSasTID[2]
MonitorElement * meClusterwidthSasTIB[4]
LocalVector localDirection() const
MonitorElement * meErrxMFTrackwidthProfileCategory2SasTID[2]
MonitorElement * meTrackangleSasTIB[4]
MonitorElement * meExpectedwidthSasTEC[5]
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
MonitorElement * mePosxMatchedTID[2]
MonitorElement * meNstpSasTIB[4]
MonitorElement * mePullxMatchedTID[2]
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")
MonitorElement * mePosxRphiTOB[6]
Sin< T >::type sin(const T &t)
float rechitsastrackanglebeta
MonitorElement * mePullTrackangleProfileSasTEC[5]
MonitorElement * meErryMatchedTOB[2]
LocalVector drift(const StripGeomDetUnit *, const MagneticField &, const SiStripLorentzAngle &)
MonitorElement * meNstpRphiTIB[4]
MonitorElement * mePosxMatchedTOB[2]
MonitorElement * meErrxMFTrackwidthProfileCategory4SasTIB[4]
MonitorElement * mePullTrackangleProfileRphiTEC[7]
MonitorElement * meErrxMFTrackwidthProfileWclus4RphiTOB[6]
MonitorElement * mePullTrackwidthProfileCategory1RphiTEC[7]
MonitorElement * mePullTrackwidthProfileCategory4SasTEC[5]
MonitorElement * meExpectedwidthRphiTOB[6]
MonitorElement * mePullTrackwidthProfileCategory2SasTEC[5]
MonitorElement * meClusterwidth
MonitorElement * mePullLFRphiTOB[6]
MonitorElement * meNstpRphiTID[3]
MonitorElement * meResxMatchedTOB[2]
MonitorElement * mePullMF
MonitorElement * meErrxMFTrackwidthProfileCategory4SasTOB[2]
MonitorElement * meResMFRphiTEC[7]
LocalPoint toLocal(const GlobalPoint &gp) const
Conversion to the R.F. of the GeomDet.
MonitorElement * meErrxMFTrackwidthProfileWclus4RphiTIB[4]
MonitorElement * meErrxMFTrackwidthProfileWclus3RphiTOB[6]
MonitorElement * meCategoryRphiTID[3]
const Bounds & bounds() const
MonitorElement * meErrxMatchedTOB[2]
MonitorElement * meTrackanglealpha
MonitorElement * meResMFTrackwidthProfileWclus4RphiTIB[4]
MonitorElement * mePullTrackangleProfileRphiTID[3]
MonitorElement * meErrxMFAngleProfile
MonitorElement * meResMFSasTID[2]
MonitorElement * meCategory
MonitorElement * meTrackangle2RphiTIB[4]
const StripTopology * topol
virtual const StripTopology & specificTopology() const
Returns a reference to the strip proxy topology.
MonitorElement * meErrxMFTrackwidthProfileCategory3RphiTEC[7]
MonitorElement * mePullLFSasTEC[5]
MonitorElement * meTrackanglebetaRphiTID[3]
const Plane & surface() const
The nominal surface of the GeomDet.
MonitorElement * meErrxMFTrackwidthProfileCategory1SasTIB[4]
MonitorElement * meTrackanglebetaSasTEC[5]
static int position[TOTALCHAMBERS][3]
MonitorElement * meResLFSasTIB[4]
MonitorElement * meTrackwidth
MonitorElement * meResMFTrackwidthProfileWclus3RphiTIB[4]
MonitorElement * meAdcRphiTEC[7]
MonitorElement * meTrackwidthSasTID[2]
MonitorElement * meResMFTrackwidthProfileWClus1
MonitorElement * mePullyMatchedTID[2]
MonitorElement * mePullTrackwidthProfileCategory4RphiTOB[6]
MonitorElement * meErrxMFTrackwidthProfileCategory1SasTOB[2]
MonitorElement * meCategorySasTEC[5]
MonitorElement * meResxMatchedTEC[5]
MonitorElement * meTrackangleRphiTOB[6]
int rechitrphiexpectedwidth
MonitorElement * mePosxRphiTEC[7]
MonitorElement * meErrxMFTrackwidthProfileWclus1RphiTOB[6]
virtual float strip(const LocalPoint &) const =0
MonitorElement * mePullTrackwidthProfileCategory3RphiTID[3]
MonitorElement * meErrxLFSasTID[2]
MonitorElement * meResMFTrackwidthProfileWclus3RphiTOB[6]
virtual float thickness() const =0
MonitorElement * meResMFTrackwidthProfileWclus2RphiTOB[6]
MonitorElement * mePosxSasTIB[4]
MonitorElement * meTrackangleRphiTIB[4]
MonitorElement * mePullTrackwidthProfileCategory4RphiTIB[4]
MonitorElement * mePullTrackwidthProfileCategory4RphiTID[3]
MonitorElement * mePosxRphiTID[3]
MonitorElement * meErrxMFClusterwidthProfileCategory1RphiTEC[7]
MonitorElement * meErrxMFClusterwidthProfileCategory1SasTEC[7]
float rechitsastrackwidth
MonitorElement * meErrxMFTrackwidthProfileCategory3RphiTOB[6]
virtual float localPitch(const LocalPoint &) const =0
MonitorElement * meResMFTrackwidthProfileWClus2
const MagneticField * magfield2_
MonitorElement * meTrackwidthRphiTOB[6]
MonitorElement * meResLFRphiTEC[7]
MonitorElement * mePullTrackwidthProfileCategory1SasTOB[2]
MonitorElement * meTrackwidthRphiTIB[4]
MonitorElement * meErrxMFTrackwidthProfileCategory3
MonitorElement * meErrxLFSasTEC[5]
MonitorElement * meErrxMFTrackwidthProfileRphiTOB[6]
MonitorElement * meErrxMFAngleProfileSasTOB[2]
MonitorElement * mePosyMatchedTIB[2]
MonitorElement * mePullTrackwidthProfileCategory4SasTID[2]
MonitorElement * meResMFSasTIB[4]
virtual MeasurementError measurementError(const LocalPoint &, const LocalError &) const =0
MonitorElement * meErrxMFTrackwidthProfileCategory3SasTEC[7]
MonitorElement * meErrxMFTrackwidthProfileSasTIB[4]
MonitorElement * meTrackwidthSasTOB[2]
MonitorElement * meAdcRphiTOB[6]
Local3DPoint localPosition() const
MonitorElement * mePullLFSasTID[2]
MonitorElement * meErryMatchedTIB[2]
MonitorElement * mePullTrackwidthProfileCategory1RphiTIB[4]
MonitorElement * meResMFTrackwidthProfileWClus3
MonitorElement * meErrxMFAngleProfileSasTIB[4]
MonitorElement * meCategoryRphiTEC[7]
MonitorElement * mePullyMatchedTOB[2]
MonitorElement * meErrxLFSasTIB[4]
MonitorElement * meResLFRphiTIB[4]
MonitorElement * mePullTrackwidthProfileCategory2SasTIB[4]
MonitorElement * mePullTrackangleProfileSasTIB[4]
MonitorElement * mePullTrackwidthProfileCategory2RphiTOB[6]
MonitorElement * mePullTrackwidthProfileCategory1SasTID[2]
MonitorElement * meErrxMF
MonitorElement * meResMFRphiTOB[6]
MonitorElement * meErrxMFRphiTID[3]
MonitorElement * meAdcSasTIB[4]
Cos< T >::type cos(const T &t)
MonitorElement * meErrxMatchedTIB[2]
MonitorElement * meTrackanglebeta
MonitorElement * mePullTrackwidthProfileCategory3RphiTEC[7]
MonitorElement * meTrackangleSasTOB[2]
MonitorElement * mePullTrackangleProfileSasTID[2]
MonitorElement * meTrackwidthSasTEC[5]
virtual const Topology & topology() const
Returns a reference to the strip proxy topology.
MonitorElement * meErrxMFTrackwidthProfileCategory4RphiTEC[7]
MonitorElement * meErrxLFRphiTIB[4]
MonitorElement * meResLFRphiTID[3]
MonitorElement * meErrxMFSasTEC[5]
Tan< T >::type tan(const T &t)
float rechitsastrackangle
MonitorElement * meResxMatchedTIB[2]
MonitorElement * meErrxLFRphiTEC[7]
MonitorElement * mePullLF
MonitorElement * mePullTrackwidthProfileCategory3RphiTOB[6]
MonitorElement * meErryMatchedTID[2]
MonitorElement * meErrxMFTrackwidthProfileSasTOB[2]
MonitorElement * meErrxMFTrackwidthProfileCategory2
virtual MeasurementPoint measurementPosition(const LocalPoint &) const =0
MonitorElement * meErrxMFTrackwidthProfileWclus2RphiTIB[4]
MonitorElement * meClusterwidthRphiTID[3]
MonitorElement * meErrxMFTrackwidthProfileCategory4SasTEC[7]
MonitorElement * meAdcSasTOB[2]
MonitorElement * meAdcSasTEC[5]
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
MonitorElement * mePullTrackwidthProfileCategory2RphiTIB[4]
MonitorElement * mePullMFRphiTIB[4]
MonitorElement * mePullTrackwidthProfileRphiTOB[6]
MonitorElement * mePosxSasTID[2]
int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum) ...
MonitorElement * mePullMFSasTID[2]
MonitorElement * meClusterwidthSasTEC[5]
virtual const GeomDet * idToDet(DetId) const
MonitorElement * meResyMatchedTID[2]
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
MonitorElement * mePullTrackwidthProfileCategory4SasTOB[2]
MonitorElement * meTrackanglebetaSasTIB[4]
MonitorElement * mePullTrackangleProfileRphiTIB[4]
std::pair< LocalPoint, LocalVector > projectHit(const PSimHit &hit, const StripGeomDetUnit *stripDet, const BoundPlane &plane)
MonitorElement * meErryMatchedTEC[5]
MonitorElement * mePullLFRphiTID[3]
float rechitrphithickness
MonitorElement * mePullTrackwidthProfileCategory2RphiTEC[7]
MonitorElement * meErrxMFTrackwidthProfileCategory1SasTEC[7]
MonitorElement * meErrxMFClusterwidthProfileCategory1
MonitorElement * meNstpRphiTEC[7]
MonitorElement * meNstpSasTID[2]
MonitorElement * meResMFRphiTIB[4]
MonitorElement * meErrxMFSasTIB[4]
LocalVector localDirection() const
Obsolete. Same as momentumAtEntry().unit(), for backward compatibility.
MonitorElement * meErrxMFAngleProfileRphiTEC[7]
MonitorElement * mePullyMatchedTEC[5]
MonitorElement * mePullLFRphiTIB[4]
float rechitrphitrackwidth
MonitorElement * meResMFTrackwidthProfileWclus2RphiTIB[4]
MonitorElement * meErrxMFSasTOB[2]
MonitorElement * meTrackangleSasTID[2]
MonitorElement * mePullTrackwidthProfileCategory1RphiTID[3]
MonitorElement * meErrxMFTrackwidthProfileCategory1
unsigned int stereo() const
stereo
MonitorElement * meErrxMFTrackwidthProfileWclus1RphiTIB[4]
MonitorElement * meErrxMFTrackwidthProfileWclus3RphiTIB[4]
MonitorElement * mePullLFSasTOB[2]
MonitorElement * meErrxMFTrackwidthProfile
float rechitrphitrackangle
MonitorElement * mePullLFSasTIB[4]
MonitorElement * meErrxMFTrackwidthProfileCategory4
MonitorElement * mePosxRphiTIB[4]
MonitorElement * meResMFTrackwidthProfileWClus22
MonitorElement * meResMFTrackwidthProfileWclus1RphiTIB[4]
MonitorElement * mePullxMatchedTOB[2]
MonitorElement * mePullTrackwidthProfileCategory4SasTIB[4]
MonitorElement * mePullLFRphiTEC[7]
MonitorElement * meErrxMFRphiTEC[7]
MonitorElement * meErrxMFTrackwidthProfileCategory3RphiTID[3]
virtual LocalError localPositionError() const =0
MonitorElement * meErrxMFAngleProfileRphiTID[3]
float rechitrphitrackanglebeta
MonitorElement * mePullMFRphiTID[3]
MonitorElement * meClusterwidthRphiTOB[6]
MonitorElement * mePullTrackwidthProfileCategory3SasTID[2]
MonitorElement * meExpectedwidthRphiTIB[4]
MonitorElement * meCategorySasTOB[2]
MonitorElement * meErrxMFTrackwidthProfileRphiTIB[4]
MonitorElement * meErrxMFTrackwidthProfileSasTEC[7]
MonitorElement * meErrxMFClusterwidthProfileCategory1SasTIB[4]
MonitorElement * meClusterwidthRphiTEC[7]
MonitorElement * meResMFRphiTID[3]
MonitorElement * mePullTrackwidthProfileCategory3RphiTIB[4]
MonitorElement * meResMFTrackwidthProfileWclus1RphiTOB[6]
MonitorElement * meErrxMFTrackwidthProfileCategory2SasTOB[2]
MonitorElement * meErrxMFTrackwidthProfileCategory2SasTEC[7]
MonitorElement * mePullTrackwidthProfileSasTID[2]
MonitorElement * meErrxLFRphiTID[3]
MonitorElement * meErrxMFTrackwidthProfileCategory4RphiTOB[6]
MonitorElement * meErrxMFTrackwidthProfileWclus2RphiTOB[6]
MonitorElement * mePosxSasTEC[5]
MonitorElement * meErrxMFTrackwidthProfileCategory3RphiTIB[4]
MonitorElement * mePullTrackwidthProfileCategory1SasTIB[4]
MonitorElement * meClusterwidthSasTID[2]
MonitorElement * meErrxMFAngleProfileSasTID[2]
MonitorElement * mePosxMatchedTEC[5]
MonitorElement * meErrxMFSasTID[2]
MonitorElement * mePosyMatchedTID[2]
MonitorElement * meErrxMFTrackwidthProfileSasTID[2]
MonitorElement * mePullTrackwidthProfileSasTEC[5]
MonitorElement * mePullxMatchedTEC[5]
virtual const GeomDetUnit * idToDetUnit(DetId) const
Return the pointer to the GeomDetUnit corresponding to a given DetId.
MonitorElement * meCategoryRphiTOB[6]
MonitorElement * mePullTrackwidthProfileSasTIB[4]
MonitorElement * meCategorySasTIB[4]
MonitorElement * meTrackanglebetaRphiTOB[6]
MonitorElement * mePullMFSasTIB[4]
MonitorElement * meErrxMFTrackwidthProfileCategory2RphiTOB[6]
MonitorElement * meErrxMFTrackwidthProfileCategory3SasTIB[4]
MonitorElement * meNstpSasTEC[5]
MonitorElement * meErrxMFTrackwidthProfileRphiTID[3]
MonitorElement * meResMFTrackwidthProfileWClus4
MonitorElement * meResMFTrackwidthProfileWClus23
MonitorElement * mePullxMatchedTIB[2]
MonitorElement * meErrxMFClusterwidthProfileCategory1RphiTOB[6]
MonitorElement * meTrackangleSasTEC[5]
MonitorElement * meErrxMFTrackwidthProfileCategory1RphiTIB[4]
MonitorElement * meResxMatchedTID[2]
MonitorElement * meExpectedwidthSasTIB[4]
MonitorElement * mePullTrackwidthProfileCategory3SasTOB[2]
MonitorElement * meErrxLFRphiTOB[6]
std::vector< PSimHit > associateHit(const TrackingRecHit &thit)
MonitorElement * meTrackangleRphiTID[3]
MonitorElement * meAdcSasTID[2]
MonitorElement * meErrxMFTrackwidthProfileCategory1RphiTEC[7]
MonitorElement * meErrxMFTrackwidthProfileCategory2RphiTEC[7]
MonitorElement * meClusterwidthRphiTIB[4]
MonitorElement * meErrxMFTrackwidthProfileCategory1RphiTOB[6]
MonitorElement * meErrxLF
MonitorElement * meErrxMFRphiTOB[6]
MonitorElement * meResLFRphiTOB[6]
MonitorElement * meTrackanglebetaRphiTIB[4]
MonitorElement * meErrxMFTrackwidthProfileWClus3
MonitorElement * meCategorySasTID[2]
MonitorElement * meNstpRphiTOB[6]
MonitorElement * meTrackanglebetaRphiTEC[7]
MonitorElement * meErrxMFTrackwidthProfileCategory4SasTID[2]
MonitorElement * meErrxMFTrackwidthProfileWClus4
MonitorElement * meResMFTrackwidthProfileWClus21
MonitorElement * meErrxMFClusterwidthProfileCategory1RphiTID[3]
MonitorElement * mePullMFRphiTOB[6]
MonitorElement * meErrxMFRphiTIB[4]
MonitorElement * meErrxMFTrackwidthProfileWClus1
MonitorElement * meResLFSasTOB[2]
MonitorElement * mePullTrackwidthProfileCategory2SasTOB[2]
MonitorElement * meClusterwidthSasTOB[2]
MonitorElement * meTrackanglebetaSasTID[2]
MonitorElement * meErrxMFAngleProfileRphiTIB[4]
virtual LocalPoint localPosition() const =0
MonitorElement * mePullTrackangleProfileSasTOB[2]
MonitorElement * mePullMFSasTEC[5]
MonitorElement * meErrxMFTrackwidthProfileCategory1RphiTID[3]
MonitorElement * meErrxMFTrackwidthProfileCategory3SasTID[2]
MonitorElement * meExpectedwidth
float rechitrphitrackangle2
MonitorElement * mePullTrackangleProfileRphiTOB[6]
MonitorElement * meErrxMFTrackwidthProfileCategory2RphiTID[3]
MonitorElement * meAdcRphiTIB[4]
MonitorElement * meErrxMFAngleProfileRphiTOB[6]
MonitorElement * meErrxLFSasTOB[2]
const GeomDetUnit * stereoDet() const
MonitorElement * meExpectedwidthRphiTID[3]
MonitorElement * mePullTrackwidthProfileCategory3SasTEC[5]
MonitorElement * meResMFSasTEC[5]
MonitorElement * mePosyMatchedTOB[2]
MonitorElement * meErrxMFTrackwidthProfileCategory2SasTIB[4]
MonitorElement * meErrxMFClusterwidthProfileCategory1SasTOB[2]
MonitorElement * meErrxMFAngleProfileSasTEC[5]
MonitorElement * meTrackanglebetaSasTOB[2]
MonitorElement * mePullTrackwidthProfileRphiTIB[4]
int rechitsasexpectedwidth
MonitorElement * mePosxSasTOB[2]
MonitorElement * meErrxMFTrackwidthProfileCategory4RphiTID[3]