#include <SUSY_HLT_DiJet_MET.h>

Inheritance diagram for SUSY_HLT_DiJet_MET:

Public Member Functions
	SUSY_HLT_DiJet_MET (const edm::ParameterSet &ps)

	~SUSY_HLT_DiJet_MET () 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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override

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

void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	bookHistos (DQMStore::IBooker &)

Private Attributes
double	etaThrJet_

double	etaThrJetTrig_

HLTConfigProvider	fHltConfig

MonitorElement *	h_caloMetvsPFMet

MonitorElement *	h_pfJet1Jet2DPhi

MonitorElement *	h_pfJet2PtTurnOn_den

MonitorElement *	h_pfJet2PtTurnOn_num

MonitorElement *	h_pfJetEta

MonitorElement *	h_pfJetPhi

MonitorElement *	h_pfJetPt

MonitorElement *	h_pfMetPhi

MonitorElement *	h_pfMetTurnOn_den

MonitorElement *	h_pfMetTurnOn_num

MonitorElement *	h_triggerJetEta

MonitorElement *	h_triggerJetPhi

MonitorElement *	h_triggerJetPt

MonitorElement *	h_triggerMet

MonitorElement *	h_triggerMetPhi

std::string	HLTProcess_

double	metCut_

double	ptThrJet_

double	ptThrJetTrig_

edm::EDGetTokenT< reco::CaloJetCollection >	theCaloJetCollection_

edm::EDGetTokenT< reco::CaloMETCollection >	theCaloMETCollection_

edm::EDGetTokenT< reco::PFJetCollection >	thePfJetCollection_

edm::EDGetTokenT< reco::PFMETCollection >	thePfMETCollection_

edm::EDGetTokenT< trigger::TriggerEvent >	theTrigSummary_

edm::InputTag	triggerFilter_

edm::InputTag	triggerJetFilter_

std::string	triggerPath_

std::string	triggerPathAuxiliaryForHadronic_

edm::EDGetTokenT< edm::TriggerResults >	triggerResults_

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 >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	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 33 of file SUSY_HLT_DiJet_MET.h.

Constructor & Destructor Documentation

◆ SUSY_HLT_DiJet_MET()

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

Definition at line 9 of file SUSY_HLT_DiJet_MET.cc.

                                                                 {
   edm::LogInfo("SUSY_HLT_DiJet_MET") << "Constructor SUSY_HLT_DiJet_MET::SUSY_HLT_DiJet_MET " << std::endl;
   // Get parameters from configuration file
   theTrigSummary_ = consumes<trigger::TriggerEvent>(ps.getParameter<edm::InputTag>("trigSummary"));
   thePfMETCollection_ = consumes<reco::PFMETCollection>(ps.getParameter<edm::InputTag>("pfMETCollection"));
   theCaloMETCollection_ = consumes<reco::CaloMETCollection>(ps.getParameter<edm::InputTag>("caloMETCollection"));
   thePfJetCollection_ = consumes<reco::PFJetCollection>(ps.getParameter<edm::InputTag>("pfJetCollection"));
   theCaloJetCollection_ = consumes<reco::CaloJetCollection>(ps.getParameter<edm::InputTag>("caloJetCollection"));
   triggerResults_ = consumes<edm::TriggerResults>(ps.getParameter<edm::InputTag>("TriggerResults"));
   HLTProcess_ = ps.getParameter<std::string>("HLTProcess");
   triggerPath_ = ps.getParameter<std::string>("TriggerPath");
   triggerPathAuxiliaryForHadronic_ = ps.getParameter<std::string>("TriggerPathAuxiliaryForHadronic");
   triggerFilter_ = ps.getParameter<edm::InputTag>("TriggerFilter");
   triggerJetFilter_ = ps.getParameter<edm::InputTag>("TriggerJetFilter");
   ptThrJetTrig_ = ps.getUntrackedParameter<double>("PtThrJetTrig");
   etaThrJetTrig_ = ps.getUntrackedParameter<double>("EtaThrJetTrig");
   ptThrJet_ = ps.getUntrackedParameter<double>("PtThrJet");
   etaThrJet_ = ps.getUntrackedParameter<double>("EtaThrJet");
   metCut_ = ps.getUntrackedParameter<double>("OfflineMetCut");
 }

References etaThrJet_, etaThrJetTrig_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLTProcess_, metCut_, ptThrJet_, ptThrJetTrig_, AlCaHLTBitMon_QueryRunRegistry::string, theCaloJetCollection_, theCaloMETCollection_, thePfJetCollection_, thePfMETCollection_, theTrigSummary_, triggerFilter_, triggerJetFilter_, triggerPath_, triggerPathAuxiliaryForHadronic_, and triggerResults_.

◆ ~SUSY_HLT_DiJet_MET()

