#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 50 of file SUSY_HLT_SingleLepton.h.

Constructor & Destructor Documentation

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

Definition at line 139 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 215 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 339 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 242 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 220 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 694 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 117 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::csvlCut_

private

Definition at line 105 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::csvmCut_

private

Definition at line 106 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::csvtCut_

private

Definition at line 107 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

HLTConfigProvider SUSY_HLT_SingleLepton::fHltConfig_

private

Definition at line 97 of file SUSY_HLT_SingleLepton.h.

Referenced by dqmBeginRun().

MonitorElement* SUSY_HLT_SingleLepton::h_btagTurnOn_den_

private

Definition at line 135 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_btagTurnOn_num_

private

Definition at line 134 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_CSVTurnOn_den_

private

Definition at line 133 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_CSVTurnOn_num_

private

Definition at line 132 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_HT_

private

Definition at line 123 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_leptonTurnOn_den_

private

Definition at line 127 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_leptonTurnOn_num_

private

Definition at line 126 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_maxCSV_

private

Definition at line 125 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_MET_

private

Definition at line 124 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_pfHTTurnOn_den_

private

Definition at line 129 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_pfHTTurnOn_num_

private

Definition at line 128 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_pfMetTurnOn_den_

private

Definition at line 131 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_pfMetTurnOn_num_

private

Definition at line 130 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_triggerLepEta_

private

Definition at line 121 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_triggerLepPhi_

private

Definition at line 122 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* SUSY_HLT_SingleLepton::h_triggerLepPt_

private

Definition at line 120 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

std::string SUSY_HLT_SingleLepton::HLTProcess_

private

Definition at line 99 of file SUSY_HLT_SingleLepton.h.

Referenced by dqmBeginRun().

double SUSY_HLT_SingleLepton::ht_threshold_

private

Definition at line 115 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::htCut_

private

Definition at line 112 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::jetEtaCut_

private

Definition at line 110 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::jetPtCut_

private

Definition at line 109 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::lep_pt_threshold_

private

Definition at line 114 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::met_threshold_

private

Definition at line 116 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

double SUSY_HLT_SingleLepton::metCut_

private

Definition at line 111 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 80 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theBeamSpotTag_

private

Definition at line 79 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 78 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theConversionCollectionTag_

private

Definition at line 77 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 65 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theElectronTag_

private

Definition at line 64 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 84 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theHLTHTTag_

private

Definition at line 83 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 88 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theHLTJetCollectionTag_

private

Definition at line 87 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 90 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theHLTJetTagCollectionTag_

private

Definition at line 89 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 86 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theHLTMETTag_

private

Definition at line 85 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 73 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theJetTagTag_

private

Definition at line 72 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theLeptonFilterTag_

private

Definition at line 82 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 67 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theMuonTag_

private

Definition at line 66 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 71 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::thePfJetTag_

private

Definition at line 70 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 69 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::thePfMETTag_

private

Definition at line 68 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 93 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theTriggerResultsTag_

private

Definition at line 92 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 95 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theTrigSummaryTag_

private

Definition at line 94 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

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

private

Definition at line 76 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

edm::InputTag SUSY_HLT_SingleLepton::theVertexCollectionTag_

private

Definition at line 75 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

std::string SUSY_HLT_SingleLepton::triggerPath_

private

Definition at line 101 of file SUSY_HLT_SingleLepton.h.

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

std::string SUSY_HLT_SingleLepton::triggerPathAuxiliary_

private

Definition at line 102 of file SUSY_HLT_SingleLepton.h.

Referenced by analyze().

std::string SUSY_HLT_SingleLepton::triggerPathLeptonAuxiliary_

private

Definition at line 103 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