CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
MuonIdDQM.cc
Go to the documentation of this file.
2 
4 {
5  inputMuonCollection_ = iConfig.getParameter<edm::InputTag>("inputMuonCollection");
6  inputDTRecSegment4DCollection_ = iConfig.getParameter<edm::InputTag>("inputDTRecSegment4DCollection");
7  inputCSCSegmentCollection_ = iConfig.getParameter<edm::InputTag>("inputCSCSegmentCollection");
8  useTrackerMuons_ = iConfig.getUntrackedParameter<bool>("useTrackerMuons");
9  useGlobalMuons_ = iConfig.getUntrackedParameter<bool>("useGlobalMuons");
10  useTrackerMuonsNotGlobalMuons_ = iConfig.getUntrackedParameter<bool>("useTrackerMuonsNotGlobalMuons");
11  useGlobalMuonsNotTrackerMuons_ = iConfig.getUntrackedParameter<bool>("useGlobalMuonsNotTrackerMuons");
12  baseFolder_ = iConfig.getUntrackedParameter<std::string>("baseFolder");
13 
14  dbe_ = 0;
15  dbe_ = edm::Service<DQMStore>().operator->();
16 }
17 
19 
20 void
22 {
23  char name[100], title[200];
24 
25  // trackerMuon == 0; globalMuon == 1; trackerMuon && !globalMuon == 2; globalMuon && !trackerMuon == 3
26  for (unsigned int i = 0; i < 4; i++) {
27  if ((i == 0 && ! useTrackerMuons_) || (i == 1 && ! useGlobalMuons_)) continue;
28  if ((i == 2 && ! useTrackerMuonsNotGlobalMuons_) || (i == 3 && ! useGlobalMuonsNotTrackerMuons_)) continue;
29  if (i == 0) dbe_->setCurrentFolder(baseFolder_+"/TrackerMuons");
30  if (i == 1) dbe_->setCurrentFolder(baseFolder_+"/GlobalMuons");
31  if (i == 2) dbe_->setCurrentFolder(baseFolder_+"/TrackerMuonsNotGlobalMuons");
32  if (i == 3) dbe_->setCurrentFolder(baseFolder_+"/GlobalMuonsNotTrackerMuons");
33 
34  hNumChambers[i] = dbe_->book1D("hNumChambers", "Number of Chambers", 17, -0.5, 16.5);
35  hNumMatches[i] = dbe_->book1D("hNumMatches", "Number of Matches", 11, -0.5, 10.5);
36  hNumChambersNoRPC[i] = dbe_->book1D("hNumChambersNoRPC", "Number of Chambers No RPC", 11, -0.5, 10.5);
37 
38  // by station
39  for(int station = 0; station < 4; ++station)
40  {
41  sprintf(name, "hDT%iNumSegments", station+1);
42  sprintf(title, "DT Station %i Number of Segments (No Arbitration)", station+1);
43  hDTNumSegments[i][station] = dbe_->book1D(name, title, 11, -0.5, 10.5);
44 
45  sprintf(name, "hDT%iDx", station+1);
46  sprintf(title, "DT Station %i Delta X", station+1);
47  hDTDx[i][station] = dbe_->book1D(name, title, 100, -100., 100.);
48 
49  sprintf(name, "hDT%iPullx", station+1);
50  sprintf(title, "DT Station %i Pull X", station+1);
51  hDTPullx[i][station] = dbe_->book1D(name, title, 100, -20., 20.);
52 
53  sprintf(name, "hDT%iDdXdZ", station+1);
54  sprintf(title, "DT Station %i Delta DxDz", station+1);
55  hDTDdXdZ[i][station] = dbe_->book1D(name, title, 100, -1., 1.);
56 
57  sprintf(name, "hDT%iPulldXdZ", station+1);
58  sprintf(title, "DT Station %i Pull DxDz", station+1);
59  hDTPulldXdZ[i][station] = dbe_->book1D(name, title, 100, -20., 20.);
60 
61  if (station < 3) {
62  sprintf(name, "hDT%iDy", station+1);
63  sprintf(title, "DT Station %i Delta Y", station+1);
64  hDTDy[i][station] = dbe_->book1D(name, title, 100, -150., 150.);
65 
66  sprintf(name, "hDT%iPully", station+1);
67  sprintf(title, "DT Station %i Pull Y", station+1);
68  hDTPully[i][station] = dbe_->book1D(name, title, 100, -20., 20.);
69 
70  sprintf(name, "hDT%iDdYdZ", station+1);
71  sprintf(title, "DT Station %i Delta DyDz", station+1);
72  hDTDdYdZ[i][station] = dbe_->book1D(name, title, 100, -2., 2.);
73 
74  sprintf(name, "hDT%iPulldYdZ", station+1);
75  sprintf(title, "DT Station %i Pull DyDz", station+1);
76  hDTPulldYdZ[i][station] = dbe_->book1D(name, title, 100, -20., 20.);
77  }
78 
79  sprintf(name, "hCSC%iNumSegments", station+1);
80  sprintf(title, "CSC Station %i Number of Segments (No Arbitration)", station+1);
81  hCSCNumSegments[i][station] = dbe_->book1D(name, title, 11, -0.5, 10.5);
82 
83  sprintf(name, "hCSC%iDx", station+1);
84  sprintf(title, "CSC Station %i Delta X", station+1);
85  hCSCDx[i][station] = dbe_->book1D(name, title, 100, -50., 50.);
86 
87  sprintf(name, "hCSC%iPullx", station+1);
88  sprintf(title, "CSC Station %i Pull X", station+1);
89  hCSCPullx[i][station] = dbe_->book1D(name, title, 100, -20., 20.);
90 
91  sprintf(name, "hCSC%iDdXdZ", station+1);
92  sprintf(title, "CSC Station %i Delta DxDz", station+1);
93  hCSCDdXdZ[i][station] = dbe_->book1D(name, title, 100, -1., 1.);
94 
95  sprintf(name, "hCSC%iPulldXdZ", station+1);
96  sprintf(title, "CSC Station %i Pull DxDz", station+1);
97  hCSCPulldXdZ[i][station] = dbe_->book1D(name, title, 100, -20., 20.);
98 
99  sprintf(name, "hCSC%iDy", station+1);
100  sprintf(title, "CSC Station %i Delta Y", station+1);
101  hCSCDy[i][station] = dbe_->book1D(name, title, 100, -50., 50.);
102 
103  sprintf(name, "hCSC%iPully", station+1);
104  sprintf(title, "CSC Station %i Pull Y", station+1);
105  hCSCPully[i][station] = dbe_->book1D(name, title, 100, -20., 20.);
106 
107  sprintf(name, "hCSC%iDdYdZ", station+1);
108  sprintf(title, "CSC Station %i Delta DyDz", station+1);
109  hCSCDdYdZ[i][station] = dbe_->book1D(name, title, 100, -1., 1.);
110 
111  sprintf(name, "hCSC%iPulldYdZ", station+1);
112  sprintf(title, "CSC Station %i Pull DyDz", station+1);
113  hCSCPulldYdZ[i][station] = dbe_->book1D(name, title, 100, -20., 20.);
114  }// station
115  }
116 
118  hSegmentIsAssociatedBool = dbe_->book1D("hSegmentIsAssociatedBool", "Segment Is Associated Boolean", 2, -0.5, 1.5);
119 }
120 
121 void
123 {
124  using namespace edm;
125  using namespace reco;
126 
131 
132  for(MuonCollection::const_iterator muon = muonCollectionH_->begin();
133  muon != muonCollectionH_->end(); ++muon)
134  {
135  // trackerMuon == 0; globalMuon == 1; trackerMuon && !globalMuon == 2; globalMuon && !trackerMuon == 3
136  for (unsigned int i = 0; i < 4; i++) {
137  if (i == 0 && (! useTrackerMuons_ || ! muon->isTrackerMuon())) continue;
138  if (i == 1 && (! useGlobalMuons_ || ! muon->isGlobalMuon())) continue;
139  if (i == 2 && (! useTrackerMuonsNotGlobalMuons_ || (! (muon->isTrackerMuon() && ! muon->isGlobalMuon())))) continue;
140  if (i == 3 && (! useGlobalMuonsNotTrackerMuons_ || (! (muon->isGlobalMuon() && ! muon->isTrackerMuon())))) continue;
141 
142  hNumChambers[i]->Fill(muon->numberOfChambers());
143  hNumMatches[i]->Fill(muon->numberOfMatches(Muon::SegmentAndTrackArbitration));
144  hNumChambersNoRPC[i]->Fill(muon->numberOfChambersNoRPC());
145 
146  // by station
147  for(int station = 0; station < 4; ++station)
148  {
149  // only fill num segments if we crossed (or nearly crossed) a chamber
150  if (muon->trackX(station+1, MuonSubdetId::DT, Muon::NoArbitration) < 900000)
151  hDTNumSegments[i][station]->Fill(muon->numberOfSegments(station+1, MuonSubdetId::DT, Muon::NoArbitration));
152  Fill(hDTDx[i][station], muon->dX(station+1, MuonSubdetId::DT, Muon::SegmentAndTrackArbitration));
153  Fill(hDTPullx[i][station], muon->pullX(station+1, MuonSubdetId::DT, Muon::SegmentAndTrackArbitration, true));
154  Fill(hDTDdXdZ[i][station], muon->dDxDz(station+1, MuonSubdetId::DT, Muon::SegmentAndTrackArbitration));
155  Fill(hDTPulldXdZ[i][station], muon->pullDxDz(station+1, MuonSubdetId::DT, Muon::SegmentAndTrackArbitration, true));
156 
157  if (station < 3) {
158  Fill(hDTDy[i][station], muon->dY(station+1, MuonSubdetId::DT, Muon::SegmentAndTrackArbitration));
159  Fill(hDTPully[i][station], muon->pullY(station+1, MuonSubdetId::DT, Muon::SegmentAndTrackArbitration, true));
160  Fill(hDTDdYdZ[i][station], muon->dDyDz(station+1, MuonSubdetId::DT, Muon::SegmentAndTrackArbitration));
161  Fill(hDTPulldYdZ[i][station], muon->pullDyDz(station+1, MuonSubdetId::DT, Muon::SegmentAndTrackArbitration, true));
162  }
163 
164  // only fill num segments if we crossed (or nearly crossed) a chamber
165  if (muon->trackX(station+1, MuonSubdetId::CSC, Muon::NoArbitration) < 900000)
166  hCSCNumSegments[i][station]->Fill(muon->numberOfSegments(station+1, MuonSubdetId::CSC, Muon::NoArbitration));
167  Fill(hCSCDx[i][station], muon->dX(station+1, MuonSubdetId::CSC, Muon::SegmentAndTrackArbitration));
168  Fill(hCSCPullx[i][station], muon->pullX(station+1, MuonSubdetId::CSC, Muon::SegmentAndTrackArbitration, true));
169  Fill(hCSCDdXdZ[i][station], muon->dDxDz(station+1, MuonSubdetId::CSC, Muon::SegmentAndTrackArbitration));
170  Fill(hCSCPulldXdZ[i][station], muon->pullDxDz(station+1, MuonSubdetId::CSC, Muon::SegmentAndTrackArbitration, true));
171  Fill(hCSCDy[i][station], muon->dY(station+1, MuonSubdetId::CSC, Muon::SegmentAndTrackArbitration));
172  Fill(hCSCPully[i][station], muon->pullY(station+1, MuonSubdetId::CSC, Muon::SegmentAndTrackArbitration, true));
173  Fill(hCSCDdYdZ[i][station], muon->dDyDz(station+1, MuonSubdetId::CSC, Muon::SegmentAndTrackArbitration));
174  Fill(hCSCPulldYdZ[i][station], muon->pullDyDz(station+1, MuonSubdetId::CSC, Muon::SegmentAndTrackArbitration, true));
175  }
176  }
177  }// muon
178 
180  segment != dtSegmentCollectionH_->end(); ++segment)
181  {
182  LocalPoint segmentLocalPosition = segment->localPosition();
183  LocalVector segmentLocalDirection = segment->localDirection();
184  LocalError segmentLocalPositionError = segment->localPositionError();
185  LocalError segmentLocalDirectionError = segment->localDirectionError();
186  bool segmentFound = false;
187 
188  for(MuonCollection::const_iterator muon = muonCollectionH_->begin();
189  muon != muonCollectionH_->end(); ++muon)
190  {
191  if (! muon->isMatchesValid())
192  continue;
193 
194  for(std::vector<MuonChamberMatch>::const_iterator chamberMatch = muon->matches().begin();
195  chamberMatch != muon->matches().end(); ++chamberMatch) {
196  for(std::vector<MuonSegmentMatch>::const_iterator segmentMatch = chamberMatch->segmentMatches.begin();
197  segmentMatch != chamberMatch->segmentMatches.end(); ++segmentMatch)
198  {
199  if (fabs(segmentMatch->x - segmentLocalPosition.x() ) < 1E-6 &&
200  fabs(segmentMatch->y - segmentLocalPosition.y() ) < 1E-6 &&
201  fabs(segmentMatch->dXdZ - segmentLocalDirection.x()/segmentLocalDirection.z()) < 1E-6 &&
202  fabs(segmentMatch->dYdZ - segmentLocalDirection.y()/segmentLocalDirection.z()) < 1E-6 &&
203  fabs(segmentMatch->xErr - sqrt(segmentLocalPositionError.xx()) ) < 1E-6 &&
204  fabs(segmentMatch->yErr - sqrt(segmentLocalPositionError.yy()) ) < 1E-6 &&
205  fabs(segmentMatch->dXdZErr - sqrt(segmentLocalDirectionError.xx()) ) < 1E-6 &&
206  fabs(segmentMatch->dYdZErr - sqrt(segmentLocalDirectionError.yy()) ) < 1E-6)
207  {
208  segmentFound = true;
209  break;
210  }
211  }// segmentMatch
212  if (segmentFound) break;
213  }// chamberMatch
214  if (segmentFound) break;
215  }// muon
216 
217  if (segmentFound)
219  else
221  }// dt segment
222 
224  segment != cscSegmentCollectionH_->end(); ++segment)
225  {
226  LocalPoint segmentLocalPosition = segment->localPosition();
227  LocalVector segmentLocalDirection = segment->localDirection();
228  LocalError segmentLocalPositionError = segment->localPositionError();
229  LocalError segmentLocalDirectionError = segment->localDirectionError();
230  bool segmentFound = false;
231 
232  for(MuonCollection::const_iterator muon = muonCollectionH_->begin();
233  muon != muonCollectionH_->end(); ++muon)
234  {
235  if (! muon->isMatchesValid())
236  continue;
237 
238  for(std::vector<MuonChamberMatch>::const_iterator chamberMatch = muon->matches().begin();
239  chamberMatch != muon->matches().end(); ++chamberMatch) {
240  for(std::vector<MuonSegmentMatch>::const_iterator segmentMatch = chamberMatch->segmentMatches.begin();
241  segmentMatch != chamberMatch->segmentMatches.end(); ++segmentMatch)
242  {
243  if (fabs(segmentMatch->x - segmentLocalPosition.x() ) < 1E-6 &&
244  fabs(segmentMatch->y - segmentLocalPosition.y() ) < 1E-6 &&
245  fabs(segmentMatch->dXdZ - segmentLocalDirection.x()/segmentLocalDirection.z()) < 1E-6 &&
246  fabs(segmentMatch->dYdZ - segmentLocalDirection.y()/segmentLocalDirection.z()) < 1E-6 &&
247  fabs(segmentMatch->xErr - sqrt(segmentLocalPositionError.xx()) ) < 1E-6 &&
248  fabs(segmentMatch->yErr - sqrt(segmentLocalPositionError.yy()) ) < 1E-6 &&
249  fabs(segmentMatch->dXdZErr - sqrt(segmentLocalDirectionError.xx()) ) < 1E-6 &&
250  fabs(segmentMatch->dYdZErr - sqrt(segmentLocalDirectionError.yy()) ) < 1E-6)
251  {
252  segmentFound = true;
253  break;
254  }
255  }// segmentMatch
256  if (segmentFound) break;
257  }// chamberMatch
258  if (segmentFound) break;
259  }// muon
260 
261  if (segmentFound)
263  else
265  }// csc segment
266 }
267 
268 void
270 
271 void MuonIdDQM::Fill(MonitorElement* me, float f) {
272  if (fabs(f) > 900000) return;
273  //if (fabs(f) < 1E-8) return;
274  me->Fill(f);
275 }
276 
277 //define this as a plug-in
bool useTrackerMuons_
Definition: MuonIdDQM.h:68
T getParameter(std::string const &) const
MonitorElement * hCSCDx[4][4]
Definition: MuonIdDQM.h:95
T getUntrackedParameter(std::string const &, T const &) const
int i
Definition: DBlmapReader.cc:9
float xx() const
Definition: LocalError.h:24
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:722
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
MonitorElement * hCSCDdYdZ[4][4]
Definition: MuonIdDQM.h:101
MuonIdDQM(const edm::ParameterSet &)
Definition: MuonIdDQM.cc:3
virtual void endJob()
Definition: MuonIdDQM.cc:269
MonitorElement * hCSCDdXdZ[4][4]
Definition: MuonIdDQM.h:97
MonitorElement * hCSCPulldXdZ[4][4]
Definition: MuonIdDQM.h:98
T y() const
Definition: PV3DBase.h:63
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: MuonIdDQM.cc:122
bool useGlobalMuonsNotTrackerMuons_
Definition: MuonIdDQM.h:71
DQMStore * dbe_
Definition: MuonIdDQM.h:62
edm::InputTag inputDTRecSegment4DCollection_
Definition: MuonIdDQM.h:66
edm::InputTag inputCSCSegmentCollection_
Definition: MuonIdDQM.h:67
bool useTrackerMuonsNotGlobalMuons_
Definition: MuonIdDQM.h:70
virtual void beginJob()
Definition: MuonIdDQM.cc:21
MonitorElement * hDTDdYdZ[4][3]
Definition: MuonIdDQM.h:92
void Fill(long long x)
MonitorElement * hDTPully[4][3]
Definition: MuonIdDQM.h:91
C::const_iterator const_iterator
constant access iterator type
Definition: RangeMap.h:45
MonitorElement * hCSCDy[4][4]
Definition: MuonIdDQM.h:99
int iEvent
Definition: GenABIO.cc:243
static const int CSC
Definition: MuonSubdetId.h:15
edm::Handle< DTRecSegment4DCollection > dtSegmentCollectionH_
Definition: MuonIdDQM.h:75
MonitorElement * hCSCPully[4][4]
Definition: MuonIdDQM.h:100
MonitorElement * hSegmentIsAssociatedBool
Definition: MuonIdDQM.h:105
float yy() const
Definition: LocalError.h:26
MonitorElement * hDTPullx[4][4]
Definition: MuonIdDQM.h:87
MonitorElement * hCSCPullx[4][4]
Definition: MuonIdDQM.h:96
T sqrt(T t)
Definition: SSEVec.h:48
T z() const
Definition: PV3DBase.h:64
edm::InputTag inputMuonCollection_
Definition: MuonIdDQM.h:65
bool useGlobalMuons_
Definition: MuonIdDQM.h:69
double f[11][100]
MonitorElement * hNumMatches[4]
Definition: MuonIdDQM.h:81
MonitorElement * hDTDdXdZ[4][4]
Definition: MuonIdDQM.h:88
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:361
edm::Handle< reco::MuonCollection > muonCollectionH_
Definition: MuonIdDQM.h:74
MonitorElement * hDTDx[4][4]
Definition: MuonIdDQM.h:86
MonitorElement * hDTPulldYdZ[4][3]
Definition: MuonIdDQM.h:93
MonitorElement * hDTDy[4][3]
Definition: MuonIdDQM.h:90
MonitorElement * hCSCPulldYdZ[4][4]
Definition: MuonIdDQM.h:102
const T & get() const
Definition: EventSetup.h:55
virtual void Fill(MonitorElement *, float)
Definition: MuonIdDQM.cc:271
MonitorElement * hDTNumSegments[4][4]
Definition: MuonIdDQM.h:85
MonitorElement * hCSCNumSegments[4][4]
Definition: MuonIdDQM.h:94
MonitorElement * hNumChambersNoRPC[4]
Definition: MuonIdDQM.h:82
edm::ESHandle< GlobalTrackingGeometry > geometry_
Definition: MuonIdDQM.h:77
std::string baseFolder_
Definition: MuonIdDQM.h:72
static const int DT
Definition: MuonSubdetId.h:14
T x() const
Definition: PV3DBase.h:62
MonitorElement * hDTPulldXdZ[4][4]
Definition: MuonIdDQM.h:89
~MuonIdDQM()
Definition: MuonIdDQM.cc:18
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:434
MonitorElement * hNumChambers[4]
Definition: MuonIdDQM.h:80
edm::Handle< CSCSegmentCollection > cscSegmentCollectionH_
Definition: MuonIdDQM.h:76