SUSY_HLT_DiJet_MET::~SUSY_HLT_DiJet_MET ( )

override

Definition at line 30 of file SUSY_HLT_DiJet_MET.cc.

                                         {
   edm::LogInfo("SUSY_HLT_DiJet_MET") << "Destructor SUSY_HLT_DiJet_MET::~SUSY_HLT_DiJet_MET " << std::endl;
 }

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 65 of file SUSY_HLT_DiJet_MET.cc.

                                                                                {
   edm::LogInfo("SUSY_HLT_DiJet_MET") << "SUSY_HLT_DiJet_MET::analyze" << std::endl;
  
   //-------------------------------
   //--- MET
   //-------------------------------
   edm::Handle<reco::PFMETCollection> pfMETCollection;
   e.getByToken(thePfMETCollection_, pfMETCollection);
   if (!pfMETCollection.isValid()) {
     edm::LogError("SUSY_HLT_DiJet_MET") << "invalid collection: PFMET"
                                         << "\n";
     return;
   }
   edm::Handle<reco::CaloMETCollection> caloMETCollection;
   e.getByToken(theCaloMETCollection_, caloMETCollection);
   if (!caloMETCollection.isValid()) {
     edm::LogError("SUSY_HLT_DiJet_MET") << "invalid collection: CaloMET"
                                         << "\n";
     return;
   }
   //-------------------------------
   //--- Jets
   //-------------------------------
   edm::Handle<reco::PFJetCollection> pfJetCollection;
   e.getByToken(thePfJetCollection_, pfJetCollection);
   if (!pfJetCollection.isValid()) {
     edm::LogError("SUSY_HLT_DiJet_MET") << "invalid collection: PFJets"
                                         << "\n";
     return;
   }
   edm::Handle<reco::CaloJetCollection> caloJetCollection;
   e.getByToken(theCaloJetCollection_, caloJetCollection);
   if (!caloJetCollection.isValid()) {
     edm::LogError("SUSY_HLT_DiJet_MET") << "invalid collection: CaloJets"
                                         << "\n";
     return;
   }
  
   //-------------------------------
   //--- Trigger
   //-------------------------------
   edm::Handle<edm::TriggerResults> hltresults;
   e.getByToken(triggerResults_, hltresults);
   if (!hltresults.isValid()) {
     edm::LogError("SUSY_HLT_DiJet_MET") << "invalid collection: TriggerResults"
                                         << "\n";
     return;
   }
   edm::Handle<trigger::TriggerEvent> triggerSummary;
   e.getByToken(theTrigSummary_, triggerSummary);
   if (!triggerSummary.isValid()) {
     edm::LogError("SUSY_HLT_DiJet_MET") << "invalid collection: TriggerSummary"
                                         << "\n";
     return;
   }
  
   // get online objects
  
   size_t filterIndex = triggerSummary->filterIndex(triggerFilter_);
   size_t jetFilterIndex = triggerSummary->filterIndex(triggerJetFilter_);
   trigger::TriggerObjectCollection triggerObjects = triggerSummary->getObjects();
  
   if (!(filterIndex >= triggerSummary->sizeFilters())) {
     const trigger::Keys &keys = triggerSummary->filterKeys(filterIndex);
     for (size_t j = 0; j < keys.size(); ++j) {
       trigger::TriggerObject foundObject = triggerObjects[keys[j]];
       h_triggerMet->Fill(foundObject.pt());
       h_triggerMetPhi->Fill(foundObject.phi());
     }
   }
  
   std::vector<float> ptJet, etaJet, phiJet;
   if (!(jetFilterIndex >= triggerSummary->sizeFilters())) {
     const trigger::Keys &keys_jetfilter = triggerSummary->filterKeys(jetFilterIndex);
     for (size_t j = 0; j < keys_jetfilter.size(); ++j) {
       trigger::TriggerObject foundObject = triggerObjects[keys_jetfilter[j]];
       h_triggerJetPt->Fill(foundObject.pt());
       h_triggerJetEta->Fill(foundObject.eta());
       h_triggerJetPhi->Fill(foundObject.phi());
       if (foundObject.pt() > ptThrJetTrig_ && fabs(foundObject.eta()) < etaThrJetTrig_) {
         ptJet.push_back(foundObject.pt());
         etaJet.push_back(foundObject.eta());
         phiJet.push_back(foundObject.phi());
       }
     }
   }
  
   bool hasFired = false;
   bool hasFiredAuxiliaryForHadronicLeg = false;
   const edm::TriggerNames &trigNames = e.triggerNames(*hltresults);
   unsigned int numTriggers = trigNames.size();
   for (unsigned int hltIndex = 0; hltIndex < numTriggers; ++hltIndex) {
     if (trigNames.triggerName(hltIndex).find(triggerPath_) != std::string::npos && hltresults->wasrun(hltIndex) &&
         hltresults->accept(hltIndex))
       hasFired = true;
     if (trigNames.triggerName(hltIndex).find(triggerPathAuxiliaryForHadronic_) != std::string::npos &&
         hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex))
       hasFiredAuxiliaryForHadronicLeg = true;
   }
  
   if (hasFiredAuxiliaryForHadronicLeg) {
     int offlineJetCounter = 0;
     int offlineCentralJets = 0;
     int nMatch = 0, jet1Index = -1, jet2Index = -1;
  
     for (reco::PFJetCollection::const_iterator i_pfjet = pfJetCollection->begin(); i_pfjet != pfJetCollection->end();
          ++i_pfjet) {
       if (i_pfjet->pt() > ptThrJet_ && fabs(i_pfjet->eta()) < etaThrJet_) {
         offlineCentralJets++;
         if (offlineCentralJets == 2 && pfMETCollection->begin()->et() > metCut_)
           h_pfJet2PtTurnOn_den->Fill(i_pfjet->pt());
  
         for (unsigned int itrigjet = 0; itrigjet < ptJet.size(); ++itrigjet) {
           if (itrigjet < 2 &&
               (sqrt((i_pfjet->phi() - phiJet.at(itrigjet)) * (i_pfjet->phi() - phiJet.at(itrigjet)) +
                     (i_pfjet->eta() - etaJet.at(itrigjet)) * (i_pfjet->eta() - etaJet.at(itrigjet))) < 0.4)) {
             nMatch++;
             if (nMatch == 1)
               jet1Index = offlineJetCounter;
             if (nMatch == 2)
               jet2Index = offlineJetCounter;
  
             if (hasFired) {
               h_pfJetPt->Fill(i_pfjet->pt());
               h_pfJetEta->Fill(i_pfjet->eta());
               h_pfJetPhi->Fill(i_pfjet->phi());
               if (offlineCentralJets == 2 && pfMETCollection->begin()->et() > metCut_)
                 h_pfJet2PtTurnOn_num->Fill(i_pfjet->pt());
             }
  
             break;
           }
         }
       }
  
       offlineJetCounter++;
     }
  
     if (hasFired) {
       h_pfMetTurnOn_num->Fill(pfMETCollection->begin()->et());
       h_pfMetPhi->Fill(pfMETCollection->begin()->phi());
       h_caloMetvsPFMet->Fill(pfMETCollection->begin()->et(), caloMETCollection->begin()->et());
  
       if (jet1Index > -1 && jet2Index > -1) {
         h_pfJet1Jet2DPhi->Fill(
             fabs(reco::deltaPhi(pfJetCollection->at(jet1Index).phi(), pfJetCollection->at(jet2Index).phi())));
       }
     }
  
     h_pfMetTurnOn_den->Fill(pfMETCollection->begin()->et());
   }
  
   ptJet.clear();
   etaJet.clear();
   phiJet.clear();
 }

