SUSY_HLT_SingleLepton Class Reference

#include <SUSY_HLT_SingleLepton.h>

Inheritance diagram for SUSY_HLT_SingleLepton:

Public Member Functions
	SUSY_HLT_SingleLepton (const edm::ParameterSet &ps)
	~SUSY_HLT_SingleLepton () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &eSetup) override
void	bookHistograms (DQMStore::IBooker &ibooker, const edm::Run &, const edm::EventSetup &) override
void	dqmBeginRun (const edm::Run &run, const edm::EventSetup &e) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
double	csv_threshold_
double	csvlCut_
double	csvmCut_
double	csvtCut_
HLTConfigProvider	fHltConfig_
MonitorElement *	h_btagTurnOn_den_
MonitorElement *	h_btagTurnOn_num_
MonitorElement *	h_CSVTurnOn_den_
MonitorElement *	h_CSVTurnOn_num_
MonitorElement *	h_HT_
MonitorElement *	h_leptonTurnOn_den_
MonitorElement *	h_leptonTurnOn_num_
MonitorElement *	h_maxCSV_
MonitorElement *	h_MET_
MonitorElement *	h_pfHTTurnOn_den_
MonitorElement *	h_pfHTTurnOn_num_
MonitorElement *	h_pfMetTurnOn_den_
MonitorElement *	h_pfMetTurnOn_num_
MonitorElement *	h_triggerLepEta_
MonitorElement *	h_triggerLepPhi_
MonitorElement *	h_triggerLepPt_
std::string	HLTProcess_
double	ht_threshold_
double	htCut_
double	jetEtaCut_
double	jetPtCut_
double	lep_pt_threshold_
double	met_threshold_
double	metCut_
edm::EDGetTokenT< reco::BeamSpot >	theBeamSpot_
edm::InputTag	theBeamSpotTag_
edm::EDGetTokenT < reco::ConversionCollection >	theConversionCollection_
edm::InputTag	theConversionCollectionTag_
edm::EDGetTokenT < reco::GsfElectronCollection >	theElectronCollection_
edm::InputTag	theElectronTag_
edm::EDGetTokenT < reco::METCollection >	theHLTHT_
edm::InputTag	theHLTHTTag_
edm::EDGetTokenT < reco::CaloJetCollection >	theHLTJetCollection_
edm::InputTag	theHLTJetCollectionTag_
edm::EDGetTokenT < reco::JetTagCollection >	theHLTJetTagCollection_
edm::InputTag	theHLTJetTagCollectionTag_
edm::EDGetTokenT < reco::METCollection >	theHLTMET_
edm::InputTag	theHLTMETTag_
edm::EDGetTokenT < reco::JetTagCollection >	theJetTagCollection_
edm::InputTag	theJetTagTag_
edm::InputTag	theLeptonFilterTag_
edm::EDGetTokenT < reco::MuonCollection >	theMuonCollection_
edm::InputTag	theMuonTag_
edm::EDGetTokenT < reco::PFJetCollection >	thePfJetCollection_
edm::InputTag	thePfJetTag_
edm::EDGetTokenT < reco::PFMETCollection >	thePfMETCollection_
edm::InputTag	thePfMETTag_
edm::EDGetTokenT < edm::TriggerResults >	theTriggerResults_
edm::InputTag	theTriggerResultsTag_
edm::EDGetTokenT < trigger::TriggerEvent >	theTrigSummary_
edm::InputTag	theTrigSummaryTag_
edm::EDGetTokenT < reco::VertexCollection >	theVertexCollection_
edm::InputTag	theVertexCollectionTag_
std::string	triggerPath_
std::string	triggerPathAuxiliary_
std::string	triggerPathLeptonAuxiliary_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr < DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 48 of file SUSY_HLT_SingleLepton.h.

Constructor & Destructor Documentation

SUSY_HLT_SingleLepton::SUSY_HLT_SingleLepton

(

const edm::ParameterSet &

ps

)

Definition at line 164 of file SUSY_HLT_SingleLepton.cc.

    : theElectronTag_(ps.getParameter<edm::InputTag>("electronCollection")),
      theElectronCollection_(consumes<reco::GsfElectronCollection>(theElectronTag_)),
      theMuonTag_(ps.getParameter<edm::InputTag>("muonCollection")),
      theMuonCollection_(consumes<reco::MuonCollection>(theMuonTag_)),
      thePfMETTag_(ps.getParameter<edm::InputTag>("pfMetCollection")),
      thePfMETCollection_(consumes<reco::PFMETCollection>(thePfMETTag_)),
      thePfJetTag_(ps.getParameter<edm::InputTag>("pfJetCollection")),
      thePfJetCollection_(consumes<reco::PFJetCollection>(thePfJetTag_)),
      theJetTagTag_(ps.getParameter<edm::InputTag>("jetTagCollection")),
      theJetTagCollection_(consumes<reco::JetTagCollection>(theJetTagTag_)),

      theVertexCollectionTag_(ps.getParameter<edm::InputTag>("vertexCollection")),
      theVertexCollection_(consumes<reco::VertexCollection>(theVertexCollectionTag_)),
      theConversionCollectionTag_(ps.getParameter<edm::InputTag>("conversionCollection")),
      theConversionCollection_(consumes<reco::ConversionCollection>(theConversionCollectionTag_)),
      theBeamSpotTag_(ps.getParameter<edm::InputTag>("beamSpot")),
      theBeamSpot_(consumes<reco::BeamSpot>(theBeamSpotTag_)),

      theLeptonFilterTag_(ps.getParameter<edm::InputTag>("leptonFilter")),
      theHLTHTTag_(ps.getParameter<edm::InputTag>("hltHt")),
      theHLTHT_(consumes<reco::METCollection>(theHLTHTTag_)),
      theHLTMETTag_(ps.getParameter<edm::InputTag>("hltMet")),
      theHLTMET_(consumes<reco::METCollection>(theHLTMETTag_)),
      theHLTJetCollectionTag_(ps.getParameter<edm::InputTag>("hltJets")),
      theHLTJetCollection_(consumes<reco::CaloJetCollection>(theHLTJetCollectionTag_)),
      theHLTJetTagCollectionTag_(ps.getParameter<edm::InputTag>("hltJetTags")),
      theHLTJetTagCollection_(consumes<reco::JetTagCollection>(theHLTJetTagCollectionTag_)),

      theTriggerResultsTag_(ps.getParameter<edm::InputTag>("triggerResults")),
      theTriggerResults_(consumes<edm::TriggerResults>(theTriggerResultsTag_)),
      theTrigSummaryTag_(ps.getParameter<edm::InputTag>("trigSummary")),
      theTrigSummary_(consumes<trigger::TriggerEvent>(theTrigSummaryTag_)),

      fHltConfig_(),

      HLTProcess_(ps.getParameter<std::string>("hltProcess")),

      triggerPath_(ps.getParameter<std::string>("triggerPath")),
      triggerPathAuxiliary_(ps.getParameter<std::string>("triggerPathAuxiliary")),
      triggerPathLeptonAuxiliary_(ps.getParameter<std::string>("triggerPathLeptonAuxiliary")),

      csvlCut_(ps.getUntrackedParameter<double>("csvlCut")),
      csvmCut_(ps.getUntrackedParameter<double>("csvmCut")),
      csvtCut_(ps.getUntrackedParameter<double>("csvtCut")),

      jetPtCut_(ps.getUntrackedParameter<double>("jetPtCut")),
      jetEtaCut_(ps.getUntrackedParameter<double>("jetEtaCut")),
      metCut_(ps.getUntrackedParameter<double>("metCut")),
      htCut_(ps.getUntrackedParameter<double>("htCut")),

      lep_pt_threshold_(ps.getUntrackedParameter<double>("leptonPtThreshold")),
      ht_threshold_(ps.getUntrackedParameter<double>("htThreshold")),
      met_threshold_(ps.getUntrackedParameter<double>("metThreshold")),
      csv_threshold_(ps.getUntrackedParameter<double>("csvThreshold")),

      h_triggerLepPt_(nullptr),
      h_triggerLepEta_(nullptr),
      h_triggerLepPhi_(nullptr),
      h_HT_(nullptr),
      h_MET_(nullptr),
      h_maxCSV_(nullptr),
      h_leptonTurnOn_num_(nullptr),
      h_leptonTurnOn_den_(nullptr),
      h_pfHTTurnOn_num_(nullptr),
      h_pfHTTurnOn_den_(nullptr),
      h_pfMetTurnOn_num_(nullptr),
      h_pfMetTurnOn_den_(nullptr),
      h_CSVTurnOn_num_(nullptr),
      h_CSVTurnOn_den_(nullptr),
      h_btagTurnOn_num_(nullptr),
      h_btagTurnOn_den_(nullptr) {
  edm::LogInfo("SUSY_HLT_SingleLepton") << "Constructor SUSY_HLT_SingleLepton::SUSY_HLT_SingleLepton\n";
}

SUSY_HLT_SingleLepton::~SUSY_HLT_SingleLepton ( )

override

Definition at line 239 of file SUSY_HLT_SingleLepton.cc.

                                              {
  edm::LogInfo("SUSY_HLT_SingleLepton") << "Destructor SUSY_HLT_SingleLepton::~SUSY_HLT_SingleLepton\n";
}

Member Function Documentation

