CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
METplusTrackMonitor.cc
Go to the documentation of this file.
1 #include "FWCore/Framework/interface/Frameworkfwd.h"
2 #include "FWCore/Framework/interface/MakerMacros.h"
3 #include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
4 #include "FWCore/MessageLogger/interface/MessageLogger.h"
5 #include "DQMServices/Core/interface/DQMStore.h"
6 #include "DQMServices/Core/interface/DQMEDAnalyzer.h"
7 #include "DQMOffline/Trigger/plugins/TriggerDQMBase.h"
8 #include "CommonTools/TriggerUtils/interface/GenericTriggerEventFlag.h"
9 #include "CommonTools/Utils/interface/StringCutObjectSelector.h"
10 #include "DataFormats/METReco/interface/CaloMET.h"
11 #include "DataFormats/METReco/interface/CaloMETCollection.h"
12 #include "DataFormats/MuonReco/interface/Muon.h"
13 #include "DataFormats/MuonReco/interface/MuonFwd.h"
14 #include "DataFormats/MuonReco/interface/MuonSelectors.h"
15 #include "DataFormats/JetReco/interface/PFJet.h"
16 #include "DataFormats/JetReco/interface/PFJetCollection.h"
17 #include "DataFormats/HLTReco/interface/TriggerObject.h"
18 #include "DataFormats/HLTReco/interface/TriggerEvent.h"
19 #include "DataFormats/VertexReco/interface/Vertex.h"
20 #include "DataFormats/VertexReco/interface/VertexFwd.h"
21 #include "DataFormats/Math/interface/deltaR.h"
22 
23 class METplusTrackMonitor : public DQMEDAnalyzer, public TriggerDQMBase {
24 public:
25  typedef dqm::reco::MonitorElement MonitorElement;
26  typedef dqm::reco::DQMStore DQMStore;
27 
28  METplusTrackMonitor(const edm::ParameterSet&);
29  ~METplusTrackMonitor() noexcept(true) override {}
30  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
31 
32 protected:
33  void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
34  void analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) override;
35 
36 private:
37  bool getHLTObj(const edm::Handle<trigger::TriggerEvent>& trigSummary,
38  const edm::InputTag& filterTag,
39  trigger::TriggerObject& obj) const;
40 
41  const std::string folderName_;
42 
43  const bool requireValidHLTPaths_;
44  bool hltPathsAreValid_;
45 
46  edm::EDGetTokenT<reco::CaloMETCollection> metToken_;
47  edm::EDGetTokenT<reco::MuonCollection> muonToken_;
48  edm::EDGetTokenT<reco::PFJetCollection> jetToken_;
49  edm::EDGetTokenT<reco::VertexCollection> vtxToken_;
50  edm::EDGetTokenT<trigger::TriggerEvent> theTrigSummary_;
51 
52  edm::InputTag hltMetTag_;
53  edm::InputTag hltMetCleanTag_;
54  edm::InputTag trackLegFilterTag_;
55 
56  std::vector<double> met_variable_binning_;
57  std::vector<double> muonPt_variable_binning_;
58 
59  MEbinning met_binning_;
60  MEbinning ls_binning_;
61  MEbinning pt_binning_;
62  MEbinning eta_binning_;
63  MEbinning phi_binning_;
64 
65  ObjME metME_variableBinning_;
66  ObjME metVsLS_;
67  ObjME metPhiME_;
68  ObjME deltaphimetj1ME_;
69  ObjME metVsHltMet_;
70  ObjME metVsHltMetClean_;
71 
72  ObjME muonPtME_variableBinning_;
73  ObjME muonPtVsLS_;
74  ObjME muonEtaME_;
75  ObjME deltaphimetmuonME_;
76  ObjME muonEtaVsPhi_;
77 
78  std::unique_ptr<GenericTriggerEventFlag> num_genTriggerEventFlag_;
79  std::unique_ptr<GenericTriggerEventFlag> den_genTriggerEventFlag_;
80 
81  StringCutObjectSelector<reco::CaloMET, true> metSelection_;
82  StringCutObjectSelector<reco::Muon, true> muonSelection_;
83  StringCutObjectSelector<reco::PFJet, true> jetSelection_;
84  StringCutObjectSelector<reco::Vertex, true> vtxSelection_;
85 
86  unsigned nmuons_;
87  unsigned njets_;
88 
89  double leadJetEtaCut_;
90 
91  bool requireLeadMatched_;
92  double maxMatchDeltaR_;
93 };
94 
95 METplusTrackMonitor::METplusTrackMonitor(const edm::ParameterSet& iConfig)
96  : folderName_(iConfig.getParameter<std::string>("FolderName")),
97  requireValidHLTPaths_(iConfig.getParameter<bool>("requireValidHLTPaths")),
98  hltPathsAreValid_(false),
99  metToken_(consumes<reco::CaloMETCollection>(iConfig.getParameter<edm::InputTag>("met"))),
100  muonToken_(consumes<reco::MuonCollection>(iConfig.getParameter<edm::InputTag>("muons"))),
101  jetToken_(consumes<reco::PFJetCollection>(iConfig.getParameter<edm::InputTag>("jets"))),
102  vtxToken_(consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("vertices"))),
103  theTrigSummary_(consumes<trigger::TriggerEvent>(iConfig.getParameter<edm::InputTag>("trigSummary"))),
104  hltMetTag_(iConfig.getParameter<edm::InputTag>("hltMetFilter")),
105  hltMetCleanTag_(iConfig.getParameter<edm::InputTag>("hltMetCleanFilter")),
106  trackLegFilterTag_(iConfig.getParameter<edm::InputTag>("trackLegFilter")),
107  met_variable_binning_(
108  iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("metBinning")),
109  muonPt_variable_binning_(
110  iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("ptBinning")),
111  met_binning_(getHistoPSet(
112  iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("metPSet"))),
113  ls_binning_(
114  getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("lsPSet"))),
115  pt_binning_(
116  getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("ptPSet"))),
117  eta_binning_(getHistoPSet(
118  iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("etaPSet"))),
119  phi_binning_(getHistoPSet(
120  iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("phiPSet"))),
121  num_genTriggerEventFlag_(new GenericTriggerEventFlag(
122  iConfig.getParameter<edm::ParameterSet>("numGenericTriggerEventPSet"), consumesCollector(), *this)),
123  den_genTriggerEventFlag_(new GenericTriggerEventFlag(
124  iConfig.getParameter<edm::ParameterSet>("denGenericTriggerEventPSet"), consumesCollector(), *this)),
125  metSelection_(iConfig.getParameter<std::string>("metSelection")),
126  muonSelection_(iConfig.getParameter<std::string>("muonSelection")),
127  jetSelection_(iConfig.getParameter<std::string>("jetSelection")),
128  vtxSelection_(iConfig.getParameter<std::string>("vtxSelection")),
129  nmuons_(iConfig.getParameter<unsigned>("nmuons")),
130  njets_(iConfig.getParameter<unsigned>("njets")),
131  leadJetEtaCut_(iConfig.getParameter<double>("leadJetEtaCut")),
132  requireLeadMatched_(iConfig.getParameter<bool>("requireLeadMatched")),
133  maxMatchDeltaR_(iConfig.getParameter<double>("maxMatchDeltaR")) {}
134 
135 void METplusTrackMonitor::bookHistograms(DQMStore::IBooker& ibooker,
136  edm::Run const& iRun,
137  edm::EventSetup const& iSetup) {
138  // Initialize the GenericTriggerEventFlag
139  if (num_genTriggerEventFlag_ && num_genTriggerEventFlag_->on()) {
140  num_genTriggerEventFlag_->initRun(iRun, iSetup);
141  }
142  if (den_genTriggerEventFlag_ && den_genTriggerEventFlag_->on()) {
143  den_genTriggerEventFlag_->initRun(iRun, iSetup);
144  }
145 
146  // check if every HLT path specified in numerator and denominator has a valid match in the HLT Menu
147  hltPathsAreValid_ = (num_genTriggerEventFlag_ && den_genTriggerEventFlag_ && num_genTriggerEventFlag_->on() &&
148  den_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->allHLTPathsAreValid() &&
149  den_genTriggerEventFlag_->allHLTPathsAreValid());
150 
151  // if valid HLT paths are required,
152  // create DQM outputs only if all paths are valid
153  if (requireValidHLTPaths_ and (not hltPathsAreValid_)) {
154  return;
155  }
156 
157  std::string histname, histtitle;
158 
159  std::string currentFolder = folderName_;
160  ibooker.setCurrentFolder(currentFolder);
161 
162  // MET leg histograms
163  histname = "met_variable";
164  histtitle = "CaloMET";
165  bookME(ibooker, metME_variableBinning_, histname, histtitle, met_variable_binning_);
166  setMETitle(metME_variableBinning_, "CaloMET [GeV]", "events / [GeV]");
167 
168  histname = "metVsLS";
169  histtitle = "CaloMET vs LS";
170  bookME(ibooker,
171  metVsLS_,
172  histname,
173  histtitle,
174  ls_binning_.nbins,
175  ls_binning_.xmin,
176  ls_binning_.xmax,
177  met_binning_.xmin,
178  met_binning_.xmax);
179  setMETitle(metVsLS_, "LS", "CaloMET [GeV]");
180 
181  histname = "metPhi";
182  histtitle = "CaloMET phi";
183  bookME(ibooker, metPhiME_, histname, histtitle, phi_binning_.nbins, phi_binning_.xmin, phi_binning_.xmax);
184  setMETitle(metPhiME_, "CaloMET #phi", "events / 0.2 rad");
185 
186  histname = "deltaphi_metjet1";
187  histtitle = "dPhi(CaloMET, jet1)";
188  bookME(ibooker, deltaphimetj1ME_, histname, histtitle, phi_binning_.nbins, 0, phi_binning_.xmax);
189  setMETitle(deltaphimetj1ME_, "#Delta#phi (CaloMET, j1)", "events / 0.1 rad");
190 
191  histname = "metVsHltMet";
192  histtitle = "CaloMET vs hltMet";
193  bookME(ibooker,
194  metVsHltMet_,
195  histname,
196  histtitle,
197  met_binning_.nbins,
198  met_binning_.xmin,
199  met_binning_.xmax,
200  met_binning_.nbins,
201  met_binning_.xmin,
202  met_binning_.xmax);
203  setMETitle(metVsHltMet_, "hltMet (online) [GeV]", "CaloMET (offline) [GeV]");
204 
205  histname = "metVsHltMetClean";
206  histtitle = "CaloMET vs hltMetClean";
207  bookME(ibooker,
208  metVsHltMetClean_,
209  histname,
210  histtitle,
211  met_binning_.nbins,
212  met_binning_.xmin,
213  met_binning_.xmax,
214  met_binning_.nbins,
215  met_binning_.xmin,
216  met_binning_.xmax);
217  setMETitle(metVsHltMetClean_, "hltMetClean (online) [GeV]", "CaloMET (offline) [GeV]");
218 
219  // Track leg histograms
220 
221  histname = "muonPt_variable";
222  histtitle = "Muon PT";
223  bookME(ibooker, muonPtME_variableBinning_, histname, histtitle, muonPt_variable_binning_);
224  setMETitle(muonPtME_variableBinning_, "Muon p_{T} [GeV]", "events / [GeV]");
225 
226  histname = "muonEta";
227  histtitle = "Muon eta";
228  bookME(ibooker, muonEtaME_, histname, histtitle, eta_binning_.nbins, eta_binning_.xmin, eta_binning_.xmax);
229  setMETitle(muonEtaME_, "Muon #eta", "events / 0.2");
230 
231  histname = "deltaphi_muonmet";
232  histtitle = "dPhi(Muon, CaloMET)";
233  bookME(ibooker, deltaphimetmuonME_, histname, histtitle, phi_binning_.nbins, 0, phi_binning_.xmax);
234  setMETitle(deltaphimetmuonME_, "#Delta#phi (Muon, CaloMET)", "events / 0.1 rad");
235 
236  histname = "muonEtaVsPhi";
237  histtitle = "Muon eta vs phi";
238  bookME(ibooker,
239  muonEtaVsPhi_,
240  histname,
241  histtitle,
242  phi_binning_.nbins,
243  phi_binning_.xmin,
244  phi_binning_.xmax,
245  eta_binning_.nbins,
246  eta_binning_.xmin,
247  eta_binning_.xmax);
248  setMETitle(muonEtaVsPhi_, "Muon #phi", "Muon #eta");
249 
250  histname = "muonPtVsLS";
251  histtitle = "Muon PT vs LS";
252  bookME(ibooker,
253  muonPtVsLS_,
254  histname,
255  histtitle,
256  ls_binning_.nbins,
257  ls_binning_.xmin,
258  ls_binning_.xmax,
259  pt_binning_.xmin,
260  pt_binning_.xmax);
261  setMETitle(muonPtVsLS_, "LS", "Muon p_{T} [GeV]");
262 }
263 
264 void METplusTrackMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) {
265  // if valid HLT paths are required,
266  // analyze event only if all paths are valid
267  if (requireValidHLTPaths_ and (not hltPathsAreValid_)) {
268  return;
269  }
270 
271  // Filter out events if Trigger Filtering is requested
272  if (den_genTriggerEventFlag_->on() && !den_genTriggerEventFlag_->accept(iEvent, iSetup))
273  return;
274 
275  edm::Handle<reco::CaloMETCollection> metHandle;
276  iEvent.getByToken(metToken_, metHandle);
277  reco::CaloMET caloMet = metHandle->front();
278  if (!metSelection_(caloMet))
279  return;
280 
281  float met = caloMet.pt();
282  float metPhi = caloMet.phi();
283 
284  edm::Handle<reco::PFJetCollection> jetsHandle;
285  iEvent.getByToken(jetToken_, jetsHandle);
286  if (jetsHandle->size() < njets_)
287  return;
288  std::vector<reco::PFJet> jets;
289  for (auto const& j : *jetsHandle) {
290  if (jetSelection_(j))
291  jets.push_back(j);
292  }
293  if (jets.size() < njets_)
294  return;
295  if (njets_ > 0 && !(jets.empty()) && fabs(jets[0].eta()) > leadJetEtaCut_)
296  return;
297  float deltaphi_metjet1 = !(jets.empty()) ? fabs(deltaPhi(caloMet.phi(), jets[0].phi())) : 10.0;
298 
299  edm::Handle<reco::VertexCollection> primaryVertices;
300  iEvent.getByToken(vtxToken_, primaryVertices);
301  if (primaryVertices->empty())
302  return;
303  const reco::Vertex* pv = nullptr;
304  for (auto const& v : *primaryVertices) {
305  if (!vtxSelection_(v))
306  continue;
307  pv = &v;
308  break;
309  }
310  if (pv == nullptr)
311  return;
312 
313  edm::Handle<reco::MuonCollection> muonHandle;
314  iEvent.getByToken(muonToken_, muonHandle);
315  if (muonHandle->size() < nmuons_)
316  return;
317  std::vector<reco::Muon> muons;
318  for (auto const& m : *muonHandle) {
319  bool passTightID =
320  muon::isTightMuon(m, *pv) &&
321  m.innerTrack()->hitPattern().trackerLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS) == 0 &&
322  m.innerTrack()->hitPattern().trackerLayersWithoutMeasurement(reco::HitPattern::TRACK_HITS) == 0;
323  if (muonSelection_(m) && passTightID)
324  muons.push_back(m);
325  }
326  if (muons.size() < nmuons_)
327  return;
328 
329  // Filling MET leg histograms (denominator)
330  metME_variableBinning_.denominator->Fill(met);
331  metPhiME_.denominator->Fill(metPhi);
332  deltaphimetj1ME_.denominator->Fill(deltaphi_metjet1);
333 
334  int ls = iEvent.id().luminosityBlock();
335  metVsLS_.denominator->Fill(ls, met);
336 
337  // Apply the selection for the MET leg numerator
338  edm::Handle<trigger::TriggerEvent> triggerSummary;
339  iEvent.getByToken(theTrigSummary_, triggerSummary);
340  if (!triggerSummary.isValid()) {
341  edm::LogError("METplusTrackMonitor") << "Invalid collection: TriggerSummary"
342  << "\n";
343  return;
344  }
345 
346  trigger::TriggerObject hltMet, hltMetClean;
347  bool passesHltMetFilter = getHLTObj(triggerSummary, hltMetTag_, hltMet);
348  bool passesHltMetCleanFilter = getHLTObj(triggerSummary, hltMetCleanTag_, hltMetClean);
349 
350  if (!passesHltMetFilter || !passesHltMetCleanFilter)
351  return;
352 
353  // Filling MET leg histograms (numerator)
354  metME_variableBinning_.numerator->Fill(met);
355  metPhiME_.numerator->Fill(metPhi);
356  deltaphimetj1ME_.numerator->Fill(deltaphi_metjet1);
357  metVsLS_.numerator->Fill(ls, met);
358  metVsHltMet_.numerator->Fill(hltMet.pt(), met);
359  metVsHltMetClean_.numerator->Fill(hltMetClean.pt(), met);
360 
361  // Filling track leg histograms (denominator)
362  double leadMuonPt = !(muons.empty()) ? muons[0].pt() : -1.0;
363  double leadMuonEta = !(muons.empty()) ? muons[0].eta() : 10.0;
364  double leadMuonPhi = !(muons.empty()) ? muons[0].phi() : 10.0;
365  float deltaphi_metmuon = !(muons.empty()) ? fabs(deltaPhi(caloMet.phi(), muons[0].phi())) : 10.0;
366 
367  muonPtME_variableBinning_.denominator->Fill(leadMuonPt);
368  muonPtVsLS_.denominator->Fill(ls, leadMuonPt);
369  muonEtaME_.denominator->Fill(leadMuonEta);
370  deltaphimetmuonME_.denominator->Fill(deltaphi_metmuon);
371  muonEtaVsPhi_.denominator->Fill(leadMuonPhi, leadMuonEta);
372 
373  // Apply the selection for the track leg numerator
374  trigger::TriggerObject isoTrk;
375  bool passesTrackLegFilter = getHLTObj(triggerSummary, trackLegFilterTag_, isoTrk);
376 
377  // require track leg filter
378  if (!passesTrackLegFilter)
379  return;
380 
381  // if requested, require lead selected muon is matched to the track leg filter object
382  if (requireLeadMatched_ && !(muons.empty()) && deltaR(muons[0], isoTrk) < maxMatchDeltaR_)
383  return;
384 
385  // require the full HLT path is fired
386  if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup))
387  return;
388 
389  // Filling track leg histograms (denominator)
390  muonPtME_variableBinning_.numerator->Fill(leadMuonPt);
391  muonPtVsLS_.numerator->Fill(ls, leadMuonPt);
392  muonEtaME_.numerator->Fill(leadMuonEta);
393  deltaphimetmuonME_.numerator->Fill(deltaphi_metmuon);
394  muonEtaVsPhi_.numerator->Fill(leadMuonPhi, leadMuonEta);
395 }
396 
397 bool METplusTrackMonitor::getHLTObj(const edm::Handle<trigger::TriggerEvent>& trigSummary,
398  const edm::InputTag& filterTag,
399  trigger::TriggerObject& obj) const {
400  double leadingPt = -1.0;
401 
402  size_t filterIndex = trigSummary->filterIndex(filterTag);
403  trigger::TriggerObjectCollection triggerObjects = trigSummary->getObjects();
404 
405  if (!(filterIndex >= trigSummary->sizeFilters())) {
406  const trigger::Keys& keys = trigSummary->filterKeys(filterIndex);
407  for (unsigned short key : keys) {
408  trigger::TriggerObject foundObject = triggerObjects[key];
409  if (foundObject.pt() > leadingPt) {
410  obj = foundObject;
411  leadingPt = obj.pt();
412  }
413  }
414  }
415 
416  return (leadingPt > 0.0);
417 }
418 
419 void METplusTrackMonitor::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
420  edm::ParameterSetDescription desc;
421  desc.add<std::string>("FolderName", "HLT/MET");
422  desc.add<bool>("requireValidHLTPaths", true);
423 
424  desc.add<edm::InputTag>("met", edm::InputTag("caloMet"));
425  desc.add<edm::InputTag>("jets", edm::InputTag("ak4PFJetsCHS"));
426  desc.add<edm::InputTag>("electrons", edm::InputTag("gedGsfElectrons"));
427  desc.add<edm::InputTag>("muons", edm::InputTag("muons"));
428  desc.add<edm::InputTag>("vertices", edm::InputTag("offlinePrimaryVertices"));
429  desc.add<edm::InputTag>("trigSummary", edm::InputTag("hltTriggerSummaryAOD"));
430  desc.add<edm::InputTag>("hltMetFilter", edm::InputTag("hltMET105", "", "HLT"));
431  desc.add<edm::InputTag>("hltMetCleanFilter", edm::InputTag("hltMETClean65", "", "HLT"));
432  desc.add<edm::InputTag>("trackLegFilter", edm::InputTag("hltTrk50Filter", "", "HLT"));
433 
434  desc.add<std::string>("metSelection", "pt > 0");
435  desc.add<std::string>("jetSelection", "pt > 0");
436  desc.add<std::string>("muonSelection", "pt > 0");
437  desc.add<std::string>("vtxSelection", "!isFake");
438  desc.add<unsigned>("njets", 0);
439  desc.add<unsigned>("nmuons", 0);
440  desc.add<double>("leadJetEtaCut", 2.4);
441  desc.add<bool>("requireLeadMatched", true);
442  desc.add<double>("maxMatchDeltaR", 0.1);
443 
444  edm::ParameterSetDescription genericTriggerEventPSet;
445  genericTriggerEventPSet.add<bool>("andOr");
446  genericTriggerEventPSet.add<edm::InputTag>("dcsInputTag", edm::InputTag("scalersRawToDigi"));
447  genericTriggerEventPSet.add<std::vector<int> >("dcsPartitions", {});
448  genericTriggerEventPSet.add<bool>("andOrDcs", false);
449  genericTriggerEventPSet.add<bool>("errorReplyDcs", true);
450  genericTriggerEventPSet.add<std::string>("dbLabel", "");
451  genericTriggerEventPSet.add<bool>("andOrHlt", true);
452  genericTriggerEventPSet.add<edm::InputTag>("hltInputTag", edm::InputTag("TriggerResults::HLT"));
453  genericTriggerEventPSet.add<std::vector<std::string> >("hltPaths", {});
454  genericTriggerEventPSet.add<std::string>("hltDBKey", "");
455  genericTriggerEventPSet.add<bool>("errorReplyHlt", false);
456  genericTriggerEventPSet.add<unsigned int>("verbosityLevel", 1);
457 
458  desc.add<edm::ParameterSetDescription>("numGenericTriggerEventPSet", genericTriggerEventPSet);
459  desc.add<edm::ParameterSetDescription>("denGenericTriggerEventPSet", genericTriggerEventPSet);
460 
461  edm::ParameterSetDescription histoPSet;
462  edm::ParameterSetDescription metPSet;
463  edm::ParameterSetDescription phiPSet;
464  edm::ParameterSetDescription ptPSet;
465  edm::ParameterSetDescription etaPSet;
466 
467  fillHistoPSetDescription(metPSet);
468  fillHistoPSetDescription(phiPSet);
469  fillHistoPSetDescription(ptPSet);
470  fillHistoPSetDescription(etaPSet);
471 
472  histoPSet.add<edm::ParameterSetDescription>("metPSet", metPSet);
473  histoPSet.add<edm::ParameterSetDescription>("phiPSet", phiPSet);
474  histoPSet.add<edm::ParameterSetDescription>("ptPSet", ptPSet);
475  histoPSet.add<edm::ParameterSetDescription>("etaPSet", etaPSet);
476 
477  std::vector<double> bins = {0., 20., 40., 60., 80., 90., 100., 110., 120., 130., 140., 150., 160.,
478  170., 180., 190., 200., 220., 240., 260., 280., 300., 350., 400., 450., 1000.};
479 
480  histoPSet.add<std::vector<double> >("metBinning", bins);
481  histoPSet.add<std::vector<double> >("ptBinning", bins);
482 
483  edm::ParameterSetDescription lsPSet;
484  fillHistoLSPSetDescription(lsPSet);
485  histoPSet.add<edm::ParameterSetDescription>("lsPSet", lsPSet);
486 
487  desc.add<edm::ParameterSetDescription>("histoPSet", histoPSet);
488 
489  descriptions.add("metPlusTrackMonitoring", desc);
490 }
491 
492 DEFINE_FWK_MODULE(METplusTrackMonitor);
METplusTrackMonitor::jetToken_
edm::EDGetTokenT< reco::PFJetCollection > jetToken_
Definition: METplusTrackMonitor.cc:48
METplusTrackMonitor::hltMetCleanTag_
edm::InputTag hltMetCleanTag_
Definition: METplusTrackMonitor.cc:53
HLT_FULL_cff.deltaR
tuple deltaR
Definition: HLT_FULL_cff.py:7048
METplusTrackMonitor::METplusTrackMonitor
METplusTrackMonitor(const edm::ParameterSet &)
Definition: METplusTrackMonitor.cc:95
DiDispStaMuonMonitor_cfi.pt
tuple pt
Definition: DiDispStaMuonMonitor_cfi.py:39
TriggerDQMBase::fillHistoLSPSetDescription
static void fillHistoLSPSetDescription(edm::ParameterSetDescription &pset)
Definition: TriggerDQMBase.cc:114
reco::LeafCandidate::pt
double pt() const final
transverse momentum
Definition: LeafCandidate.h:146
METplusTrackMonitor::getHLTObj
bool getHLTObj(const edm::Handle< trigger::TriggerEvent > &trigSummary, const edm::InputTag &filterTag, trigger::TriggerObject &obj) const
Definition: METplusTrackMonitor.cc:397
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
CaloMET_cfi.caloMet
tuple caloMet
____________________________________________________________________________||
Definition: CaloMET_cfi.py:4
eta_binning_
eta_binning_(getHistoPSet(iConfig.getParameter< edm::ParameterSet >("histoPSet").getParameter< edm::ParameterSet >("etaPSet")))
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
HLT_FULL_cff.primaryVertices
tuple primaryVertices
Definition: HLT_FULL_cff.py:52244
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
den_genTriggerEventFlag_
den_genTriggerEventFlag_(new GenericTriggerEventFlag(iConfig.getParameter< edm::ParameterSet >("denGenericTriggerEventPSet"), consumesCollector(),*this))
METplusTrackMonitor::MonitorElement
dqm::reco::MonitorElement MonitorElement
Definition: METplusTrackMonitor.cc:25
eostools.ls
def ls
Definition: eostools.py:349
GenericTriggerEventFlag
Provides a code based selection for trigger and DCS information in order to have no failing filters i...
Definition: GenericTriggerEventFlag.h:42
METplusTrackMonitor::folderName_
const std::string folderName_
Definition: METplusTrackMonitor.cc:41
TriggerDQMBase::fillHistoPSetDescription
static void fillHistoPSetDescription(edm::ParameterSetDescription &pset)
Definition: TriggerDQMBase.cc:108
ls_binning_
ls_binning_(getHistoPSet(iConfig.getParameter< edm::ParameterSet >("histoPSet").getParameter< edm::ParameterSet >("lsPSet")))
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:111
l1t::VertexCollection
std::vector< Vertex > VertexCollection
Definition: Vertex.h:12
edm::EventID::luminosityBlock
LuminosityBlockNumber_t luminosityBlock() const
Definition: EventID.h:39
srCondWrite_cfg.deltaPhi
tuple deltaPhi
Definition: srCondWrite_cfg.py:54
METplusTrackMonitor::trackLegFilterTag_
edm::InputTag trackLegFilterTag_
Definition: METplusTrackMonitor.cc:54
TriggerDQMBase::setMETitle
void setMETitle(ObjME &me, const std::string &titleX, const std::string &titleY)
Definition: TriggerDQMBase.cc:3
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
METplusTrackMonitor::DQMStore
dqm::reco::DQMStore DQMStore
Definition: METplusTrackMonitor.cc:26
phi_binning_
phi_binning_(getHistoPSet(iConfig.getParameter< edm::ParameterSet >("histoPSet").getParameter< edm::ParameterSet >("phiPSet")))
reco::MuonCollection
std::vector< Muon > MuonCollection
collection of Muon objects
Definition: MuonFwd.h:9
METplusTrackMonitor::muonPt_variable_binning_
std::vector< double > muonPt_variable_binning_
Definition: METplusTrackMonitor.cc:57
METplusTrackMonitor::muonToken_
edm::EDGetTokenT< reco::MuonCollection > muonToken_
Definition: METplusTrackMonitor.cc:47
METplusTrackMonitor::analyze
void analyze(edm::Event const &iEvent, edm::EventSetup const &iSetup) override
Definition: METplusTrackMonitor.cc:264
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:290
CaloMETCollection
Collection of Calo MET.
trigger::TriggerObject
Single trigger physics object (e.g., an isolated muon)
Definition: TriggerObject.h:21
METplusTrackMonitor::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: METplusTrackMonitor.cc:419
pt_binning_
pt_binning_(getHistoPSet(iConfig.getParameter< edm::ParameterSet >("histoPSet").getParameter< edm::ParameterSet >("ptPSet")))
visualization-live-secondInstance_cfg.m
tuple m
Definition: visualization-live-secondInstance_cfg.py:79
iEvent
int iEvent
Definition: GenABIO.cc:224
METplusTrackMonitor::den_genTriggerEventFlag_
std::unique_ptr< GenericTriggerEventFlag > den_genTriggerEventFlag_
Definition: METplusTrackMonitor.cc:79
nmuons_
nmuons_(iConfig.getParameter< unsigned int >("nmuons"))
METplusTrackMonitor::met_variable_binning_
std::vector< double > met_variable_binning_
Definition: METplusTrackMonitor.cc:56
fwrapper::jets
vector< PseudoJet > jets
Definition: fastjetfortran_madfks.cc:44
TriggerDQMBase::ObjME::denominator
MonitorElement * denominator
Definition: TriggerDQMBase.h:28
TriggerDQMBase::ObjME::numerator
MonitorElement * numerator
Definition: TriggerDQMBase.h:27
submitPVResolutionJobs.key
tuple key
prepare the HTCondor submission files and eventually submit them
Definition: submitPVResolutionJobs.py:398
njets_
njets_(iConfig.getParameter< unsigned int >("njets"))
metSelection_
metSelection_(iConfig.getParameter< std::string >("metSelection"))
DQMOffline_LumiMontiroring_cff.lsPSet
tuple lsPSet
Definition: DQMOffline_LumiMontiroring_cff.py:14
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
METplusTrackMonitor::num_genTriggerEventFlag_
std::unique_ptr< GenericTriggerEventFlag > num_genTriggerEventFlag_
Definition: METplusTrackMonitor.cc:78
metToken_
metToken_(consumes< reco::PFMETCollection >(iConfig.getParameter< edm::InputTag >("met")))
met_variable_binning_
met_variable_binning_(iConfig.getParameter< edm::ParameterSet >("histoPSet").getParameter< std::vector< double > >("metBinning"))
METplusTrackMonitor::metSelection_
StringCutObjectSelector< reco::CaloMET, true > metSelection_
Definition: METplusTrackMonitor.cc:81
num_genTriggerEventFlag_
num_genTriggerEventFlag_(new GenericTriggerEventFlag(iConfig.getParameter< edm::ParameterSet >("numGenericTriggerEventPSet"), consumesCollector(),*this))
met_binning_
met_binning_(getHistoPSet(iConfig.getParameter< edm::ParameterSet >("histoPSet").getParameter< edm::ParameterSet >("metPSet")))
trigger::TriggerObjectCollection
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:75
METplusTrackMonitor::jetSelection_
StringCutObjectSelector< reco::PFJet, true > jetSelection_
Definition: METplusTrackMonitor.cc:83
vtxSelection_
vtxSelection_(iConfig.getParameter< std::string >("vertexSelection"))
METplusTrackMonitor::muonSelection_
StringCutObjectSelector< reco::Muon, true > muonSelection_
Definition: METplusTrackMonitor.cc:82
METplusTrackMonitor::vtxToken_
edm::EDGetTokenT< reco::VertexCollection > vtxToken_
Definition: METplusTrackMonitor.cc:49
trigger::Keys
std::vector< size_type > Keys
Definition: TriggerTypeDefs.h:19
METplusTrackMonitor::vtxSelection_
StringCutObjectSelector< reco::Vertex, true > vtxSelection_
Definition: METplusTrackMonitor.cc:84
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
METplusTrackMonitor::metToken_
edm::EDGetTokenT< reco::CaloMETCollection > metToken_
Definition: METplusTrackMonitor.cc:46
TriggerDQMBase::bookME
void bookME(DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const uint nbins, const double xmin, const double xmax, const bool bookDen=true)
Definition: TriggerDQMBase.cc:15
METplusTrackMonitor::~METplusTrackMonitor
~METplusTrackMonitor() noexcept(true) override
Definition: METplusTrackMonitor.cc:29
HLT_FULL_cff.InputTag
tuple InputTag
Definition: HLT_FULL_cff.py:86956
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:59
reco::PFJetCollection
std::vector< PFJet > PFJetCollection
collection of PFJet objects
Definition: PFJetCollection.h:14
patZpeak.muons
tuple muons
Definition: patZpeak.py:39
muon::isTightMuon
bool isTightMuon(const reco::Muon &, const reco::Vertex &)
Definition: MuonSelectors.cc:880
jetSelection_
jetSelection_(iConfig.getParameter< std::string >("jetSelection"))
reco::LeafCandidate::phi
double phi() const final
momentum azimuthal angle
Definition: LeafCandidate.h:148
METplusTrackMonitor::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: METplusTrackMonitor.cc:135
METplusTrackMonitor::theTrigSummary_
edm::EDGetTokenT< trigger::TriggerEvent > theTrigSummary_
Definition: METplusTrackMonitor.cc:50
edm::Run
Definition: Run.h:45
Generated for CMSSW Reference Manual by doxygen 1.8.5