CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EwkMuDQM.cc
Go to the documentation of this file.
1 #include "DQM/Physics/src/EwkMuDQM.h"
2 
3 #include "FWCore/MessageLogger/interface/MessageLogger.h"
4 #include "FWCore/ParameterSet/interface/ParameterSet.h"
5 #include "DataFormats/Common/interface/Handle.h"
6 
7 #include "FWCore/ServiceRegistry/interface/Service.h"
8 #include "CommonTools/UtilAlgos/interface/TFileService.h"
9 
10 #include "DQMServices/Core/interface/DQMStore.h"
11 #include "DQMServices/Core/interface/MonitorElement.h"
12 
13 #include "DataFormats/TrackReco/interface/Track.h"
14 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
15 #include "DataFormats/VertexReco/interface/Vertex.h"
16 
17 #include "DataFormats/MuonReco/interface/Muon.h"
18 #include "DataFormats/MuonReco/interface/MuonSelectors.h"
19 #include "DataFormats/METReco/interface/MET.h"
20 #include "DataFormats/JetReco/interface/Jet.h"
21 #include "DataFormats/EgammaCandidates/interface/Photon.h"
22 
23 #include "DataFormats/GeometryVector/interface/Phi.h"
24 
25 #include "FWCore/Common/interface/TriggerNames.h"
26 #include "FWCore/Framework/interface/Event.h"
27 #include "DataFormats/Common/interface/TriggerResults.h"
28 
29 #include "DataFormats/Common/interface/View.h"
30 
31 #include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
32 
33 using namespace edm;
34 using namespace std;
35 using namespace reco;
36 
37 EwkMuDQM::EwkMuDQM(const ParameterSet& cfg)
38  :
39  // Input collections
40  metTag_(cfg.getUntrackedParameter<edm::InputTag>("METTag",
41  edm::InputTag("pfmet"))),
42  jetTag_(cfg.getUntrackedParameter<edm::InputTag>(
43  "JetTag", edm::InputTag("ak4PFJets"))),
44  trigTag_(consumes<edm::TriggerResults>(
45  cfg.getUntrackedParameter<edm::InputTag>(
46  "TrigTag", edm::InputTag("TriggerResults::HLT")))),
47  muonTag_(consumes<edm::View<reco::Muon> >(
48  cfg.getUntrackedParameter<edm::InputTag>("MuonTag",
49  edm::InputTag("muons")))),
50  metToken_(consumes<edm::View<reco::MET> >(
51  cfg.getUntrackedParameter<edm::InputTag>("METTag",
52  edm::InputTag("pfmet")))),
53  jetToken_(consumes<edm::View<reco::Jet> >(
54  cfg.getUntrackedParameter<edm::InputTag>(
55  "JetTag", edm::InputTag("ak4PFJets")))),
56  phoTag_(consumes<edm::View<reco::Photon> >(
57  cfg.getUntrackedParameter<edm::InputTag>("phoTag",
58  edm::InputTag("photons")))),
59  vertexTag_(consumes<edm::View<reco::Vertex> >(
60  cfg.getUntrackedParameter<edm::InputTag>(
61  "VertexTag", edm::InputTag("offlinePrimaryVertices")))),
62  beamSpotTag_(
63  consumes<reco::BeamSpot>(cfg.getUntrackedParameter<edm::InputTag>(
64  "beamSpotTag", edm::InputTag("offlineBeamSpot")))),
65  trigPathNames_(cfg.getUntrackedParameter<std::vector<std::string> >(
66  "TrigPathNames")),
67 
68  // Muon quality cuts
69  isAlsoTrackerMuon_(cfg.getUntrackedParameter<bool>(
70  "IsAlsoTrackerMuon", true)), // Glb muon also tracker muon
71  dxyCut_(
72  cfg.getUntrackedParameter<double>("DxyCut", 0.2)), // dxy < 0.2 cm
73  normalizedChi2Cut_(cfg.getUntrackedParameter<double>(
74  "NormalizedChi2Cut", 10.)), // chi2/ndof (of global fit) <10.0
75  trackerHitsCut_(cfg.getUntrackedParameter<int>("TrackerHitsCut",
76  11)), // Tracker Hits >10
77  pixelHitsCut_(
78  cfg.getUntrackedParameter<int>("PixelHitsCut", 1)), // Pixel Hits >0
79  muonHitsCut_(cfg.getUntrackedParameter<int>("MuonHitsCut",
80  1)), // Valid Muon Hits >0
81  nMatchesCut_(cfg.getUntrackedParameter<int>(
82  "NMatchesCut", 2)), // At least 2 Chambers with matches
83 
84  // W-boson cuts
85  isRelativeIso_(cfg.getUntrackedParameter<bool>("IsRelativeIso", true)),
86  isCombinedIso_(cfg.getUntrackedParameter<bool>("IsCombinedIso", false)),
87  isoCut03_(cfg.getUntrackedParameter<double>("IsoCut03", 0.1)),
88  acopCut_(cfg.getUntrackedParameter<double>("AcopCut", 999.)),
89  metMin_(cfg.getUntrackedParameter<double>("MetMin", -999999.)),
90  metMax_(cfg.getUntrackedParameter<double>("MetMax", 999999.)),
91  mtMin_(cfg.getUntrackedParameter<double>("MtMin", 50.)),
92  mtMax_(cfg.getUntrackedParameter<double>("MtMax", 200.)),
93  ptCut_(cfg.getUntrackedParameter<double>("PtCut", 20.)),
94  etaCut_(cfg.getUntrackedParameter<double>("EtaCut", 2.4)),
95 
96  // Z rejection
97  ptThrForZ1_(cfg.getUntrackedParameter<double>("PtThrForZ1", 20.)),
98  ptThrForZ2_(cfg.getUntrackedParameter<double>("PtThrForZ2", 10.)),
99 
100  // Z selection
101  dimuonMassMin_(cfg.getUntrackedParameter<double>("dimuonMassMin", 80.)),
102  dimuonMassMax_(cfg.getUntrackedParameter<double>("dimuonMassMax", 120.)),
103 
104  // Top rejection
105  eJetMin_(cfg.getUntrackedParameter<double>("EJetMin", 999999.)),
106  nJetMax_(cfg.getUntrackedParameter<int>("NJetMax", 999999)),
107 
108  // Photon cuts
109  ptThrForPhoton_(cfg.getUntrackedParameter<double>("ptThrForPhoton", 5.)),
110  nPhoMax_(cfg.getUntrackedParameter<int>("nPhoMax", 999999)) {
111  isValidHltConfig_ = false;
112 
113 }
114 
115 void EwkMuDQM::dqmBeginRun(const Run& iRun, const EventSetup& iSet) {
116  nall = 0;
117  nsel = 0;
118  nz = 0;
119 
120  nrec = 0;
121  niso = 0;
122  nhlt = 0;
123  nmet = 0;
124 
125  // passed as parameter to HLTConfigProvider::init(), not yet used
126  bool isConfigChanged = false;
127  // isValidHltConfig_ used to short-circuit analyze() in case of problems
128  isValidHltConfig_ =
129  hltConfigProvider_.init(iRun, iSet, "HLT", isConfigChanged);
130 }
131 
132 void EwkMuDQM::bookHistograms(DQMStore::IBooker & ibooker,
133  edm::Run const &, edm::EventSetup const &) {
134 
135  ibooker.setCurrentFolder("Physics/EwkMuDQM");
136 
137  char chtitle[256] = "";
138 
139  pt_before_ = ibooker.book1D("PT_BEFORECUTS",
140  "Muon transverse momentum (global muon) [GeV]", 100, 0., 100.);
141  pt_after_ = ibooker.book1D("PT_AFTERWCUTS",
142  "Muon transverse momentum (global muon) [GeV]", 100, 0., 100.);
143 
144  eta_before_ = ibooker.book1D("ETA_BEFORECUTS",
145  "Muon pseudo-rapidity", 50, -2.5, 2.5);
146  eta_after_ = ibooker.book1D("ETA_AFTERWCUTS",
147  "Muon pseudo-rapidity", 50, -2.5, 2.5);
148 
149  dxy_before_ = ibooker.book1D("DXY_BEFORECUTS",
150  "Muon transverse distance to beam spot [cm]", 100, -0.5, 0.5);
151  dxy_after_ = ibooker.book1D("DXY_AFTERWCUTS",
152  "Muon transverse distance to beam spot [cm]", 100, -0.5, 0.5);
153 
154  goodewkmuon_before_ = ibooker.book1D("GOODEWKMUON_BEFORECUTS",
155  "Quality-muon flag", 2, -0.5, 1.5);
156  goodewkmuon_after_ = ibooker.book1D("GOODEWKMUON_AFTERWCUTS",
157  "Quality-muon flag", 2, -0.5, 1.5);
158 
159  if (isRelativeIso_) {
160  if (isCombinedIso_) {
161  iso_before_ = ibooker.book1D("ISO_BEFORECUTS",
162  "Relative (combined) isolation variable", 100, 0., 1.);
163  iso_after_ = ibooker.book1D("ISO_AFTERWCUTS",
164  "Relative (combined) isolation variable", 100, 0., 1.);
165  } else {
166  iso_before_ = ibooker.book1D("ISO_BEFORECUTS",
167  "Relative (tracker) isolation variable", 100, 0., 1.);
168  iso_after_ = ibooker.book1D("ISO_AFTERWCUTS",
169  "Relative (tracker) isolation variable", 100, 0., 1.);
170  }
171  } else {
172  if (isCombinedIso_) {
173  iso_before_ = ibooker.book1D("ISO_BEFORECUTS",
174  "Absolute (combined) isolation variable [GeV]", 100, 0., 20.);
175  iso_after_ = ibooker.book1D("ISO_AFTERWCUTS",
176  "Absolute (combined) isolation variable [GeV]", 100, 0., 20.);
177  } else {
178  iso_before_ = ibooker.book1D("ISO_BEFORECUTS",
179  "Absolute (tracker) isolation variable [GeV]", 100, 0., 20.);
180  iso_after_ = ibooker.book1D("ISO_AFTERWCUTS",
181  "Absolute (tracker) isolation variable [GeV]", 100, 0., 20.);
182  }
183  }
184 
185  trig_before_ = ibooker.book1D("TRIG_BEFORECUTS",
186  "Trigger response (boolean of muon triggers)", 2, -0.5, 1.5);
187  trig_after_ = ibooker.book1D("TRIG_AFTERWCUTS",
188  "Trigger response (boolean of muon triggers)", 2, -0.5, 1.5);
189 
190  snprintf(chtitle, 255, "Transverse mass (%s) [GeV]", metTag_.label().data());
191  mt_before_ = ibooker.book1D("MT_BEFORECUTS", chtitle, 150, 0., 300.);
192  mt_after_ = ibooker.book1D("MT_AFTERWCUTS", chtitle, 150, 0., 300.);
193 
194  snprintf(chtitle, 255, "Missing transverse energy (%s) [GeV]",
195  metTag_.label().data());
196  met_before_ = ibooker.book1D("MET_BEFORECUTS", chtitle, 100, 0., 200.);
197  met_after_ = ibooker.book1D("MET_AFTERWCUTS", chtitle, 100, 0., 200.);
198  met_afterZ_ = ibooker.book1D("MET_AFTERZCUTS", chtitle, 100, 0., 200.);
199 
200  snprintf(chtitle, 255, "MU-MET (%s) acoplanarity", metTag_.label().data());
201  acop_before_ = ibooker.book1D("ACOP_BEFORECUTS", chtitle, 50, 0., M_PI);
202  acop_after_ = ibooker.book1D("ACOP_AFTERWCUTS", chtitle, 50, 0., M_PI);
203 
204  snprintf(chtitle, 255, "Z selection: muons above %.2f GeV", ptThrForZ1_);
205  n_zselPt1thr_ = ibooker.book1D("NZSELPT1THR", chtitle, 10, -0.5, 9.5);
206  snprintf(chtitle, 255, "Z selection: muons above %.2f GeV", ptThrForZ2_);
207  n_zselPt2thr_ = ibooker.book1D("NZSELPT2THR", chtitle, 10, -0.5, 9.5);
208 
209  snprintf(chtitle, 255, "Number of jets (%s) above %.2f GeV",
210  jetTag_.label().data(), eJetMin_);
211  njets_before_ = ibooker.book1D("NJETS_BEFORECUTS", chtitle, 16, -0.5, 15.5);
212  njets_after_ = ibooker.book1D("NJETS_AFTERWCUTS", chtitle, 16, -0.5, 15.5);
213  njets_afterZ_ = ibooker.book1D("NJETS_AFTERZCUTS", chtitle, 16, -0.5, 15.5);
214 
215  leadingjet_pt_before_ = ibooker.book1D("LEADINGJET_PT_BEFORECUTS",
216  "Leading Jet transverse momentum", 300, 0., 300.);
217  leadingjet_pt_after_ = ibooker.book1D("LEADINGJET_PT_AFTERWCUTS",
218  "Leading Jet transverse momentum", 300, 0., 300.);
219  leadingjet_pt_afterZ_ = ibooker.book1D("LEADINGJET_PT_AFTERZCUTS",
220  "Leading Jet transverse momentum", 300, 0., 300.);
221 
222  leadingjet_eta_before_ = ibooker.book1D("LEADINGJET_ETA_BEFORECUTS",
223  "Leading Jet pseudo-rapidity", 50, -2.5, 2.5);
224  leadingjet_eta_after_ = ibooker.book1D("LEADINGJET_ETA_AFTERWCUTS",
225  "Leading Jet pseudo-rapidity", 50, -2.5, 2.5);
226  leadingjet_eta_afterZ_ = ibooker.book1D("LEADINGJET_ETA_AFTERZCUTS",
227  "Leading Jet pseudo-rapidity", 50, -2.5, 2.5);
228 
229  ptDiffPM_before_ = ibooker.book1D("PTDIFFPM_BEFORE_CUTS",
230  "pt(Muon+)-pt(Muon-) after Z cuts [GeV]", 200, -100., 100.);
231  ptDiffPM_afterZ_ = ibooker.book1D("PTDIFFPM_AFTERZ_CUTS",
232  "pt(Muon+)-pt(Muon-) after Z cuts [GeV]", 200, -100., 100.);
233 
236  pt1_afterZ_ = ibooker.book1D("PT1_AFTERZCUTS",
237  "Muon transverse momentum (global muon) [GeV]", 100, 0., 100.);
238  eta1_afterZ_ = ibooker.book1D("ETA1_AFTERZCUTS",
239  "Muon pseudo-rapidity", 50, -2.5, 2.5);
240  dxy1_afterZ_ = ibooker.book1D("DXY1_AFTERZCUTS",
241  "Muon transverse distance to beam spot [cm]", 100, -0.5, 0.5);
242  goodewkmuon1_afterZ_ = ibooker.book1D("GOODEWKMUON1_AFTERZCUTS",
243  "Quality-muon flag", 2, -0.5, 1.5);
244 
245  if (isRelativeIso_) {
246  if (isCombinedIso_) {
247  iso1_afterZ_ = ibooker.book1D("ISO1_AFTERZCUTS",
248  "Relative (combined) isolation variable", 100, 0., 1.);
249  iso2_afterZ_ = ibooker.book1D("ISO2_AFTERZCUTS",
250  "Relative (combined) isolation variable", 100, 0., 1.);
251  } else {
252  iso1_afterZ_ = ibooker.book1D("ISO1_AFTERZCUTS",
253  "Relative (tracker) isolation variable", 100, 0., 1.);
254  iso2_afterZ_ = ibooker.book1D("ISO2_AFTERZCUTS",
255  "Relative (tracker) isolation variable", 100, 0., 1.);
256  }
257  } else {
258  if (isCombinedIso_) {
259  iso1_afterZ_ = ibooker.book1D("ISO1_AFTERZCUTS",
260  "Absolute (combined) isolation variable [GeV]", 100, 0., 20.);
261  iso2_afterZ_ = ibooker.book1D("ISO2_AFTERZCUTS",
262  "Absolute (combined) isolation variable [GeV]", 100, 0., 20.);
263  } else {
264  iso1_afterZ_ = ibooker.book1D("ISO1_AFTERZCUTS",
265  "Absolute (tracker) isolation variable [GeV]", 100, 0., 20.);
266  iso2_afterZ_ = ibooker.book1D("ISO2_AFTERZCUTS",
267  "Absolute (tracker) isolation variable [GeV]", 100, 0., 20.);
268  }
269  }
270 
271  pt2_afterZ_ = ibooker.book1D("PT2_AFTERZCUTS",
272  "Muon transverse momentum (global muon) [GeV]", 100, 0., 100.);
273  eta2_afterZ_ = ibooker.book1D("ETA2_AFTERZCUTS",
274  "Muon pseudo-rapidity", 50, -2.5, 2.5);
275  dxy2_afterZ_ = ibooker.book1D("DXY2_AFTERZCUTS",
276  "Muon transverse distance to beam spot [cm]", 100, -0.5, 0.5);
277  goodewkmuon2_afterZ_ = ibooker.book1D("GOODEWKMUON2_AFTERZCUTS",
278  "Quality-muon flag", 2, -0.5, 1.5);
279  ztrig_afterZ_ = ibooker.book1D("ZTRIG_AFTERZCUTS",
280  "Trigger response (boolean of muon triggers)", 2, -0.5, 1.5);
281  dimuonmass_before_ = ibooker.book1D("DIMUONMASS_BEFORECUTS",
282  "DiMuonMass (2 globals)", 100, 0, 200);
283  dimuonmass_afterZ_ = ibooker.book1D("DIMUONMASS_AFTERZCUTS",
284  "DiMuonMass (2 globals)", 100, 0, 200);
285  npvs_before_ = ibooker.book1D("NPVs_BEFORECUTS",
286  "Number of Valid Primary Vertices", 51, -0.5, 50.5);
287  npvs_after_ = ibooker.book1D("NPVs_AFTERWCUTS",
288  "Number of Valid Primary Vertices", 51, -0.5, 50.5);
289  npvs_afterZ_ = ibooker.book1D("NPVs_AFTERZCUTS",
290  "Number of Valid Primary Vertices", 51, -0.5, 50.5);
291  muoncharge_before_ = ibooker.book1D("MUONCHARGE_BEFORECUTS",
292  "Muon Charge", 3, -1.5, 1.5);
293  muoncharge_after_ = ibooker.book1D("MUONCHARGE_AFTERWCUTS",
294  "Muon Charge", 3, -1.5, 1.5);
295  muoncharge_afterZ_ = ibooker.book1D("MUONCHARGE_AFTERZCUTS",
296  "Muon Charge", 3, -1.5, 1.5);
297 
298  // Adding these to replace the NZ ones (more useful, since they are more
299  // general?)
300  nmuons_ = ibooker.book1D("NMuons",
301  "Number of muons in the event", 10, -0.5, 9.5);
302  ngoodmuons_ = ibooker.book1D("NGoodMuons",
303  "Number of muons passing the quality criteria", 10, -0.5, 9.5);
304 
305  nph_ = ibooker.book1D("nph", "Number of photons in the event", 20, 0., 20.);
306  phPt_ = ibooker.book1D("phPt", "Photon transverse momentum [GeV]", 100, 0.,1000.);
307  snprintf(chtitle, 255, "Photon pseudorapidity (pT>%4.1f)", ptThrForPhoton_);
308  phEta_ = ibooker.book1D("phEta", chtitle, 100, -2.5, 2.5);
309 }
310 
311 void EwkMuDQM::endRun(const Run& r, const EventSetup& iSet) {}
312 
313 void EwkMuDQM::analyze(const Event& ev, const EventSetup& iSet) {
314 
315  // Muon collection
316  Handle<View<Muon> > muonCollection;
317  if (!ev.getByToken(muonTag_, muonCollection)) {
318  // LogWarning("") << ">>> Muon collection does not exist !!!";
319  return;
320  }
321  unsigned int muonCollectionSize = muonCollection->size();
322 
323  // Beam spot
324  Handle<reco::BeamSpot> beamSpotHandle;
325  if (!ev.getByToken(beamSpotTag_, beamSpotHandle)) {
326  // LogWarning("") << ">>> No beam spot found !!!";
327  return;
328  }
329 
330  // Loop to reject/control Z->mumu is done separately
331  unsigned int nmuonsForZ1 = 0;
332  unsigned int nmuonsForZ2 = 0;
333  bool cosmic = false;
334  for (unsigned int i = 0; i < muonCollectionSize; i++) {
335  const Muon& mu = muonCollection->at(i);
336  if (!mu.isGlobalMuon()) continue;
337  double pt = mu.pt();
338  double dxy = mu.innerTrack()->dxy(beamSpotHandle->position());
339 
340  if (fabs(dxy) > 1) {
341  cosmic = true;
342  break;
343  }
344 
345  if (pt > ptThrForZ1_) nmuonsForZ1++;
346  if (pt > ptThrForZ2_) nmuonsForZ2++;
347 
348  for (unsigned int j = i + 1; j < muonCollectionSize; j++) {
349  const Muon& mu2 = muonCollection->at(j);
350  if (mu2.isGlobalMuon() && (mu.charge() * mu2.charge() == -1)) {
351  const math::XYZTLorentzVector ZRecoGlb(
352  mu.px() + mu2.px(), mu.py() + mu2.py(), mu.pz() + mu2.pz(),
353  mu.p() + mu2.p());
354  dimuonmass_before_->Fill(ZRecoGlb.mass());
355  if (mu.charge() > 0) {
356  ptDiffPM_before_->Fill(mu.pt() - mu2.pt());
357  } else {
358  ptDiffPM_before_->Fill(mu2.pt() - mu.pt());
359  }
360  }
361  }
362  }
363  if (cosmic) return;
364 
365  LogTrace("") << "> Z rejection: muons above " << ptThrForZ1_
366  << " [GeV]: " << nmuonsForZ1;
367  LogTrace("") << "> Z rejection: muons above " << ptThrForZ2_
368  << " [GeV]: " << nmuonsForZ2;
369 
370  // MET
371  Handle<View<MET> > metCollection;
372  if (!ev.getByToken(metToken_, metCollection)) {
373  // LogWarning("") << ">>> MET collection does not exist !!!";
374  return;
375  }
376  const MET& met = metCollection->at(0);
377  double met_et = met.pt();
378  LogTrace("") << ">>> MET, MET_px, MET_py: " << met_et << ", " << met.px()
379  << ", " << met.py() << " [GeV]";
380  met_before_->Fill(met_et);
381 
382  // Vertices in the event
383  Handle<View<reco::Vertex> > vertexCollection;
384  if (!ev.getByToken(vertexTag_, vertexCollection)) {
385  LogError("") << ">>> Vertex collection does not exist !!!";
386  return;
387  }
388  unsigned int vertexCollectionSize = vertexCollection->size();
389 
390  int nvvertex = 0;
391  for (unsigned int i = 0; i < vertexCollectionSize; i++) {
392  const Vertex& vertex = vertexCollection->at(i);
393  if (vertex.isValid()) nvvertex++;
394  }
395 
396  npvs_before_->Fill(nvvertex);
397 
398  bool trigger_fired = false;
399  Handle<TriggerResults> triggerResults;
400  if (!ev.getByToken(trigTag_, triggerResults)) {
401  // LogWarning("") << ">>> TRIGGER collection does not exist !!!";
402  return;
403  }
404  const edm::TriggerNames& trigNames = ev.triggerNames(*triggerResults);
405  // LogWarning("")<<"Loop over triggers";
406 
407  /* change faulty logic of triggering
408  for (unsigned int i=0; i<triggerResults->size(); i++)
409  {
410  const std::string trigName = trigNames.triggerName(i);
411 
412  bool found=false;
413  for(unsigned int index=0; index<trigPathNames_.size() && found==false;
414  index++) {
415  size_t trigPath = trigName.find(trigPathNames_[index]); // 0 if
416  found, pos if not
417  if (trigPath==0) found=true;
418  }
419  if(!found) {continue;}
420 
421  bool prescaled=false;
422  for (unsigned int ps= 0; ps< hltConfigProvider_.prescaleSize();
423  ps++){
424  const unsigned int prescaleValue =
425  hltConfigProvider_.prescaleValue(ps, trigName) ;
426  if (prescaleValue != 1) prescaled =true;
427  }
428 
429  if( triggerResults->accept(i) && !prescaled){ trigger_fired=true;}
430  // LogWarning("")<<"TrigNo: "<<i<<" "<<found<<"
431  "<<trigName<<" ---> FIRED";}
432  }
433  */
434 
435  // get the prescale set for this event
436  const int prescaleSet = hltConfigProvider_.prescaleSet(ev, iSet);
437  if (prescaleSet == -1) {
438  LogTrace("") << "Failed to determine prescaleSet\n";
439  // std::cout << "Failed to determine prescaleSet. Check the GlobalTag in
440  // cfg\n";
441  return;
442  }
443 
444  for (unsigned int i = 0;
445  (i < triggerResults->size()) && (trigger_fired == false); i++) {
446  // skip trigger, if it did not fire
447  if (!triggerResults->accept(i)) continue;
448 
449  // skip trigger, if it is not on our list
450  bool found = false;
451  const std::string trigName = trigNames.triggerName(i);
452  for (unsigned int index = 0;
453  index < trigPathNames_.size() && found == false; index++) {
454  if (trigName.find(trigPathNames_.at(index)) == 0) found = true;
455  }
456  if (!found) continue;
457 
458  // skip trigger, if it is prescaled
459  if (prescaleSet != -1) {
460  if (hltConfigProvider_.prescaleValue(prescaleSet, trigName) != 1)
461  continue;
462  } else {
463  // prescaleSet is not known.
464  // This branch is not needed, if prescaleSet=-1 forces to skip event
465  int prescaled = 0;
466  for (unsigned int ps = 0;
467  !prescaled && (ps < hltConfigProvider_.prescaleSize()); ++ps) {
468  if (hltConfigProvider_.prescaleValue(ps, trigName) != 1) {
469  prescaled = 1;
470  }
471  }
472  if (prescaled) {
473  // std::cout << "trigger prescaled\n";
474  continue;
475  }
476  }
477 
478  // std::cout << "found unprescaled trigger that fired: " << trigName <<
479  // "\n";
480  trigger_fired = true;
481  }
482  // if (trigger_fired) std::cout << "\n\tGot Trigger\n";
483 
484  trig_before_->Fill(trigger_fired);
485 
486  // Jet collection
487  Handle<View<Jet> > jetCollection;
488  if (!ev.getByToken(jetToken_, jetCollection)) {
489  // LogError("") << ">>> JET collection does not exist !!!";
490  return;
491  }
492  unsigned int jetCollectionSize = jetCollection->size();
493  int njets = 0;
494  int LEADJET = -1;
495  double max_pt = 0;
496  for (unsigned int i = 0; i < jetCollectionSize; i++) {
497  const Jet& jet = jetCollection->at(i);
498  double minDistance = 99999; // This is in order to use PFJets
499  for (unsigned int j = 0; j < muonCollectionSize; j++) {
500  const Muon& mu = muonCollection->at(j);
501  double distance = sqrt((mu.eta() - jet.eta()) * (mu.eta() - jet.eta()) +
502  (mu.phi() - jet.phi()) * (mu.phi() - jet.phi()));
503  if (minDistance > distance) minDistance = distance;
504  }
505  if (minDistance < 0.3)
506  continue; // 0.3 is the isolation cone around the muon
507  if (jet.et() > max_pt) {
508  LEADJET = i;
509  max_pt = jet.et();
510  }
511  if (jet.et() > eJetMin_) {
512  njets++;
513  }
514  }
515 
516  LogTrace("") << ">>> Total number of jets: " << jetCollectionSize;
517  LogTrace("") << ">>> Number of jets above " << eJetMin_
518  << " [GeV]: " << njets;
519  njets_before_->Fill(njets);
520  double lead_jet_pt = -1;
521  double lead_jet_eta = -100;
522  if (LEADJET != -1) {
523  const Jet& leadJet = jetCollection->at(LEADJET);
524  leadingjet_pt_before_->Fill(leadJet.pt());
525  leadingjet_eta_before_->Fill(leadJet.eta());
526  lead_jet_pt = leadJet.pt();
527  lead_jet_eta = leadJet.eta();
528  }
529  // Photon Collection
530  Handle<View<Photon> > photonCollection;
531  if (!ev.getByToken(phoTag_, photonCollection)) {
532  // LogError("")
533  return;
534  }
535  unsigned int ngam = 0;
536 
537  for (unsigned int i = 0; i < photonCollection->size(); i++) {
538  const Photon& ph = photonCollection->at(i);
539  double photonPt = ph.pt();
540  if (photonPt > ptThrForPhoton_) {
541  ngam++;
542  phEta_->Fill(ph.eta());
543  }
544  phPt_->Fill(photonPt);
545  }
546  nph_->Fill(ngam);
547  LogTrace("") << " >>> N photons " << ngam << std::endl;
548 
549  nmuons_->Fill(muonCollectionSize);
550 
551  // Start counting
552  nall++;
553 
554  // Histograms per event should be done only once, so keep track of them
555  bool hlt_hist_done = false;
556  bool zhlt_hist_done = false;
557  bool zjets_hist_done = false;
558  bool zfullsel_hist_done = false;
559  bool met_hist_done = false;
560  bool njets_hist_done = false;
561  bool wfullsel_hist_done = false;
562 
563  // Central W->mu nu selection criteria
564  const int NFLAGS = 11;
565  bool muon_sel[NFLAGS];
566  const int NFLAGSZ = 13;
567  bool zmuon_sel[NFLAGSZ];
568  bool muon4Z = false;
569 
570  double number_of_muons = 0;
571  double number_of_goodMuons = 0;
572 
573  for (unsigned int i = 0; i < muonCollectionSize; i++) {
574  for (int j = 0; j < NFLAGS; ++j) {
575  muon_sel[j] = false;
576  }
577 
578  number_of_muons++;
579 
580  const Muon& mu = muonCollection->at(i);
581  if (!mu.isGlobalMuon()) continue;
582  if (mu.globalTrack().isNull()) continue;
583  if (mu.innerTrack().isNull()) continue;
584 
585  LogTrace("") << "> Wsel: processing muon number " << i << "...";
586  reco::TrackRef gm = mu.globalTrack();
587  reco::TrackRef tk = mu.innerTrack();
588 
589  // Pt,eta cuts
590  double pt = mu.pt();
591  double eta = mu.eta();
592  LogTrace("") << "\t... pt, eta: " << pt << " [GeV], " << eta;
593  ;
594  if (pt > ptCut_) muon_sel[0] = true;
595  if (fabs(eta) < etaCut_) muon_sel[1] = true;
596 
597  double charge = mu.charge();
598 
599  // d0, chi2, nhits quality cuts
600  double dxy = gm->dxy(beamSpotHandle->position());
601  double normalizedChi2 = gm->normalizedChi2();
602  double trackerHits = tk->hitPattern().numberOfValidTrackerHits();
603  int pixelHits = tk->hitPattern().numberOfValidPixelHits();
604  int muonHits = gm->hitPattern().numberOfValidMuonHits();
605  int nMatches = mu.numberOfMatches();
606 
607  LogTrace("") << "\t... dxy, normalizedChi2, trackerHits, isTrackerMuon?: "
608  << dxy << " [cm], " << normalizedChi2 << ", " << trackerHits
609  << ", " << mu.isTrackerMuon();
610  if (fabs(dxy) < dxyCut_) muon_sel[2] = true;
611 
612  bool quality = true;
613 
614  if (normalizedChi2 > normalizedChi2Cut_) quality = false;
615  if (trackerHits < trackerHitsCut_) quality = false;
616  if (pixelHits < pixelHitsCut_) quality = false;
617  if (muonHits < muonHitsCut_) quality = false;
618  ;
619  if (!mu.isTrackerMuon()) quality = false;
620  if (nMatches < nMatchesCut_) quality = false;
621  muon_sel[3] = quality;
622  if (quality) number_of_goodMuons++;
623 
624  pt_before_->Fill(pt);
625  eta_before_->Fill(eta);
626  dxy_before_->Fill(dxy);
627  muoncharge_before_->Fill(charge);
628  goodewkmuon_before_->Fill(quality);
629 
630  // Charge asymmetry
631  // if (quality) {
632  // if (charge>0) ptPlus_before_->Fill(pt);
633  // if (charge<0) ptMinus_before_->Fill(pt);
634  //}
635 
636  // Isolation cuts
637  double isovar = mu.isolationR03().sumPt;
638  if (isCombinedIso_) {
639  isovar += mu.isolationR03().emEt;
640  isovar += mu.isolationR03().hadEt;
641  }
642  if (isRelativeIso_) isovar /= pt;
643  if (isovar < isoCut03_) muon_sel[4] = true;
644 
645  LogTrace("") << "\t... isolation value" << isovar << ", isolated? "
646  << muon_sel[6];
647  iso_before_->Fill(isovar);
648 
649  // HLT (not mtched to muon for the time being)
650  if (trigger_fired) muon_sel[5] = true;
651 
652  // For Z:
653  if (pt > ptThrForZ1_ && fabs(eta) < etaCut_ && fabs(dxy) < dxyCut_ &&
654  quality && trigger_fired && isovar < isoCut03_) {
655  muon4Z = true;
656  }
657 
658  // MET/MT cuts
659  double w_et = met_et + mu.pt();
660  double w_px = met.px() + mu.px();
661  double w_py = met.py() + mu.py();
662 
663  double massT = w_et * w_et - w_px * w_px - w_py * w_py;
664  massT = (massT > 0) ? sqrt(massT) : 0;
665 
666  LogTrace("") << "\t... W mass, W_et, W_px, W_py: " << massT << ", " << w_et
667  << ", " << w_px << ", " << w_py << " [GeV]";
668  if (massT > mtMin_ && massT < mtMax_) muon_sel[6] = true;
669  mt_before_->Fill(massT);
670  if (met_et > metMin_ && met_et < metMax_) muon_sel[7] = true;
671 
672  // Acoplanarity cuts
673  Geom::Phi<double> deltaphi(mu.phi() - atan2(met.py(), met.px()));
674  double acop = deltaphi.value();
675  if (acop < 0) acop = -acop;
676  acop = M_PI - acop;
677  LogTrace("") << "\t... acoplanarity: " << acop;
678  if (acop < acopCut_) muon_sel[8] = true;
679  acop_before_->Fill(acop);
680 
681  // Remaining flags (from global event information)
682  if (nmuonsForZ1 < 1 || nmuonsForZ2 < 2) muon_sel[9] = true;
683  if (njets <= nJetMax_) muon_sel[10] = true;
684 
685  // Collect necessary flags "per muon"
686  int flags_passed = 0;
687  for (int j = 0; j < NFLAGS; ++j) {
688  if (muon_sel[j]) flags_passed += 1;
689  }
690 
691  // Do N-1 histograms now (and only once for global event quantities)
692  if (flags_passed >= (NFLAGS - 1)) {
693  if (!muon_sel[0] || flags_passed == NFLAGS) pt_after_->Fill(pt);
694  if (!muon_sel[1] || flags_passed == NFLAGS) eta_after_->Fill(eta);
695  if (!muon_sel[2] || flags_passed == NFLAGS) dxy_after_->Fill(dxy);
696  if (!muon_sel[3] || flags_passed == NFLAGS)
697  goodewkmuon_after_->Fill(quality);
698  if (!muon_sel[4] || flags_passed == NFLAGS) iso_after_->Fill(isovar);
699  if (!muon_sel[5] || flags_passed == NFLAGS)
700  if (!hlt_hist_done) trig_after_->Fill(trigger_fired);
701  hlt_hist_done = true;
702  if (!muon_sel[6] || flags_passed == NFLAGS) mt_after_->Fill(massT);
703  if (!muon_sel[7] || flags_passed == NFLAGS)
704  if (!met_hist_done) met_after_->Fill(met_et);
705  met_hist_done = true;
706  if (!muon_sel[8] || flags_passed == NFLAGS) acop_after_->Fill(acop);
707  // no action here for muon_sel[9]
708  if (!muon_sel[10] || flags_passed == NFLAGS) {
709  if (!njets_hist_done) {
710  njets_after_->Fill(njets);
711  leadingjet_pt_after_->Fill(lead_jet_pt);
712  leadingjet_eta_after_->Fill(lead_jet_eta);
713  }
714  njets_hist_done = true;
715  }
716  if (flags_passed == NFLAGS) {
717  if (!wfullsel_hist_done) {
718  npvs_after_->Fill(nvvertex);
719  muoncharge_after_->Fill(charge);
720  // if (charge>0) ptPlus_afterW_->Fill(pt);
721  // if (charge<0) ptMinus_afterW_->Fill(pt);
722  }
723  wfullsel_hist_done = true;
724  }
725  }
726 
727  // The cases in which the event is rejected as a Z are considered
728  // independently:
729  if (muon4Z && !muon_sel[9]) {
730 
731  // Plots for 2 muons
732  for (unsigned int j = i + 1; j < muonCollectionSize; j++) {
733 
734  for (int ij = 0; ij < NFLAGSZ; ++ij) {
735  zmuon_sel[ij] = false;
736  }
737 
738  for (int ji = 0; ji < 5; ++ji) {
739  zmuon_sel[ji] = muon_sel[ji];
740  }
741 
742  const Muon& mu2 = muonCollection->at(j);
743  if (!mu2.isGlobalMuon()) continue;
744  if (mu2.charge() * charge != -1) continue;
745  reco::TrackRef gm2 = mu2.globalTrack();
746  reco::TrackRef tk2 = mu2.innerTrack();
747  double pt2 = mu2.pt();
748  if (pt2 > ptThrForZ2_) zmuon_sel[5] = true;
749  double eta2 = mu2.eta();
750  if (fabs(eta2) < etaCut_) zmuon_sel[6] = true;
751  double dxy2 = gm2->dxy(beamSpotHandle->position());
752  if (fabs(dxy2) < dxyCut_) zmuon_sel[7] = true;
753  double normalizedChi22 = gm2->normalizedChi2();
754  double trackerHits2 = tk2->hitPattern().numberOfValidTrackerHits();
755  int pixelHits2 = tk2->hitPattern().numberOfValidPixelHits();
756  int muonHits2 = gm2->hitPattern().numberOfValidMuonHits();
757  int nMatches2 = mu2.numberOfMatches();
758  bool quality2 = true;
759  if (normalizedChi22 > normalizedChi2Cut_) quality2 = false;
760  if (trackerHits2 < trackerHitsCut_) quality2 = false;
761  if (pixelHits2 < pixelHitsCut_) quality2 = false;
762  if (muonHits2 < muonHitsCut_) quality2 = false;
763  if (!mu2.isTrackerMuon()) quality2 = false;
764  if (nMatches2 < nMatchesCut_) quality2 = false;
765  zmuon_sel[8] = quality2;
766  double isovar2 = mu2.isolationR03().sumPt;
767  if (isCombinedIso_) {
768  isovar2 += mu2.isolationR03().emEt;
769  isovar2 += mu2.isolationR03().hadEt;
770  }
771  if (isRelativeIso_) isovar2 /= pt2;
772  if (isovar2 < isoCut03_) zmuon_sel[9] = true;
773  if (trigger_fired) zmuon_sel[10] = true;
774  const math::XYZTLorentzVector ZRecoGlb(
775  mu.px() + mu2.px(), mu.py() + mu2.py(), mu.pz() + mu2.pz(),
776  mu.p() + mu2.p());
777  if (ZRecoGlb.mass() > dimuonMassMin_ &&
778  ZRecoGlb.mass() < dimuonMassMax_)
779  zmuon_sel[11] = true;
780 
781  // jet flag
782  if (njets <= nJetMax_) zmuon_sel[12] = true;
783 
784  // start filling histos: N-1 plots
785  int flags_passed_z = 0;
786 
787  for (int jj = 0; jj < NFLAGSZ; ++jj) {
788  if (zmuon_sel[jj]) ++flags_passed_z;
789  }
790 
791  if (flags_passed_z >= (NFLAGSZ - 1)) {
792  if (!zmuon_sel[0] || flags_passed_z == NFLAGSZ) {
793  pt1_afterZ_->Fill(pt);
794  }
795  if (!zmuon_sel[1] || flags_passed_z == NFLAGSZ) {
796  eta1_afterZ_->Fill(eta);
797  }
798  if (!zmuon_sel[2] || flags_passed_z == NFLAGSZ) {
799  dxy1_afterZ_->Fill(dxy);
800  }
801  if (!zmuon_sel[3] || flags_passed_z == NFLAGSZ) {
802  goodewkmuon1_afterZ_->Fill(quality);
803  }
804  if (!zmuon_sel[4] || flags_passed_z == NFLAGSZ) {
805  iso1_afterZ_->Fill(isovar);
806  }
807  if (!zmuon_sel[5] || flags_passed_z == NFLAGSZ) {
808  pt2_afterZ_->Fill(pt2);
809  }
810  if (!zmuon_sel[6] || flags_passed_z == NFLAGSZ) {
811  eta2_afterZ_->Fill(eta2);
812  }
813  if (!zmuon_sel[7] || flags_passed_z == NFLAGSZ) {
814  dxy2_afterZ_->Fill(dxy2);
815  }
816  if (!zmuon_sel[8] || flags_passed_z == NFLAGSZ) {
817  goodewkmuon2_afterZ_->Fill(quality2);
818  }
819  if (!zmuon_sel[9] || flags_passed_z == NFLAGSZ) {
820  iso2_afterZ_->Fill(isovar2);
821  }
822  if (!zmuon_sel[10] || flags_passed_z == NFLAGSZ) {
823  if (!zhlt_hist_done) ztrig_afterZ_->Fill(trigger_fired);
824  zhlt_hist_done = true;
825  }
826  if (!zmuon_sel[11] || flags_passed_z == NFLAGSZ) {
827  dimuonmass_afterZ_->Fill(ZRecoGlb.mass());
828  }
829  if (!zmuon_sel[12] || flags_passed_z == NFLAGSZ) {
830  if (!zjets_hist_done) {
831  njets_afterZ_->Fill(njets);
832  leadingjet_pt_afterZ_->Fill(lead_jet_pt);
833  leadingjet_eta_afterZ_->Fill(lead_jet_eta);
834  }
835  zjets_hist_done = true;
836  }
837  if (flags_passed_z == NFLAGSZ) {
838  met_afterZ_->Fill(met_et);
839  if (!zfullsel_hist_done) {
840  npvs_afterZ_->Fill(nvvertex);
841  muoncharge_afterZ_->Fill(charge);
842  if (charge > 0) {
843  // ptPlus_afterZ_->Fill(mu.pt());
844  // ptMinus_afterZ_->Fill(mu2.pt());
845  ptDiffPM_afterZ_->Fill(mu.pt() - mu2.pt());
846  } else {
847  // ptPlus_afterZ_->Fill(mu2.pt());
848  // ptMinus_afterZ_->Fill(mu.pt());
849  ptDiffPM_afterZ_->Fill(mu2.pt() - mu.pt());
850  }
851  }
852  zfullsel_hist_done = true;
853  }
854  }
855  }
856  }
857  }
858 
859  if (zfullsel_hist_done) {
860  // here was a Z candidate
861  n_zselPt1thr_->Fill(nmuonsForZ1);
862  n_zselPt2thr_->Fill(nmuonsForZ2);
863  }
864 
865  // nmuons_->Fill(number_of_muons);
866  // nmuons_->Fill(muonCollectionSize);
867  ngoodmuons_->Fill(number_of_goodMuons);
868 
869  return;
870 }
871 
872 // Local Variables:
873 // show-trailing-whitespace: t
874 // truncate-lines: t
875 // End:
EwkMuDQM::met_before_
MonitorElement * met_before_
Definition: EwkMuDQM.h:112
reco::MuonIsolation::hadEt
float hadEt
hcal sum-Et
Definition: MuonIsolation.h:9
i
int i
Definition: DBlmapReader.cc:9
reco::LeafCandidate::p
virtual double p() const
magnitude of momentum vector
Definition: LeafCandidate.h:157
EwkMuDQM::ptDiffPM_afterZ_
MonitorElement * ptDiffPM_afterZ_
Definition: EwkMuDQM.h:134
EwkMuDQM::etaCut_
double etaCut_
Definition: EwkMuDQM.h:73
edm::Event::triggerNames
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
Definition: Event.cc:199
EwkMuDQM::dimuonmass_afterZ_
MonitorElement * dimuonmass_afterZ_
Definition: EwkMuDQM.h:154
EwkMuDQM::nall
unsigned int nall
Definition: EwkMuDQM.h:90
EwkMuDQM::nMatchesCut_
int nMatchesCut_
Definition: EwkMuDQM.h:60
EwkMuDQM::phEta_
MonitorElement * phEta_
Definition: EwkMuDQM.h:172
EwkMuDQM::EwkMuDQM
EwkMuDQM(const edm::ParameterSet &)
Definition: EwkMuDQM.cc:37
reco::MuonIsolation::sumPt
float sumPt
sum-pt of tracks
Definition: MuonIsolation.h:7
EwkMuDQM::njets_afterZ_
MonitorElement * njets_afterZ_
Definition: EwkMuDQM.h:119
reco::LeafCandidate::pt
virtual float pt() const
transverse momentum
Definition: LeafCandidate.h:178
EwkMuDQM::nrec
unsigned int nrec
Definition: EwkMuDQM.h:91
EcalCondDBWriter_cfi.cosmic
string cosmic
Definition: EcalCondDBWriter_cfi.py:43
EwkMuDQM::iso2_afterZ_
MonitorElement * iso2_afterZ_
Definition: EwkMuDQM.h:146
reco::LeafCandidate::et
virtual double et() const
transverse energy
Definition: LeafCandidate.h:161
EwkMuDQM::mtMax_
double mtMax_
Definition: EwkMuDQM.h:70
EwkMuDQM::npvs_before_
MonitorElement * npvs_before_
Definition: EwkMuDQM.h:156
EwkMuDQM::ptThrForPhoton_
double ptThrForPhoton_
Definition: EwkMuDQM.h:84
reco::Muon::innerTrack
virtual TrackRef innerTrack() const
Definition: Muon.h:48
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:446
EwkMuDQM::dxy1_afterZ_
MonitorElement * dxy1_afterZ_
Definition: EwkMuDQM.h:139
EwkMuDQM::njets_before_
MonitorElement * njets_before_
Definition: EwkMuDQM.h:117
EwkMuDQM::dxy_before_
MonitorElement * dxy_before_
Definition: EwkMuDQM.h:102
EwkMuDQM::ztrig_afterZ_
MonitorElement * ztrig_afterZ_
Definition: EwkMuDQM.h:152
EwkMuDQM::trig_before_
MonitorElement * trig_before_
Definition: EwkMuDQM.h:108
EwkMuDQM::trig_after_
MonitorElement * trig_after_
Definition: EwkMuDQM.h:109
EwkMuDQM::leadingjet_eta_afterZ_
MonitorElement * leadingjet_eta_afterZ_
Definition: EwkMuDQM.h:125
EwkMuDQM::metMax_
double metMax_
Definition: EwkMuDQM.h:68
reco::LeafCandidate::phi
virtual float phi() const
momentum azimuthal angle
Definition: LeafCandidate.h:180
reco::Muon::isTrackerMuon
bool isTrackerMuon() const
Definition: Muon.h:219
EwkMuDQM::muoncharge_before_
MonitorElement * muoncharge_before_
Definition: EwkMuDQM.h:160
reco::Jet
Base class for all types of Jets.
Definition: Jet.h:20
reco::Vertex::isValid
bool isValid() const
Tells whether the vertex is valid.
Definition: Vertex.h:60
EwkMuDQM::beamSpotTag_
edm::EDGetTokenT< reco::BeamSpot > beamSpotTag_
Definition: EwkMuDQM.h:51
EwkMuDQM::eta2_afterZ_
MonitorElement * eta2_afterZ_
Definition: EwkMuDQM.h:143
EwkMuDQM::trigTag_
edm::EDGetTokenT< edm::TriggerResults > trigTag_
Definition: EwkMuDQM.h:45
EwkMuDQM::ptCut_
double ptCut_
Definition: EwkMuDQM.h:72
EwkMuDQM::isoCut03_
double isoCut03_
Definition: EwkMuDQM.h:64
EwkMuDQM::npvs_afterZ_
MonitorElement * npvs_afterZ_
Definition: EwkMuDQM.h:158
reco::Muon::isGlobalMuon
bool isGlobalMuon() const
Definition: Muon.h:218
EwkMuDQM::leadingjet_eta_before_
MonitorElement * leadingjet_eta_before_
Definition: EwkMuDQM.h:123
ev
bool ev
Definition: Hydjet2Hadronizer.cc:90
EwkMuDQM::trigPathNames_
std::vector< std::string > trigPathNames_
Definition: EwkMuDQM.h:52
EwkMuDQM::eJetMin_
double eJetMin_
Definition: EwkMuDQM.h:81
eta
T eta() const
Definition: Basic3DVectorLD.h:177
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:107
EwkMuDQM::acop_before_
MonitorElement * acop_before_
Definition: EwkMuDQM.h:114
EwkMuDQM::leadingjet_pt_after_
MonitorElement * leadingjet_pt_after_
Definition: EwkMuDQM.h:121
GoodVertex_cfg.vertexCollection
tuple vertexCollection
Definition: GoodVertex_cfg.py:61
EwkMuDQM::goodewkmuon_before_
MonitorElement * goodewkmuon_before_
Definition: EwkMuDQM.h:104
EwkMuDQM::eta_after_
MonitorElement * eta_after_
Definition: EwkMuDQM.h:101
EwkMuDQM::muoncharge_afterZ_
MonitorElement * muoncharge_afterZ_
Definition: EwkMuDQM.h:162
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
EwkMuDQM::pt_before_
MonitorElement * pt_before_
Definition: EwkMuDQM.h:98
EwkMuDQM::nhlt
unsigned int nhlt
Definition: EwkMuDQM.h:93
EwkMuDQM::dxy2_afterZ_
MonitorElement * dxy2_afterZ_
Definition: EwkMuDQM.h:144
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:164
EwkMuDQM::leadingjet_pt_before_
MonitorElement * leadingjet_pt_before_
Definition: EwkMuDQM.h:120
EwkMuDQM::goodewkmuon1_afterZ_
MonitorElement * goodewkmuon1_afterZ_
Definition: EwkMuDQM.h:140
EwkMuDQM::iso_before_
MonitorElement * iso_before_
Definition: EwkMuDQM.h:106
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
EwkMuDQM::npvs_after_
MonitorElement * npvs_after_
Definition: EwkMuDQM.h:157
EwkMuDQM::goodewkmuon2_afterZ_
MonitorElement * goodewkmuon2_afterZ_
Definition: EwkMuDQM.h:145
HLTConfigProvider::prescaleValue
unsigned int prescaleValue(unsigned int set, const std::string &trigger) const
HLT prescale value in specific prescale set for a specific trigger path.
Definition: HLTConfigProvider.h:232
EwkMuDQM::njets_after_
MonitorElement * njets_after_
Definition: EwkMuDQM.h:118
EwkMuDQM::endRun
virtual void endRun(const edm::Run &, const edm::EventSetup &)
Definition: EwkMuDQM.cc:311
EwkMuDQM::ngoodmuons_
MonitorElement * ngoodmuons_
Definition: EwkMuDQM.h:165
EwkMuDQM::n_zselPt2thr_
MonitorElement * n_zselPt2thr_
Definition: EwkMuDQM.h:150
EwkMuDQM::metMin_
double metMin_
Definition: EwkMuDQM.h:67
EwkMuDQM::pt2_afterZ_
MonitorElement * pt2_afterZ_
Definition: EwkMuDQM.h:142
EwkMuDQM::ptThrForZ2_
double ptThrForZ2_
Definition: EwkMuDQM.h:76
EwkMuDQM::nJetMax_
int nJetMax_
Definition: EwkMuDQM.h:82
reco::MET
Definition: MET.h:42
reco::LeafCandidate::eta
virtual float eta() const
momentum pseudorapidity
Definition: LeafCandidate.h:184
EwkMuDQM::hltConfigProvider_
HLTConfigProvider hltConfigProvider_
Definition: EwkMuDQM.h:88
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
EwkMuDQM::nsel
unsigned int nsel
Definition: EwkMuDQM.h:95
reco::MuonIsolation::emEt
float emEt
ecal sum-Et
Definition: MuonIsolation.h:8
EwkMuDQM::met_afterZ_
MonitorElement * met_afterZ_
Definition: EwkMuDQM.h:136
EwkMuDQM::dxyCut_
double dxyCut_
Definition: EwkMuDQM.h:55
reco::LeafCandidate::charge
virtual int charge() const
electric charge
Definition: LeafCandidate.h:140
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
EwkMuDQM::leadingjet_eta_after_
MonitorElement * leadingjet_eta_after_
Definition: EwkMuDQM.h:124
j
int j
Definition: DBlmapReader.cc:9
EwkMuDQM::dimuonMassMin_
double dimuonMassMin_
Definition: EwkMuDQM.h:78
RPCpg::mu
const int mu
Definition: Constants.h:22
EwkMuDQM::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: EwkMuDQM.cc:132
EwkMuDQM::isCombinedIso_
bool isCombinedIso_
Definition: EwkMuDQM.h:63
EwkMuDQM::dimuonMassMax_
double dimuonMassMax_
Definition: EwkMuDQM.h:79
triggerResults
static std::string const triggerResults
Definition: EdmProvDump.cc:41
EwkMuDQM::mt_after_
MonitorElement * mt_after_
Definition: EwkMuDQM.h:111
EwkMuDQM::goodewkmuon_after_
MonitorElement * goodewkmuon_after_
Definition: EwkMuDQM.h:105
EwkMuDQM::muoncharge_after_
MonitorElement * muoncharge_after_
Definition: EwkMuDQM.h:161
EwkMuDQM::mt_before_
MonitorElement * mt_before_
Definition: EwkMuDQM.h:110
edm::Ref::isNull
bool isNull() const
Checks for null.
Definition: Ref.h:247
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
M_PI
#define M_PI
Definition: readTestVector.cc:25
EwkMuDQM::eta_before_
MonitorElement * eta_before_
Definition: EwkMuDQM.h:100
EwkMuDQM::iso1_afterZ_
MonitorElement * iso1_afterZ_
Definition: EwkMuDQM.h:141
EwkMuDQM::isRelativeIso_
bool isRelativeIso_
Definition: EwkMuDQM.h:62
Geom::Phi::value
T value() const
Explicit access to value in case implicit conversion not OK.
Definition: Phi.h:38
trigNames
static const char *const trigNames[]
Definition: EcalDumpRaw.cc:74
EwkMuDQM::jetTag_
edm::InputTag jetTag_
Definition: EwkMuDQM.h:44
EwkMuDQM::pixelHitsCut_
int pixelHitsCut_
Definition: EwkMuDQM.h:58
EwkMuDQM::vertexTag_
edm::EDGetTokenT< edm::View< reco::Vertex > > vertexTag_
Definition: EwkMuDQM.h:50
EwkMuDQM::nmuons_
MonitorElement * nmuons_
Definition: EwkMuDQM.h:164
EwkMuDQM::normalizedChi2Cut_
double normalizedChi2Cut_
Definition: EwkMuDQM.h:56
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:274
EwkMuDQM::metTag_
edm::InputTag metTag_
Definition: EwkMuDQM.h:43
HLTConfigProvider::prescaleSet
int prescaleSet(const edm::Event &iEvent, const edm::EventSetup &iSetup) const
Definition: HLTConfigProvider.cc:252
EwkMuDQM::pt_after_
MonitorElement * pt_after_
Definition: EwkMuDQM.h:99
EwkMuDQM::iso_after_
MonitorElement * iso_after_
Definition: EwkMuDQM.h:107
reco::LeafCandidate::px
virtual double px() const
x coordinate of momentum vector
Definition: LeafCandidate.h:172
EwkMuDQM::dqmBeginRun
virtual void dqmBeginRun(const edm::Run &, const edm::EventSetup &)
Definition: EwkMuDQM.cc:115
edm::TriggerNames::triggerName
std::string const & triggerName(unsigned int index) const
Definition: TriggerNames.cc:27
HLTConfigProvider::init
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d&#39;tor
Definition: HLTConfigProvider.cc:42
EwkMuDQM::eta1_afterZ_
MonitorElement * eta1_afterZ_
Definition: EwkMuDQM.h:138
EwkMuDQM::ptDiffPM_before_
MonitorElement * ptDiffPM_before_
Definition: EwkMuDQM.h:129
reco::LeafCandidate::pz
virtual double pz() const
z coordinate of momentum vector
Definition: LeafCandidate.h:176
reco::Muon::numberOfMatches
int numberOfMatches(ArbitrationType type=SegmentAndTrackArbitration) const
get number of chambers with matched segments
EwkMuDQM::pt1_afterZ_
MonitorElement * pt1_afterZ_
Definition: EwkMuDQM.h:137
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:42
EwkMuDQM::nph_
MonitorElement * nph_
Definition: EwkMuDQM.h:168
EwkMuDQM::trackerHitsCut_
int trackerHitsCut_
Definition: EwkMuDQM.h:57
EwkMuDQM::mtMin_
double mtMin_
Definition: EwkMuDQM.h:69
EwkMuDQM::nz
unsigned int nz
Definition: EwkMuDQM.h:96
EwkMuDQM::phPt_
MonitorElement * phPt_
Definition: EwkMuDQM.h:170
EwkMuDQM::ptThrForZ1_
double ptThrForZ1_
Definition: EwkMuDQM.h:75
EwkMuDQM::jetToken_
edm::EDGetTokenT< edm::View< reco::Jet > > jetToken_
Definition: EwkMuDQM.h:48
noPileUpMEtUtilities::leadJet
reco::MVAMEtJetInfo leadJet(const reco::MVAMEtJetInfoCollection &)
Definition: noPileUpMEtUtilities.cc:111
EwkMuDQM::dxy_after_
MonitorElement * dxy_after_
Definition: EwkMuDQM.h:103
HLTConfigProvider::prescaleSize
unsigned int prescaleSize() const
Definition: HLTConfigProvider.h:228
EwkMuDQM::muonHitsCut_
int muonHitsCut_
Definition: EwkMuDQM.h:59
caloMETBenchmarkGeneric_cfi.met
string met
Definition: caloMETBenchmarkGeneric_cfi.py:3
EwkMuDQM::met_after_
MonitorElement * met_after_
Definition: EwkMuDQM.h:113
edm::false
volatile std::atomic< bool > shutdown_flag false
Definition: UnixSignalHandlers.cc:22
EwkMuDQM::muonTag_
edm::EDGetTokenT< edm::View< reco::Muon > > muonTag_
Definition: EwkMuDQM.h:46
EwkMuDQM::n_zselPt1thr_
MonitorElement * n_zselPt1thr_
Definition: EwkMuDQM.h:149
EwkMuDQM::niso
unsigned int niso
Definition: EwkMuDQM.h:92
EwkMuDQM::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: EwkMuDQM.cc:313
EwkMuDQM::metToken_
edm::EDGetTokenT< edm::View< reco::MET > > metToken_
Definition: EwkMuDQM.h:47
EwkMuDQM::phoTag_
edm::EDGetTokenT< edm::View< reco::Photon > > phoTag_
Definition: EwkMuDQM.h:49
reco::LeafCandidate::py
virtual double py() const
y coordinate of momentum vector
Definition: LeafCandidate.h:174
EwkMuDQM::dimuonmass_before_
MonitorElement * dimuonmass_before_
Definition: EwkMuDQM.h:153
reco::Muon::isolationR03
const MuonIsolation & isolationR03() const
Definition: Muon.h:158
edm::Run
Definition: Run.h:41
EwkMuDQM::acopCut_
double acopCut_
Definition: EwkMuDQM.h:66
EwkMuDQM::nmet
unsigned int nmet
Definition: EwkMuDQM.h:94
EwkMuDQM::isValidHltConfig_
bool isValidHltConfig_
Definition: EwkMuDQM.h:87
reco::Muon::globalTrack
virtual TrackRef globalTrack() const
reference to Track reconstructed in both tracked and muon detector
Definition: Muon.h:54
EwkMuDQM::leadingjet_pt_afterZ_
MonitorElement * leadingjet_pt_afterZ_
Definition: EwkMuDQM.h:122
EwkMuDQM::acop_after_
MonitorElement * acop_after_
Definition: EwkMuDQM.h:115
Generated for CMSSW Reference Manual by doxygen 1.8.5