References pdwgLeptonRecoSkim_cfi::caloJetCollection, ewkDQM_cfi::caloMETCollection, reco::deltaPhi(), MillePedeFileConverter_cfg::e, trigger::TriggerObject::eta(), etaThrJet_, etaThrJetTrig_, dqm::impl::MonitorElement::Fill(), h_caloMetvsPFMet, h_pfJet1Jet2DPhi, h_pfJet2PtTurnOn_den, h_pfJet2PtTurnOn_num, h_pfJetEta, h_pfJetPhi, h_pfJetPt, h_pfMetPhi, h_pfMetTurnOn_den, h_pfMetTurnOn_num, h_triggerJetEta, h_triggerJetPhi, h_triggerJetPt, h_triggerMet, h_triggerMetPhi, HLTBitAnalyser_cfi::hltresults, dqmiolumiharvest::j, relativeConstraints::keys, metCut_, ExoticaDQM_cfi::pfJetCollection, B2GDQM_cfi::pfMETCollection, trigger::TriggerObject::phi(), trigger::TriggerObject::pt(), ptThrJet_, ptThrJetTrig_, mathSSE::sqrt(), theCaloJetCollection_, theCaloMETCollection_, thePfJetCollection_, thePfMETCollection_, theTrigSummary_, triggerFilter_, triggerJetFilter_, triggerMatchMonitor_cfi::triggerObjects, triggerPath_, triggerPathAuxiliaryForHadronic_, triggerResults_, HLTMuonOfflineAnalyzer_cfi::triggerSummary, and trigNames.

◆ bookHistograms()