void SUSY_HLT_SingleLepton::analyze ( const edm::Event &  e, const edm::EventSetup &  eSetup  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 330 of file SUSY_HLT_SingleLepton.cc.

References align::BeamSpot, csv_threshold_, csvlCut_, csvmCut_, csvtCut_, metsig::electron, trigger::TriggerObject::eta(), dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), h_btagTurnOn_den_, h_btagTurnOn_num_, h_CSVTurnOn_den_, h_CSVTurnOn_num_, h_HT_, h_leptonTurnOn_den_, h_leptonTurnOn_num_, h_maxCSV_, h_MET_, h_pfHTTurnOn_den_, h_pfHTTurnOn_num_, h_pfMetTurnOn_den_, h_pfMetTurnOn_num_, h_triggerLepEta_, h_triggerLepPhi_, h_triggerLepPt_, ht_threshold_, htCut_, edm::EventBase::isRealData(), edm::HandleBase::isValid(), metsig::jet, jetEtaCut_, jetPtCut_, submitPVResolutionJobs::key, edm::InputTag::label(), lep_pt_threshold_, met_threshold_, metCut_, HLT_FULL_cff::muon, mergeVDriftHistosByStation::name, HLT_FULL_cff::pfJetCollection, pfMET_cfi::pfMET, trigger::TriggerObject::phi(), trigger::TriggerObject::pt(), edm::TriggerNames::size(), AlCaHLTBitMon_QueryRunRegistry::string, theBeamSpot_, theBeamSpotTag_, theConversionCollection_, theConversionCollectionTag_, theElectronCollection_, theElectronTag_, theHLTHT_, theHLTHTTag_, theHLTJetCollection_, theHLTJetCollectionTag_, theHLTJetTagCollection_, theHLTJetTagCollectionTag_, theHLTMET_, theHLTMETTag_, theJetTagCollection_, theJetTagTag_, theLeptonFilterTag_, theMuonCollection_, theMuonTag_, thePfJetCollection_, thePfJetTag_, thePfMETCollection_, thePfMETTag_, theTriggerResults_, theTriggerResultsTag_, theTrigSummary_, theTrigSummaryTag_, theVertexCollection_, theVertexCollectionTag_, edm::TriggerNames::triggerName(), edm::Event::triggerNames(), triggerPath_, triggerPathAuxiliary_, triggerPathLeptonAuxiliary_, and trigNames.

                                                                                  {
  edm::LogInfo("SUSY_HLT_SingleLepton") << "SUSY_HLT_SingleLepton::analyze\n";

  // HLT HT
  edm::Handle<reco::METCollection> HLTHT;
  if (!theHLTHTTag_.label().empty()) {
    e.getByToken(theHLTHT_, HLTHT);
    if (!HLTHT.isValid()) {
      edm::LogInfo("SUSY_HLT_SingleLepton") << "Invalid METCollection: " << theHLTHTTag_.label() << '\n';
    }
  }

  // HLT MET
  edm::Handle<reco::METCollection> HLTMET;
  if (!theHLTMETTag_.label().empty()) {
    e.getByToken(theHLTMET_, HLTMET);
    if (!HLTMET.isValid()) {
      edm::LogInfo("SUSY_HLT_SingleLepton") << "Invalid METCollection: " << theHLTMETTag_.label() << '\n';
    }
  }

  // HLT Jets
  edm::Handle<reco::CaloJetCollection> HLTJetCollection;
  if (!theHLTJetCollectionTag_.label().empty()) {
    e.getByToken(theHLTJetCollection_, HLTJetCollection);
    if (!HLTJetCollection.isValid()) {
      edm::LogInfo("SUSY_HLT_SingleLepton") << "Invalid CaloJetCollection: " << theHLTJetCollectionTag_.label() << '\n';
    }
  }

  // HLT Jet Tags
  edm::Handle<reco::JetTagCollection> HLTJetTagCollection;
  if (!theHLTJetTagCollectionTag_.label().empty()) {
    e.getByToken(theHLTJetTagCollection_, HLTJetTagCollection);
    if (!HLTJetTagCollection.isValid()) {
      edm::LogInfo("SUSY_HLT_SingleLepton")
          << "Invalid JetTagCollection: " << theHLTJetTagCollectionTag_.label() << '\n';
    }
  }

  // Vertex
  edm::Handle<reco::VertexCollection> VertexCollection;
  if (!theVertexCollectionTag_.label().empty()) {
    e.getByToken(theVertexCollection_, VertexCollection);
    if (!VertexCollection.isValid()) {
      edm::LogWarning("SUSY_HLT_SingleLepton")
          << "Invalid VertexCollection: " << theVertexCollectionTag_.label() << '\n';
    }
  }

  // Conversions
  edm::Handle<reco::ConversionCollection> ConversionCollection;
  if (!theConversionCollectionTag_.label().empty()) {
    e.getByToken(theConversionCollection_, ConversionCollection);
    if (!ConversionCollection.isValid()) {
      edm::LogWarning("SUSY_HLT_SingleLepton")
          << "Invalid ConversionCollection: " << theConversionCollectionTag_.label() << '\n';
    }
  }

  // Beam Spot
  edm::Handle<reco::BeamSpot> BeamSpot;
  if (!theBeamSpotTag_.label().empty()) {
    e.getByToken(theBeamSpot_, BeamSpot);
    if (!BeamSpot.isValid()) {
      edm::LogWarning("SUSY_HLT_SingleLepton") << "Invalid BeamSpot: " << theBeamSpotTag_.label() << '\n';
    }
  }

  // MET
  edm::Handle<reco::PFMETCollection> pfMETCollection;
  if (!thePfMETTag_.label().empty()) {
    e.getByToken(thePfMETCollection_, pfMETCollection);
    if (!pfMETCollection.isValid()) {
      edm::LogWarning("SUSY_HLT_SingleLepton") << "Invalid PFMETCollection: " << thePfMETTag_.label() << '\n';
    }
  }

  // Jets
  edm::Handle<reco::PFJetCollection> pfJetCollection;
  if (!thePfJetTag_.label().empty()) {
    e.getByToken(thePfJetCollection_, pfJetCollection);
    if (!pfJetCollection.isValid()) {
      edm::LogWarning("SUSY_HLT_SingleLepton") << "Invalid PFJetCollection: " << thePfJetTag_.label() << '\n';
    }
  }

  // b-tags
  edm::Handle<reco::JetTagCollection> jetTagCollection;
  if (!theJetTagTag_.label().empty()) {
    e.getByToken(theJetTagCollection_, jetTagCollection);
    if (!jetTagCollection.isValid()) {
      edm::LogWarning("SUSY_HLT_SingleLepton") << "Invalid JetTagCollection: " << theJetTagTag_.label() << '\n';
    }
  }

  // Electron
  edm::Handle<reco::GsfElectronCollection> ElectronCollection;
  if (!theElectronTag_.label().empty()) {
    e.getByToken(theElectronCollection_, ElectronCollection);
    if (!ElectronCollection.isValid()) {
      edm::LogWarning("SUSY_HLT_SingleLepton") << "Invalid GsfElectronCollection: " << theElectronTag_.label() << '\n';
    }
  }

  // Muon
  edm::Handle<reco::MuonCollection> MuonCollection;
  if (!theMuonTag_.label().empty()) {
    e.getByToken(theMuonCollection_, MuonCollection);
    if (!MuonCollection.isValid()) {
      edm::LogWarning("SUSY_HLT_SingleLepton") << "Invalid MuonCollection: " << theMuonTag_.label() << '\n';
    }
  }

  // Trigger
  edm::Handle<edm::TriggerResults> hltresults;
  if (!theTriggerResultsTag_.label().empty()) {
    e.getByToken(theTriggerResults_, hltresults);
    if (!hltresults.isValid()) {
      edm::LogWarning("SUSY_HLT_SingleLepton") << "Invalid TriggerResults: " << theTriggerResultsTag_.label() << '\n';
    }
  }
  edm::Handle<trigger::TriggerEvent> triggerSummary;
  if (!theTrigSummaryTag_.label().empty()) {
    e.getByToken(theTrigSummary_, triggerSummary);
    if (!triggerSummary.isValid()) {
      edm::LogWarning("SUSY_HLT_SingleLepton") << "Invalid TriggerEvent: " << theTrigSummaryTag_.label() << '\n';
    }
  }

  // Get online leptons
  std::vector<float> ptLepton, etaLepton, phiLepton;
  if (triggerSummary.isValid()) {
    // Leptons
    size_t filterIndex = triggerSummary->filterIndex(theLeptonFilterTag_);
    trigger::TriggerObjectCollection triggerObjects = triggerSummary->getObjects();
    if (!(filterIndex >= triggerSummary->sizeFilters())) {
      size_t ilep = 0, num_keys = triggerSummary->filterKeys(filterIndex).size();
      ptLepton.resize(num_keys);
      etaLepton.resize(num_keys);
      phiLepton.resize(num_keys);
      for (const auto &key : triggerSummary->filterKeys(filterIndex)) {
        const trigger::TriggerObject &foundObject = triggerObjects[key];

        if (h_triggerLepPt_)
          h_triggerLepPt_->Fill(foundObject.pt());
        if (h_triggerLepEta_)
          h_triggerLepEta_->Fill(foundObject.eta());
        if (h_triggerLepPhi_)
          h_triggerLepPhi_->Fill(foundObject.phi());

        ptLepton.at(ilep) = foundObject.pt();
        etaLepton.at(ilep) = foundObject.eta();
        phiLepton.at(ilep) = foundObject.phi();
        ++ilep;
      }
    }
  }

  // Get online ht and met
  const float hlt_ht = ((HLTHT.isValid() && !HLTHT->empty()) ? HLTHT->front().sumEt() : -1.0);
  if (h_HT_)
    h_HT_->Fill(hlt_ht);
  const float hlt_met = ((HLTMET.isValid() && !HLTMET->empty()) ? HLTMET->front().pt() : -1.0);
  if (h_MET_)
    h_MET_->Fill(hlt_met);

  // Get online csv and fill plot
  float hlt_csv = -1.0;
  if (HLTJetCollection.isValid() && HLTJetTagCollection.isValid()) {
    for (const auto &jet : *HLTJetTagCollection) {
      if (jet.second > hlt_csv)
        hlt_csv = jet.second;
    }
  }
  if (h_maxCSV_)
    h_maxCSV_->Fill(hlt_csv);

  // Test whether main and auxilliary triggers fired
  bool hasFired = false;
  bool hasFiredAuxiliary = false;
  bool hasFiredLeptonAuxiliary = false;
  if (hltresults.isValid()) {
    const edm::TriggerNames &trigNames = e.triggerNames(*hltresults);
    for (unsigned int hltIndex = 0; hltIndex < trigNames.size(); ++hltIndex) {
      if (hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex)) {
        const std::string &name = trigNames.triggerName(hltIndex);
        if (Contains(name, triggerPath_))
          hasFired = true;
        if (Contains(name, triggerPathAuxiliary_))
          hasFiredAuxiliary = true;
        if (Contains(name, triggerPathLeptonAuxiliary_))
          hasFiredLeptonAuxiliary = true;
      }
    }
  }

  // Get offline HT
  double pfHT = -1.0;
  if (pfJetCollection.isValid()) {
    pfHT = 0.0;
    for (const auto &pfjet : *pfJetCollection) {
      if (pfjet.pt() < jetPtCut_)
        continue;
      if (fabs(pfjet.eta()) > jetEtaCut_)
        continue;
      pfHT += pfjet.pt();
    }
  }

  // Get offline MET
  double pfMET = -1.0;
  if (pfMETCollection.isValid() && !pfMETCollection->empty()) {
    pfMET = pfMETCollection->front().et();
  }

  // Get offline b-tagging info
  float maxCSV = -1.0;
  unsigned num_csvl = 0;
  unsigned num_csvm = 0;
  unsigned num_csvt = 0;
  if (jetTagCollection.isValid()) {
    for (const auto &jet : *jetTagCollection) {
      const float CSV = jet.second;
      if (jet.first->pt() > jetPtCut_) {
        if (CSV > maxCSV) {
          maxCSV = CSV;
        }
        if (CSV > csvlCut_) {
          ++num_csvl;
          if (CSV > csvmCut_) {
            ++num_csvm;
            if (CSV > csvtCut_) {
              ++num_csvt;
            }
          }
        }
      }
    }
  }
  if (h_maxCSV_)
    h_maxCSV_->Fill(maxCSV);

  // Fill lepton pt efficiency plot
  double lep_max_pt = -1.0;
  if (VertexCollection.isValid() && !VertexCollection->empty()) {  // for quality checks
    // Try to find a reco electron
    if (ElectronCollection.isValid() && ConversionCollection.isValid() && BeamSpot.isValid()) {
      for (const auto &electron : *ElectronCollection) {
        if (IsGood(electron, VertexCollection->front().position(), BeamSpot->position(), ConversionCollection)) {
          if (electron.pt() > lep_max_pt)
            lep_max_pt = electron.pt();
        }
      }
    }

    // Try to find a reco muon
    if (MuonCollection.isValid()) {
      for (const auto &muon : *MuonCollection) {
        if (IsGood(muon, VertexCollection->front().position())) {
          if (muon.pt() > lep_max_pt) {
            lep_max_pt = muon.pt();
          }
        }
      }
    }
  }

  const bool lep_plateau = lep_max_pt > lep_pt_threshold_ || lep_pt_threshold_ < 0.0;
  const bool ht_plateau = pfHT > ht_threshold_ || ht_threshold_ < 0.0;
  const bool met_plateau = pfMET > met_threshold_ || met_threshold_ < 0.0;
  const bool csv_plateau = maxCSV > csv_threshold_ || csv_threshold_ < 0.0;

  // Fill lepton turn-on histograms
  if (hasFiredLeptonAuxiliary || triggerPathLeptonAuxiliary_.empty() || !e.isRealData()) {
    // Fill histograms using highest pt reco lepton
    if (ht_plateau && met_plateau && csv_plateau && (pfMET > metCut_ || metCut_ < 0.0) &&
        (pfHT > htCut_ || htCut_ < 0.0)) {
      if (h_leptonTurnOn_den_)
        h_leptonTurnOn_den_->Fill(lep_max_pt);
      if (h_leptonTurnOn_num_ && hasFired)
        h_leptonTurnOn_num_->Fill(lep_max_pt);
    }
  }

  // Fill remaining turn-on histograms
  if (hasFiredAuxiliary || triggerPathAuxiliary_.empty() || !e.isRealData()) {
    // Fill HT efficiency plot
    if (lep_plateau && met_plateau && csv_plateau) {
      if (h_pfHTTurnOn_den_)
        h_pfHTTurnOn_den_->Fill(pfHT);
      if (h_pfHTTurnOn_num_ && hasFired)
        h_pfHTTurnOn_num_->Fill(pfHT);
    }

    // Fill MET efficiency plot
    if (lep_plateau && ht_plateau && csv_plateau) {
      if (h_pfMetTurnOn_den_)
        h_pfMetTurnOn_den_->Fill(pfMET);
      if (h_pfMetTurnOn_num_ && hasFired)
        h_pfMetTurnOn_num_->Fill(pfMET);
    }

    // Fill CSV efficiency plot
    if (lep_plateau && ht_plateau && met_plateau) {
      if (h_CSVTurnOn_den_)
        h_CSVTurnOn_den_->Fill(maxCSV);
      if (h_CSVTurnOn_num_ && hasFired)
        h_CSVTurnOn_num_->Fill(maxCSV);

      if (h_btagTurnOn_den_) {
        switch (num_csvl) {
          default:
            h_btagTurnOn_den_->Fill(4);
            [[fallthrough]];
          case 3:
            h_btagTurnOn_den_->Fill(3);
            [[fallthrough]];
          case 2:
            h_btagTurnOn_den_->Fill(2);
            [[fallthrough]];
          case 1:
            h_btagTurnOn_den_->Fill(1);
            [[fallthrough]];
          case 0:
            h_btagTurnOn_den_->Fill(0);
        }
        switch (num_csvm) {
          default:
            h_btagTurnOn_den_->Fill(8);
            [[fallthrough]];
          case 3:
            h_btagTurnOn_den_->Fill(7);
            [[fallthrough]];
          case 2:
            h_btagTurnOn_den_->Fill(6);
            [[fallthrough]];
          case 1:
            h_btagTurnOn_den_->Fill(5);
            [[fallthrough]];
          case 0:
            break;  // Don't double count in the no tag bin
        }
        switch (num_csvt) {
          default:
            h_btagTurnOn_den_->Fill(12);
            [[fallthrough]];
          case 3:
            h_btagTurnOn_den_->Fill(11);
            [[fallthrough]];
          case 2:
            h_btagTurnOn_den_->Fill(10);
            [[fallthrough]];
          case 1:
            h_btagTurnOn_den_->Fill(9);
            [[fallthrough]];
          case 0:
            break;  // Don't double count in the no tag bin
        }
      }
      if (h_btagTurnOn_num_ && hasFired) {
        switch (num_csvl) {
          default:
            h_btagTurnOn_num_->Fill(4);
            [[fallthrough]];
          case 3:
            h_btagTurnOn_num_->Fill(3);
            [[fallthrough]];
          case 2:
            h_btagTurnOn_num_->Fill(2);
            [[fallthrough]];
          case 1:
            h_btagTurnOn_num_->Fill(1);
            [[fallthrough]];
          case 0:
            h_btagTurnOn_num_->Fill(0);
        }
        switch (num_csvm) {
          default:
            h_btagTurnOn_num_->Fill(8);
            [[fallthrough]];
          case 3:
            h_btagTurnOn_num_->Fill(7);
            [[fallthrough]];
          case 2:
            h_btagTurnOn_num_->Fill(6);
            [[fallthrough]];
          case 1:
            h_btagTurnOn_num_->Fill(5);
            [[fallthrough]];
          case 0:
            break;  // Don't double count in the no tag bin
        }
        switch (num_csvt) {
          default:
            h_btagTurnOn_num_->Fill(12);
            [[fallthrough]];
          case 3:
            h_btagTurnOn_num_->Fill(11);
            [[fallthrough]];
          case 2:
            h_btagTurnOn_num_->Fill(10);
            [[fallthrough]];
          case 1:
            h_btagTurnOn_num_->Fill(9);
            [[fallthrough]];
          case 0:
            break;  // Don't double count in the no tag bin
        }
      }
    }
  }
}

