diff --git a/DQMOffline/Trigger/plugins/DiJetMonitor.cc b/DQMOffline/Trigger/plugins/DiJetMonitor.cc index 78f050ef78fe9..b82d3ee7c4ef0 100644 --- a/DQMOffline/Trigger/plugins/DiJetMonitor.cc +++ b/DQMOffline/Trigger/plugins/DiJetMonitor.cc @@ -67,6 +67,11 @@ class DiJetMonitor : public DQMEDAnalyzer, public TriggerDQMBase { ObjME jetpt1ME_; ObjME jetpt2ME_; + ObjME jetPhi1ME_; + ObjME jetPhi2ME_; + ObjME jetEta1ME_; + ObjME jetEta2ME_; + ObjME jetphiTagME_; ObjME jetptAvgaME_; ObjME jetptAvgaThrME_; ObjME jetptAvgbME_; @@ -87,7 +92,7 @@ class DiJetMonitor : public DQMEDAnalyzer, public TriggerDQMBase { // Define Phi Bin // const double DiJet_MAX_PHI = 3.2; - //unsigned int DiJet_N_PHI = 64; + // unsigned int DiJet_N_PHI = 64; unsigned int DiJet_N_PHI = 32; MEbinning dijet_phi_binning{DiJet_N_PHI, -DiJet_MAX_PHI, DiJet_MAX_PHI}; // Define Eta Bin // @@ -161,6 +166,30 @@ void DiJetMonitor::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iR bookME(ibooker, jetpt2ME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax); setMETitle(jetpt2ME_, "Pt_2 [GeV]", "events"); + histname = "jetphi1"; + histtitle = "leading Jet Phi"; + bookME( + ibooker, jetPhi1ME_, histname, histtitle, dijet_phi_binning.nbins, dijet_phi_binning.xmin, dijet_phi_binning.xmax); + setMETitle(jetPhi1ME_, "Jet_Phi_1", "events"); + + histname = "jetphi2"; + histtitle = "second leading Jet Phi"; + bookME( + ibooker, jetPhi2ME_, histname, histtitle, dijet_phi_binning.nbins, dijet_phi_binning.xmin, dijet_phi_binning.xmax); + setMETitle(jetPhi2ME_, "Jet_Phi_2", "events"); + + histname = "jeteta1"; + histtitle = "leading Jet Eta"; + bookME( + ibooker, jetEta1ME_, histname, histtitle, dijet_eta_binning.nbins, dijet_eta_binning.xmin, dijet_eta_binning.xmax); + setMETitle(jetEta1ME_, "Jet_Eta_1", "events"); + + histname = "jeteta2"; + histtitle = "second leading Jet Eta"; + bookME( + ibooker, jetEta2ME_, histname, histtitle, dijet_eta_binning.nbins, dijet_eta_binning.xmin, dijet_eta_binning.xmax); + setMETitle(jetEta2ME_, "Jet_Eta_2", "events"); + histname = "jetptAvgB"; histtitle = "Pt average before offline selection"; bookME( @@ -261,6 +290,17 @@ void DiJetMonitor::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iR dijet_phi_binning.xmin, dijet_phi_binning.xmax); setMETitle(jetphiPrbME_, "Phi_probe #phi", "events"); + + histname = "jetphiTag"; + histtitle = "Tag Jet phi"; + bookME(ibooker, + jetphiTagME_, + histname, + histtitle, + dijet_phi_binning.nbins, + dijet_phi_binning.xmin, + dijet_phi_binning.xmax); + setMETitle(jetphiTagME_, "Phi_tag #phi", "events"); } void DiJetMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) { @@ -312,6 +352,10 @@ void DiJetMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSet jetpt1ME_.denominator->Fill(pt_1); jetpt2ME_.denominator->Fill(pt_2); + jetPhi1ME_.denominator->Fill(phi_1); + jetPhi2ME_.denominator->Fill(phi_2); + jetEta1ME_.denominator->Fill(eta_1); + jetEta2ME_.denominator->Fill(eta_2); jetptAvgbME_.denominator->Fill(pt_avg_b); if (dijet_selection(eta_1, phi_1, eta_2, phi_2, pt_1, pt_2, tag_id, probe_id, iEvent.id().event())) { @@ -326,6 +370,7 @@ void DiJetMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSet jetetaTagME_.denominator->Fill(eta_1); jetptAsyME_.denominator->Fill(pt_asy); jetphiPrbME_.denominator->Fill(phi_2); + jetphiTagME_.denominator->Fill(phi_1); jetAsyEtaME_.denominator->Fill(pt_asy, eta_2); jetEtaPhiME_.denominator->Fill(eta_2, phi_2); } @@ -340,6 +385,7 @@ void DiJetMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSet jetetaTagME_.denominator->Fill(eta_2); jetptAsyME_.denominator->Fill(pt_asy); jetphiPrbME_.denominator->Fill(phi_1); + jetphiTagME_.denominator->Fill(phi_2); jetAsyEtaME_.denominator->Fill(pt_asy, eta_1); jetEtaPhiME_.denominator->Fill(eta_1, phi_1); } @@ -349,6 +395,10 @@ void DiJetMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSet jetpt1ME_.numerator->Fill(pt_1); jetpt2ME_.numerator->Fill(pt_2); + jetPhi1ME_.numerator->Fill(phi_1); + jetPhi2ME_.numerator->Fill(phi_2); + jetEta1ME_.numerator->Fill(eta_1); + jetEta2ME_.numerator->Fill(eta_2); jetptAvgbME_.numerator->Fill(pt_avg_b); if (tag_id == 0 && probe_id == 1) { @@ -362,6 +412,7 @@ void DiJetMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSet jetetaTagME_.numerator->Fill(eta_1); jetptAsyME_.numerator->Fill(pt_asy); jetphiPrbME_.numerator->Fill(phi_2); + jetphiTagME_.numerator->Fill(phi_1); jetAsyEtaME_.numerator->Fill(pt_asy, eta_2); jetEtaPhiME_.numerator->Fill(eta_2, phi_2); } @@ -376,6 +427,7 @@ void DiJetMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSet jetetaTagME_.numerator->Fill(eta_2); jetptAsyME_.numerator->Fill(pt_asy); jetphiPrbME_.numerator->Fill(phi_1); + jetphiTagME_.numerator->Fill(phi_2); jetAsyEtaME_.numerator->Fill(pt_asy, eta_1); jetEtaPhiME_.numerator->Fill(eta_1, phi_1); } @@ -436,6 +488,7 @@ void DiJetMonitor::fillDescriptions(edm::ConfigurationDescriptions& descriptions edm::ParameterSetDescription genericTriggerEventPSet; GenericTriggerEventFlag::fillPSetDescription(genericTriggerEventPSet); + desc.add("numGenericTriggerEventPSet", genericTriggerEventPSet); desc.add("denGenericTriggerEventPSet", genericTriggerEventPSet); diff --git a/DQMOffline/Trigger/plugins/HTMonitor.cc b/DQMOffline/Trigger/plugins/HTMonitor.cc index 32ee2dda234b5..1c3d41f807816 100644 --- a/DQMOffline/Trigger/plugins/HTMonitor.cc +++ b/DQMOffline/Trigger/plugins/HTMonitor.cc @@ -53,6 +53,7 @@ class HTMonitor : public DQMEDAnalyzer, public TriggerDQMBase { const bool requireValidHLTPaths_; bool hltPathsAreValid_; + const bool enableFullMonitoring_; edm::InputTag metInputTag_; edm::InputTag jetInputTag_; @@ -67,13 +68,21 @@ class HTMonitor : public DQMEDAnalyzer, public TriggerDQMBase { edm::EDGetTokenT vtxToken_; std::vector ht_variable_binning_; + std::vector jetptBinning_; MEbinning ht_binning_; MEbinning ls_binning_; ObjME qME_variableBinning_; ObjME htVsLS_; - ObjME deltaphimetj1ME_; - ObjME deltaphij1j2ME_; + ObjME METPhiME_; + ObjME jetpt1ME_; + ObjME jetpt2ME_; + ObjME phij1ME_; + ObjME phij2ME_; + ObjME etaj1ME_; + ObjME etaj2ME_; + ObjME nJetsME_; //number of jets passing jet selection (Pt>0) + ObjME nJetsHTME_; //number of jets passing jet selection HT (Pt>30 & eta<2.5) std::unique_ptr num_genTriggerEventFlag_; std::unique_ptr den_genTriggerEventFlag_; @@ -91,17 +100,30 @@ class HTMonitor : public DQMEDAnalyzer, public TriggerDQMBase { static constexpr double MAXedge_PHI = 3.2; static constexpr int Nbin_PHI = 64; static constexpr MEbinning phi_binning_{Nbin_PHI, -MAXedge_PHI, MAXedge_PHI}; - + // Define Eta Bining + static constexpr double MAX_ETA = 5.0; + static constexpr int N_ETA = 20; + static constexpr MEbinning eta_binning{N_ETA, -MAX_ETA, MAX_ETA}; + //Define nJets Binning HT selection Pt>30 && eta<2.4 + static constexpr int MIN_NJETS_HT = 0; + static constexpr int MAX_NJETS_HT = 30; + static constexpr int N_BIN_NJETS_HT = 30; + static constexpr MEbinning nJets_HT_binning{N_BIN_NJETS_HT, MIN_NJETS_HT, MAX_NJETS_HT}; + //Define nJets Binning general selection Pt>0 + static constexpr int MIN_NJETS = 0; + static constexpr int MAX_NJETS = 200; + static constexpr int N_BIN_NJETS = 200; + static constexpr MEbinning nJets_binning{N_BIN_NJETS, MIN_NJETS, MAX_NJETS}; + + std::string quantity_; bool warningWasPrinted_; - - enum quant { HT, MJJ, SOFTDROP }; - quant quantity_; }; HTMonitor::HTMonitor(const edm::ParameterSet& iConfig) : folderName_(iConfig.getParameter("FolderName")), requireValidHLTPaths_(iConfig.getParameter("requireValidHLTPaths")), hltPathsAreValid_(false), + enableFullMonitoring_(iConfig.getParameter("enableFullMonitoring")), metInputTag_(iConfig.getParameter("met")), jetInputTag_(iConfig.getParameter("jets")), eleInputTag_(iConfig.getParameter("electrons")), @@ -114,6 +136,8 @@ HTMonitor::HTMonitor(const edm::ParameterSet& iConfig) vtxToken_(mayConsume(vtxInputTag_)), ht_variable_binning_( iConfig.getParameter("histoPSet").getParameter >("htBinning")), + jetptBinning_( + iConfig.getParameter("histoPSet").getParameter >("jetptBinning")), ht_binning_( getHistoPSet(iConfig.getParameter("histoPSet").getParameter("htPSet"))), ls_binning_( @@ -131,20 +155,8 @@ HTMonitor::HTMonitor(const edm::ParameterSet& iConfig) nelectrons_(iConfig.getParameter("nelectrons")), nmuons_(iConfig.getParameter("nmuons")), dEtaCut_(iConfig.getParameter("dEtaCut")), - warningWasPrinted_(false) { - /* mia: THIS CODE SHOULD BE DELETED !!!! */ - string quantity = iConfig.getParameter("quantity"); - if (quantity == "HT") { - quantity_ = HT; - } else if (quantity == "Mjj") { - quantity_ = MJJ; - } else if (quantity == "softdrop") { - quantity_ = SOFTDROP; - } else { - throw cms::Exception("quantity not defined") - << "the quantity '" << quantity << "' is undefined. Please check your config!" << std::endl; - } -} + quantity_(iConfig.getParameter("quantity")), + warningWasPrinted_(false) {} HTMonitor::~HTMonitor() throw() { if (num_genTriggerEventFlag_) { @@ -179,55 +191,87 @@ void HTMonitor::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iRun, std::string currentFolder = folderName_; ibooker.setCurrentFolder(currentFolder); - - switch (quantity_) { - case HT: { - histname = "ht_variable"; - histtitle = "HT"; - bookME(ibooker, qME_variableBinning_, histname, histtitle, ht_variable_binning_); - setMETitle(qME_variableBinning_, "HT [GeV]", "events / [GeV]"); - - histname = "htVsLS"; - histtitle = "HT vs LS"; - bookME(ibooker, - htVsLS_, - histname, - histtitle, - ls_binning_.nbins, - ls_binning_.xmin, - ls_binning_.xmax, - ht_binning_.xmin, - ht_binning_.xmax); - setMETitle(htVsLS_, "LS", "HT [GeV]"); - - histname = "deltaphi_metjet1"; - histtitle = "DPHI_METJ1"; - bookME(ibooker, deltaphimetj1ME_, histname, histtitle, phi_binning_.nbins, phi_binning_.xmin, phi_binning_.xmax); - setMETitle(deltaphimetj1ME_, "delta phi (met, j1)", "events / 0.1 rad"); - - histname = "deltaphi_jet1jet2"; - histtitle = "DPHI_J1J2"; - bookME(ibooker, deltaphij1j2ME_, histname, histtitle, phi_binning_.nbins, phi_binning_.xmin, phi_binning_.xmax); - setMETitle(deltaphij1j2ME_, "delta phi (j1, j2)", "events / 0.1 rad"); - break; - } - - case MJJ: { - histname = "mjj_variable"; - histtitle = "Mjj"; - bookME(ibooker, qME_variableBinning_, histname, histtitle, ht_variable_binning_); - setMETitle(qME_variableBinning_, "Mjj [GeV]", "events / [GeV]"); - break; + if (quantity_ == "HT") { + histname = "ht_variable"; + histtitle = "HT"; + bookME(ibooker, qME_variableBinning_, histname, histtitle, ht_variable_binning_); + setMETitle(qME_variableBinning_, "HT [GeV]", "events / [GeV]"); + + histname = "htVsLS"; + histtitle = "HT vs LS"; + bookME(ibooker, + htVsLS_, + histname, + histtitle, + ls_binning_.nbins, + ls_binning_.xmin, + ls_binning_.xmax, + ht_binning_.xmin, + ht_binning_.xmax); + setMETitle(htVsLS_, "LS", "HT [GeV]"); + + if (!enableFullMonitoring_) { + return; } - case SOFTDROP: { - histname = "softdrop_variable"; - histtitle = "softdropmass"; - bookME(ibooker, qME_variableBinning_, histname, histtitle, ht_variable_binning_); - setMETitle(qME_variableBinning_, "leading jet softdropmass [GeV]", "events / [GeV]"); - break; - } - } + histname = "METPhi"; + histtitle = "METPhi"; + bookME(ibooker, METPhiME_, histname, histtitle, phi_binning_.nbins, phi_binning_.xmin, phi_binning_.xmax); + setMETitle(METPhiME_, "MET phi", "events / 0.1 rad"); + + histname = "jetPt1"; + histtitle = "leading Jet Pt"; + bookME(ibooker, jetpt1ME_, histname, histtitle, jetptBinning_); + setMETitle(jetpt1ME_, "Pt_1 [GeV]", "events"); + + histname = "jetPt2"; + histtitle = "second leading Jet Pt"; + bookME(ibooker, jetpt2ME_, histname, histtitle, jetptBinning_); + setMETitle(jetpt2ME_, "Pt_2 [GeV]", "events"); + + histname = "jetEta1"; + histtitle = "leading Jet eta"; + bookME(ibooker, etaj1ME_, histname, histtitle, eta_binning.nbins, eta_binning.xmin, eta_binning.xmax); + setMETitle(etaj1ME_, "Jet_1 #eta", "events"); + + histname = "jetEta2"; + histtitle = "subleading Jet eta"; + bookME(ibooker, etaj2ME_, histname, histtitle, eta_binning.nbins, eta_binning.xmin, eta_binning.xmax); + setMETitle(etaj2ME_, "Jet_2 #eta", "events"); + + histname = "jetPhi1"; + histtitle = "leading Jet phi"; + bookME(ibooker, phij1ME_, histname, histtitle, phi_binning_.nbins, phi_binning_.xmin, phi_binning_.xmax); + setMETitle(phij1ME_, "Jet_1 #phi", "events / 0.1 rad"); + + histname = "jetPhi2"; + histtitle = "subleading Jet phi"; + bookME(ibooker, phij2ME_, histname, histtitle, phi_binning_.nbins, phi_binning_.xmin, phi_binning_.xmax); + setMETitle(phij2ME_, "Jet_2 #phi", "events / 0.1 rad"); + + histname = "nJets"; + histtitle = "number of Jets"; + bookME(ibooker, nJetsME_, histname, histtitle, nJets_binning.nbins, nJets_binning.xmin, nJets_binning.xmax); + setMETitle(nJetsME_, "number of Jets", "events"); + + histname = "nJetsHT"; + histtitle = "number of Jets HT"; + bookME( + ibooker, nJetsHTME_, histname, histtitle, nJets_HT_binning.nbins, nJets_HT_binning.xmin, nJets_HT_binning.xmax); + setMETitle(nJetsHTME_, "number of Jets HT", "events"); + } //end if -- HT quantity + else if (quantity_ == "Mjj") { + histname = "mjj_variable"; + histtitle = "Mjj"; + bookME(ibooker, qME_variableBinning_, histname, histtitle, ht_variable_binning_); + setMETitle(qME_variableBinning_, "Mjj [GeV]", "events / [GeV]"); + } //end if -- Mjj quantity + else if (quantity_ == "softdrop") { + histname = "softdrop_variable"; + histtitle = "softdropmass"; + bookME(ibooker, qME_variableBinning_, histname, histtitle, ht_variable_binning_); + setMETitle(qME_variableBinning_, "leading jet softdropmass [GeV]", "events / [GeV]"); + } //end if -- softdrop quantity } void HTMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) { @@ -277,13 +321,14 @@ void HTMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) if (jets.size() < njets_) return; - float deltaPhi_met_j1 = 10.0; - float deltaPhi_j1_j2 = 10.0; - - if (!jets.empty()) - deltaPhi_met_j1 = fabs(deltaPhi(pfmet.phi(), jets[0].phi())); - if (jets.size() >= 2) - deltaPhi_j1_j2 = fabs(deltaPhi(jets[0].phi(), jets[1].phi())); + const float metphi = pfmet.phi(); + const int nJetsSel = jets.size(); //jetSelection: PT>0 + const float Pt_J1 = !jets.empty() ? jets[0].pt() : -10.0; + const float Pt_J2 = jets.size() >= 2 ? jets[1].pt() : -10.0; + const float Phi_J1 = !jets.empty() ? jets[0].phi() : -10.0; + const float Phi_J2 = jets.size() >= 2 ? jets[1].phi() : -10.0; + const float Eta_J1 = !jets.empty() ? jets[0].p4().eta() : -10.0; + const float Eta_J2 = jets.size() >= 2 ? jets[1].p4().eta() : -10.0; std::vector electrons; edm::Handle eleHandle; @@ -368,82 +413,98 @@ void HTMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) } // fill histograms - switch (quantity_) { - case HT: { - float ht = 0.0; - for (auto const& j : *jetHandle) { - if (jetSelection_HT_(j)) - ht += j.pt(); + if (quantity_ == "HT") { + float ht = 0.0; + int nJets_HT = 0; + for (auto const& j : *jetHandle) { + if (jetSelection_HT_(j)) { + ht += j.pt(); + nJets_HT = nJets_HT + 1; //===== jetSelection HT: Pt>30 & eta<2.5 } + } - // filling histograms (denominator) - qME_variableBinning_.denominator->Fill(ht); - - deltaphimetj1ME_.denominator->Fill(deltaPhi_met_j1); - deltaphij1j2ME_.denominator->Fill(deltaPhi_j1_j2); - - int ls = iEvent.id().luminosityBlock(); - htVsLS_.denominator->Fill(ls, ht); + // filling histograms (denominator) + qME_variableBinning_.denominator->Fill(ht); - // applying selection for numerator - if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup)) - return; + int ls = iEvent.id().luminosityBlock(); + htVsLS_.denominator->Fill(ls, ht); - // filling histograms (num_genTriggerEventFlag_) - qME_variableBinning_.numerator->Fill(ht); + if (enableFullMonitoring_) { //===check the flag - htVsLS_.numerator->Fill(ls, ht); - deltaphimetj1ME_.numerator->Fill(deltaPhi_met_j1); - deltaphij1j2ME_.numerator->Fill(deltaPhi_j1_j2); - break; + METPhiME_.denominator->Fill(metphi); + jetpt1ME_.denominator->Fill(Pt_J1); + jetpt2ME_.denominator->Fill(Pt_J2); + phij1ME_.denominator->Fill(Phi_J1); + phij2ME_.denominator->Fill(Phi_J2); + etaj1ME_.denominator->Fill(Eta_J1); + etaj2ME_.denominator->Fill(Eta_J2); + nJetsME_.denominator->Fill(nJetsSel); + nJetsHTME_.denominator->Fill(nJets_HT); } - case MJJ: { - if (jets.size() < 2) - return; + // applying selection for numerator + if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup)) + return; - // deltaEta cut - if (fabs(jets[0].p4().Eta() - jets[1].p4().Eta()) >= dEtaCut_) - return; - float mjj = (jets[0].p4() + jets[1].p4()).M(); + // filling histograms (num_genTriggerEventFlag_) + qME_variableBinning_.numerator->Fill(ht); + + htVsLS_.numerator->Fill(ls, ht); + + if (enableFullMonitoring_) { //===check the flag + METPhiME_.numerator->Fill(metphi); + jetpt1ME_.numerator->Fill(Pt_J1); + jetpt2ME_.numerator->Fill(Pt_J2); + phij1ME_.numerator->Fill(Phi_J1); + phij2ME_.numerator->Fill(Phi_J2); + etaj1ME_.numerator->Fill(Eta_J1); + etaj2ME_.numerator->Fill(Eta_J2); + nJetsME_.numerator->Fill(nJetsSel); + nJetsHTME_.numerator->Fill(nJets_HT); + } + } //end if -- HT + else if (quantity_ == "Mjj") { + if (jets.size() < 2) + return; - qME_variableBinning_.denominator->Fill(mjj); + // deltaEta cut + if (fabs(jets[0].p4().Eta() - jets[1].p4().Eta()) >= dEtaCut_) + return; + float mjj = (jets[0].p4() + jets[1].p4()).M(); - // applying selection for numerator - if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup)) - return; + qME_variableBinning_.denominator->Fill(mjj); - qME_variableBinning_.numerator->Fill(mjj); - break; - } + // applying selection for numerator + if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup)) + return; - case SOFTDROP: { - if (jets.size() < 2) - return; + qME_variableBinning_.numerator->Fill(mjj); + } // end if -- Mjj + else if (quantity_ == "softdrop") { + if (jets.size() < 2) + return; - // deltaEta cut - if (fabs(jets[0].p4().Eta() - jets[1].p4().Eta()) >= dEtaCut_) - return; + // deltaEta cut + if (fabs(jets[0].p4().Eta() - jets[1].p4().Eta()) >= dEtaCut_) + return; - float softdrop = jets[0].p4().M(); + float softdrop = jets[0].p4().M(); - qME_variableBinning_.denominator->Fill(softdrop); + qME_variableBinning_.denominator->Fill(softdrop); - // applying selection for numerator - if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup)) - return; + // applying selection for numerator + if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup)) + return; - qME_variableBinning_.numerator->Fill(softdrop); - break; - } - } + qME_variableBinning_.numerator->Fill(softdrop); + } //end if -- softdrop } void HTMonitor::fillDescriptions(edm::ConfigurationDescriptions& descriptions) { edm::ParameterSetDescription desc; desc.add("FolderName", "HLT/HT"); desc.add("requireValidHLTPaths", true); - + desc.add("enableFullMonitoring", false); desc.add("met", edm::InputTag("pfMet")); desc.add("jets", edm::InputTag("ak4PFJetsCHS")); desc.add("electrons", edm::InputTag("gedGsfElectrons")); @@ -461,6 +522,7 @@ void HTMonitor::fillDescriptions(edm::ConfigurationDescriptions& descriptions) { edm::ParameterSetDescription genericTriggerEventPSet; GenericTriggerEventFlag::fillPSetDescription(genericTriggerEventPSet); + desc.add("numGenericTriggerEventPSet", genericTriggerEventPSet); desc.add("denGenericTriggerEventPSet", genericTriggerEventPSet); @@ -473,7 +535,10 @@ void HTMonitor::fillDescriptions(edm::ConfigurationDescriptions& descriptions) { 350., 400., 450., 500., 550., 600., 650., 700., 750., 800., 850., 900., 950., 1000., 1050., 1100., 1200., 1300., 1400., 1500., 2000., 2500.}; histoPSet.add >("htBinning", bins); - + std::vector bins_2 = { + 0., 20., 40., 60., 80., 90., 100., 110., 120., 130., 140., 150., 160., + 170., 180., 190., 200., 220., 240., 260., 280., 300., 350., 400., 450., 1000.}; // Jet pT Binning + histoPSet.add >("jetptBinning", bins_2); edm::ParameterSetDescription lsPSet; fillHistoLSPSetDescription(lsPSet); histoPSet.add("lsPSet", lsPSet); diff --git a/DQMOffline/Trigger/plugins/JetMonitor.cc b/DQMOffline/Trigger/plugins/JetMonitor.cc index 4fb8852ee2201..c8e1727344a25 100644 --- a/DQMOffline/Trigger/plugins/JetMonitor.cc +++ b/DQMOffline/Trigger/plugins/JetMonitor.cc @@ -34,9 +34,6 @@ class JetMonitor : public DQMEDAnalyzer, public TriggerDQMBase { bool isEndCapP(double eta); bool isEndCapM(double eta); bool isForward(double eta); - bool isHEP17(double eta, double phi); - bool isHEM17(double eta, double phi); - bool isHEP18(double eta, double phi); // -0.87< Phi < -1.22 void bookMESub(DQMStore::IBooker&, ObjME* a_me, @@ -70,6 +67,8 @@ class JetMonitor : public DQMEDAnalyzer, public TriggerDQMBase { bool isPFJetTrig; bool isCaloJetTrig; + const bool enableFullMonitoring_; + edm::EDGetTokenT > jetSrc_; std::unique_ptr num_genTriggerEventFlag_; @@ -85,28 +84,14 @@ class JetMonitor : public DQMEDAnalyzer, public TriggerDQMBase { ObjME a_ME_HF[7]; ObjME a_ME_HE_p[7]; ObjME a_ME_HE_m[7]; - ObjME a_ME_HEM17[7]; - ObjME a_ME_HEP17[7]; - ObjME a_ME_HEP18[7]; - - ObjME jetHEP17_AbsEtaVsPhi_; - ObjME jetHEM17_AbsEtaVsPhi_; - ObjME jetHEP17_AbsEta_; - ObjME jetHEM17_AbsEta_; std::vector v_jetpt; std::vector v_jeteta; std::vector v_jetphi; // (mia) not optimal, we should make use of variable binning which reflects the detector ! - MEbinning jet_phi_binning_{32, -3.2, 3.2}; - MEbinning jet_eta_binning_{20, -5, 5}; - - MEbinning eta_binning_hep17_{9, 1.3, 3.0}; - MEbinning eta_binning_hem17_{9, -3.0, -1.3}; - - MEbinning phi_binning_hep17_{7, -0.87, -0.52}; - MEbinning phi_binning_hep18_{7, -0.52, -0.17}; + MEbinning jet_phi_binning_{64, -3.2, 3.2}; + MEbinning jet_eta_binning_{50, -5, 5}; }; JetMonitor::JetMonitor(const edm::ParameterSet& iConfig) @@ -116,6 +101,7 @@ JetMonitor::JetMonitor(const edm::ParameterSet& iConfig) ptcut_(iConfig.getParameter("ptcut")), isPFJetTrig(iConfig.getParameter("ispfjettrg")), isCaloJetTrig(iConfig.getParameter("iscalojettrg")), + enableFullMonitoring_(iConfig.getParameter("enableFullMonitoring")), jetSrc_(mayConsume >(iConfig.getParameter("jetSrc"))), num_genTriggerEventFlag_(new GenericTriggerEventFlag( iConfig.getParameter("numGenericTriggerEventPSet"), consumesCollector(), *this)), @@ -208,6 +194,12 @@ void JetMonitor::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iRun true, true, false); + + //check the flag + if (!enableFullMonitoring_) { + return; + } + bookMESub(ibooker, a_ME_HE_p, sizeof(a_ME_HE_p) / sizeof(a_ME_HE_p[0]), @@ -216,7 +208,7 @@ void JetMonitor::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iRun "HE_p", "(HE+)", true, - true, + false, true, false); bookMESub(ibooker, @@ -227,76 +219,9 @@ void JetMonitor::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iRun "HE_m", "(HE-)", true, - true, - true, - false); - bookMESub(ibooker, - a_ME_HEP17, - sizeof(a_ME_HEP17) / sizeof(a_ME_HEP17[0]), - hist_obtag, - histtitle_obtag, - "HEP17", - "(HEP17)", - true, false, - false, - false); - bookMESub(ibooker, - a_ME_HEM17, - sizeof(a_ME_HEM17) / sizeof(a_ME_HEM17[0]), - hist_obtag, - histtitle_obtag, - "HEM17", - "(HEM17)", true, - false, - false, - false); - bookMESub(ibooker, - a_ME_HEP18, - sizeof(a_ME_HEP18) / sizeof(a_ME_HEP18[0]), - hist_obtag, - histtitle_obtag, - "HEP18", - "(HEP18)", - false, - false, - false, false); - - /* - WE WOULD NEED TURNON CURVES TO BE COMPARED NOT JUST THE ZOOM OF A 2D MAP !!! - - histname = hist_obtag +"AbsEtaVsPhi_HEP17"; histtitle = histtitle_obtag + " |eta| Vs phi (HEP17) "; - bookME(ibooker,jetHEP17_AbsEtaVsPhi_,histname,histtitle, eta_binning_hep17_.nbins, eta_binning_hep17_.xmin, eta_binning_hep17_.xmax, phi_binning_hep17_.nbins,phi_binning_hep17_.xmin,phi_binning_hep17_.xmax); - setMETitle(jetHEP17_AbsEtaVsPhi_,histtitle_obtag + " |#eta|","#phi"); - - histname = hist_obtag +"AbsEtaVsPhi_HEM17"; histtitle = histtitle_obtag + " |eta| Vs phi (HEM17) "; - bookME(ibooker,jetHEM17_AbsEtaVsPhi_,histname,histtitle, eta_binning_hep17_.nbins, eta_binning_hep17_.xmin, eta_binning_hep17_.xmax, phi_binning_hep17_.nbins,phi_binning_hep17_.xmin,phi_binning_hep17_.xmax); - setMETitle(jetHEM17_AbsEtaVsPhi_,histtitle_obtag + " |#eta|","#phi"); - */ - - histname = hist_obtag + "abseta_HEP17"; - histtitle = histtitle_obtag + " |#eta| (HEP17) "; - bookME(ibooker, - jetHEP17_AbsEta_, - histname, - histtitle, - eta_binning_hep17_.nbins, - eta_binning_hep17_.xmin, - eta_binning_hep17_.xmax); - setMETitle(jetHEP17_AbsEta_, histtitle_obtag + " |#eta|", "events / |#eta|"); - - histname = hist_obtag + "abseta_HEM17"; - histtitle = histtitle_obtag + " |eta| (HEM17) "; - bookME(ibooker, - jetHEM17_AbsEta_, - histname, - histtitle, - eta_binning_hep17_.nbins, - eta_binning_hep17_.xmin, - eta_binning_hep17_.xmax); - setMETitle(jetHEM17_AbsEta_, histtitle_obtag + " |#eta|", "events / |#eta|"); } void JetMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) { @@ -325,7 +250,6 @@ void JetMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup for (edm::View::const_iterator ibegin = offjets->begin(), iend = offjets->end(), ijet = ibegin; ijet != iend; ++ijet) { - //if (ijet->pt()< 20) {continue;} if (ijet->pt() < ptcut_) { continue; } @@ -346,35 +270,18 @@ void JetMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup FillME(a_ME_HB, jetpt_, jetphi_, jeteta_, ls, "denominator", true, true, true, false); } else if (isEndCapP(jeteta_)) { FillME(a_ME_HE, jetpt_, jetphi_, jeteta_, ls, "denominator", true, true, true, false); - FillME(a_ME_HE_p, jetpt_, jetphi_, jeteta_, ls, "denominator", true, true, true, false); + if (enableFullMonitoring_) { + FillME(a_ME_HE_p, jetpt_, jetphi_, jeteta_, ls, "denominator", true, false, true, false); + } } else if (isEndCapM(jeteta_)) { FillME(a_ME_HE, jetpt_, jetphi_, jeteta_, ls, "denominator", true, true, true, false); - FillME(a_ME_HE_m, jetpt_, jetphi_, jeteta_, ls, "denominator", true, true, true, false); + if (enableFullMonitoring_) { + FillME(a_ME_HE_m, jetpt_, jetphi_, jeteta_, ls, "denominator", true, false, true, false); + } } else if (isForward(jeteta_)) { FillME(a_ME_HF, jetpt_, jetphi_, jeteta_, ls, "denominator", true, true, true, false); } - if (isHEP17(jeteta_, jetphi_)) { - FillME(a_ME_HEP17, - jetpt_, - jetphi_, - jeteta_, - ls, - "denominator", - true, - false, - false, - false); // doPhi, doEta, doEtaPhi, doVsLS - jetHEP17_AbsEta_.denominator->Fill(abs(jeteta_)); - } else if (isHEM17(jeteta_, jetphi_)) { - FillME( - a_ME_HEM17, jetpt_, jetphi_, jeteta_, ls, "denominator", true, false, false, false); // doPhi, doEta, doEtaPhi - jetHEM17_AbsEta_.denominator->Fill(abs(jeteta_)); - } else if (isHEP18(jeteta_, jetphi_)) { - FillME( - a_ME_HEP18, jetpt_, jetphi_, jeteta_, ls, "denominator", false, false, false, false); // doPhi, doEta, doEtaPhi - } - if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup)) return; // Require Numerator // @@ -383,50 +290,17 @@ void JetMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup FillME(a_ME_HB, jetpt_, jetphi_, jeteta_, ls, "numerator", true, true, true, false); } else if (isEndCapP(jeteta_)) { FillME(a_ME_HE, jetpt_, jetphi_, jeteta_, ls, "numerator", true, true, true, false); - FillME(a_ME_HE_p, jetpt_, jetphi_, jeteta_, ls, "numerator", true, true, true, false); + if (enableFullMonitoring_) { + FillME(a_ME_HE_p, jetpt_, jetphi_, jeteta_, ls, "numerator", true, false, true, false); + } } else if (isEndCapM(jeteta_)) { FillME(a_ME_HE, jetpt_, jetphi_, jeteta_, ls, "numerator", true, true, true, false); - FillME(a_ME_HE_m, jetpt_, jetphi_, jeteta_, ls, "numerator", true, true, true, false); + if (enableFullMonitoring_) { + FillME(a_ME_HE_m, jetpt_, jetphi_, jeteta_, ls, "numerator", true, false, true, false); + } } else if (isForward(jeteta_)) { FillME(a_ME_HF, jetpt_, jetphi_, jeteta_, ls, "numerator", true, true, true, false); } - - if (isHEP17(jeteta_, jetphi_)) { - FillME(a_ME_HEP17, - jetpt_, - jetphi_, - jeteta_, - ls, - "numerator", - true, - false, - false, - false); // doPhi, doEta, doEtaPhi, doVsLS - jetHEP17_AbsEta_.numerator->Fill(abs(jeteta_)); - } else if (isHEM17(jeteta_, jetphi_)) { - FillME(a_ME_HEM17, - jetpt_, - jetphi_, - jeteta_, - ls, - "numerator", - true, - false, - false, - false); // doPhi, doEta, doEtaPhi, doVsLS - jetHEM17_AbsEta_.numerator->Fill(abs(jeteta_)); - } else if (isHEP18(jeteta_, jetphi_)) { - FillME(a_ME_HEP18, - jetpt_, - jetphi_, - jeteta_, - ls, - "numerator", - false, - false, - false, - false); // doPhi, doEta, doEtaPhi, doVsLS - } } bool JetMonitor::isBarrel(double eta) { @@ -446,13 +320,12 @@ bool JetMonitor::isEndCapM(double eta) { /// For Hcal Endcap Plus Area bool JetMonitor::isEndCapP(double eta) { bool output = false; - //if ( eta<=3.0 && eta >1.3) output=true; if (fabs(eta) <= 3.0 && fabs(eta) > 1.3 && (eta > 0)) output = true; // (mia) this magic number should come from some file in CMSSW !!! return output; } -/// For Hcal Forward Plus Area +/// For Hcal Forward Area bool JetMonitor::isForward(double eta) { bool output = false; if (fabs(eta) > 3.0) @@ -460,35 +333,6 @@ bool JetMonitor::isForward(double eta) { return output; } -/// For Hcal HEP17 Area -bool JetMonitor::isHEP17(double eta, double phi) { - bool output = false; - // phi -0.87 to -0.52 - if (fabs(eta) <= 3.0 && fabs(eta) > 1.3 && (eta > 0) && phi > -0.87 && phi <= -0.52) { - output = true; - } // (mia) this magic number should come from some file in CMSSW !!! - return output; -} - -/// For Hcal HEM17 Area -bool JetMonitor::isHEM17(double eta, double phi) { - bool output = false; - if (fabs(eta) <= 3.0 && fabs(eta) > 1.3 && (eta < 0) && phi > -0.87 && phi <= -0.52) { - output = true; - } // (mia) this magic number should come from some file in CMSSW !!! - return output; -} - -/// For Hcal HEP18 Area -bool JetMonitor::isHEP18(double eta, double phi) { - bool output = false; - // phi -0.87 to -0.52 - if (fabs(eta) <= 3.0 && fabs(eta) > 1.3 && (eta > 0) && phi > -0.52 && phi <= -0.17) { - output = true; - } // (mia) this magic number should come from some file in CMSSW !!! - return output; -} - void JetMonitor::FillME(ObjME* a_me, const double pt_, const double phi_, @@ -512,7 +356,7 @@ void JetMonitor::FillME(ObjME* a_me, a_me[4].denominator->Fill(eta_); // eta if (doEtaPhi) a_me[5].denominator->Fill(eta_, phi_); // eta vs phi - if (doEta) + if (doEtaPhi) a_me[6].denominator->Fill(eta_, pt_); // eta vs pT } else if (DenoOrNume == "numerator") { a_me[0].numerator->Fill(pt_); // pt @@ -522,10 +366,10 @@ void JetMonitor::FillME(ObjME* a_me, if (doPhi) a_me[3].numerator->Fill(phi_); // phi if (doEta) - a_me[4].numerator->Fill(eta_); // eat + a_me[4].numerator->Fill(eta_); // eta if (doEtaPhi) a_me[5].numerator->Fill(eta_, phi_); // eta vs phi - if (doEta) + if (doEtaPhi) a_me[6].numerator->Fill(eta_, pt_); // eta vs pT } else { edm::LogWarning("JetMonitor") << "CHECK OUT denu option in FillME !!! DenoOrNume ? : " << DenoOrNume << std::endl; @@ -555,33 +399,6 @@ void JetMonitor::bookMESub(DQMStore::IBooker& Ibooker, double maxbin_eta = jet_eta_binning_.xmax; double minbin_eta = jet_eta_binning_.xmin; - if (h_subOptName == "HEP17") { - nbin_phi = phi_binning_hep17_.nbins; - maxbin_phi = phi_binning_hep17_.xmax; - minbin_phi = phi_binning_hep17_.xmin; - - nbin_eta = eta_binning_hep17_.nbins; - maxbin_eta = eta_binning_hep17_.xmax; - minbin_eta = eta_binning_hep17_.xmin; - } - if (h_subOptName == "HEM17") { - nbin_phi = phi_binning_hep17_.nbins; - maxbin_phi = phi_binning_hep17_.xmax; - minbin_phi = phi_binning_hep17_.xmin; - - nbin_eta = eta_binning_hem17_.nbins; - maxbin_eta = eta_binning_hem17_.xmax; - minbin_eta = eta_binning_hem17_.xmin; - } - if (h_subOptName == "HEP18") { - nbin_phi = phi_binning_hep18_.nbins; - maxbin_phi = phi_binning_hep18_.xmax; - minbin_phi = phi_binning_hep18_.xmin; - - nbin_eta = eta_binning_hep17_.nbins; - maxbin_eta = eta_binning_hep17_.xmax; - minbin_eta = eta_binning_hep17_.xmin; - } hName = h_Name + "pT" + hSubN; hTitle = h_Title + " pT " + hSubT; bookME(Ibooker, a_me[0], hName, hTitle, jetpt_binning_.nbins, jetpt_binning_.xmin, jetpt_binning_.xmax); @@ -618,7 +435,7 @@ void JetMonitor::bookMESub(DQMStore::IBooker& Ibooker, hName = h_Name + "eta" + hSubN; hTitle = h_Title + " eta " + hSubT; bookME(Ibooker, a_me[4], hName, hTitle, nbin_eta, minbin_eta, maxbin_eta); - setMETitle(a_me[4], h_Title + " #eta", "events / #eta"); + setMETitle(a_me[4], h_Title + " #eta", "events"); } if (doEtaPhi) { @@ -628,7 +445,7 @@ void JetMonitor::bookMESub(DQMStore::IBooker& Ibooker, setMETitle(a_me[5], h_Title + " #eta", "#phi"); } - if (doEta) { + if (doEtaPhi) { hName = h_Name + "EtaVspT" + hSubN; hTitle = h_Title + " eta Vs pT " + hSubT; bookME(Ibooker, @@ -655,6 +472,8 @@ void JetMonitor::fillDescriptions(edm::ConfigurationDescriptions& descriptions) desc.add("ispfjettrg", true); desc.add("iscalojettrg", false); + desc.add("enableFullMonitoring", true); + edm::ParameterSetDescription genericTriggerEventPSet; GenericTriggerEventFlag::fillPSetDescription(genericTriggerEventPSet); diff --git a/DQMOffline/Trigger/plugins/METMonitor.cc b/DQMOffline/Trigger/plugins/METMonitor.cc index f3237162f6c36..f98a8b836c4f6 100644 --- a/DQMOffline/Trigger/plugins/METMonitor.cc +++ b/DQMOffline/Trigger/plugins/METMonitor.cc @@ -42,6 +42,7 @@ class METMonitor : public DQMEDAnalyzer, public TriggerDQMBase { const bool requireValidHLTPaths_; bool hltPathsAreValid_; + const bool enableFullMonitoring_; edm::InputTag metInputTag_; edm::InputTag jetInputTag_; @@ -56,6 +57,7 @@ class METMonitor : public DQMEDAnalyzer, public TriggerDQMBase { edm::EDGetTokenT vtxToken_; std::vector met_variable_binning_; + std::vector jetptBinning_; MEbinning met_binning_; MEbinning ls_binning_; @@ -65,6 +67,16 @@ class METMonitor : public DQMEDAnalyzer, public TriggerDQMBase { ObjME metPhiME_; ObjME deltaphimetj1ME_; ObjME deltaphij1j2ME_; + ObjME phij1ME_; + ObjME phij2ME_; + ObjME jetPhiME_; + ObjME etaj1ME_; + ObjME etaj2ME_; + ObjME jetEtaME_; + ObjME jetpt1ME_; + ObjME jetpt2ME_; + ObjME jetPtME_; + ObjME nJetsME_; std::unique_ptr num_genTriggerEventFlag_; std::unique_ptr den_genTriggerEventFlag_; @@ -81,6 +93,15 @@ class METMonitor : public DQMEDAnalyzer, public TriggerDQMBase { static constexpr double MAX_PHI = 3.2; static constexpr int N_PHI = 64; static constexpr MEbinning phi_binning_{N_PHI, -MAX_PHI, MAX_PHI}; + // Define Eta Bining + static constexpr double MAX_ETA = 5.0; + static constexpr int N_ETA = 50; + static constexpr MEbinning eta_binning{N_ETA, -MAX_ETA, MAX_ETA}; + //Define nJets Binning general selection Pt>0 + static constexpr int MIN_NJETS = 0; + static constexpr int MAX_NJETS = 200; + static constexpr int N_BIN_NJETS = 200; + static constexpr MEbinning nJets_binning{N_BIN_NJETS, MIN_NJETS, MAX_NJETS}; std::vector warningPrinted4token_; }; @@ -89,6 +110,7 @@ METMonitor::METMonitor(const edm::ParameterSet& iConfig) : folderName_(iConfig.getParameter("FolderName")), requireValidHLTPaths_(iConfig.getParameter("requireValidHLTPaths")), hltPathsAreValid_(false), + enableFullMonitoring_(iConfig.getParameter("enableFullMonitoring")), metInputTag_(iConfig.getParameter("met")), jetInputTag_(iConfig.getParameter("jets")), eleInputTag_(iConfig.getParameter("electrons")), @@ -101,6 +123,8 @@ METMonitor::METMonitor(const edm::ParameterSet& iConfig) vtxToken_(mayConsume(vtxInputTag_)), met_variable_binning_( iConfig.getParameter("histoPSet").getParameter >("metBinning")), + jetptBinning_( + iConfig.getParameter("histoPSet").getParameter >("jetptBinning")), met_binning_(getHistoPSet( iConfig.getParameter("histoPSet").getParameter("metPSet"))), ls_binning_( @@ -195,6 +219,61 @@ void METMonitor::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iRun histtitle = "PFMET phi"; bookME(ibooker, metPhiME_, histname, histtitle, phi_binning_.nbins, phi_binning_.xmin, phi_binning_.xmax); setMETitle(metPhiME_, "PF MET #phi", "events / 0.1 rad"); + + histname = "jetEta"; + histtitle = "Jet eta"; + bookME(ibooker, jetEtaME_, histname, histtitle, eta_binning.nbins, eta_binning.xmin, eta_binning.xmax); + setMETitle(jetEtaME_, "Jet #eta", "events"); + + histname = "jetPhi"; + histtitle = "Jet phi"; + bookME(ibooker, jetPhiME_, histname, histtitle, phi_binning_.nbins, phi_binning_.xmin, phi_binning_.xmax); + setMETitle(jetPhiME_, "Jet #phi", "events / 0.1 rad"); + + histname = "jetPt"; + histtitle = "Jet Pt"; + bookME(ibooker, jetPtME_, histname, histtitle, jetptBinning_); + setMETitle(jetPtME_, "jet Pt[GeV]", "events"); + + histname = "NJets"; + histtitle = "number of Jets"; + bookME(ibooker, nJetsME_, histname, histtitle, nJets_binning.nbins, nJets_binning.xmin, nJets_binning.xmax); + setMETitle(nJetsME_, "number of Jets", "events"); + + //check the flag + if (!enableFullMonitoring_) { + return; + } + + histname = "jetEta_1"; + histtitle = "leading Jet eta"; + bookME(ibooker, etaj1ME_, histname, histtitle, eta_binning.nbins, eta_binning.xmin, eta_binning.xmax); + setMETitle(etaj1ME_, "Jet_1 #eta", "events"); + + histname = "jetEta_2"; + histtitle = "subleading Jet eta"; + bookME(ibooker, etaj2ME_, histname, histtitle, eta_binning.nbins, eta_binning.xmin, eta_binning.xmax); + setMETitle(etaj2ME_, "Jet_2 #eta", "events"); + + histname = "jetPhi_1"; + histtitle = "leading Jet phi"; + bookME(ibooker, phij1ME_, histname, histtitle, phi_binning_.nbins, phi_binning_.xmin, phi_binning_.xmax); + setMETitle(phij1ME_, "Jet_1 #phi", "events / 0.1 rad"); + + histname = "jetPhi_2"; + histtitle = "subleading Jet phi"; + bookME(ibooker, phij2ME_, histname, histtitle, phi_binning_.nbins, phi_binning_.xmin, phi_binning_.xmax); + setMETitle(phij2ME_, "Jet_2 #phi", "events / 0.1 rad"); + + histname = "jetPt_1"; + histtitle = "leading Jet Pt"; + bookME(ibooker, jetpt1ME_, histname, histtitle, jetptBinning_); + setMETitle(jetpt1ME_, "Pt_1 [GeV]", "events"); + + histname = "jetPt_2"; + histtitle = "second leading Jet Pt"; + bookME(ibooker, jetpt2ME_, histname, histtitle, jetptBinning_); + setMETitle(jetpt2ME_, "Pt_2 [GeV]", "events"); } void METMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) { @@ -222,8 +301,8 @@ void METMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup if (!metSelection_(pfmet)) return; - float met = pfmet.pt(); - float phi = pfmet.phi(); + const float met = pfmet.pt(); + const float phi = pfmet.phi(); std::vector jets; edm::Handle jetHandle; @@ -249,13 +328,15 @@ void METMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup if (!jetInputTag_.label().empty()) return; } - float deltaPhi_met_j1 = 10.0; - float deltaPhi_j1_j2 = 10.0; - - if (!jets.empty()) - deltaPhi_met_j1 = fabs(deltaPhi(pfmet.phi(), jets[0].phi())); - if (jets.size() >= 2) - deltaPhi_j1_j2 = fabs(deltaPhi(jets[0].phi(), jets[1].phi())); + const float deltaPhi_met_j1 = !jets.empty() ? fabs(deltaPhi(pfmet.phi(), jets[0].phi())) : -10.0; + const float deltaPhi_j1_j2 = jets.size() >= 2 ? fabs(deltaPhi(jets[0].phi(), jets[1].phi())) : -10.0; + const int nJetsSel = jets.size(); + const float Pt_J1 = !jets.empty() ? jets[0].pt() : -10.; + const float Pt_J2 = jets.size() >= 2 ? jets[1].pt() : -10.0; + const float Phi_J1 = !jets.empty() ? jets[0].phi() : -10.0; + const float Phi_J2 = jets.size() >= 2 ? jets[1].phi() : -10.0; + const float Eta_J1 = !jets.empty() ? jets[0].p4().eta() : -10.0; + const float Eta_J2 = jets.size() >= 2 ? jets[1].p4().eta() : -10.0; std::vector electrons; edm::Handle eleHandle; @@ -347,6 +428,22 @@ void METMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup const int ls = iEvent.id().luminosityBlock(); metVsLS_.denominator->Fill(ls, met); + + nJetsME_.denominator->Fill(nJetsSel); + for (auto const& jet : jets) { + jetPtME_.denominator->Fill(jet.pt()); + jetPhiME_.denominator->Fill(jet.phi()); + jetEtaME_.denominator->Fill(jet.p4().eta()); + } + if (enableFullMonitoring_) { //===check the flag + jetpt1ME_.denominator->Fill(Pt_J1); + jetpt2ME_.denominator->Fill(Pt_J2); + phij1ME_.denominator->Fill(Phi_J1); + phij2ME_.denominator->Fill(Phi_J2); + etaj1ME_.denominator->Fill(Eta_J1); + etaj2ME_.denominator->Fill(Eta_J2); + } + // applying selection for numerator if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup)) return; @@ -358,12 +455,29 @@ void METMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup metVsLS_.numerator->Fill(ls, met); deltaphimetj1ME_.numerator->Fill(deltaPhi_met_j1); deltaphij1j2ME_.numerator->Fill(deltaPhi_j1_j2); + + nJetsME_.numerator->Fill(nJetsSel); + for (auto const& jet : jets) { + jetPtME_.numerator->Fill(jet.pt()); + jetPhiME_.numerator->Fill(jet.phi()); + jetEtaME_.numerator->Fill(jet.p4().eta()); + } + if (enableFullMonitoring_) { //===check the flag + jetpt1ME_.numerator->Fill(Pt_J1); + jetpt2ME_.numerator->Fill(Pt_J2); + phij1ME_.numerator->Fill(Phi_J1); + phij2ME_.numerator->Fill(Phi_J2); + etaj1ME_.numerator->Fill(Eta_J1); + etaj2ME_.numerator->Fill(Eta_J2); + } } void METMonitor::fillDescriptions(edm::ConfigurationDescriptions& descriptions) { edm::ParameterSetDescription desc; desc.add("FolderName", "HLT/MET"); desc.add("requireValidHLTPaths", true); + //========== flag to enable more or less outputs in each trigger ===== + desc.add("enableFullMonitoring", false); desc.add("met", edm::InputTag("pfMet")); desc.add("jets", edm::InputTag("ak4PFJetsCHS")); @@ -391,6 +505,10 @@ void METMonitor::fillDescriptions(edm::ConfigurationDescriptions& descriptions) std::vector bins = {0., 20., 40., 60., 80., 90., 100., 110., 120., 130., 140., 150., 160., 170., 180., 190., 200., 220., 240., 260., 280., 300., 350., 400., 450., 1000.}; histoPSet.add >("metBinning", bins); + std::vector bins_ = { + 0., 20., 40., 60., 80., 90., 100., 110., 120., 130., 140., 150., 160., + 170., 180., 190., 200., 220., 240., 260., 280., 300., 350., 400., 450., 1000.}; // Jet pT Binning + histoPSet.add >("jetptBinning", bins_); edm::ParameterSetDescription lsPSet; fillHistoLSPSetDescription(lsPSet); diff --git a/DQMOffline/Trigger/python/HTMonitor_cff.py b/DQMOffline/Trigger/python/HTMonitor_cff.py index b0bac4b38dc8e..1e03db0401ff2 100644 --- a/DQMOffline/Trigger/python/HTMonitor_cff.py +++ b/DQMOffline/Trigger/python/HTMonitor_cff.py @@ -3,28 +3,19 @@ from DQMOffline.Trigger.HTMonitor_cfi import hltHTmonitoring -# HLT_PFMETNoMu90_PFMHTNoMu90_IDTight -PFMETNoMu90_PFMHTNoMu90_HTmonitoring = hltHTmonitoring.clone( - FolderName = 'HLT/JME/MET/PFMETNoMu90/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETNoMu90_PFMHTNoMu90_IDTight_v*"]) -) # HLT_PFMETNoMu120_PFMHTNoMu120_IDTight PFMETNoMu120_PFMHTNoMu120_HTmonitoring = hltHTmonitoring.clone( FolderName = 'HLT/JME/MET/PFMETNoMu120/', + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETNoMu120_PFMHTNoMu120_IDTight_v*"]) ) -# HLT_MET200 -MET200_HTmonitoring = hltHTmonitoring.clone( - FolderName = 'HLT/JME/MET/MET200/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_MET200_v*"]) -) - # HLT_MonoCentralPFJet80_PFMETNoMu110_PFMHTNoMu110_IDTight MonoCentralPFJet80_PFMETNoMu110_PFMHTNoMu110_HTmonitoring = hltHTmonitoring.clone( FolderName = 'HLT/EXO/MET/MonoCentralPFJet80_PFMETNoMu110/', jetSelection = "pt > 100 && eta < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_MonoCentralPFJet80_PFMETNoMu110_PFMHTNoMu110_IDTight_v*"]) ) # HLT_MonoCentralPFJet80_PFMETNoMu120_PFMHTNoMu120_IDTight @@ -32,6 +23,7 @@ FolderName = 'HLT/EXO/MET/MonoCentralPFJet80_PFMETNoMu120/', jetSelection = "pt > 100 && eta < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_MonoCentralPFJet80_PFMETNoMu120_PFMHTNoMu120_IDTight_v*"]) ) # HLT_MonoCentralPFJet80_PFMETNoMu130_PFMHTNoMu130_IDTight @@ -39,6 +31,7 @@ FolderName = 'HLT/EXO/MET/MonoCentralPFJet80_PFMETNoMu130/', jetSelection = "pt > 100 && eta < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_MonoCentralPFJet80_PFMETNoMu130_PFMHTNoMu130_IDTight_v*"]) ) # HLT_MonoCentralPFJet80_PFMETNoMu140_PFMHTNoMu140_IDTight @@ -46,19 +39,15 @@ FolderName = 'HLT/EXO/MET/MonoCentralPFJet80_PFMETNoMu140/', jetSelection = "pt > 100 && eta < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_MonoCentralPFJet80_PFMETNoMu140_PFMHTNoMu140_IDTight_v*"]) ) -# HLT_MonoCentralPFJet80_PFMETNoMu90_PFMHTNoMu90_IDTight -MonoCentralPFJet80_PFMETNoMu90_PFMHTNoMu90_HTmonitoring = hltHTmonitoring.clone( - FolderName = 'HLT/EXO/MET/MonoCentralPFJet80_PFMETNoMu90/', - jetSelection_HT = "pt > 30 && eta < 2.5", - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_MonoCentralPFJet80_PFMETNoMu110_PFMHTNoMu110_IDTight_v*"]) -) # HLT_PFHT350MinPFJet15 PFHT350MinPFJet15_HTmonitoring = hltHTmonitoring.clone( FolderName = 'HLT/HT/HLT_PFHT350MinPFJet15/', jetSelection = "pt > 15", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT350MinPFJet15_v*"]) ) # HLT_PFHT500_PFMET100_PFMHT100_IDTight @@ -67,6 +56,7 @@ metSelection = "pt > 200", jetSelection = "pt > 0", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT500_PFMET100_PFMHT100_IDTight_v*"]) ) # HLT_PFHT500_PFMET110_PFMHT110_IDTight @@ -75,6 +65,7 @@ metSelection = "pt > 210", jetSelection = "pt > 0", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT500_PFMET110_PFMHT110_IDTight_v*"]) ) # HLT_PFHT700_PFMET85_PFMHT85_IDTight @@ -83,6 +74,7 @@ metSelection = "pt > 185", jetSelection = "pt > 0", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT700_PFMET85_PFMHT85_IDTight_v*"]) ) # HLT_PFHT700_PFMET95_PFMHT95_IDTight @@ -91,6 +83,7 @@ metSelection = "pt > 195", jetSelection = "pt > 0", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT700_PFMET95_PFMHT95_IDTight_v*"]) ) # HLT_PFHT800_PFMET75_PFMHT75_IDTight @@ -99,6 +92,7 @@ metSelection = "pt > 175", jetSelection = "pt > 0", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT800_PFMET75_PFMHT75_IDTight_v*"]) ) # HLT_PFHT800_PFMET85_PFMHT85_IDTight @@ -107,6 +101,7 @@ metSelection = "pt > 185", jetSelection = "pt > 0", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT800_PFMET85_PFMHT85_IDTight_v*"]) ) # HLT_PFHT1050_v7 @@ -114,24 +109,28 @@ FolderName = 'HLT/JME/HT/PFHT1050/', jetSelection = "pt > 0", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = True, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT1050_v*"]) ) #HLT_PFHT890 PFHT890_HTmonitoring = hltHTmonitoring.clone( FolderName = 'HLT/JME/HT/PFHT890/', jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = True, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT890_v*"]) ) #HLT_PFHT780 PFHT780_HTmonitoring = hltHTmonitoring.clone( FolderName = 'HLT/JME/HT/PFHT780/', jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = True, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT780_v*"]) ) #HLT_PFHT680 PFHT680_HTmonitoring = hltHTmonitoring.clone( FolderName = 'HLT/JME/HT/PFHT680/', jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = True, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT680_v*"]) ) @@ -139,36 +138,42 @@ PFHT590_HTmonitoring = hltHTmonitoring.clone( FolderName = 'HLT/JME/HT/PFHT590/', jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = True, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT590_v*"]) ) #HLT_PFHT510 PFHT510_HTmonitoring = hltHTmonitoring.clone( FolderName = 'HLT/JME/HT/PFHT510/', jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = True, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT510_v*"]) ) #HLT_PFHT430 PFHT430_HTmonitoring = hltHTmonitoring.clone( FolderName = 'HLT/JME/HT/PFHT430/', jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = True, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT430_v*"]) ) #HLT_PFHT370 PFHT370_HTmonitoring = hltHTmonitoring.clone( FolderName = 'HLT/JME/HT/PFHT370/', jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = True, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT370_v*"]) ) #HLT_PFHT250 PFHT250_HTmonitoring = hltHTmonitoring.clone( FolderName = 'HLT/JME/HT/PFHT250/', jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = True, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT250_v*"]) ) #HLT_PFHT180 PFHT180_HTmonitoring = hltHTmonitoring.clone( FolderName = 'HLT/JME/HT/PFHT180/', jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = True, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT180_v*"]) ) # HLT_PFMETTypeOne110_PFMHT110_IDTight @@ -176,6 +181,7 @@ FolderName = 'HLT/EXO/MET/PFMETTypeOne110/', jetSelection = "pt > 100 && eta < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETTypeOne110_PFMHT110_IDTight_v*"]) ) # HLT_PFMETTypeOne120_PFMHT120_IDTight @@ -183,6 +189,7 @@ FolderName = 'HLT/EXO/MET/PFMETTypeOne120/', jetSelection = "pt > 100 && eta < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PPFMETTypeOne120_PFMHT120_IDTight_v*"]) ) # HLT_PFMETTypeOne130_PFMHT130_IDTight @@ -190,6 +197,7 @@ FolderName = 'HLT/EXO/MET/PFMETTypeOne130/', jetSelection = "pt > 100 && eta < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETTypeOne130_PFMHT130_IDTight_v*"]) ) # HLT_PFMETTypeOne140_PFMHT140_IDTight @@ -197,6 +205,7 @@ FolderName = 'HLT/EXO/MET/PFMETTypeOne140/', jetSelection = "pt > 100 && eta < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1", jetSelection_HT = "pt > 30 && eta < 2.5", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETTypeOne140_PFMHT140_IDTight_v*"]) ) # HLT_PFMET120_PFMHT120_IDTight_PFHT60 @@ -204,6 +213,7 @@ FolderName = 'HLT/EXO/MET/PFMET120_PFHT60/', metSelection = "pt > 220", jetSelection = "pt > 70 && abs(eta) < 2.4 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET120_PFMHT120_IDTight_PFHT60_v*"]) ) # HLT_PFMETNoMu120_PFMHTNoMu120_IDTight_PFHT60 @@ -211,6 +221,7 @@ FolderName = 'HLT/EXO/MET/PFMETNoMu120_PFHT60/', metSelection = "pt > 220", jetSelection = "pt > 70 && abs(eta) < 2.4 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETNoMu120_PFMHTNoMu120_IDTight_PFHT60_v*"]) ) # HLT_PFMETTypeOne120_PFMHT120_IDTight_PFHT60 @@ -218,14 +229,12 @@ FolderName = 'HLT/EXO/MET/PFMETTypeOne120_PFHT60/', metSelection = "pt > 220", jetSelection = "pt > 70 && abs(eta) < 2.4 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1", + enableFullMonitoring = False, numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETTypeOne120_PFMHT120_IDTight_PFHT60_v*"]) ) exoHLTHTmonitoring = cms.Sequence( - PFMETNoMu90_PFMHTNoMu90_HTmonitoring - + PFMETNoMu120_PFMHTNoMu120_HTmonitoring - + MET200_HTmonitoring + PFMETNoMu120_PFMHTNoMu120_HTmonitoring + MonoCentralPFJet80_PFMETNoMu120_PFMHTNoMu120_HTmonitoring - + MonoCentralPFJet80_PFMETNoMu90_PFMHTNoMu90_HTmonitoring + PFHT350MinPFJet15_HTmonitoring + PFHT500_PFMET100_PFMHT100_HTmonitoring + PFHT500_PFMET110_PFMHT110_HTmonitoring diff --git a/DQMOffline/Trigger/python/HTMonitoring_Client_cff.py b/DQMOffline/Trigger/python/HTMonitoring_Client_cff.py index 07890344e485c..01937d2918ec7 100644 --- a/DQMOffline/Trigger/python/HTMonitoring_Client_cff.py +++ b/DQMOffline/Trigger/python/HTMonitoring_Client_cff.py @@ -6,13 +6,20 @@ verbose = cms.untracked.uint32(0), # Set to 2 for all messages resolution = cms.vstring(), efficiency = cms.vstring( - "effic_ht 'HT turnON; PF HT [GeV]; efficiency' ht_numerator ht_denominator", - "effic_ht_variable 'HT turnON; PF HT [GeV]; efficiency' ht_variable_numerator ht_variable_denominator", - "effic_deltaphimetj1 'DELTAPHI turnON; DELTA PHI (PFMET, PFJET1); efficiency' deltaphimetj1_numerator deltaphimetj1_denominator", - "effic_deltaphij1j2 'DELTAPHI turnON; DELTA PHI (PFJET1, PFJET2); efficiency' deltaphij1j2_numerator deltaphij1j2_denominator" + "effic_ht 'HT turnON; PF HT [GeV]; efficiency' ht_numerator ht_denominator", + "effic_ht_variable 'HT turnON; PF HT [GeV]; efficiency' ht_variable_numerator ht_variable_denominator", + "effic_METPhi 'METPhi turnON; MET #phi; efficiency' METPhi_numerator METPhi_denominator", + "effic_jetPt1 'JET_PT1 turnON; LEADING JET PT [GeV]; efficiency' jetPt1_numerator jetPt1_denominator", + "effic_jetPt2 'JET_PT2 turnON; SUBLEADING JET PT [GeV]; efficiency' jetPt2_numerator jetPt2_denominator", + "effic_jetEta1 'JET_ETA1 turnON; LEADING JET #eta; efficiency' jetEta1_numerator jetEta1_denominator", + "effic_jetEta2 'JET_ETA2 turnON; SUBLEADING JET #eta; efficiency' jetEta2_numerator jetEta2_denominator", + "effic_jetPhi1 'JET_PHI1 turnON; LEADING JET #phi; efficiency' jetPhi1_numerator jetPhi1_denominator", + "effic_jetPhi2 'JET_PHI2 turnON; SUBLEADING JET #phi; efficiency' jetPhi2_numerator jetPhi2_denominator", + "effic_nJets 'nJets; nJets; efficiency' nJets_numerator nJets_denominator", + "effic_nJets_HT 'nJets (Pt>30 && eta<2.5); nJets_; efficiency' nJetsHT_numerator nJetsHT_denominator" ), efficiencyProfile = cms.untracked.vstring( - "effic_ht_vs_LS 'HT efficiency vs LS; LS; PF HT efficiency' htVsLS_numerator htVsLS_denominator" + "effic_ht_vs_LS 'HT efficiency vs LS; LS; PF HT efficiency' htVsLS_numerator htVsLS_denominator" ), ) diff --git a/DQMOffline/Trigger/python/JetMETPromptMonitoring_Client_cff.py b/DQMOffline/Trigger/python/JetMETPromptMonitoring_Client_cff.py index 4838387fc16c7..57719d61878c4 100644 --- a/DQMOffline/Trigger/python/JetMETPromptMonitoring_Client_cff.py +++ b/DQMOffline/Trigger/python/JetMETPromptMonitoring_Client_cff.py @@ -14,39 +14,22 @@ "effic_pfjetpT_HB_pTThresh 'Jet pT turnON (HB); PFJet(pT) [GeV]; efficiency' pfjetpT_pTThresh_HB_numerator pfjetpT_pTThresh_HB_denominator", "effic_pfjetphi_HB 'Jet efficiency vs #phi (HB); PF Jet #phi [rad]; efficiency' pfjetphi_HB_numerator pfjetphi_HB_denominator", "effic_pfjeteta_HB 'Jet efficiency vs #eta (HB); PF Jet #eta; efficiency' pfjeteta_HB_numerator pfjeteta_HB_denominator", - + ## HE "effic_pfjetpT_HE 'Jet pT turnON (HE); PFJet(pT) [GeV]; efficiency' pfjetpT_HE_numerator pfjetpT_HE_denominator", "effic_pfjetpT_HE_pTThresh 'Jet pT turnON (HE); PFJet(pT) [GeV]; efficiency' pfjetpT_pTThresh_HE_numerator pfjetpT_pTThresh_HE_denominator", "effic_pfjetphi_HE 'Jet efficiency vs #phi (HE); PF Jet #phi [rad]; efficiency' pfjetphi_HE_numerator pfjetphi_HE_denominator", "effic_pfjeteta_HE 'Jet efficiency vs #eta (HE); PF Jet #eta; efficiency' pfjeteta_HE_numerator pfjeteta_HE_denominator", - + ## HE_p "effic_pfjetpT_HE_p 'Jet pT turnON (HEP); PFJet(pT) [GeV]; efficiency' pfjetpT_HE_p_numerator pfjetpT_HE_p_denominator", "effic_pfjetpT_HE_p_pTThresh 'Jet pT turnON (HEP); PFJet(pT) [GeV]; efficiency' pfjetpT_pTThresh_HE_p_numerator pfjetpT_pTThresh_HE_p_denominator", "effic_pfjetphi_HE_p 'Jet efficiency vs #phi (HEP); PF Jet #phi [rad]; efficiency' pfjetphi_HE_p_numerator pfjetphi_HE_p_denominator", "effic_pfjeteta_HE_p 'Jet efficiency vs #eta (HEP); PF Jet #eta; efficiency' pfjeteta_HE_p_numerator pfjeteta_HE_p_denominator", - + ## HE_m "effic_pfjetpT_HE_m 'Jet pT turnON (HEM); PFJet(pT) [GeV]; efficiency' pfjetpT_HE_m_numerator pfjetpT_HE_m_denominator", "effic_pfjetpT_HE_m_pTThresh 'Jet pT turnON (HEM); PFJet(pT) [GeV]; efficiency' pfjetpT_pTThresh_HE_m_numerator pfjetpT_pTThresh_HE_m_denominator", "effic_pfjetphi_HE_m 'Jet efficiency vs #phi (HEM); PF Jet #phi [rad]; efficiency' pfjetphi_HE_m_numerator pfjetphi_HE_m_denominator", - "effic_pfjeteta_HE_m 'Jet efficiency vs #eta (HEM); PF Jet #eta; efficiency' pfjeteta_HE_m_numerator pfjeteta_HE_m_denominator", - - "effic_pfjetpT_HEP17 'Jet pT turnON (HEP17); PFJet(pT) [GeV]; efficiency' pfjetpT_HEP17_numerator pfjetpT_HEP17_denominator", - "effic_pfjetpT_HEP17_pTThresh 'Jet pT turnON (HEP17); PFJet(pT) [GeV]; efficiency' pfjetpT_pTThresh_HEP17_numerator pfjetpT_pTThresh_HEP17_denominator", - "effic_pfjetphi_HEP17 'Jet efficiency vs #phi (HEP17); PF Jet #phi [rad]; efficiency' pfjetphi_HEP17_numerator pfjetphi_HEP17_denominator", -# "effic_pfjeteta_HEP17 'Jet efficiency vs #eta (HEP17); PF Jet #eta; efficiency' pfjeteta_HEP17_numerator pfjeteta_HEP17_denominator", - "effic_pfjetabseta_HEP17 'Jet efficiency vs |#eta| (HEP17); PF Jet |#eta|; efficiency' pfjetabseta_HEP17_numerator pfjetabseta_HEP17_denominator", - - "effic_pfjetpT_HEM17 'Jet pT turnON (HEM17); PFJet(pT) [GeV]; efficiency' pfjetpT_HEM17_numerator pfjetpT_HEM17_denominator", - "effic_pfjetpT_HEM17_pTThresh 'Jet pT turnON (HEM17); PFJet(pT) [GeV]; efficiency' pfjetpT_pTThresh_HEM17_numerator pfjetpT_pTThresh_HEM17_denominator", - "effic_pfjetphi_HEM17 'Jet efficiency vs #phi (HEM17); PF Jet #phi [rad]; efficiency' pfjetphi_HEM17_numerator pfjetphi_HEM17_denominator", -# "effic_pfjeteta_HEM17 'Jet efficiency vs #eta (HEM17); PF Jet #eta; efficiency' pfjeteta_HEM17_numerator pfjeteta_HEM17_denominator", - "effic_pfjetabseta_HEM17 'Jet efficiency vs |#eta| (HEM17); PF Jet |#eta|; efficiency' pfjetabseta_HEM17_numerator pfjetabseta_HEM17_denominator", - - "effic_pfjetpT_HEP18 'Jet pT turnON (HEP18); PFJet(pT) [GeV]; efficiency' pfjetpT_HEP18_numerator pfjetpT_HEP18_denominator", - "effic_pfjetpT_HEP18_pTThresh 'Jet pT turnON (HEP18); PFJet(pT) [GeV]; efficiency' pfjetpT_pTThresh_HEP18_numerator pfjetpT_pTThresh_HEP18_denominator", -# "effic_pfjetphi_HEP18 'Jet efficiency vs #phi (HEP18); PF Jet #phi [rad]; efficiency' pfjetphi_HEP18_numerator pfjetphi_HEP18_denominator", -# "effic_pfjeteta_HEP18 'Jet efficiency vs #eta (HEP18); PF Jet #eta; efficiency' pfjeteta_HEP18_numerator pfjeteta_HEP18_denominator", - + "effic_pfjeteta_HE_m 'Jet efficiency vs #eta (HEM); PF Jet #eta; efficiency' pfjeteta_HE_m_numerator pfjeteta_HE_m_denominator", + ## HF "effic_pfjetpT_HF 'Jet pT turnON (HF); PFJet(pT) [GeV]; efficiency' pfjetpT_HF_numerator pfjetpT_HF_denominator", "effic_pfjetpT_HF_pTThresh 'Jet pT turnON (HF); PFJet(pT) [GeV]; efficiency' pfjetpT_pTThresh_HF_numerator pfjetpT_pTThresh_HF_denominator", "effic_pfjetphi_HF 'Jet efficiency vs #phi (HF); PF Jet #phi [rad]; efficiency' pfjetphi_HF_numerator pfjetphi_HF_denominator", @@ -56,21 +39,23 @@ "effic_pfjetEtaVsPhi_HB 'Jet efficiency vs #eta and #phi(HB); PF Jet #eta; #phi' pfjetEtaVsPhi_HB_numerator pfjetEtaVsPhi_HB_denominator", "effic_pfjetEtaVsPhi_HE 'Jet efficiency vs #eta and #phi(HE); PF Jet #eta; #phi' pfjetEtaVsPhi_HE_numerator pfjetEtaVsPhi_HE_denominator", "effic_pfjetEtaVsPhi_HF 'Jet efficiency vs #eta and #phi(HF); PF Jet #eta; #phi' pfjetEtaVsPhi_HF_numerator pfjetEtaVsPhi_HF_denominator", -# "effic_pfjetEtaVsPhi_HEM17 'Jet efficiency vs #eta and #phi(HEM17); PF Jet #eta; #phi' pfjetEtaVsPhi_HEM17_numerator pfjetEtaVsPhi_HEM17_denominator", -# "effic_pfjetEtaVsPhi_HEP17 'Jet efficiency vs #eta and #phi(HEP17); PF Jet #eta; #phi' pfjetEtaVsPhi_HEP17_numerator pfjetEtaVsPhi_HEP17_denominator", -# "effic_pfjetEtaVsPhi_HEP18 'Jet efficiency vs #eta and #phi(HEP18); PF Jet #eta; #phi' pfjetEtaVsPhi_HEP18_numerator pfjetEtaVsPhi_HEP18_denominator", + "effic_pfjetEtaVsPhi_HE_p 'Jet efficiency vs #eta and #phi(HE_p); PF Jet #eta; #phi' pfjetEtaVsPhi_HE_p_numerator pfjetEtaVsPhi_HE_p_denominator", "effic_pfjetEtaVsPhi_HE_m 'Jet efficiency vs #eta and #phi(HE_m); PF Jet #eta; #phi' pfjetEtaVsPhi_HE_m_numerator pfjetEtaVsPhi_HE_m_denominator", -# "effic_pfjetAbsEtaVsPhi_HEM17 'Jet efficiency vs |#eta| and #phi(HEM17); PF Jet |#eta|; #phi' pfjetAbsEtaVsPhi_HEM17_numerator pfjetAbsEtaVsPhi_HEM17_denominator", -# "effic_pfjetAbsEtaVsPhi_HEP17 'Jet efficiency vs |#eta| and #phi(HEP17); PF Jet |#eta|; #phi' pfjetAbsEtaVsPhi_HEP17_numerator pfjetAbsEtaVsPhi_HEP17_denominator", + + "effic_pfjetEtaVspT 'Jet efficiency #eta vs Pt; PF Jet #eta; Pt' pfjetEtaVspT_numerator pfjetEtaVspT_denominator", + "effic_pfjetEtaVspT_HB 'Jet efficiency #eta vs Pt(HB); PF Jet #eta; Pt' pfjetEtaVspT_HB_numerator pfjetEtaVspT_HB_denominator", + "effic_pfjetEtaVspT_HE 'Jet efficiency #eta vs Pt(HE); PF Jet #eta; Pt' pfjetEtaVspT_HE_numerator pfjetEtaVspT_HE_denominator", + "effic_pfjetEtaVspT_HF 'Jet efficiency #eta vs Pt(HF); PF Jet #eta; Pt' pfjetEtaVspT_HF_numerator pfjetEtaVspT_HF_denominator", + + "effic_pfjetEtaVspT_HE_p 'Jet efficiency #eta vs Pt(HE_p); PF Jet #eta; Pt' pfjetEtaVspT_HE_p_numerator pfjetEtaVspT_HE_p_denominator", + "effic_pfjetEtaVspT_HE_m 'Jet efficiency #eta vs Pt(HE_m); PF Jet #eta; Pt' pfjetEtaVspT_HE_m_numerator pfjetEtaVspT_HE_m_denominator" ), efficiencyProfile = cms.untracked.vstring( "effic_pfjetpT_vs_LS 'JET efficiency vs LS; LS; PF JET efficiency' pfjetpTVsLS_numerator pfjetpTVsLS_denominator", # "effic_pfjetpT_HBvs_LS 'JET efficiency vs LS; LS; PF JET efficiency' pfjetpTVsLS_HB_numerator pfjetpTVsLS_HB_denominator", # "effic_pfjetpT_HEvs_LS 'JET efficiency vs LS; LS; PF JET efficiency' pfjetpTVsLS_HE_numerator pfjetpTVsLS_HE_denominator", # "effic_pfjetpT_HFvs_LS 'JET efficiency vs LS; LS; PF JET efficiency' pfjetpTVsLS_HF_numerator pfjetpTVsLS_HF_denominator", -# "effic_pfjetpT_HEP17vs_LS 'JET efficiency vs LS; LS; PF JET efficiency' pfjetpTVsLS_HEP17_numerator pfjetpTVsLS_HEP17_denominator", -# "effic_pfjetpT_HEM17vs_LS 'JET efficiency vs LS; LS; PF JET efficiency' pfjetpTVsLS_HEM17_numerator pfjetpTVsLS_HEM17_denominator", # "effic_pfjetpT_HE_mvs_LS 'JET efficiency vs LS; LS; PF JET efficiency' pfjetpTVsLS_HE_m_numerator pfjetpTVsLS_HE_m_denominator", # "effic_pfjetpT_HE_pvs_LS 'JET efficiency vs LS; LS; PF JET efficiency' pfjetpTVsLS_HE_p_numerator pfjetpTVsLS_HE_p_denominator", ), @@ -107,40 +92,27 @@ "effic_calojetphi_HE_m 'Jet efficiency vs #phi (HEM); Calo Jet #phi [rad]; efficiency' calojetphi_HE_m_numerator calojetphi_HE_m_denominator", "effic_calojeteta_HE_m 'Jet efficiency vs #eta (HEM); Calo Jet #eta; efficiency' calojeteta_HE_m_numerator calojeteta_HE_m_denominator", - "effic_calojetpT_HEP17 'Jet pT turnON (HEP17); CaloJet(pT) [GeV]; efficiency' calojetpT_HEP17_numerator calojetpT_HEP17_denominator", - "effic_calojetpT_HEP17_pTThresh 'Jet pT turnON (HEP17); CaloJet(pT) [GeV]; efficiency' calojetpT_pTThresh_HEP17_numerator calojetpT_pTThresh_HEP17_denominator", -# "effic_calojetphi_HEP17 'Jet efficiency vs #phi (HEP17); Calo Jet #phi [rad]; efficiency' calojetphi_HEP17_numerator calojetphi_HEP17_denominator", -# "effic_calojeteta_HEP17 'Jet efficiency vs #eta (HEP17); Calo Jet #eta; efficiency' calojeteta_HEP17_numerator calojeteta_HEP17_denominator", -# "effic_calojetabseta_HEP17 'Jet efficiency vs |#eta| (HEP17); Calo Jet |#eta|; efficiency' calojetabseta_HEP17_numerator calojetabseta_HEP17_denominator", - - "effic_calojetpT_HEP18 'Jet pT turnON (HEP18); CaloJet(pT) [GeV]; efficiency' calojetpT_HEP18_numerator calojetpT_HEP18_denominator", - "effic_calojetpT_HEP18_pTThresh 'Jet pT turnON (HEP18); CaloJet(pT) [GeV]; efficiency' calojetpT_pTThresh_HEP18_numerator calojetpT_pTThresh_HEP18_denominator", -# "effic_calojetphi_HEP18 'Jet efficiency vs #phi (HEP18); Calo Jet phi [rad]; efficiency' calojetphi_HEP18_numerator calojetphi_HEP18_denominator", -# "effic_calojeteta_HEP18 'Jet efficiency vs #eta (HEP18); Calo Jet #eta; efficiency' calojeteta_HEP18_numerator calojeteta_HEP18_denominator", - - "effic_calojetpT_HEM17 'Jet pT turnON (HEM17); CaloJet(pT) [GeV]; efficiency' calojetpT_HEM17_numerator calojetpT_HEM17_denominator", - "effic_calojetpT_HEM17_pTThresh 'Jet pT turnON (HEM17); CaloJet(pT) [GeV]; efficiency' calojetpT_pTThresh_HEM17_numerator calojetpT_pTThresh_HEM17_denominator", -# "effic_calojetphi_HEM17 'Jet efficiency vs #phi (HEM17); Calo Jet #phi; efficiency' calojetphi_HEM17_numerator calojetphi_HEM17_denominator", -# "effic_calojeteta_HEM17 'Jet efficiency vs #eta (HEM17); Calo Jet #eta; efficiency' calojeteta_HEM17_numerator calojeteta_HEM17_denominator", -# "effic_calojetabseta_HEM17 'Jet efficiency vs #eta (HEM17); Calo Jet |#eta|; efficiency' calojetabseta_HEM17_numerator calojetabseta_HEM17_denominator", - "effic_calojetpT_HF 'Jet pT turnON; CaloJet(pT) [GeV]; efficiency' calojetpT_HF_numerator calojetpT_HF_denominator", "effic_calojetpT_HF_pTThresh 'Jet pT turnON; CaloJet(pT) [GeV]; efficiency' calojetpT_pTThresh_HF_numerator calojetpT_pTThresh_HF_denominator", "effic_calojetphi_HF 'Jet efficiency vs #phi; Calo Jet #phi [rad]; efficiency' calojetphi_HF_numerator calojetphi_HF_denominator", "effic_calojeteta_HF 'Jet efficiency vs #eta; Calo Jet #eta; efficiency' calojeteta_HF_numerator calojeteta_HF_denominator", ## 2D Eff - "effic_calojetEtaVsPhi 'Jet efficiency vs #eta and #phi; PF Jet #eta; #phi' calojetEtaVsPhi_numerator calojetEtaVsPhi_denominator", - "effic_calojetEtaVsPhi_HB 'Jet efficiency vs #eta and #phi(HB); PF Jet #eta; #phi' calojetEtaVsPhi_HB_numerator calojetEtaVsPhi_HB_denominator", - "effic_calojetEtaVsPhi_HE 'Jet efficiency vs #eta and #phi(HE); PF Jet #eta; #phi' calojetEtaVsPhi_HE_numerator calojetEtaVsPhi_HE_denominator", - "effic_calojetEtaVsPhi_HF 'Jet efficiency vs #eta and #phi(HF); PF Jet #eta; #phi' calojetEtaVsPhi_HF_numerator calojetEtaVsPhi_HF_denominator", -# "effic_calojetEtaVsPhi_HEM17 'Jet efficiency vs #eta and #phi(HEM17); PF Jet #eta; #phi' calojetEtaVsPhi_HEM17_numerator calojetEtaVsPhi_HEM17_denominator", -# "effic_calojetEtaVsPhi_HEP17 'Jet efficiency vs #eta and #phi(HEP17); PF Jet #eta; #phi' calojetEtaVsPhi_HEP17_numerator calojetEtaVsPhi_HEP17_denominator", -# "effic_calojetEtaVsPhi_HEP18 'Jet efficiency vs #eta and #phi(HEP18); PF Jet #eta; #phi' calojetEtaVsPhi_HEP18_numerator calojetEtaVsPhi_HEP18_denominator", - "effic_calojetEtaVsPhi_HE_p 'Jet efficiency vs #eta and #phi(HE_p); PF Jet #eta; #phi' calojetEtaVsPhi_HE_p_numerator calojetEtaVsPhi_HE_p_denominator", - "effic_calojetEtaVsPhi_HE_m 'Jet efficiency vs #eta and #phi(HE_m); PF Jet #eta; #phi' calojetEtaVsPhi_HE_m_numerator calojetEtaVsPhi_HE_m_denominator", -# "effic_calojetAbsEtaVsPhi_HEM17 'Jet efficiency vs |#eta| and #phi(HEM17); PF Jet |#eta|; #phi' calojetAbsEtaVsPhi_HEM17_numerator calojetAbsEtaVsPhi_HEM17_denominator", -# "effic_calojetAbsEtaVsPhi_HEP17 'Jet efficiency vs |#eta| and #phi(HEP17); PF Jet |#eta|; #phi' calojetAbsEtaVsPhi_HEP17_numerator calojetAbsEtaVsPhi_HEP17_denominator", + "effic_calojetEtaVsPhi 'Jet efficiency vs #eta and #phi; Calo Jet #eta; #phi' calojetEtaVsPhi_numerator calojetEtaVsPhi_denominator", + "effic_calojetEtaVsPhi_HB 'Jet efficiency vs #eta and #phi(HB); Calo Jet #eta; #phi' calojetEtaVsPhi_HB_numerator calojetEtaVsPhi_HB_denominator", + "effic_calojetEtaVsPhi_HE 'Jet efficiency vs #eta and #phi(HE); Calo Jet #eta; #phi' calojetEtaVsPhi_HE_numerator calojetEtaVsPhi_HE_denominator", + "effic_calojetEtaVsPhi_HF 'Jet efficiency vs #eta and #phi(HF); Calo Jet #eta; #phi' calojetEtaVsPhi_HF_numerator calojetEtaVsPhi_HF_denominator", + "effic_calojetEtaVsPhi_HE_p 'Jet efficiency vs #eta and #phi(HE_p); Calo Jet #eta; #phi' calojetEtaVsPhi_HE_p_numerator calojetEtaVsPhi_HE_p_denominator", + "effic_calojetEtaVsPhi_HE_m 'Jet efficiency vs #eta and #phi(HE_m); Calo Jet #eta; #phi' calojetEtaVsPhi_HE_m_numerator calojetEtaVsPhi_HE_m_denominator", + + "effic_calojetEtaVspT 'Jet efficiency #eta vs Pt; Calo Jet #eta; Pt' calojetEtaVspT_numerator calojetEtaVspT_denominator", + "effic_calojetEtaVspT_HB 'Jet efficiency #eta vs Pt(HB); Calo Jet #eta; Pt' calojetEtaVspT_HB_numerator calojetEtaVspT_HB_denominator", + "effic_calojetEtaVspT_HE 'Jet efficiency #eta vs Pt(HE); Calo Jet #eta; Pt' calojetEtaVspT_HE_numerator calojetEtaVspT_HE_denominator", + "effic_calojetEtaVspT_HF 'Jet efficiency #eta vs Pt(HF); Calo Jet #eta; Pt' calojetEtaVspT_HF_numerator calojetEtaVspT_HF_denominator", + + "effic_calojetEtaVspT_HE_p 'Jet efficiency #eta vs Pt(HE_p); Calo Jet #eta; Pt' calojetEtaVspT_HE_p_numerator calojetEtaVspT_HE_p_denominator", + "effic_calojetEtaVspT_HE_m 'Jet efficiency #eta vs Pt(HE_m); Calo Jet #eta; Pt' calojetEtaVspT_HE_m_numerator calojetEtaVspT_HE_m_denominator" + ), efficiencyProfile = cms.untracked.vstring( "effic_calojetpT_vs_LS 'JET efficiency vs LS; LS; Calo Jet efficiency' calojetpTVsLS_numerator calojetpTVsLS_denominator", @@ -148,46 +120,15 @@ # "effic_calojetpT_HBvs_LS 'JET efficiency vs LS; LS; Calo JET efficiency' calojetpTVsLS_HB_numerator calojetpTVsLS_HB_denominator", # "effic_calojetpT_HEvs_LS 'JET efficiency vs LS; LS; Calo JET efficiency' calojetpTVsLS_HE_numerator calojetpTVsLS_HE_denominator", # "effic_calojetpT_HFvs_LS 'JET efficiency vs LS; LS; Calo JET efficiency' calojetpTVsLS_HF_numerator calojetpTVsLS_HF_denominator", -# "effic_calojetpT_HEP17vs_LS 'JET efficiency vs LS; LS; Calo JET efficiency' calojetpTVsLS_HEP17_numerator calojetpTVsLS_HEP17_denominator", -# "effic_calojetpT_HEM17vs_LS 'JET efficiency vs LS; LS; Calo JET efficiency' calojetpTVsLS_HEM17_numerator calojetpTVsLS_HEM17_denominator", # "effic_calojetpT_HE_mvs_LS 'JET efficiency vs LS; LS; Calo JET efficiency' calojetpTVsLS_HE_m_numerator calojetpTVsLS_HE_m_denominator", # "effic_calojetpT_HE_pvs_LS 'JET efficiency vs LS; LS; Calo JET efficiency' calojetpTVsLS_HE_p_numerator calojetpTVsLS_HE_p_denominator", ), ) -pfjetRatio = DQMEDHarvester("DQMGenericClient", - subDirs = cms.untracked.vstring("HLT/JME/*"), - verbose = cms.untracked.uint32(0), # Set to 2 for all messages - resolution = cms.vstring(), - efficiency = cms.vstring( - "ratio_pfjetpT_HEP17VSHEM17 'HEP17/HEM17 vs pT; PFJet(pT) [GeV]; Ratio' effic_pfjetpT_HEP17 effic_pfjetpT_HEM17 simpleratio" , - "ratio_pfjetpT_pTTresh_HEP17VSHEM17 'HEP17/HEM17 vs pT; PFJet(pT) [GeV]; Ratio' effic_pfjetpT_HEP17_pTThresh effic_pfjetpT_HEM17_pTThresh simpleratio" , - "ratio_pfjetphi_HEP17VSHEM17 'HEP17/HEM17 vs #phi; PFJet #phi [GeV]; Ratio' effic_pfjetphi_HEP17 effic_pfjetphi_HEM17 simpleratio", - "ratio_pfjeteta_HEP17VSHEM17 'HEP17/HEM17 vs |#eta|; PFJet(|#eta|) ; Ratio' effic_pfjetabseta_HEP17 effic_pfjetabseta_HEM17 simpleratio", - "ratio_pfjetpT_HEP17VSHEP18 'HEP17/HEP18 vs pT; PFJet(pT) [GeV]; Ratio' effic_pfjetpT_HEP17 effic_pfjetpT_HEP18 simpleratio" , - "ratio_pfjetpT_pTTresh_HEP17VSHEP18 'HEP17/HEP18 vs pT; PFJet(pT) [GeV]; Ratio' effic_pfjetpT_HEP17_pTThresh effic_pfjetpT_HEP18_pTThresh simpleratio" , - ) -) -calojetRatio = DQMEDHarvester("DQMGenericClient", - subDirs = cms.untracked.vstring("HLT/JME/*"), - verbose = cms.untracked.uint32(0), # Set to 2 for all messages - resolution = cms.vstring(), - efficiency = cms.vstring( - "ratio_calojetpT_HEP17VSHEM17 'HEP17/HEM17 vs pT; CaloJet(pT) [GeV]; Ratio' effic_calojetpT_HEP17 effic_calojetpT_HEM17 simpleratio" , - "ratio_calojetpT_pTTresh_HEP17VSHEM17 'HEP17/HEM17 vs pT; CaloJet(pT) [GeV]; Ratio' effic_calojetpT_HEP17_pTThresh effic_calojetpT_HEM17_pTThresh simpleratio" , - "ratio_calojetphi_HEP17VSHEM17 'HEP17/HEM17 vs #phi; CaloJet #phi [GeV]; Ratio' effic_calojetphi_HEP17 effic_calojetphi_HEM17 simpleratio", - "ratio_calojeteta_HEP17VSHEM17 'HEP17/HEM17 vs |#eta|; CaloJet(|#eta|) ; Ratio' effic_calojetabseta_HEP17 effic_calojetabseta_HEM17 simpleratio", - - "ratio_calojetpT_HEP17VSHEP18 'HEP17/HEP18 vs pT; CaloJet(pT) [GeV]; Ratio' effic_calojetpT_HEP17 effic_calojetpT_HEP18 simpleratio" , - "ratio_calojetpT_pTTresh_HEP17VSHEP18 'HEP17/HEP18 vs pT; CaloJet(pT) [GeV]; Ratio' effic_calojetpT_HEP17_pTThresh effic_calojetpT_HEP18_pTThresh simpleratio" , - ) -) JetMetPromClient = cms.Sequence( pfjetEfficiency *calojetEfficiency - *pfjetRatio - *calojetRatio ) diff --git a/DQMOffline/Trigger/python/METMonitor_cff.py b/DQMOffline/Trigger/python/METMonitor_cff.py index fcdedacfb5ccc..54f23e2f15aba 100644 --- a/DQMOffline/Trigger/python/METMonitor_cff.py +++ b/DQMOffline/Trigger/python/METMonitor_cff.py @@ -6,270 +6,238 @@ PFMET110_PFMHT110_IDTight_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/PFMET110/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET110_PFMHT110_IDTight_v*"]), + enableFullMonitoring = True, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) #HLT_PFMET120_PFMHT120_IDTight PFMET120_PFMHT120_IDTight_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/PFMET120', + FolderName = 'HLT/JME/MET/PFMET120/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET120_PFMHT120_IDTight_v*"]), + enableFullMonitoring = True, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_PFMET130_PFMHT130_IDTight PFMET130_PFMHT130_IDTight_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/PFMET130/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET130_PFMHT130_IDTight_v*"]), + enableFullMonitoring = True, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_PFMET140_PFMHT140_IDTight PFMET140_PFMHT140_IDTight_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/PFMET140/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET140_PFMHT140_IDTight_v*"]), + enableFullMonitoring = True, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_PFMETTypeOne110_PFMHT110_IDTight PFMETTypeOne110_PFMHT110_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/PFMETTypeOne110/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETTypeOne110_PFMHT110_IDTight_v*"]), + enableFullMonitoring = True, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_PFMETTypeOne120_PFMHT120_IDTight PFMETTypeOne120_PFMHT120_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/PFMETTypeOne120/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PPFMETTypeOne120_PFMHT120_IDTight_v*"]), + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETTypeOne120_PFMHT120_IDTight_v*"]), + enableFullMonitoring = True, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_PFMETTypeOne130_PFMHT130_IDTight PFMETTypeOne130_PFMHT130_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/PFMETTypeOne130/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETTypeOne130_PFMHT130_IDTight_v*"]), + enableFullMonitoring = True, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_PFMETTypeOne140_PFMHT140_IDTight PFMETTypeOne140_PFMHT140_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/PFMETTypeOne140/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETTypeOne140_PFMHT140_IDTight_v*"]), - jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" -) -# HLT_PFMETNoMu90_PFMHTNoMu90_IDTight -PFMETNoMu90_PFMHTNoMu90_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/PFMETNoMu90/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETNoMu90_PFMHTNoMu90_IDTight_v*"]), + enableFullMonitoring = True, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_PFMETNoMu110_PFMHTNoMu110_IDTight PFMETNoMu110_PFMHTNoMu110_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/PFMETNoMu110/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETNoMu110_PFMHTNoMu110_IDTight_v*"]), + enableFullMonitoring = True, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_PFMETNoMu120_PFMHTNoMu120_IDTight PFMETNoMu120_PFMHTNoMu120_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/PFMETNoMu120/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETNoMu120_PFMHTNoMu120_IDTight_v*"]), + enableFullMonitoring = True, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_PFMETNoMu130_PFMHTNoMu130_IDTight PFMETNoMu130_PFMHTNoMu130_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/PFMETNoMu130/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETNoMu130_PFMHTNoMu130_IDTight_v*"]), + enableFullMonitoring = True, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_PFMETNoMu140_PFMHTNoMu140_IDTight PFMETNoMu140_PFMHTNoMu140_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/PFMETNoMu140/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETNoMu140_PFMHTNoMu140_IDTight_v*"]), + enableFullMonitoring = True, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) -# HLT_MET200 -MET200_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/MET200/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_MET200_v*"]) -) # HLT_MonoCentralPFJet80_PFMETNoMu110_PFMHTNoMu110_IDTight MonoCentralPFJet80_PFMETNoMu110_PFMHTNoMu110_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/EXO/MET/MonoCentralPFJet80_PFMETNoMu110/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_MonoCentralPFJet80_PFMETNoMu110_PFMHTNoMu110_IDTight"]), + enableFullMonitoring = False, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_MonoCentralPFJet80_PFMETNoMu120_PFMHTNoMu120_IDTight MonoCentralPFJet80_PFMETNoMu120_PFMHTNoMu120_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/EXO/MET/MonoCentralPFJet80_PFMETNoMu120/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_MonoCentralPFJet80_PFMETNoMu120_PFMHTNoMu120_IDTight_v*"]), + enableFullMonitoring = False, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_MonoCentralPFJet80_PFMETNoMu130_PFMHTNoMu130_IDTight MonoCentralPFJet80_PFMETNoMu130_PFMHTNoMu130_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/EXO/MET/MonoCentralPFJet80_PFMETNoMu130/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_MonoCentralPFJet80_PFMETNoMu130_PFMHTNoMu130_IDTight_v*"]), + enableFullMonitoring = False, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_MonoCentralPFJet80_PFMETNoMu140_PFMHTNoMu140_IDTight MonoCentralPFJet80_PFMETNoMu140_PFMHTNoMu140_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/EXO/MET/MonoCentralPFJet80_PFMETNoMu140/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_MonoCentralPFJet80_PFMETNoMu140_PFMHTNoMu140_IDTight_v*"]), - jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" -) -# HLT_MonoCentralPFJet80_PFMETNoMu90_PFMHTNoMu90_IDTight -MonoCentralPFJet80_PFMETNoMu90_PFMHTNoMu90_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/EXO/MET/MonoCentralPFJet80_PFMETNoMu90/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_MonoCentralPFJet80_PFMETNoMu110_PFMHTNoMu110_IDTight_v*"]), + enableFullMonitoring = False, jetSelection = "pt > 100 && abs(eta) < 2.5 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_PFHT500_PFMET100_PFMHT100_IDTight PFHT500_PFMET100_PFMHT100_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/EXO/MET/PFHT500_PFMET100_PFMHT100/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT500_PFMET100_PFMHT100_IDTight_v*"]) + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT500_PFMET100_PFMHT100_IDTight_v*"]), + enableFullMonitoring = False ) # HLT_PFHT500_PFMET110_PFMHT110_IDTight PFHT500_PFMET110_PFMHT110_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/EXO/MET/PFHT500_PFMET110_PFMHT110/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT500_PFMET110_PFMHT110_IDTight_v*"]) + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT500_PFMET110_PFMHT110_IDTight_v*"]), + enableFullMonitoring = False, ) # HLT_PFHT700_PFMET85_PFMHT85_IDTight PFHT700_PFMET85_PFMHT85_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/EXO/MET/PFHT700_PFMET85_PFMHT85/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT700_PFMET85_PFMHT85_IDTight_v*"]) + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT700_PFMET85_PFMHT85_IDTight_v*"]), + enableFullMonitoring = False, ) # HLT_PFHT700_PFMET95_PFMHT95_IDTight PFHT700_PFMET95_PFMHT95_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/EXO/MET/PFHT700_PFMET95_PFMHT95/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT700_PFMET95_PFMHT95_IDTight_v*"]) + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT700_PFMET95_PFMHT95_IDTight_v*"]), + enableFullMonitoring = False, ) # HLT_PFHT800_PFMET75_PFMHT75_IDTight PFHT800_PFMET75_PFMHT75_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/EXO/MET/PFHT800_PFMET75_PFMHT75/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT800_PFMET75_PFMHT75_IDTight_v*"]) + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT800_PFMET75_PFMHT75_IDTight_v*"]), + enableFullMonitoring = False, ) # HLT_PFHT800_PFMET85_PFMHT85_IDTight PFHT800_PFMET85_PFMHT85_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/EXO/MET/PFHT800_PFMET85_PFMHT85/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT800_PFMET85_PFMHT85_IDTight_v*"]) + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT800_PFMET85_PFMHT85_IDTight_v*"]), + enableFullMonitoring = False, ) # HLT_PFMET120_PFMHT120_IDTight_PFHT60 PFMET120_PFMHT120_IDTight_PFHT60_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/EXO/MET/PFMET120_PFHT60/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET120_PFMHT120_IDTight_PFHT60_v*"]), + enableFullMonitoring = False, jetSelection = "pt > 70 && abs(eta) < 2.4 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) # HLT_PFMETNoMu120_PFMHTNoMu120_IDTight_PFHT60 PFMETNoMu120_PFMHTNoMu120_IDTight_PFHT60_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/EXO/MET/PFMETNoMu120_PFHT60/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETNoMu120_PFMHTNoMu120_IDTight_PFHT60_v*"]), + enableFullMonitoring = False, jetSelection = cms.string("pt > 70 && abs(eta) < 2.4 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1") ) # HLT_PFMETTypeOne120_PFMHT120_IDTight_PFHT60 PFMETTypeOne120_PFMHT120_IDTight_PFHT60_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/EXO/MET/PFMETTypeOne120_PFHT60/', + FolderName = 'HLT/EXO/MET/PFMETTypeOne120_PFMHT120_PFHT60/', numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETTypeOne120_PFMHT120_IDTight_PFHT60_v*"]), + enableFullMonitoring = False, jetSelection = "pt > 70 && abs(eta) < 2.4 && neutralHadronEnergyFraction < 0.8 && chargedHadronEnergyFraction > 0.1" ) ## Add Pure MET Trigger ## # HLT_L1ETMHadSeeds_v L1ETMHadSeeds_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/L1ETMHadSeeds/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_L1ETMHadSeeds_v*"]) + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_L1ETMHadSeeds_v*"]), + enableFullMonitoring = True ) # HLT_CaloMHT90_v CaloMHT90_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/CaloMHT/CaloMHT90/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMHT90_v*"]) -) -# HLT_CaloMET70_HBHECleaned_v -CaloMET70_HBHECleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/Calo/CaloMET70_HBHECleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET70_HBHECleaned_v*"]) -) -# HLT_CaloMET80_HBHECleaned_v -CaloMET80_HBHECleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/Calo/CaloMET80_HBHECleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET80_HBHECleaned_v*"]) + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMHT90_v*"]), + enableFullMonitoring = False ) # HLT_CaloMET80_NotCleaned_v CaloMET80_NotCleaned_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/Calo/CaloMET80_NotCleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET80_NotCleaned_v*"]) -) -# HLT_CaloMET90_HBHECleaned_v -CaloMET90_HBHECleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/Calo/CaloMET90_HBHECleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET90_HBHECleaned_v*"]) + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET80_NotCleaned_v*"]), + enableFullMonitoring = True ) # HLT_CaloMET90_NotCleaned_v CaloMET90_NotCleaned_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/Calo/CaloMET90_NotCleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET90_NotCleaned_v*"]) -) -# HLT_CaloMET100_HBHECleaned_v -CaloMET100_HBHECleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/Calo/CaloMET100_HBHECleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET100_HBHECleaned_v*"]) + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET90_NotCleaned_v*"]), + enableFullMonitoring = True ) # HLT_CaloMET100_NotCleaned_v CaloMET100_NotCleaned_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/Calo/CaloMET100_NotCleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET100_NotCleaned_v*"]) + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET100_NotCleaned_v*"]), + enableFullMonitoring = True ) # HLT_CaloMET110_NotCleaned_v CaloMET110_NotCleaned_METmonitoring = hltMETmonitoring.clone( FolderName = 'HLT/JME/MET/Calo/CaloMET110_NotCleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET110_NotCleaned_v*"]) -) -# HLT_CaloMET250_HBHECleaned_v -CaloMET250_HBHECleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/Calo/CaloMET250_HBHECleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET250_HBHECleaned_v*"]) -) -# HLT_CaloMET250_NotCleaned_v -CaloMET250_NotCleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/Calo/CaloMET250_NotCleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET250_NotCleaned_v*"]) -) -# HLT_CaloMET300_HBHECleaned_v -CaloMET300_HBHECleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/Calo/CaloMET300_HBHECleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET300_HBHECleaned_v*"]) -) -# HLT_CaloMET300_HBHECleaned_v -CaloMET350_HBHECleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/Calo/CaloMET350_HBHECleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET350_HBHECleaned_v*"]) -) -# HLT_CaloMET300_HBHECleaned_v -PFMET200_HBHE_BeamHaloCleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/PFMET200_HBHE_BeamHaloCleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET200_HBHE_BeamHaloCleaned_v*"]) -) -# HLT_CaloMET300_HBHECleaned_v -PFMET200_HBHECleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/PFMET200_HBHECleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET200_HBHECleaned_v*"]) -) -# HLT_CaloMET300_HBHECleaned_v -PFMET200_NotCleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/PFMET200_NotCleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET200_NotCleaned_v*"]) -) -# HLT_CaloMET300_HBHECleaned_v -PFMET250_HBHECleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/PFMET250_HBHECleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET250_HBHECleaned_v*"]) + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_CaloMET110_NotCleaned_v*"]), + enableFullMonitoring = True ) -# HLT_CaloMET300_HBHECleaned_v -PFMET300_HBHECleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/PFMET300_HBHECleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET300_HBHECleaned_v*"]) +# HLT_PFMET200_BeamHaloCleaned_v +PFMET200_BeamHaloCleaned_METmonitoring = hltMETmonitoring.clone( + FolderName = 'HLT/JME/MET/PFMET200_BeamHaloCleaned/', + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET200_BeamHaloCleaned_v*"]), + enableFullMonitoring = True ) -# HLT_CaloMET300_HBHECleaned_v -PFMETTypeOne200_HBHE_BeamHaloCleaned_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/JME/MET/PFMETTypeOne200_HBHE_BeamHaloCleaned/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETTypeOne200_HBHE_BeamHaloCleaned_v*"]) -) -# HLT_L1MET_DTClusterNoMB1S50_v -L1MET_DTClusterNoMB1S50_METmonitoring = hltMETmonitoring.clone( - FolderName = 'HLT/EXO/MET/L1MET_DTClusterNoMB1S50/', - numGenericTriggerEventPSet = dict(hltPaths = ["HLT_L1MET_DTClusterNoMB1S50_v*"]) +# HLT_PFMET200_NotCleaned_v +PFMET200_NotCleaned_METmonitoring = hltMETmonitoring.clone( + FolderName = 'HLT/JME/MET/PFMET200_NotCleaned/', + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET200_NotCleaned_v*"]), + enableFullMonitoring = True +) +# HLT_PFMET250_NotCleaned_v +PFMET250_NotCleaned_METmonitoring = hltMETmonitoring.clone( + FolderName = 'HLT/JME/MET/PFMET250_NotCleaned/', + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET250_NotCleaned_v*"]), + enableFullMonitoring = True +) +# HLT_PFMET300_NotCleaned_v +PFMET300_NotCleaned_METmonitoring = hltMETmonitoring.clone( + FolderName = 'HLT/JME/MET/PFMET300_NotCleaned/', + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMET300_NotCleaned_v*"]), + enableFullMonitoring = True +) +# HLT_PFMETTypeOne200_BeamHaloCleaned_v +PFMETTypeOne200_BeamHaloCleaned_METmonitoring = hltMETmonitoring.clone( + FolderName = 'HLT/JME/MET/PFMETTypeOne200_BeamHaloCleaned/', + numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFMETTypeOne200_BeamHaloCleaned_v*"]), + enableFullMonitoring = True ) - exoHLTMETmonitoring = cms.Sequence( PFMET110_PFMHT110_IDTight_METmonitoring + PFMET120_PFMHT120_IDTight_METmonitoring @@ -293,31 +261,18 @@ + PFHT700_PFMET95_PFMHT95_METmonitoring + PFHT800_PFMET75_PFMHT75_METmonitoring + PFHT800_PFMET85_PFMHT85_METmonitoring - + PFMETNoMu90_PFMHTNoMu90_METmonitoring - + MET200_METmonitoring - + MonoCentralPFJet80_PFMETNoMu90_PFMHTNoMu90_METmonitoring + PFMET120_PFMHT120_IDTight_PFHT60_METmonitoring + PFMETNoMu120_PFMHTNoMu120_IDTight_PFHT60_METmonitoring + PFMETTypeOne120_PFMHT120_IDTight_PFHT60_METmonitoring + L1ETMHadSeeds_METmonitoring + CaloMHT90_METmonitoring - + CaloMET70_HBHECleaned_METmonitoring - + CaloMET80_HBHECleaned_METmonitoring + CaloMET80_NotCleaned_METmonitoring - + CaloMET90_HBHECleaned_METmonitoring + CaloMET90_NotCleaned_METmonitoring - + CaloMET100_HBHECleaned_METmonitoring + CaloMET100_NotCleaned_METmonitoring + CaloMET110_NotCleaned_METmonitoring - + CaloMET250_HBHECleaned_METmonitoring - + CaloMET250_NotCleaned_METmonitoring - + CaloMET300_HBHECleaned_METmonitoring - + CaloMET350_HBHECleaned_METmonitoring - + PFMET200_HBHE_BeamHaloCleaned_METmonitoring - + PFMET200_HBHECleaned_METmonitoring + + PFMET200_BeamHaloCleaned_METmonitoring + PFMET200_NotCleaned_METmonitoring - + PFMET250_HBHECleaned_METmonitoring - + PFMET300_HBHECleaned_METmonitoring - + PFMETTypeOne200_HBHE_BeamHaloCleaned_METmonitoring - + L1MET_DTClusterNoMB1S50_METmonitoring + + PFMET250_NotCleaned_METmonitoring + + PFMET300_NotCleaned_METmonitoring + + PFMETTypeOne200_BeamHaloCleaned_METmonitoring ) diff --git a/DQMOffline/Trigger/python/METMonitor_cfi.py b/DQMOffline/Trigger/python/METMonitor_cfi.py index 312896594947e..492d4e9828d26 100644 --- a/DQMOffline/Trigger/python/METMonitor_cfi.py +++ b/DQMOffline/Trigger/python/METMonitor_cfi.py @@ -37,8 +37,8 @@ andOrDcs = False, errorReplyDcs = True, verbosityLevel = 1, - #hltPaths = ["HLT_IsoMu27_v*"], - hltPaths = []) + hltPaths = ["HLT_IsoMu27_v*"]) + # hltPaths = []) ) from Configuration.Eras.Modifier_stage2L1Trigger_cff import stage2L1Trigger diff --git a/DQMOffline/Trigger/python/METMonitoring_Client_cff.py b/DQMOffline/Trigger/python/METMonitoring_Client_cff.py index 61f88ee0f6d09..060261b859875 100644 --- a/DQMOffline/Trigger/python/METMonitoring_Client_cff.py +++ b/DQMOffline/Trigger/python/METMonitoring_Client_cff.py @@ -6,23 +6,67 @@ verbose = cms.untracked.uint32(0), # Set to 2 for all messages resolution = cms.vstring(), efficiency = cms.vstring( - "effic_met 'MET turnON; PF MET [GeV]; efficiency' met_numerator met_denominator", - "effic_met_variable 'MET turnON; PF MET [GeV]; efficiency' met_variable_numerator met_variable_denominator", - "effic_metPhi 'MET efficiency vs phi; PF MET phi [rad]; efficiency' metPhi_numerator metPhi_denominator", - "effic_ht 'HT turnON; PF HT [GeV]; efficiency' ht_numerator ht_denominator", - "effic_ht_variable 'HT turnON; PF HT [GeV]; efficiency' ht_variable_numerator ht_variable_denominator", - "effic_deltaphimetj1 'DELTAPHI turnON; DELTA PHI (PFMET, PFJET1); efficiency' deltaphimetj1_numerator deltaphimetj1_denominator", - "effic_deltaphij1j2 'DELTAPHI turnON; DELTA PHI (PFJET1, PFJET2); efficiency' deltaphij1j2_numerator deltaphij1j2_denominator" - + "effic_met 'MET turnON; PF MET [GeV]; efficiency' met_numerator met_denominator", + "effic_met_variable 'MET turnON; PF MET [GeV]; efficiency' met_variable_numerator met_variable_denominator", + "effic_metPhi 'MET efficiency vs phi; PF MET phi [rad]; efficiency' metPhi_numerator metPhi_denominator", + "effic_ht 'HT turnON; PF HT [GeV]; efficiency' ht_numerator ht_denominator", + "effic_ht_variable 'HT turnON; PF HT [GeV]; efficiency' ht_variable_numerator ht_variable_denominator", + "effic_deltaphimetj1 'DELTAPHI turnON; DELTA PHI (PFMET, PFJET1); efficiency' deltaphi_metjet1_numerator deltaphi_metjet1_denominator", + "effic_deltaphij1j2 'DELTAPHI turnON; DELTA PHI (PFJET1, PFJET2); efficiency' deltaphi_jet1jet2_numerator deltaphi_jet1jet2_denominator", + #============================================================================================================================================== + "effic_jetPhi_1 'JET_PHI_1; LEADING JET #phi; efficiency' jetPhi_1_numerator jetPhi_1_denominator", + "effic_jetPhi_2 'JET_PHI_2; SUBLEADING JET #phi; efficiency' jetPhi_2_numerator jetPhi_2_denominator", + "effic_jetEta_1 'JET_ETA1; LEADING JET #eta; efficiency' jetEta_1_numerator jetEta_1_denominator", + "effic_jetEta_2 'JET_ETA2; SUBLEADING JET #eta; efficiency' jetEta_2_numerator jetEta_2_denominator", + "effic_jetPt_1 'JET_PT1 turnON; LEADING JET PT [GeV]; efficiency' jetPt_1_numerator jetPt_1_denominator", + "effic_jetPt_2 'JET_PT2 turnON; SUBLEADING JET PT [GeV]; efficiency' jetPt_2_numerator jetPt_2_denominator", + "effic_NJets 'nJets; nJets; efficiency' NJets_numerator NJets_denominator", + "effic_jetPhi_inlc 'JET_PHI; JET #phi; efficiency' jetPhi_numerator jetPhi_denominator", + "effic_jetEta_incl 'JET_ETA; JET #eta; efficiency' jetEta_numerator jetEta_denominator", + "effic_jetPt_incl 'JET_PT turnON; JET PT [GeV]; efficiency' jetPt_numerator jetPt_denominator" + ), efficiencyProfile = cms.untracked.vstring( - "effic_met_vs_LS 'MET efficiency vs LS; LS; PF MET efficiency' metVsLS_numerator metVsLS_denominator", - "effic_ht_vs_LS 'HT efficiency vs LS; LS; PF HT efficiency' htVsLS_numerator htVsLS_denominator" + "effic_met_vs_LS 'MET efficiency vs LS; LS; PF MET efficiency' metVsLS_numerator metVsLS_denominator", + "effic_ht_vs_LS 'HT efficiency vs LS; LS; PF HT efficiency' htVsLS_numerator htVsLS_denominator" ), ) +caloMHTEfficiency = DQMEDHarvester("DQMGenericClient", + subDirs = cms.untracked.vstring("HLT/JME/CaloMHT/*"), + verbose = cms.untracked.uint32(0), # Set to 2 for all messages + resolution = cms.vstring(), + efficiency = cms.vstring( + "effic_met 'MET turnON; PF MET [GeV]; efficiency' met_numerator met_denominator", + "effic_met_variable 'MET turnON; PF MET [GeV]; efficiency' met_variable_numerator met_variable_denominator", + "effic_metPhi 'MET efficiency vs phi; PF MET phi [rad]; efficiency' metPhi_numerator metPhi_denominator", + "effic_ht 'HT turnON; PF HT [GeV]; efficiency' ht_numerator ht_denominator", + "effic_ht_variable 'HT turnON; PF HT [GeV]; efficiency' ht_variable_numerator ht_variable_denominator", + "effic_deltaphimetj1 'DELTAPHI turnON; DELTA PHI (PFMET, PFJET1); efficiency' deltaphi_metjet1_numerator deltaphi_metjet1_denominator", + "effic_deltaphij1j2 'DELTAPHI turnON; DELTA PHI (PFJET1, PFJET2); efficiency' deltaphi_jet1jet2_numerator deltaphi_jet1jet2_denominator", + #============================================================================================================================================== + "effic_jetPhi_1 'JET_PHI_1; LEADING JET #phi; efficiency' jetPhi_1_numerator jetPhi_1_denominator", + "effic_jetPhi_2 'JET_PHI_2; SUBLEADING JET #phi; efficiency' jetPhi_2_numerator jetPhi_2_denominator", + "effic_jetEta_1 'JET_ETA1; LEADING JET #eta; efficiency' jetEta_1_numerator jetEta_1_denominator", + "effic_jetEta_2 'JET_ETA2; SUBLEADING JET #eta; efficiency' jetEta_2_numerator jetEta_2_denominator", + "effic_jetPt_1 'JET_PT1 turnON; LEADING JET PT [GeV]; efficiency' jetPt_1_numerator jetPt_1_denominator", + "effic_jetPt_2 'JET_PT2 turnON; SUBLEADING JET PT [GeV]; efficiency' jetPt_2_numerator jetPt_2_denominator", + "effic_NJets 'nJets; nJets; efficiency' NJets_numerator NJets_denominator", + "effic_jetPhi_inlc 'JET_PHI; JET #phi; efficiency' jetPhi_numerator jetPhi_denominator", + "effic_jetEta_incl 'JET_ETA; JET #eta; efficiency' jetEta_numerator jetEta_denominator", + "effic_jetPt_incl 'JET_PT turnON; JET PT [GeV]; efficiency' jetPt_numerator jetPt_denominator" + + ), + efficiencyProfile = cms.untracked.vstring( + "effic_met_vs_LS 'MET efficiency vs LS; LS; PF MET efficiency' metVsLS_numerator metVsLS_denominator", + "effic_ht_vs_LS 'HT efficiency vs LS; LS; PF HT efficiency' htVsLS_numerator htVsLS_denominator" + ), + +) + metClient = cms.Sequence( metEfficiency + +caloMHTEfficiency ) diff --git a/HLTriggerOffline/JetMET/python/Validation/SingleJetValidation_cfi.py b/HLTriggerOffline/JetMET/python/Validation/SingleJetValidation_cfi.py index 123c1d5596952..fb0d00511031c 100644 --- a/HLTriggerOffline/JetMET/python/Validation/SingleJetValidation_cfi.py +++ b/HLTriggerOffline/JetMET/python/Validation/SingleJetValidation_cfi.py @@ -10,7 +10,7 @@ triggerEventObject = cms.untracked.InputTag("hltTriggerSummaryRAW","","HLT"), DQMFolder = cms.untracked.string(foldernm), PatternJetTrg = cms.untracked.string("HLT_PF(NoPU)?Jet([0-9])+(_v[0-9]+)?$"), - PatternMetTrg = cms.untracked.string("HLT_PF(ch)?MET([0-9])+(_HBHECleaned+)+(_v[0-9]+)?$"), + PatternMetTrg = cms.untracked.string("HLT_+[Calo|PF]+MET([0-9])+[_NotCleaned|_BeamHaloCleaned]+(_v[0-9]+)?$"), PatternMuTrg = cms.untracked.string("HLT_Mu([0-9])+(_v[0-9]+)?$"), LogFileName = cms.untracked.string('JetMETSingleJetValidation.log'), PFJetAlgorithm = cms.untracked.InputTag(pfjetcoll), diff --git a/HLTriggerOffline/JetMET/src/JetMETDQMPostProcessor.cc b/HLTriggerOffline/JetMET/src/JetMETDQMPostProcessor.cc index 29e0e7e44b6a2..728348b224a53 100644 --- a/HLTriggerOffline/JetMET/src/JetMETDQMPostProcessor.cc +++ b/HLTriggerOffline/JetMET/src/JetMETDQMPostProcessor.cc @@ -102,7 +102,7 @@ void JetMETDQMPostProcessor::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGe igetter, "_meGenJetPtTrg", "_meGenJetPtTrgLow", - "_meTurnOngJetPt", + "_meTurnOngJetPtLow", "Gen Jet PtLow", "Gen Jet Pt Turn-On Data"); dividehistos(ibooker, @@ -116,7 +116,7 @@ void JetMETDQMPostProcessor::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGe igetter, "_meGenJetEtaTrg", "_meGenJetEtaTrgLow", - "_meTurnOngJetEta", + "_meTurnOngJetEtaLow", "Gen Jet EtaLow", "Gen Jet Eta Turn-On Data"); dividehistos(ibooker, @@ -130,7 +130,7 @@ void JetMETDQMPostProcessor::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGe igetter, "_meGenJetPhiTrg", "_meGenJetPhiTrgLow", - "_meTurnOngJetPhi", + "_meTurnOngJetPhiLow", "Gen Jet PhiLow", "Gen Jet Phi Turn-On Data"); @@ -146,7 +146,7 @@ void JetMETDQMPostProcessor::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGe igetter, "_meHLTJetPtTrg", "_meHLTJetPtTrgLow", - "_meTurnOnhJetPt", + "_meTurnOnhJetPtLow", "HLT Jet PtLow", "HLT Jet Pt Turn-On Data"); dividehistos(ibooker, @@ -160,7 +160,7 @@ void JetMETDQMPostProcessor::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGe igetter, "_meHLTJetEtaTrg", "_meHLTJetEtaTrgLow", - "_meTurnOnhJetEta", + "_meTurnOnhJetEtaLow", "HLT Jet EtaLow", "HLT Jet Eta Turn-On Data"); dividehistos(ibooker, @@ -174,7 +174,7 @@ void JetMETDQMPostProcessor::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGe igetter, "_meHLTJetPhiTrg", "_meHLTJetPhiTrgLow", - "_meTurnOnhJetPhi", + "_meTurnOnhJetPhiLow", "HLT Jet PhiLow", "HLT Jet Phi Turn-On Data"); }