void SUSY_HLT_DiJet_MET::bookHistograms	(	DQMStore::IBooker &	ibooker_,
		edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 59 of file SUSY_HLT_DiJet_MET.cc.

                                                                                                         {
   edm::LogInfo("SUSY_HLT_DiJet_MET") << "SUSY_HLT_DiJet_MET::bookHistograms" << std::endl;
   // book at beginRun
   bookHistos(ibooker_);
 }

References bookHistos().

◆ bookHistos()

void SUSY_HLT_DiJet_MET::bookHistos ( DQMStore::IBooker & ibooker_ )

private

Definition at line 222 of file SUSY_HLT_DiJet_MET.cc.

                                                              {
   ibooker_.cd();
   ibooker_.setCurrentFolder("HLT/SUSYBSM/" + triggerPath_);
  
   // offline quantities
   h_pfMetPhi = ibooker_.book1D("pfMetPhi", "PF MET Phi", 20, -3.5, 3.5);
   h_pfJetPt = ibooker_.book1D("pfJetPt", "PF Jet p_{T} (trigger-matched jets, |#eta| < 2.4); GeV", 20, 0.0, 500.0);
   h_pfJetEta = ibooker_.book1D("pfJetEta", "PF Jet #eta (trigger-matched jets, |#eta| < 2.4)", 20, -3.0, 3.0);
   h_pfJetPhi = ibooker_.book1D("pfJetPhi", "PF Jet #phi (trigger-matched jets, |#eta| < 2.4)", 20, -3.5, 3.5);
   h_pfJet1Jet2DPhi =
       ibooker_.book1D("pfJet1Jet2DPhi", "|#Delta#phi| between two leading trigger-matched jets", 20, 0.0, 3.5);
   h_caloMetvsPFMet = ibooker_.book2D("caloMetvsPFMet", "Calo MET vs PF MET; GeV; GeV", 25, 0.0, 500.0, 25, 0.0, 500.0);
  
   // online quantities
   h_triggerMet = ibooker_.book1D("triggerMet", "Trigger MET; GeV", 20, 0.0, 500.0);
   h_triggerMetPhi = ibooker_.book1D("triggerMetPhi", "Trigger MET Phi", 20, -3.5, 3.5);
   h_triggerJetPt = ibooker_.book1D("triggerJetPt", "Trigger Jet p_{T}; GeV", 20, 0.0, 500.0);
   h_triggerJetEta = ibooker_.book1D("triggerJetEta", "Trigger Jet Eta", 20, -3.0, 3.0);
   h_triggerJetPhi = ibooker_.book1D("triggerJetPhi", "Trigger Jet Phi", 20, -3.5, 3.5);
  
   // num and den hists to be divided in harvesting step to make turn on curves
   h_pfMetTurnOn_num = ibooker_.book1D("pfMetTurnOn_num", "PF MET", 20, 0.0, 500.0);
   h_pfMetTurnOn_den = ibooker_.book1D("pfMetTurnOn_den", "PF MET Turn On Denominator", 20, 0.0, 500.0);
   h_pfJet2PtTurnOn_num =
       ibooker_.book1D("pfJet2PtTurnOn_num", "PF Jet2 Pt (NCentralPFJets >= 2, PFMET > 250)", 20, 0.0, 500.0);
   h_pfJet2PtTurnOn_den = ibooker_.book1D("pfJet2PtTurnOn_den", "PF Jet2 Pt Turn On Denominator", 20, 0.0, 500.0);
  
   ibooker_.cd();
 }

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::NavigatorBase::cd(), h_caloMetvsPFMet, h_pfJet1Jet2DPhi, h_pfJet2PtTurnOn_den, h_pfJet2PtTurnOn_num, h_pfJetEta, h_pfJetPhi, h_pfJetPt, h_pfMetPhi, h_pfMetTurnOn_den, h_pfMetTurnOn_num, h_triggerJetEta, h_triggerJetPhi, h_triggerJetPt, h_triggerMet, h_triggerMetPhi, dqm::implementation::NavigatorBase::setCurrentFolder(), and triggerPath_.

Referenced by bookHistograms().

◆ dqmBeginRun()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 34 of file SUSY_HLT_DiJet_MET.cc.

                                                                               {
   bool changed;
  
   if (!fHltConfig.init(run, e, HLTProcess_, changed)) {
     edm::LogError("SUSY_HLT_DiJet_MET") << "Initialization of HLTConfigProvider failed!!";
     return;
   }
  
   bool pathFound = false;
   const std::vector<std::string> allTrigNames = fHltConfig.triggerNames();
   for (size_t j = 0; j < allTrigNames.size(); ++j) {
     if (allTrigNames[j].find(triggerPath_) != std::string::npos) {
       pathFound = true;
     }
   }
  
   if (!pathFound) {
     LogDebug("SUSY_HLT_DiJet_MET") << "Path not found"
                                    << "\n";
     return;
   }
  
   edm::LogInfo("SUSY_HLT_DiJet_MET") << "SUSY_HLT_DiJet_MET::beginRun" << std::endl;
 }