#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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

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

void	endRun (const edm::Run &run, const edm::EventSetup &eSetup) override

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_

Detailed Description

Definition at line 49 of file SUSY_HLT_SingleLepton.h.

Constructor & Destructor Documentation

SUSY_HLT_SingleLepton::SUSY_HLT_SingleLepton ( const edm::ParameterSet & ps )

Definition at line 165 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 240 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
	)

overrideprotected

Definition at line 331 of file SUSY_HLT_SingleLepton.cc.

References edm::HLTGlobalStatus::accept(), align::BeamSpot, csv_threshold_, csvlCut_, csvmCut_, csvtCut_, metsig::electron, trigger::TriggerObject::eta(), MonitorElement::Fill(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), edm::Event::getByToken(), trigger::TriggerEvent::getObjects(), 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_, crabWrapper::key, edm::InputTag::label(), lep_pt_threshold_, met_threshold_, metCut_, dataset::name, pfMETBenchmark_cfi::pfMET, trigger::TriggerObject::phi(), reco::BeamSpot::position(), trigger::TriggerObject::pt(), edm::TriggerNames::size(), trigger::TriggerEvent::sizeFilters(), 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(), TriggerAnalyzer::triggerObjects, triggerPath_, triggerPathAuxiliary_, triggerPathLeptonAuxiliary_, trigNames, and edm::HLTGlobalStatus::wasrun().

                                                                                   {
   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);
           case 3:
             h_btagTurnOn_den_->Fill(3);
           case 2:
             h_btagTurnOn_den_->Fill(2);
           case 1:
             h_btagTurnOn_den_->Fill(1);
           case 0:
             h_btagTurnOn_den_->Fill(0);
         }
         switch (num_csvm) {
           default:
             h_btagTurnOn_den_->Fill(8);
           case 3:
             h_btagTurnOn_den_->Fill(7);
           case 2:
             h_btagTurnOn_den_->Fill(6);
           case 1:
             h_btagTurnOn_den_->Fill(5);
           case 0:
             break;  // Don't double count in the no tag bin
         }
         switch (num_csvt) {
           default:
             h_btagTurnOn_den_->Fill(12);
           case 3:
             h_btagTurnOn_den_->Fill(11);
           case 2:
             h_btagTurnOn_den_->Fill(10);
           case 1:
             h_btagTurnOn_den_->Fill(9);
           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);
           case 3:
             h_btagTurnOn_num_->Fill(3);
           case 2:
             h_btagTurnOn_num_->Fill(2);
           case 1:
             h_btagTurnOn_num_->Fill(1);
           case 0:
             h_btagTurnOn_num_->Fill(0);
         }
         switch (num_csvm) {
           default:
             h_btagTurnOn_num_->Fill(8);
           case 3:
             h_btagTurnOn_num_->Fill(7);
           case 2:
             h_btagTurnOn_num_->Fill(6);
           case 1:
             h_btagTurnOn_num_->Fill(5);
           case 0:
             break;  // Don't double count in the no tag bin
         }
         switch (num_csvt) {
           default:
             h_btagTurnOn_num_->Fill(12);
           case 3:
             h_btagTurnOn_num_->Fill(11);
           case 2:
             h_btagTurnOn_num_->Fill(10);
           case 1:
             h_btagTurnOn_num_->Fill(9);
           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 &
	)

overrideprotected

Definition at line 266 of file SUSY_HLT_SingleLepton.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::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(), DQMStore::IBooker::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
	)

overrideprotected

Definition at line 244 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";
 }

void SUSY_HLT_SingleLepton::endRun	(	const edm::Run &	run,
		const edm::EventSetup &	eSetup
	)

overrideprotected

Definition at line 721 of file SUSY_HLT_SingleLepton.cc.

References DEFINE_FWK_MODULE.

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

Member Data Documentation

double SUSY_HLT_SingleLepton::csv_threshold_

private

Definition at line 115 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::csvlCut_

private

Definition at line 103 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::csvmCut_

private

Definition at line 104 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::csvtCut_

private

Definition at line 105 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

HLTConfigProvider SUSY_HLT_SingleLepton::fHltConfig_

private

Definition at line 95 of file SUSY_HLT_SingleLepton.h.

Referenced by dqmBeginRun().

MonitorElement* SUSY_HLT_SingleLepton::h_btagTurnOn_den_

private

Definition at line 133 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_btagTurnOn_num_

private

Definition at line 132 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_CSVTurnOn_den_

private

Definition at line 131 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_CSVTurnOn_num_

private

Definition at line 130 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_HT_

private

Definition at line 121 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_leptonTurnOn_den_

private

Definition at line 125 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_leptonTurnOn_num_

private

Definition at line 124 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_maxCSV_

private

Definition at line 123 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_MET_

private

Definition at line 122 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_pfHTTurnOn_den_

private

Definition at line 127 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_pfHTTurnOn_num_

private

Definition at line 126 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_pfMetTurnOn_den_

private

Definition at line 129 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_pfMetTurnOn_num_

private

Definition at line 128 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_triggerLepEta_

private

Definition at line 119 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_triggerLepPhi_

private

Definition at line 120 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_triggerLepPt_

private

Definition at line 118 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

std::string SUSY_HLT_SingleLepton::HLTProcess_

private

Definition at line 97 of file SUSY_HLT_SingleLepton.h.

Referenced by dqmBeginRun().

double SUSY_HLT_SingleLepton::ht_threshold_

private

Definition at line 113 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::htCut_

private

Definition at line 110 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::jetEtaCut_

private

Definition at line 108 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::jetPtCut_

private

Definition at line 107 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::lep_pt_threshold_

private

Definition at line 112 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::met_threshold_

private

Definition at line 114 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::metCut_

private

Definition at line 109 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 78 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theBeamSpotTag_

private

Definition at line 77 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 76 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theConversionCollectionTag_

private

Definition at line 75 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 63 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theElectronTag_

private

Definition at line 62 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 82 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theHLTHTTag_

private

Definition at line 81 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 86 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theHLTJetCollectionTag_

private

Definition at line 85 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 88 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theHLTJetTagCollectionTag_

private

Definition at line 87 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 84 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theHLTMETTag_

private

Definition at line 83 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 71 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theJetTagTag_

private

Definition at line 70 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theLeptonFilterTag_

private

Definition at line 80 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 65 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theMuonTag_

private

Definition at line 64 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 69 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::thePfJetTag_

private

Definition at line 68 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 67 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::thePfMETTag_

private

Definition at line 66 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 91 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theTriggerResultsTag_

private

Definition at line 90 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 93 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theTrigSummaryTag_

private

Definition at line 92 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 74 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theVertexCollectionTag_

private

Definition at line 73 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

std::string SUSY_HLT_SingleLepton::triggerPath_

private

Definition at line 99 of file SUSY_HLT_SingleLepton.h.

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

std::string SUSY_HLT_SingleLepton::triggerPathAuxiliary_

private

Definition at line 100 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

std::string SUSY_HLT_SingleLepton::triggerPathLeptonAuxiliary_

private

Definition at line 101 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation