CMS 3D CMS Logo
AlignmentMonitorSegmentDifferences.cc
Go to the documentation of this file.
1 /*
2  * Package: CommonAlignmentProducer
3  * Class : AlignmentMonitorSegmentDifferences
4  *
5  * Original Author: Jim Pivarski
6  * Created: Mon Nov 12 13:30:14 CST 2007
7  *
8  * $Id: AlignmentMonitorSegmentDifferences.cc,v 1.5 2011/04/15 23:09:38 khotilov Exp $
9  */
10 
11 #include "Alignment/CommonAlignmentMonitor/interface/AlignmentMonitorPluginFactory.h"
12 #include "Alignment/CommonAlignmentMonitor/interface/AlignmentMonitorBase.h"
13 #include "Alignment/MuonAlignmentAlgorithms/interface/MuonResidualsFromTrack.h"
14 #include "Alignment/MuonAlignmentAlgorithms/interface/MuonResidualsPositionFitter.h"
15 #include "Alignment/MuonAlignmentAlgorithms/interface/MuonResidualsAngleFitter.h"
16 #include "Alignment/MuonAlignmentAlgorithms/interface/MuonResidualsTwoBin.h"
17 
18 #include "FWCore/ParameterSet/interface/ParameterSet.h"
19 #include "FWCore/Utilities/interface/InputTag.h"
20 #include "FWCore/MessageLogger/interface/MessageLogger.h"
21 #include "FWCore/Framework/interface/EventSetup.h"
22 #include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"
23 #include "Geometry/Records/interface/GlobalTrackingGeometryRecord.h"
24 
25 #include <sstream>
26 
27 #include "TrackingTools/Records/interface/TrackingComponentsRecord.h"
28 #include "TrackingTools/Records/interface/DetIdAssociatorRecord.h"
29 #include "TrackingTools/GeomPropagators/interface/Propagator.h"
30 #include "TrackingTools/TrackAssociator/interface/DetIdAssociator.h"
31 #include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
32 #include "MagneticField/Engine/interface/MagneticField.h"
33 
34 class AlignmentMonitorSegmentDifferences : public AlignmentMonitorBase {
35 public:
36  AlignmentMonitorSegmentDifferences(const edm::ParameterSet &cfg);
37  ~AlignmentMonitorSegmentDifferences() override {}
38 
39  void book() override;
40 
41  void event(const edm::Event &iEvent,
42  const edm::EventSetup &iSetup,
43  const ConstTrajTrackPairCollection &iTrajTracks) override;
44  void processMuonResidualsFromTrack(MuonResidualsFromTrack &mrft);
45 
46 private:
47  edm::InputTag m_muonCollectionTag;
48  double m_minTrackPt;
49  double m_minTrackP;
50  double m_maxDxy;
51  int m_minTrackerHits;
52  double m_maxTrackerRedChi2;
53  bool m_allowTIDTEC;
54  bool m_minNCrossedChambers;
55  int m_minDT13Hits;
56  int m_minDT2Hits;
57  int m_minCSCHits;
58  bool m_doDT;
59  bool m_doCSC;
60 
61  // wheel, sector, stationdiff
62  TProfile *m_dt13_resid[5][12][3];
63  TProfile *m_dt13_slope[5][12][3];
64  TProfile *m_dt2_resid[5][12][2];
65  TProfile *m_dt2_slope[5][12][2];
66  TH1F *m_posdt13_resid[5][12][3];
67  TH1F *m_posdt13_slope[5][12][3];
68  TH1F *m_posdt2_resid[5][12][2];
69  TH1F *m_posdt2_slope[5][12][2];
70  TH1F *m_negdt13_resid[5][12][3];
71  TH1F *m_negdt13_slope[5][12][3];
72  TH1F *m_negdt2_resid[5][12][2];
73  TH1F *m_negdt2_slope[5][12][2];
74 
75  // endcap, chamber, stationdiff
76  TProfile *m_cscouter_resid[2][36][2];
77  TProfile *m_cscouter_slope[2][36][2];
78  TProfile *m_cscinner_resid[2][18][3];
79  TProfile *m_cscinner_slope[2][18][3];
80  TH1F *m_poscscouter_resid[2][36][2];
81  TH1F *m_poscscouter_slope[2][36][2];
82  TH1F *m_poscscinner_resid[2][18][3];
83  TH1F *m_poscscinner_slope[2][18][3];
84  TH1F *m_negcscouter_resid[2][36][2];
85  TH1F *m_negcscouter_slope[2][36][2];
86  TH1F *m_negcscinner_resid[2][18][3];
87  TH1F *m_negcscinner_slope[2][18][3];
88 
89  // cross-system segdiffs: // [endcap][dtsector]
90  TH1F *m_x_pos_dt1_csc1_resid[2][12]; // tracks going through DT W+-2 St1 and CSC St1 Ring3
91  TH1F *m_x_pos_dt1_csc2_resid[2][12]; // tracks going through DT W+-2 St1 and CSC St2 Ring2
92  TH1F *m_x_pos_dt2_csc1_resid[2][12]; // tracks going through DT W+-2 St2 and CSC St1 Ring3
93  TH1F *m_x_neg_dt1_csc1_resid[2][12];
94  TH1F *m_x_neg_dt1_csc2_resid[2][12];
95  TH1F *m_x_neg_dt2_csc1_resid[2][12];
96 };
97 
98 //
99 // constants, enums and typedefs
100 //
101 
102 //
103 // static data member definitions
104 //
105 
106 //
107 // member functions
108 //
109 
110 AlignmentMonitorSegmentDifferences::AlignmentMonitorSegmentDifferences(const edm::ParameterSet &cfg)
111  : AlignmentMonitorBase(cfg, "AlignmentMonitorSegmentDifferences"),
112  m_muonCollectionTag(cfg.getParameter<edm::InputTag>("muonCollectionTag")),
113  m_minTrackPt(cfg.getParameter<double>("minTrackPt")),
114  m_minTrackP(cfg.getParameter<double>("minTrackP")),
115  m_maxDxy(cfg.getParameter<double>("maxDxy")),
116  m_minTrackerHits(cfg.getParameter<int>("minTrackerHits")),
117  m_maxTrackerRedChi2(cfg.getParameter<double>("maxTrackerRedChi2")),
118  m_allowTIDTEC(cfg.getParameter<bool>("allowTIDTEC")),
119  m_minNCrossedChambers(cfg.getParameter<int>("minNCrossedChambers")),
120  m_minDT13Hits(cfg.getParameter<int>("minDT13Hits")),
121  m_minDT2Hits(cfg.getParameter<int>("minDT2Hits")),
122  m_minCSCHits(cfg.getParameter<int>("minCSCHits")),
123  m_doDT(cfg.getParameter<bool>("doDT")),
124  m_doCSC(cfg.getParameter<bool>("doCSC")) {}
125 
126 void AlignmentMonitorSegmentDifferences::book() {
127  char name[225], pos[228], neg[228];
128 
129  double max_curv = 1. / m_minTrackPt;
130 
131  if (m_doDT)
132  for (int wheel = -2; wheel <= +2; wheel++) {
133  char wheel_label[][2] = {"A", "B", "C", "D", "E"};
134  for (int sector = 1; sector <= 12; sector++) {
135  char wheel_sector[50];
136  sprintf(wheel_sector, "%s_%02d", wheel_label[wheel + 2], sector);
137 
138  int nb = 100;
139  double wnd = 25.;
140 
141  sprintf(name, "dt13_resid_%s_12", wheel_sector);
142  sprintf(pos, "pos%s", name);
143  sprintf(neg, "neg%s", name);
144  m_dt13_resid[wheel + 2][sector - 1][0] =
145  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
146  m_posdt13_resid[wheel + 2][sector - 1][0] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
147  m_negdt13_resid[wheel + 2][sector - 1][0] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
148 
149  sprintf(name, "dt13_resid_%s_23", wheel_sector);
150  sprintf(pos, "pos%s", name);
151  sprintf(neg, "neg%s", name);
152  m_dt13_resid[wheel + 2][sector - 1][1] =
153  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
154  m_posdt13_resid[wheel + 2][sector - 1][1] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
155  m_negdt13_resid[wheel + 2][sector - 1][1] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
156 
157  sprintf(name, "dt13_resid_%s_34", wheel_sector);
158  sprintf(pos, "pos%s", name);
159  sprintf(neg, "neg%s", name);
160  m_dt13_resid[wheel + 2][sector - 1][2] =
161  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
162  m_posdt13_resid[wheel + 2][sector - 1][2] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
163  m_negdt13_resid[wheel + 2][sector - 1][2] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
164 
165  sprintf(name, "dt2_resid_%s_12", wheel_sector);
166  sprintf(pos, "pos%s", name);
167  sprintf(neg, "neg%s", name);
168  m_dt2_resid[wheel + 2][sector - 1][0] =
169  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -200., 200., " ");
170  m_posdt2_resid[wheel + 2][sector - 1][0] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
171  m_negdt2_resid[wheel + 2][sector - 1][0] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
172 
173  sprintf(name, "dt2_resid_%s_23", wheel_sector);
174  sprintf(pos, "pos%s", name);
175  sprintf(neg, "neg%s", name);
176  m_dt2_resid[wheel + 2][sector - 1][1] =
177  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -200., 200., " ");
178  m_posdt2_resid[wheel + 2][sector - 1][1] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
179  m_negdt2_resid[wheel + 2][sector - 1][1] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
180 
181  sprintf(name, "dt13_slope_%s_12", wheel_sector);
182  sprintf(pos, "pos%s", name);
183  sprintf(neg, "neg%s", name);
184  m_dt13_slope[wheel + 2][sector - 1][0] =
185  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
186  m_posdt13_slope[wheel + 2][sector - 1][0] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
187  m_negdt13_slope[wheel + 2][sector - 1][0] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
188 
189  sprintf(name, "dt13_slope_%s_23", wheel_sector);
190  sprintf(pos, "pos%s", name);
191  sprintf(neg, "neg%s", name);
192  m_dt13_slope[wheel + 2][sector - 1][1] =
193  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
194  m_posdt13_slope[wheel + 2][sector - 1][1] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
195  m_negdt13_slope[wheel + 2][sector - 1][1] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
196 
197  sprintf(name, "dt13_slope_%s_34", wheel_sector);
198  sprintf(pos, "pos%s", name);
199  sprintf(neg, "neg%s", name);
200  m_dt13_slope[wheel + 2][sector - 1][2] =
201  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
202  m_posdt13_slope[wheel + 2][sector - 1][2] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
203  m_negdt13_slope[wheel + 2][sector - 1][2] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
204 
205  sprintf(name, "dt2_slope_%s_12", wheel_sector);
206  sprintf(pos, "pos%s", name);
207  sprintf(neg, "neg%s", name);
208  m_dt2_slope[wheel + 2][sector - 1][0] =
209  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -1000., 1000., " ");
210  m_posdt2_slope[wheel + 2][sector - 1][0] = book1D("/iterN/", pos, pos, nb, -100., 100.);
211  m_negdt2_slope[wheel + 2][sector - 1][0] = book1D("/iterN/", neg, neg, nb, -100., 100.);
212 
213  sprintf(name, "dt2_slope_%s_23", wheel_sector);
214  sprintf(pos, "pos%s", name);
215  sprintf(neg, "neg%s", name);
216  m_dt2_slope[wheel + 2][sector - 1][1] =
217  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -1000., 1000., " ");
218  m_posdt2_slope[wheel + 2][sector - 1][1] = book1D("/iterN/", pos, pos, nb, -100., 100.);
219  m_negdt2_slope[wheel + 2][sector - 1][1] = book1D("/iterN/", neg, neg, nb, -100., 100.);
220  }
221  }
222 
223  if (m_doCSC)
224  for (int endcap = 1; endcap <= 2; endcap++) {
225  std::string endcapletter;
226  if (endcap == 1)
227  endcapletter = "p";
228  else if (endcap == 2)
229  endcapletter = "m";
230 
231  for (int chamber = 1; chamber <= 36; chamber++) {
232  char ec_chamber[50];
233  sprintf(ec_chamber, "%s_%02d", endcapletter.c_str(), chamber);
234 
235  int nb = 100;
236  double wnd = 60.;
237 
238  sprintf(name, "cscouter_resid_%s_12", ec_chamber);
239  sprintf(pos, "pos%s", name);
240  sprintf(neg, "neg%s", name);
241  m_cscouter_resid[endcap - 1][chamber - 1][0] =
242  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
243  m_poscscouter_resid[endcap - 1][chamber - 1][0] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
244  m_negcscouter_resid[endcap - 1][chamber - 1][0] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
245 
246  sprintf(name, "cscouter_resid_%s_23", ec_chamber);
247  sprintf(pos, "pos%s", name);
248  sprintf(neg, "neg%s", name);
249  m_cscouter_resid[endcap - 1][chamber - 1][1] =
250  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
251  m_poscscouter_resid[endcap - 1][chamber - 1][1] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
252  m_negcscouter_resid[endcap - 1][chamber - 1][1] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
253 
254  sprintf(name, "cscouter_slope_%s_12", ec_chamber);
255  sprintf(pos, "pos%s", name);
256  sprintf(neg, "neg%s", name);
257  m_cscouter_slope[endcap - 1][chamber - 1][0] =
258  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
259  m_poscscouter_slope[endcap - 1][chamber - 1][0] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
260  m_negcscouter_slope[endcap - 1][chamber - 1][0] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
261 
262  sprintf(name, "cscouter_slope_%s_23", ec_chamber);
263  sprintf(pos, "pos%s", name);
264  sprintf(neg, "neg%s", name);
265  m_cscouter_slope[endcap - 1][chamber - 1][1] =
266  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
267  m_poscscouter_slope[endcap - 1][chamber - 1][1] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
268  m_negcscouter_slope[endcap - 1][chamber - 1][1] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
269  }
270 
271  for (int chamber = 1; chamber <= 18; chamber++) {
272  char ec_chamber[50];
273  sprintf(ec_chamber, "%s_%02d", endcapletter.c_str(), chamber);
274 
275  int nb = 100;
276  double wnd = 40.;
277 
278  sprintf(name, "cscinner_resid_%s_12", ec_chamber);
279  sprintf(pos, "pos%s", name);
280  sprintf(neg, "neg%s", name);
281  m_cscinner_resid[endcap - 1][chamber - 1][0] =
282  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
283  m_poscscinner_resid[endcap - 1][chamber - 1][0] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
284  m_negcscinner_resid[endcap - 1][chamber - 1][0] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
285 
286  sprintf(name, "cscinner_resid_%s_23", ec_chamber);
287  sprintf(pos, "pos%s", name);
288  sprintf(neg, "neg%s", name);
289  m_cscinner_resid[endcap - 1][chamber - 1][1] =
290  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
291  m_poscscinner_resid[endcap - 1][chamber - 1][1] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
292  m_negcscinner_resid[endcap - 1][chamber - 1][1] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
293 
294  sprintf(name, "cscinner_resid_%s_34", ec_chamber);
295  sprintf(pos, "pos%s", name);
296  sprintf(neg, "neg%s", name);
297  m_cscinner_resid[endcap - 1][chamber - 1][2] =
298  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
299  m_poscscinner_resid[endcap - 1][chamber - 1][2] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
300  m_negcscinner_resid[endcap - 1][chamber - 1][2] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
301 
302  sprintf(name, "cscinner_slope_%s_12", ec_chamber);
303  sprintf(pos, "pos%s", name);
304  sprintf(neg, "neg%s", name);
305  m_cscinner_slope[endcap - 1][chamber - 1][0] =
306  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
307  m_poscscinner_slope[endcap - 1][chamber - 1][0] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
308  m_negcscinner_slope[endcap - 1][chamber - 1][0] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
309 
310  sprintf(name, "cscinner_slope_%s_23", ec_chamber);
311  sprintf(pos, "pos%s", name);
312  sprintf(neg, "neg%s", name);
313  m_cscinner_slope[endcap - 1][chamber - 1][1] =
314  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
315  m_poscscinner_slope[endcap - 1][chamber - 1][1] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
316  m_negcscinner_slope[endcap - 1][chamber - 1][1] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
317 
318  sprintf(name, "cscinner_slope_%s_34", ec_chamber);
319  sprintf(pos, "pos%s", name);
320  sprintf(neg, "neg%s", name);
321  m_cscinner_slope[endcap - 1][chamber - 1][2] =
322  bookProfile("/iterN/", name, name, 20, -max_curv, max_curv, 1, -100., 100., " ");
323  m_poscscinner_slope[endcap - 1][chamber - 1][2] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
324  m_negcscinner_slope[endcap - 1][chamber - 1][2] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
325  }
326  }
327 
328  // cross-system
329  for (int e = 1; e <= 2; e++)
330  for (int s = 1; s <= 12; s++) {
331  char endcap_sector[50];
332  if (e == 1)
333  sprintf(endcap_sector, "Wp2S%02d", s);
334  if (e == 2)
335  sprintf(endcap_sector, "Wm2S%02d", s);
336 
337  int nb = 200;
338  double wnd = 100.;
339 
340  sprintf(pos, "pos_x_dt1_csc1_%s", endcap_sector);
341  sprintf(neg, "neg_x_dt1_csc1_%s", endcap_sector);
342  m_x_pos_dt1_csc1_resid[e - 1][s - 1] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
343  m_x_neg_dt1_csc1_resid[e - 1][s - 1] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
344 
345  sprintf(pos, "pos_x_dt1_csc2_%s", endcap_sector);
346  sprintf(neg, "neg_x_dt1_csc2_%s", endcap_sector);
347  m_x_pos_dt1_csc2_resid[e - 1][s - 1] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
348  m_x_neg_dt1_csc2_resid[e - 1][s - 1] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
349 
350  sprintf(pos, "pos_x_dt2_csc1_%s", endcap_sector);
351  sprintf(neg, "neg_x_dt2_csc1_%s", endcap_sector);
352  m_x_pos_dt2_csc1_resid[e - 1][s - 1] = book1D("/iterN/", pos, pos, nb, -wnd, wnd);
353  m_x_neg_dt2_csc1_resid[e - 1][s - 1] = book1D("/iterN/", neg, neg, nb, -wnd, wnd);
354  }
355 }
356 
357 void AlignmentMonitorSegmentDifferences::event(const edm::Event &iEvent,
358  const edm::EventSetup &iSetup,
359  const ConstTrajTrackPairCollection &trajtracks) {
360  edm::ESHandle<GlobalTrackingGeometry> globalGeometry;
361  iSetup.get<GlobalTrackingGeometryRecord>().get(globalGeometry);
362 
363  edm::Handle<reco::BeamSpot> beamSpot;
364  iEvent.getByLabel(m_beamSpotTag, beamSpot);
365 
366  edm::ESHandle<DetIdAssociator> muonDetIdAssociator_;
367  iSetup.get<DetIdAssociatorRecord>().get("MuonDetIdAssociator", muonDetIdAssociator_);
368 
369  edm::ESHandle<Propagator> prop;
370  iSetup.get<TrackingComponentsRecord>().get("SteppingHelixPropagatorAny", prop);
371 
372  edm::ESHandle<MagneticField> magneticField;
373  iSetup.get<IdealMagneticFieldRecord>().get(magneticField);
374 
375  if (m_muonCollectionTag.label().empty()) // use trajectories
376  {
377  for (ConstTrajTrackPairCollection::const_iterator trajtrack = trajtracks.begin(); trajtrack != trajtracks.end();
378  ++trajtrack) {
379  const Trajectory *traj = (*trajtrack).first;
380  const reco::Track *track = (*trajtrack).second;
381 
382  if (track->pt() > m_minTrackPt && track->p() > m_minTrackP && fabs(track->dxy(beamSpot->position())) < m_maxDxy) {
383  MuonResidualsFromTrack muonResidualsFromTrack(
384  iSetup, magneticField, globalGeometry, muonDetIdAssociator_, prop, traj, track, pNavigator(), 1000.);
385  processMuonResidualsFromTrack(muonResidualsFromTrack);
386  }
387  } // end loop over tracks
388  } else {
389  edm::Handle<reco::MuonCollection> muons;
390  iEvent.getByLabel(m_muonCollectionTag, muons);
391 
392  for (reco::MuonCollection::const_iterator muon = muons->begin(); muon != muons->end(); ++muon) {
393  if (!(muon->isTrackerMuon() && muon->innerTrack().isNonnull()))
394  continue;
395 
396  if (m_minTrackPt < muon->pt() && m_minTrackP < muon->p() &&
397  fabs(muon->innerTrack()->dxy(beamSpot->position())) < m_maxDxy) {
398  MuonResidualsFromTrack muonResidualsFromTrack(globalGeometry, &(*muon), pNavigator(), 100.);
399  processMuonResidualsFromTrack(muonResidualsFromTrack);
400  }
401  }
402  }
403 }
404 
405 void AlignmentMonitorSegmentDifferences::processMuonResidualsFromTrack(MuonResidualsFromTrack &mrft) {
406  if (mrft.trackerNumHits() < m_minTrackerHits)
407  return;
408  if (!m_allowTIDTEC && mrft.contains_TIDTEC())
409  return;
410  if (mrft.normalizedChi2() > m_maxTrackerRedChi2)
411  return;
412 
413  int nMuChambers = 0;
414  std::vector<DetId> chamberIds = mrft.chamberIds();
415  for (unsigned ch = 0; ch < chamberIds.size(); ch++)
416  if (chamberIds[ch].det() == DetId::Muon)
417  nMuChambers++;
418  if (nMuChambers < m_minNCrossedChambers)
419  return;
420 
421  double qoverpt = (mrft.getTrack()->charge() > 0 ? 1. : -1.) / mrft.getTrack()->pt();
422  double qoverpz = 0.;
423  if (fabs(mrft.getTrack()->pz()) > 0.01)
424  qoverpz = mrft.getTrack()->charge() / fabs(mrft.getTrack()->pz());
425 
426  for (std::vector<DetId>::const_iterator chamberId = chamberIds.begin(); chamberId != chamberIds.end(); ++chamberId) {
427  if (chamberId->det() != DetId::Muon)
428  continue;
429 
430  // **************** DT ****************
431  if (m_doDT && chamberId->subdetId() == MuonSubdetId::DT) {
432  MuonChamberResidual *dt13 = mrft.chamberResidual(*chamberId, MuonChamberResidual::kDT13);
433  MuonChamberResidual *dt2 = mrft.chamberResidual(*chamberId, MuonChamberResidual::kDT2);
434 
435  if (dt13 != nullptr && dt13->numHits() >= m_minDT13Hits) {
436  DTChamberId thisid(chamberId->rawId());
437  for (std::vector<DetId>::const_iterator otherId = chamberIds.begin(); otherId != chamberIds.end(); ++otherId) {
438  if (otherId->det() == DetId::Muon && otherId->subdetId() == MuonSubdetId::DT) {
439  DTChamberId thatid(otherId->rawId());
440  if (thisid.rawId() != thatid.rawId() && thisid.wheel() == thatid.wheel() &&
441  thisid.sector() == thatid.sector()) {
442  MuonChamberResidual *dt13other = mrft.chamberResidual(*otherId, MuonChamberResidual::kDT13);
443  if (dt13other != nullptr && dt13other->numHits() >= m_minDT13Hits) {
444  double slopediff = 1000. * (dt13->global_resslope() - dt13other->global_resslope());
445  //double length = dt13->chamberAlignable()->surface().toGlobal(align::LocalPoint(0,0,0)).perp() -
446  // dt13other->chamberAlignable()->surface().toGlobal(align::LocalPoint(0,0,0)).perp();
447  //double residdiff = 10. * (dt13->global_residual() + length*dt13->global_resslope() - dt13other->global_residual());
448  double residdiff = 10. * (dt13->global_residual() - dt13other->global_residual());
449 
450  int st = 0;
451  if (thatid.station() - thisid.station() == 1)
452  st = thisid.station();
453  if (st > 0) {
454  m_dt13_resid[thisid.wheel() + 2][thisid.sector() - 1][st - 1]->Fill(qoverpt, residdiff);
455  m_dt13_slope[thisid.wheel() + 2][thisid.sector() - 1][st - 1]->Fill(qoverpt, slopediff);
456  if (qoverpt > 0) {
457  m_posdt13_resid[thisid.wheel() + 2][thisid.sector() - 1][st - 1]->Fill(residdiff);
458  m_posdt13_slope[thisid.wheel() + 2][thisid.sector() - 1][st - 1]->Fill(slopediff);
459  } else {
460  m_negdt13_resid[thisid.wheel() + 2][thisid.sector() - 1][st - 1]->Fill(residdiff);
461  m_negdt13_slope[thisid.wheel() + 2][thisid.sector() - 1][st - 1]->Fill(slopediff);
462  }
463  }
464  } // end other numhits
465  } // end this near other
466  } // end other is DT
467 
468  // cross-system: other is CSC
469  // only do it for DT stubs in W+-2 St1&2:
470  if (!(abs(thisid.wheel()) == 2 && (thisid.station() == 1 || thisid.station() == 2)))
471  continue;
472  if (otherId->det() == DetId::Muon && otherId->subdetId() == MuonSubdetId::CSC) {
473  CSCDetId thatid(otherId->rawId());
474  //only do it for CSC stubs in St1R3 or St2R2:
475  if (!((thatid.station() == 1 && thatid.ring() == 3) || (thatid.station() == 2 && thatid.ring() == 2)))
476  continue;
477 
478  MuonChamberResidual *cscother = mrft.chamberResidual(*otherId, MuonChamberResidual::kCSC);
479  if (cscother != nullptr && cscother->numHits() >= m_minCSCHits) {
480  // scale to adjust the csc residual size to be comparabe to dt's one
481  double csc_scale = dt13->chamberAlignable()
482  ->surface()
483  .toGlobal(align::LocalPoint(dt13->trackx(), dt13->tracky(), 0))
484  .perp() /
485  cscother->chamberAlignable()
486  ->surface()
487  .toGlobal(align::LocalPoint(cscother->trackx(), cscother->tracky(), 0))
488  .perp();
489  double residdiff = 10. * (dt13->global_residual() - cscother->global_residual() * csc_scale);
490  if (thisid.station() == 1 && thatid.station() == 1) {
491  if (qoverpt > 0)
492  m_x_pos_dt1_csc1_resid[thatid.endcap() - 1][thisid.sector() - 1]->Fill(residdiff);
493  else
494  m_x_neg_dt1_csc1_resid[thatid.endcap() - 1][thisid.sector() - 1]->Fill(residdiff);
495  } else if (thisid.station() == 1 && thatid.station() == 2) {
496  if (qoverpt > 0)
497  m_x_pos_dt1_csc2_resid[thatid.endcap() - 1][thisid.sector() - 1]->Fill(residdiff);
498  else
499  m_x_neg_dt1_csc2_resid[thatid.endcap() - 1][thisid.sector() - 1]->Fill(residdiff);
500  } else if (thisid.station() == 2 && thatid.station() == 1) {
501  if (qoverpt > 0)
502  m_x_pos_dt2_csc1_resid[thatid.endcap() - 1][thisid.sector() - 1]->Fill(residdiff);
503  else
504  m_x_neg_dt2_csc1_resid[thatid.endcap() - 1][thisid.sector() - 1]->Fill(residdiff);
505  }
506  }
507  } // end other is CSC
508  } // end loop over other
509  } // end if DT13
510 
511  // z-direction
512  if (dt2 != nullptr && dt2->numHits() >= m_minDT2Hits && (dt2->chi2() / double(dt2->ndof())) < 2.0) {
513  DTChamberId thisid(chamberId->rawId());
514  for (std::vector<DetId>::const_iterator otherId = chamberIds.begin(); otherId != chamberIds.end(); ++otherId) {
515  if (otherId->det() == DetId::Muon && otherId->subdetId() == MuonSubdetId::DT) {
516  DTChamberId thatid(otherId->rawId());
517  if (thisid.rawId() != thatid.rawId() && thisid.wheel() == thatid.wheel() &&
518  thisid.sector() == thatid.sector()) {
519  MuonChamberResidual *dt2other = mrft.chamberResidual(*otherId, MuonChamberResidual::kDT2);
520  if (dt2other != nullptr && dt2other->numHits() >= m_minDT2Hits) {
521  double slopediff = 1000. * (dt2->global_resslope() - dt2other->global_resslope());
522  //double length = dt2->chamberAlignable()->surface().toGlobal(align::LocalPoint(0,0,0)).perp() -
523  // dt2other->chamberAlignable()->surface().toGlobal(align::LocalPoint(0,0,0)).perp();
524  //double residdiff = 10. * (dt2->global_residual() + length*dt2->global_resslope() - dt2other->global_residual());
525  double residdiff = 10. * (dt2->global_residual() - dt2other->global_residual());
526 
527  int st = 0;
528  if (thatid.station() - thisid.station() == 1)
529  st = thisid.station();
530  if (st > 0) {
531  m_dt2_resid[thisid.wheel() + 2][thisid.sector() - 1][st - 1]->Fill(qoverpt, residdiff);
532  m_dt2_slope[thisid.wheel() + 2][thisid.sector() - 1][st - 1]->Fill(qoverpt, slopediff);
533  if (qoverpt > 0) {
534  m_posdt2_resid[thisid.wheel() + 2][thisid.sector() - 1][st - 1]->Fill(residdiff);
535  m_posdt2_slope[thisid.wheel() + 2][thisid.sector() - 1][st - 1]->Fill(slopediff);
536  } else {
537  m_negdt2_resid[thisid.wheel() + 2][thisid.sector() - 1][st - 1]->Fill(residdiff);
538  m_negdt2_slope[thisid.wheel() + 2][thisid.sector() - 1][st - 1]->Fill(slopediff);
539  }
540  }
541  } // end other numhits
542  } // end this near other
543  } // end other is DT
544  } // end loop over other
545  } // end if DT2
546  } // end if DT
547 
548  // **************** CSC ****************
549  else if (m_doCSC && chamberId->subdetId() == MuonSubdetId::CSC) {
550  MuonChamberResidual *csc = mrft.chamberResidual(*chamberId, MuonChamberResidual::kCSC);
551  if (csc->numHits() >= m_minCSCHits) {
552  CSCDetId thisid(chamberId->rawId());
553  for (std::vector<DetId>::const_iterator otherId = chamberIds.begin(); otherId != chamberIds.end(); ++otherId) {
554  if (otherId->det() == DetId::Muon && otherId->subdetId() == MuonSubdetId::CSC) {
555  CSCDetId thatid(otherId->rawId());
556  if (thisid.rawId() != thatid.rawId() && thisid.endcap() == thatid.endcap()) {
557  MuonChamberResidual *cscother = mrft.chamberResidual(*otherId, MuonChamberResidual::kCSC);
558  if (cscother != nullptr && cscother->numHits() >= m_minCSCHits) {
559  double slopediff = 1000. * (csc->global_resslope() - cscother->global_resslope());
560  //double length = csc->chamberAlignable()->surface().toGlobal(align::LocalPoint(0,0,0)).z() -
561  // cscother->chamberAlignable()->surface().toGlobal(align::LocalPoint(0,0,0)).z();
562  //double residdiff = 10. * (csc->global_residual() + length*csc->global_resslope() - cscother->global_residual());
563  double residdiff = 10. * (csc->global_residual() - cscother->global_residual());
564 
565  int thischamber = thisid.chamber();
566  int thisring = thisid.ring();
567  if (thisid.station() == 1 && (thisring == 1 || thisring == 4)) {
568  thischamber = (thischamber - 1) / 2 + 1;
569  thisring = 1;
570  }
571 
572  if (thisring == thatid.ring() && thischamber == thatid.chamber()) {
573  bool inner = (thisring == 1);
574  bool outer = (thisring == 2);
575  int st = 0;
576  if (thatid.station() - thisid.station() == 1 && (inner || thisid.station() < 3))
577  st = thisid.station();
578 
579  if (outer && st > 0) {
580  m_cscouter_resid[thisid.endcap() - 1][thischamber - 1][st - 1]->Fill(qoverpz, residdiff);
581  m_cscouter_slope[thisid.endcap() - 1][thischamber - 1][st - 1]->Fill(qoverpz, slopediff);
582  if (qoverpz > 0) {
583  m_poscscouter_resid[thisid.endcap() - 1][thischamber - 1][st - 1]->Fill(residdiff);
584  m_poscscouter_slope[thisid.endcap() - 1][thischamber - 1][st - 1]->Fill(slopediff);
585  } else {
586  m_negcscouter_resid[thisid.endcap() - 1][thischamber - 1][st - 1]->Fill(residdiff);
587  m_negcscouter_slope[thisid.endcap() - 1][thischamber - 1][st - 1]->Fill(slopediff);
588  }
589  }
590  if (inner && st > 0) {
591  m_cscinner_resid[thisid.endcap() - 1][thischamber - 1][st - 1]->Fill(qoverpz, residdiff);
592  m_cscinner_slope[thisid.endcap() - 1][thischamber - 1][st - 1]->Fill(qoverpz, slopediff);
593  if (qoverpz > 0) {
594  m_poscscinner_resid[thisid.endcap() - 1][thischamber - 1][st - 1]->Fill(residdiff);
595  m_poscscinner_slope[thisid.endcap() - 1][thischamber - 1][st - 1]->Fill(slopediff);
596  } else {
597  m_negcscinner_resid[thisid.endcap() - 1][thischamber - 1][st - 1]->Fill(residdiff);
598  m_negcscinner_slope[thisid.endcap() - 1][thischamber - 1][st - 1]->Fill(slopediff);
599  }
600  }
601  } // end of same ring&chamber
602  } // end other min numhits
603  } // end this near other
604  } // end other is CSC
605  } // end loop over other
606 
607  } // end if this min numhits
608  } // end if CSC
609 
610  } // end loop over chamberIds
611 }
612 
613 DEFINE_EDM_PLUGIN(AlignmentMonitorPluginFactory,
614  AlignmentMonitorSegmentDifferences,
615  "AlignmentMonitorSegmentDifferences");
AlignmentMonitorSegmentDifferences::m_minTrackP
double m_minTrackP
Definition: AlignmentMonitorSegmentDifferences.cc:49
Propagator.h
PDWG_BPHSkim_cff.muons
muons
Definition: PDWG_BPHSkim_cff.py:47
MuonSubdetId::CSC
static constexpr int CSC
Definition: MuonSubdetId.h:12
AlignmentMonitorSegmentDifferences::m_x_pos_dt1_csc2_resid
TH1F * m_x_pos_dt1_csc2_resid[2][12]
Definition: AlignmentMonitorSegmentDifferences.cc:91
AlignmentMonitorSegmentDifferences::m_minTrackPt
double m_minTrackPt
Definition: AlignmentMonitorSegmentDifferences.cc:48
electrons_cff.bool
bool
Definition: electrons_cff.py:366
pwdgSkimBPark_cfi.beamSpot
beamSpot
Definition: pwdgSkimBPark_cfi.py:5
HLT_FULL_cff.track
track
Definition: HLT_FULL_cff.py:11713
MessageLogger.h
AlignmentMonitorSegmentDifferences::m_posdt2_resid
TH1F * m_posdt2_resid[5][12][2]
Definition: AlignmentMonitorSegmentDifferences.cc:68
MuonResidualsAngleFitter.h
muon
Definition: MuonCocktails.h:17
AlignmentMonitorSegmentDifferences::m_negdt13_slope
TH1F * m_negdt13_slope[5][12][3]
Definition: AlignmentMonitorSegmentDifferences.cc:71
DiDispStaMuonMonitor_cfi.pt
pt
Definition: DiDispStaMuonMonitor_cfi.py:39
MuonResidualsPositionFitter.h
edm
HLT enums.
Definition: AlignableModifier.h:19
AlignmentMonitorSegmentDifferences::m_poscscouter_slope
TH1F * m_poscscouter_slope[2][36][2]
Definition: AlignmentMonitorSegmentDifferences.cc:81
pos
Definition: PixelAliasList.h:18
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89285
MuonResidualsFromTrack::chamberResidual
MuonChamberResidual * chamberResidual(DetId chamberId, int type)
Definition: MuonResidualsFromTrack.cc:735
MuonResidualsFromTrack::trackerNumHits
int trackerNumHits() const
Definition: MuonResidualsFromTrack.h:78
MuonResidualsTwoBin.h
GlobalTrackingGeometryRecord
Definition: GlobalTrackingGeometryRecord.h:17
AlignmentMonitorSegmentDifferences::m_x_neg_dt2_csc1_resid
TH1F * m_x_neg_dt2_csc1_resid[2][12]
Definition: AlignmentMonitorSegmentDifferences.cc:95
AlignmentMonitorPluginFactory
AlignmentMonitorSegmentDifferences::m_negdt2_resid
TH1F * m_negdt2_resid[5][12][2]
Definition: AlignmentMonitorSegmentDifferences.cc:72
AlignmentMonitorSegmentDifferences::m_negcscouter_slope
TH1F * m_negcscouter_slope[2][36][2]
Definition: AlignmentMonitorSegmentDifferences.cc:85
makeMuonMisalignmentScenario.endcap
endcap
Definition: makeMuonMisalignmentScenario.py:320
HLT_FULL_cff.magneticField
magneticField
Definition: HLT_FULL_cff.py:348
edm::Handle< reco::BeamSpot >
AlignmentMonitorSegmentDifferences::m_maxDxy
double m_maxDxy
Definition: AlignmentMonitorSegmentDifferences.cc:50
AlignmentMonitorBase::pNavigator
AlignableNavigator * pNavigator()
Definition: AlignmentMonitorBase.h:114
AlignmentMonitorSegmentDifferences::m_cscouter_slope
TProfile * m_cscouter_slope[2][36][2]
Definition: AlignmentMonitorSegmentDifferences.cc:77
AlignmentMonitorSegmentDifferences::m_poscscinner_slope
TH1F * m_poscscinner_slope[2][18][3]
Definition: AlignmentMonitorSegmentDifferences.cc:83
MuonChamberResidual::global_residual
double global_residual() const
Definition: MuonChamberResidual.cc:49
AlignmentMonitorSegmentDifferences::book
void book() override
Book or retrieve histograms; MUST be reimplemented.
Definition: AlignmentMonitorSegmentDifferences.cc:126
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:637
IdealMagneticFieldRecord
Definition: IdealMagneticFieldRecord.h:11
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
AlignmentMonitorSegmentDifferences::m_minTrackerHits
int m_minTrackerHits
Definition: AlignmentMonitorSegmentDifferences.cc:51
csc
Definition: L1Track.h:19
AlignmentMonitorSegmentDifferences::m_posdt2_slope
TH1F * m_posdt2_slope[5][12][2]
Definition: AlignmentMonitorSegmentDifferences.cc:69
alignCSCRings.s
s
Definition: alignCSCRings.py:92
AlignmentMonitorSegmentDifferences::m_muonCollectionTag
edm::InputTag m_muonCollectionTag
Definition: AlignmentMonitorSegmentDifferences.cc:47
AlignmentMonitorSegmentDifferences::m_negcscinner_slope
TH1F * m_negcscinner_slope[2][18][3]
Definition: AlignmentMonitorSegmentDifferences.cc:87
edm::EventSetup::get
T get() const
Definition: EventSetup.h:87
AlignmentMonitorSegmentDifferences::m_poscscinner_resid
TH1F * m_poscscinner_resid[2][18][3]
Definition: AlignmentMonitorSegmentDifferences.cc:82
AlignmentMonitorSegmentDifferences::m_negcscinner_resid
TH1F * m_negcscinner_resid[2][18][3]
Definition: AlignmentMonitorSegmentDifferences.cc:86
AlignmentMonitorSegmentDifferences::m_minDT2Hits
int m_minDT2Hits
Definition: AlignmentMonitorSegmentDifferences.cc:56
AlignmentMonitorSegmentDifferences::m_cscinner_resid
TProfile * m_cscinner_resid[2][18][3]
Definition: AlignmentMonitorSegmentDifferences.cc:78
HLT_FULL_cff.muon
muon
Definition: HLT_FULL_cff.py:11710
SurfaceOrientation::inner
Definition: Surface.h:19
AlignmentMonitorSegmentDifferences::m_maxTrackerRedChi2
double m_maxTrackerRedChi2
Definition: AlignmentMonitorSegmentDifferences.cc:52
reco::Track
Definition: Track.h:27
IdealMagneticFieldRecord.h
MuonResidualsFromTrack
Definition: MuonResidualsFromTrack.h:52
edm::ESHandle< GlobalTrackingGeometry >
MuonResidualsFromTrack::getTrack
const reco::Track * getTrack()
Definition: MuonResidualsFromTrack.h:75
MuonChamberResidual::chi2
double chi2() const
Definition: MuonChamberResidual.h:59
MuonResidualsFromTrack::chamberIds
const std::vector< DetId > chamberIds() const
Definition: MuonResidualsFromTrack.h:86
reco::TrackBase::charge
int charge() const
track electric charge
Definition: TrackBase.h:596
AlignableSurface::toGlobal
align::GlobalPoints toGlobal(const align::LocalPoints &) const
Return in global coord given a set of local points.
Definition: AlignableSurface.cc:15
MuonSubdetId::DT
static constexpr int DT
Definition: MuonSubdetId.h:11
Alignable::surface
const AlignableSurface & surface() const
Return the Surface (global position and orientation) of the object.
Definition: Alignable.h:132
Point3DBase
Definition: Point3DBase.h:10
MuonResidualsFromTrack::normalizedChi2
double normalizedChi2() const
Definition: MuonResidualsFromTrack.cc:729
DEFINE_EDM_PLUGIN
#define DEFINE_EDM_PLUGIN(factory, type, name)
Definition: PluginFactory.h:124
AlignmentMonitorPluginFactory.h
GlobalTrackingGeometryRecord.h
AlignmentMonitorSegmentDifferences::m_dt2_resid
TProfile * m_dt2_resid[5][12][2]
Definition: AlignmentMonitorSegmentDifferences.cc:64
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
AlignmentMonitorSegmentDifferences::m_x_pos_dt1_csc1_resid
TH1F * m_x_pos_dt1_csc1_resid[2][12]
Definition: AlignmentMonitorSegmentDifferences.cc:90
AlignmentMonitorSegmentDifferences::m_posdt13_slope
TH1F * m_posdt13_slope[5][12][3]
Definition: AlignmentMonitorSegmentDifferences.cc:67
AlignmentMonitorSegmentDifferences::m_minCSCHits
int m_minCSCHits
Definition: AlignmentMonitorSegmentDifferences.cc:57
AlignmentMonitorSegmentDifferences::m_dt2_slope
TProfile * m_dt2_slope[5][12][2]
Definition: AlignmentMonitorSegmentDifferences.cc:65
MuonChamberResidual::kDT2
Definition: MuonChamberResidual.h:28
AlignmentMonitorSegmentDifferences::m_minDT13Hits
int m_minDT13Hits
Definition: AlignmentMonitorSegmentDifferences.cc:55
edm::ParameterSet
Definition: ParameterSet.h:47
AlignmentMonitorBase::bookProfile
TProfile * bookProfile(std::string dir, std::string name, std::string title, int nchX, double lowX, double highX, int nchY=1, double lowY=0., double highY=0., const char *option="s")
Definition: AlignmentMonitorBase.cc:112
MuonChamberResidual::numHits
int numHits() const
Definition: MuonChamberResidual.h:52
MuonChamberResidual::chamberAlignable
AlignableDetOrUnitPtr chamberAlignable() const
Definition: MuonChamberResidual.h:50
AlignmentMonitorBase.h
AlignmentMonitorSegmentDifferences::m_dt13_resid
TProfile * m_dt13_resid[5][12][3]
Definition: AlignmentMonitorSegmentDifferences.cc:62
MuonChamberResidual::trackx
double trackx() const
Definition: MuonChamberResidual.h:64
CSCDetId
Definition: CSCDetId.h:26
AlignmentMonitorSegmentDifferences::m_doCSC
bool m_doCSC
Definition: AlignmentMonitorSegmentDifferences.cc:59
AlignmentMonitorBase
Definition: AlignmentMonitorBase.h:42
makeMuonMisalignmentScenario.wheel
wheel
Definition: makeMuonMisalignmentScenario.py:319
createfilelist.int
int
Definition: createfilelist.py:10
MuonChamberResidual
Definition: MuonChamberResidual.h:26
iEvent
int iEvent
Definition: GenABIO.cc:224
MuonResidualsFromTrack::contains_TIDTEC
bool contains_TIDTEC() const
Definition: MuonResidualsFromTrack.h:84
AlignmentMonitorBase::m_beamSpotTag
const edm::InputTag m_beamSpotTag
Definition: AlignmentMonitorBase.h:116
DetIdAssociatorRecord.h
MagneticField.h
edm::EventSetup
Definition: EventSetup.h:58
MuonChamberResidual::kCSC
Definition: MuonChamberResidual.h:28
MuonResidualsFromTrack.h
get
#define get
AlignmentMonitorSegmentDifferences::m_x_neg_dt1_csc2_resid
TH1F * m_x_neg_dt1_csc2_resid[2][12]
Definition: AlignmentMonitorSegmentDifferences.cc:94
DetIdAssociatorRecord
Definition: DetIdAssociatorRecord.h:13
InputTag.h
AlignmentMonitorSegmentDifferences::~AlignmentMonitorSegmentDifferences
~AlignmentMonitorSegmentDifferences() override
Definition: AlignmentMonitorSegmentDifferences.cc:37
looper.cfg
cfg
Definition: looper.py:297
AlignmentMonitorSegmentDifferences::m_x_neg_dt1_csc1_resid
TH1F * m_x_neg_dt1_csc1_resid[2][12]
Definition: AlignmentMonitorSegmentDifferences.cc:93
AlignmentMonitorBase::book1D
TH1F * book1D(std::string dir, std::string name, std::string title, int nchX, double lowX, double highX)
Definition: AlignmentMonitorBase.cc:107
AlignmentMonitorSegmentDifferences::m_negcscouter_resid
TH1F * m_negcscouter_resid[2][36][2]
Definition: AlignmentMonitorSegmentDifferences.cc:84
MuonChamberResidual::tracky
double tracky() const
Definition: MuonChamberResidual.h:65
AlignmentMonitorSegmentDifferences::m_posdt13_resid
TH1F * m_posdt13_resid[5][12][3]
Definition: AlignmentMonitorSegmentDifferences.cc:66
Trajectory
Definition: Trajectory.h:38
DetIdAssociator.h
TrackingComponentsRecord.h
AlignmentMonitorSegmentDifferences::m_cscouter_resid
TProfile * m_cscouter_resid[2][36][2]
Definition: AlignmentMonitorSegmentDifferences.cc:76
AlignmentMonitorSegmentDifferences::AlignmentMonitorSegmentDifferences
AlignmentMonitorSegmentDifferences(const edm::ParameterSet &cfg)
Definition: AlignmentMonitorSegmentDifferences.cc:110
AlignmentMonitorSegmentDifferences
Definition: AlignmentMonitorSegmentDifferences.cc:34
AlignmentMonitorSegmentDifferences::m_cscinner_slope
TProfile * m_cscinner_slope[2][18][3]
Definition: AlignmentMonitorSegmentDifferences.cc:79
relativeConstraints.chamber
chamber
Definition: relativeConstraints.py:53
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
EventSetup.h
AlignmentMonitorSegmentDifferences::m_negdt2_slope
TH1F * m_negdt2_slope[5][12][2]
Definition: AlignmentMonitorSegmentDifferences.cc:73
AlignmentMonitorSegmentDifferences::event
void event(const edm::Event &iEvent, const edm::EventSetup &iSetup, const ConstTrajTrackPairCollection &iTrajTracks) override
Called for each event (by "run()"): may be reimplemented.
Definition: AlignmentMonitorSegmentDifferences.cc:357
AlignmentMonitorSegmentDifferences::m_doDT
bool m_doDT
Definition: AlignmentMonitorSegmentDifferences.cc:58
AlignmentMonitorSegmentDifferences::m_dt13_slope
TProfile * m_dt13_slope[5][12][3]
Definition: AlignmentMonitorSegmentDifferences.cc:63
reco::TrackBase::pz
double pz() const
z coordinate of momentum vector
Definition: TrackBase.h:646
AlignmentMonitorBase::ConstTrajTrackPairCollection
std::vector< ConstTrajTrackPair > ConstTrajTrackPairCollection
Definition: AlignmentMonitorBase.h:45
DetId::Muon
Definition: DetId.h:26
MuonChamberResidual::ndof
int ndof() const
Definition: MuonChamberResidual.h:60
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
AlignmentMonitorSegmentDifferences::m_x_pos_dt2_csc1_resid
TH1F * m_x_pos_dt2_csc1_resid[2][12]
Definition: AlignmentMonitorSegmentDifferences.cc:92
DTChamberId
Definition: DTChamberId.h:14
ParameterSet.h
AlignmentMonitorSegmentDifferences::m_allowTIDTEC
bool m_allowTIDTEC
Definition: AlignmentMonitorSegmentDifferences.cc:53
SurfaceOrientation::outer
Definition: Surface.h:19
AlignmentMonitorSegmentDifferences::processMuonResidualsFromTrack
void processMuonResidualsFromTrack(MuonResidualsFromTrack &mrft)
Definition: AlignmentMonitorSegmentDifferences.cc:405
GlobalTrackingGeometry.h
MuonChamberResidual::kDT13
Definition: MuonChamberResidual.h:28
edm::Event
Definition: Event.h:73
AlignmentMonitorSegmentDifferences::m_poscscouter_resid
TH1F * m_poscscouter_resid[2][36][2]
Definition: AlignmentMonitorSegmentDifferences.cc:80
MuonChamberResidual::global_resslope
double global_resslope() const
Definition: MuonChamberResidual.cc:51
AlignmentMonitorSegmentDifferences::m_negdt13_resid
TH1F * m_negdt13_resid[5][12][3]
Definition: AlignmentMonitorSegmentDifferences.cc:70
edm::InputTag
Definition: InputTag.h:15
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
AlignmentMonitorSegmentDifferences::m_minNCrossedChambers
bool m_minNCrossedChambers
Definition: AlignmentMonitorSegmentDifferences.cc:54
TrackingComponentsRecord
Definition: TrackingComponentsRecord.h:12
Generated for CMSSW Reference Manual by doxygen 1.8.16