void SUSY_HLT_SingleLepton::bookHistograms ( DQMStore::IBooker &  ibooker, const edm::Run &  , const edm::EventSetup &    )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 265 of file SUSY_HLT_SingleLepton.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::NavigatorBase::cd(), h_btagTurnOn_den_, h_btagTurnOn_num_, h_CSVTurnOn_den_, h_CSVTurnOn_num_, h_HT_, h_leptonTurnOn_den_, h_leptonTurnOn_num_, h_maxCSV_, h_MET_, h_pfHTTurnOn_den_, h_pfHTTurnOn_num_, h_pfMetTurnOn_den_, h_pfMetTurnOn_num_, h_triggerLepEta_, h_triggerLepPhi_, h_triggerLepPt_, edm::InputTag::label(), dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, theElectronTag_, theHLTHTTag_, theHLTJetCollectionTag_, theHLTJetTagCollectionTag_, theHLTMETTag_, theMuonTag_, and triggerPath_.

                                                                                                          {
  edm::LogInfo("SUSY_HLT_SingleLepton") << "SUSY_HLT_SingleLepton::bookHistograms\n";
  // book at beginRun
  ibooker.cd();
  ibooker.setCurrentFolder("HLT/SUSYBSM/" + triggerPath_);

  bool is_mu = false;
  bool is_ele = false;
  if (theElectronTag_.label().empty() && !theMuonTag_.label().empty()) {
    is_mu = true;
  } else if (!theElectronTag_.label().empty() && theMuonTag_.label().empty()) {
    is_ele = true;
  }
  std::string lepton = "lepton", Lepton = "Lepton";
  if (is_mu && !is_ele) {
    lepton = "muon";
    Lepton = "Muon";
  } else if (is_ele && !is_mu) {
    lepton = "electron";
    Lepton = "Electron";
  }

  // online quantities
  h_triggerLepPt_ = ibooker.book1D("triggerLepPt", (";" + Lepton + " p_{T} [GeV];").c_str(), 20, 0.0, 500.0);
  h_triggerLepEta_ = ibooker.book1D("triggerLepEta", (";" + Lepton + " #eta;").c_str(), 20, -3.0, 3.0);
  h_triggerLepPhi_ = ibooker.book1D("triggerLepPhi", (";" + Lepton + " #phi;").c_str(), 20, -3.5, 3.5);

  if (!theHLTHTTag_.label().empty()) {
    h_HT_ = ibooker.book1D("HT", ";HLT HT [GeV];", 40, 0.0, 1000.0);
  }

  if (!theHLTMETTag_.label().empty()) {
    h_MET_ = ibooker.book1D("MET", ";HLT MET [GeV];", 40, 0.0, 1000.0);
  }

  if (!theHLTJetCollectionTag_.label().empty() && !theHLTJetTagCollectionTag_.label().empty()) {
    h_maxCSV_ = ibooker.book1D("maxCSV", ";Max HLT CSV;", 20, 0.0, 1.0);
  }

  // num and den hists to be divided in harvesting step to make turn on curves
  h_leptonTurnOn_num_ =
      ibooker.book1D("leptonTurnOn_num", ("Numerator;Offline " + lepton + " p_{T} [GeV];").c_str(), 30, 0.0, 150);
  h_leptonTurnOn_den_ =
      ibooker.book1D("leptonTurnOn_den", ("Denominator;Offline " + lepton + " p_{T} [GeV];").c_str(), 30, 0.0, 150.0);
  h_pfHTTurnOn_num_ = ibooker.book1D("pfHTTurnOn_num", "Numerator;Offline H_{T} [GeV];", 30, 0.0, 1500.0);
  h_pfHTTurnOn_den_ = ibooker.book1D("pfHTTurnOn_den", "Denominator;Offline H_{T} [GeV];", 30, 0.0, 1500.0);

  if (!theHLTMETTag_.label().empty()) {
    h_pfMetTurnOn_num_ = ibooker.book1D("pfMetTurnOn_num", "Numerator;Offline MET [GeV];", 20, 0.0, 500.0);
    h_pfMetTurnOn_den_ = ibooker.book1D("pfMetTurnOn_den", "Denominator;Offline MET [GeV];", 20, 0.0, 500.0);
  }

  if (!theHLTJetCollectionTag_.label().empty() && !theHLTJetTagCollectionTag_.label().empty()) {
    h_CSVTurnOn_num_ = ibooker.book1D("CSVTurnOn_num", "Numerator;Offline Max CSV Discriminant;", 20, 0.0, 1.0);
    h_CSVTurnOn_den_ = ibooker.book1D("CSVTurnOn_den", "Denominator;Offline Max CSV Discriminant;", 20, 0.0, 1.0);

    h_btagTurnOn_num_ = ibooker.book1D("btagTurnOn_num", "Numerator;Offline CSV Requirement;", 13, -0.5, 12.5);
    h_btagTurnOn_den_ = ibooker.book1D("btagTurnOn_den", "Denominator;Offline CSV Requirements;", 13, -0.5, 12.5);

    SetBinLabels(h_btagTurnOn_num_);
    SetBinLabels(h_btagTurnOn_den_);
  }
  ibooker.cd();
}

void SUSY_HLT_SingleLepton::dqmBeginRun ( const edm::Run &  run, const edm::EventSetup &  e  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 243 of file SUSY_HLT_SingleLepton.cc.

References fHltConfig_, HLTProcess_, HLTConfigProvider::init(), LogDebug, HLTConfigProvider::triggerNames(), and triggerPath_.

                                                                                 {
  bool changed;

  if (!fHltConfig_.init(run, e, HLTProcess_, changed)) {
    edm::LogError("SUSY_HLT_SingleLepton") << "Initialization of HLTConfigProvider failed!!\n";
    return;
  }

  bool pathFound = false;
  for (const auto &trig_name : fHltConfig_.triggerNames()) {
    if (Contains(trig_name, triggerPath_))
      pathFound = true;
  }

  if (!pathFound) {
    LogDebug("SUSY_HLT_SingleLepton") << "Path not found: " << triggerPath_ << '\n';
    return;
  }

  edm::LogInfo("SUSY_HLT_SingleLepton") << "SUSY_HLT_SingleLepton::beginRun\n";
}

Member Data Documentation

double SUSY_HLT_SingleLepton::csv_threshold_

private

Definition at line 113 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::csvlCut_

private

Definition at line 101 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::csvmCut_

private

Definition at line 102 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::csvtCut_

private

Definition at line 103 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

HLTConfigProvider SUSY_HLT_SingleLepton::fHltConfig_

private

Definition at line 93 of file SUSY_HLT_SingleLepton.h.

Referenced by dqmBeginRun().

MonitorElement* SUSY_HLT_SingleLepton::h_btagTurnOn_den_

private

Definition at line 131 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_btagTurnOn_num_

private

Definition at line 130 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_CSVTurnOn_den_

private

Definition at line 129 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_CSVTurnOn_num_

private

Definition at line 128 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_HT_

private

Definition at line 119 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_leptonTurnOn_den_

private

Definition at line 123 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_leptonTurnOn_num_

private

Definition at line 122 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_maxCSV_

private

Definition at line 121 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_MET_

private

Definition at line 120 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_pfHTTurnOn_den_

private

Definition at line 125 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_pfHTTurnOn_num_

private

Definition at line 124 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_pfMetTurnOn_den_

private

Definition at line 127 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_pfMetTurnOn_num_

private

Definition at line 126 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_triggerLepEta_

private

Definition at line 117 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_triggerLepPhi_

private

Definition at line 118 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_triggerLepPt_

private

Definition at line 116 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

std::string SUSY_HLT_SingleLepton::HLTProcess_

private

Definition at line 95 of file SUSY_HLT_SingleLepton.h.

Referenced by dqmBeginRun().

double SUSY_HLT_SingleLepton::ht_threshold_

private

Definition at line 111 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::htCut_

private

Definition at line 108 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::jetEtaCut_

private

Definition at line 106 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::jetPtCut_

private

Definition at line 105 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::lep_pt_threshold_

private

Definition at line 110 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::met_threshold_

private

Definition at line 112 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::metCut_

private

Definition at line 107 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::EDGetTokenT<reco::BeamSpot> SUSY_HLT_SingleLepton::theBeamSpot_

private

Definition at line 76 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theBeamSpotTag_

private

Definition at line 75 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::EDGetTokenT<reco::ConversionCollection> SUSY_HLT_SingleLepton::theConversionCollection_

private

Definition at line 74 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theConversionCollectionTag_

private

Definition at line 73 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::EDGetTokenT<reco::GsfElectronCollection> SUSY_HLT_SingleLepton::theElectronCollection_

private

Definition at line 61 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theElectronTag_

private

Definition at line 60 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<reco::METCollection> SUSY_HLT_SingleLepton::theHLTHT_

private

Definition at line 80 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theHLTHTTag_

private

Definition at line 79 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<reco::CaloJetCollection> SUSY_HLT_SingleLepton::theHLTJetCollection_

private

Definition at line 84 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theHLTJetCollectionTag_

private

Definition at line 83 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<reco::JetTagCollection> SUSY_HLT_SingleLepton::theHLTJetTagCollection_

private

Definition at line 86 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theHLTJetTagCollectionTag_

private

Definition at line 85 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<reco::METCollection> SUSY_HLT_SingleLepton::theHLTMET_

private

Definition at line 82 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theHLTMETTag_

private

Definition at line 81 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<reco::JetTagCollection> SUSY_HLT_SingleLepton::theJetTagCollection_

private

Definition at line 69 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theJetTagTag_

private

Definition at line 68 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theLeptonFilterTag_

private

Definition at line 78 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::EDGetTokenT<reco::MuonCollection> SUSY_HLT_SingleLepton::theMuonCollection_

private

Definition at line 63 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theMuonTag_

private

Definition at line 62 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<reco::PFJetCollection> SUSY_HLT_SingleLepton::thePfJetCollection_

private

Definition at line 67 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::thePfJetTag_

private

Definition at line 66 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::EDGetTokenT<reco::PFMETCollection> SUSY_HLT_SingleLepton::thePfMETCollection_

private

Definition at line 65 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::thePfMETTag_

private

Definition at line 64 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::EDGetTokenT<edm::TriggerResults> SUSY_HLT_SingleLepton::theTriggerResults_

private

Definition at line 89 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theTriggerResultsTag_

private

Definition at line 88 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::EDGetTokenT<trigger::TriggerEvent> SUSY_HLT_SingleLepton::theTrigSummary_

private

Definition at line 91 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theTrigSummaryTag_

private

Definition at line 90 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::EDGetTokenT<reco::VertexCollection> SUSY_HLT_SingleLepton::theVertexCollection_

private

Definition at line 72 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theVertexCollectionTag_

private

Definition at line 71 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

std::string SUSY_HLT_SingleLepton::triggerPath_

private

Definition at line 97 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), bookHistograms(), and dqmBeginRun().

std::string SUSY_HLT_SingleLepton::triggerPathAuxiliary_

private

Definition at line 98 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

std::string SUSY_HLT_SingleLepton::triggerPathLeptonAuxiliary_

private

Definition at line 99 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().