#include <SUSY_HLT_InclusiveHT.h>

Inheritance diagram for SUSY_HLT_InclusiveHT:

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

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

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

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

Private Attributes
double	etaThrJet_

MonitorElement *	h_caloHT

MonitorElement *	h_caloJetEta

MonitorElement *	h_caloJetPhi

MonitorElement *	h_caloJetPt

MonitorElement *	h_pfHT

MonitorElement *	h_pfHTTurnOn_den

MonitorElement *	h_pfHTTurnOn_num

MonitorElement *	h_pfJetEta

MonitorElement *	h_pfJetPhi

MonitorElement *	h_pfJetPt

MonitorElement *	h_pfMet

MonitorElement *	h_pfMetPhi

MonitorElement *	h_pfMetTurnOn_den

MonitorElement *	h_pfMetTurnOn_num

MonitorElement *	h_triggerHT

MonitorElement *	h_triggerJetEta

MonitorElement *	h_triggerJetPhi

MonitorElement *	h_triggerJetPt

MonitorElement *	h_triggerMetPhi

MonitorElement *	h_triggerMetPt

double	ptThrJet_

edm::EDGetTokenT< reco::CaloJetCollection >	theCaloJetCollection_

edm::EDGetTokenT< reco::PFJetCollection >	thePfJetCollection_

edm::EDGetTokenT< reco::PFMETCollection >	thePfMETCollection_

edm::EDGetTokenT< trigger::TriggerEvent >	theTrigSummary_

edm::InputTag	triggerFilter_

std::string	triggerPath_

std::string	triggerPathAuxiliaryForHadronic_

edm::EDGetTokenT< edm::TriggerResults >	triggerResults_

Detailed Description

Definition at line 27 of file SUSY_HLT_InclusiveHT.h.

Constructor & Destructor Documentation

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

Definition at line 8 of file SUSY_HLT_InclusiveHT.cc.

References etaThrJet_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ptThrJet_, AlCaHLTBitMon_QueryRunRegistry::string, theCaloJetCollection_, thePfJetCollection_, thePfMETCollection_, theTrigSummary_, triggerFilter_, triggerPath_, triggerPathAuxiliaryForHadronic_, and triggerResults_.

                                                                     {
   edm::LogInfo("SUSY_HLT_InclusiveHT") << "Constructor SUSY_HLT_InclusiveHT::SUSY_HLT_InclusiveHT " << 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"));
   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"));
   triggerPath_ = ps.getParameter<std::string>("TriggerPath");
   triggerPathAuxiliaryForHadronic_ = ps.getParameter<std::string>("TriggerPathAuxiliaryForHadronic");
   triggerFilter_ = ps.getParameter<edm::InputTag>("TriggerFilter");
   ptThrJet_ = ps.getUntrackedParameter<double>("PtThrJet");
   etaThrJet_ = ps.getUntrackedParameter<double>("EtaThrJet");
 }

SUSY_HLT_InclusiveHT::~SUSY_HLT_InclusiveHT ( )

override

Definition at line 23 of file SUSY_HLT_InclusiveHT.cc.

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

Member Function Documentation

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

overrideprotected

Definition at line 37 of file SUSY_HLT_InclusiveHT.cc.

References edm::HLTGlobalStatus::accept(), etaThrJet_, MonitorElement::Fill(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), edm::Event::getByToken(), trigger::TriggerEvent::getObjects(), h_caloHT, h_caloJetEta, h_caloJetPhi, h_caloJetPt, h_pfHT, h_pfHTTurnOn_den, h_pfHTTurnOn_num, h_pfJetEta, h_pfJetPhi, h_pfJetPt, h_pfMet, h_pfMetPhi, h_pfMetTurnOn_den, h_pfMetTurnOn_num, h_triggerHT, h_triggerMetPhi, h_triggerMetPt, trigger::TriggerObject::id(), edm::EventBase::isRealData(), edm::HandleBase::isValid(), relativeConstraints::keys, trigger::TriggerObject::phi(), trigger::TriggerObject::pt(), ptThrJet_, edm::TriggerNames::size(), trigger::TriggerEvent::sizeFilters(), theCaloJetCollection_, thePfJetCollection_, thePfMETCollection_, theTrigSummary_, triggerFilter_, edm::TriggerNames::triggerName(), edm::Event::triggerNames(), TriggerAnalyzer::triggerObjects, triggerPath_, triggerPathAuxiliaryForHadronic_, triggerResults_, trigNames, and edm::HLTGlobalStatus::wasrun().

                                                                                  {
   edm::LogInfo("SUSY_HLT_InclusiveHT") << "SUSY_HLT_InclusiveHT::analyze" << std::endl;
 
   //-------------------------------
   //--- MET
   //-------------------------------
   edm::Handle<reco::PFMETCollection> pfMETCollection;
   e.getByToken(thePfMETCollection_, pfMETCollection);
   if (!pfMETCollection.isValid()) {
     edm::LogError("SUSY_HLT_InclusiveHT") << "invalid collection: PFMET"
                                           << "\n";
     return;
   }
   //-------------------------------
   //--- Jets
   //-------------------------------
   edm::Handle<reco::PFJetCollection> pfJetCollection;
   e.getByToken(thePfJetCollection_, pfJetCollection);
   if (!pfJetCollection.isValid()) {
     edm::LogError("SUSY_HLT_InclusiveHT") << "invalid collection: PFJets"
                                           << "\n";
     return;
   }
   edm::Handle<reco::CaloJetCollection> caloJetCollection;
   e.getByToken(theCaloJetCollection_, caloJetCollection);
   if (!caloJetCollection.isValid()) {
     edm::LogError("SUSY_HLT_InclusiveHT") << "invalid collection: CaloJets"
                                           << "\n";
     return;
   }
 
   // check what is in the menu
   edm::Handle<edm::TriggerResults> hltresults;
   e.getByToken(triggerResults_, hltresults);
   if (!hltresults.isValid()) {
     edm::LogError("SUSY_HLT_InclusiveHT") << "invalid collection: TriggerResults"
                                           << "\n";
     return;
   }
 
   //-------------------------------
   //--- Trigger
   //-------------------------------
   edm::Handle<trigger::TriggerEvent> triggerSummary;
   e.getByToken(theTrigSummary_, triggerSummary);
   if (!triggerSummary.isValid()) {
     edm::LogError("SUSY_HLT_InclusiveHT") << "invalid collection: TriggerSummary"
                                           << "\n";
     return;
   }
 
   // get online objects
   size_t filterIndex = triggerSummary->filterIndex(triggerFilter_);
   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]];
       // if(foundObject.id() == 85 && foundObject.pt() > 40.0 &&
       // fabs(foundObject.eta()) < 3.0){
       //  h_triggerJetPt->Fill(foundObject.pt());
       //  h_triggerJetEta->Fill(foundObject.eta());
       //  h_triggerJetPhi->Fill(foundObject.phi());
       //}
       if (foundObject.id() == 87) {
         h_triggerMetPt->Fill(foundObject.pt());
         h_triggerMetPhi->Fill(foundObject.phi());
       }
       if (foundObject.id() == 89) {
         h_triggerHT->Fill(foundObject.pt());
       }
     }
   }
 
   bool hasFired = false, 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 || !e.isRealData()) {
     float caloHT = 0.0;
     float pfHT = 0.0;
     for (reco::PFJetCollection::const_iterator i_pfjet = pfJetCollection->begin(); i_pfjet != pfJetCollection->end();
          ++i_pfjet) {
       if (i_pfjet->pt() < ptThrJet_)
         continue;
       if (fabs(i_pfjet->eta()) > etaThrJet_)
         continue;
       pfHT += i_pfjet->pt();
     }
 
     if (hasFired) {
       for (reco::CaloJetCollection::const_iterator i_calojet = caloJetCollection->begin();
            i_calojet != caloJetCollection->end();
            ++i_calojet) {
         if (i_calojet->pt() < ptThrJet_)
           continue;
         if (fabs(i_calojet->eta()) > etaThrJet_)
           continue;
         h_caloJetPt->Fill(i_calojet->pt());
         h_caloJetEta->Fill(i_calojet->eta());
         h_caloJetPhi->Fill(i_calojet->phi());
         caloHT += i_calojet->pt();
       }
       for (reco::PFJetCollection::const_iterator i_pfjet = pfJetCollection->begin(); i_pfjet != pfJetCollection->end();
            ++i_pfjet) {
         if (i_pfjet->pt() < ptThrJet_)
           continue;
         if (fabs(i_pfjet->eta()) > etaThrJet_)
           continue;
         h_pfJetPt->Fill(i_pfjet->pt());
         h_pfJetEta->Fill(i_pfjet->eta());
         h_pfJetPhi->Fill(i_pfjet->phi());
       }
       h_pfMet->Fill(pfMETCollection->begin()->et());
       h_pfMetPhi->Fill(pfMETCollection->begin()->phi());
       h_pfHT->Fill(pfHT);
       h_caloHT->Fill(caloHT);
 
       h_pfMetTurnOn_num->Fill(pfMETCollection->begin()->et());
       h_pfHTTurnOn_num->Fill(pfHT);
     }
     // fill denominator histograms for all events, used for turn on curves
     h_pfMetTurnOn_den->Fill(pfMETCollection->begin()->et());
     h_pfHTTurnOn_den->Fill(pfHT);
   }
 }

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

overrideprotected

Definition at line 31 of file SUSY_HLT_InclusiveHT.cc.

References bookHistos().

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

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

private

Definition at line 176 of file SUSY_HLT_InclusiveHT.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::cd(), DEFINE_FWK_MODULE, h_caloHT, h_caloJetEta, h_caloJetPhi, h_caloJetPt, h_pfHT, h_pfHTTurnOn_den, h_pfHTTurnOn_num, h_pfJetEta, h_pfJetPhi, h_pfJetPt, h_pfMet, h_pfMetPhi, h_pfMetTurnOn_den, h_pfMetTurnOn_num, h_triggerHT, h_triggerMetPhi, h_triggerMetPt, DQMStore::IBooker::setCurrentFolder(), and triggerPath_.

Referenced by bookHistograms().

                                                                {
   ibooker_.cd();
   ibooker_.setCurrentFolder("HLT/SUSYBSM/" + triggerPath_);
 
   // offline quantities
   h_pfMet = ibooker_.book1D("pfMet", "PF Missing E_{T}; GeV", 20, 0.0, 500.0);
   h_pfMetPhi = ibooker_.book1D("pfMetPhi", "PF MET Phi", 20, -3.5, 3.5);
   h_pfHT = ibooker_.book1D("pfHT", "PF H_{T}; GeV", 30, 0.0, 1500.0);
   h_caloHT = ibooker_.book1D("caloHT", "Calo H_{T}; GeV", 30, 0.0, 1500.0);
   h_pfJetPt = ibooker_.book1D("pfJetPt", "PFJet P_{T}; GeV", 20, 0.0, 500.0);
   h_pfJetEta = ibooker_.book1D("pfJetEta", "PFJet Eta", 20, -3.0, 3.0);
   h_pfJetPhi = ibooker_.book1D("pfJetPhi", "PFJet Phi", 20, -3.5, 3.5);
   h_caloJetPt = ibooker_.book1D("caloJetPt", "CaloJet P_{T}; GeV", 20, 0.0, 500.0);
   h_caloJetEta = ibooker_.book1D("caloJetEta", "CaloJet Eta", 20, -3.0, 3.0);
   h_caloJetPhi = ibooker_.book1D("caloJetPhi", "CaloJet Phi", 20, -3.5, 3.5);
 
   // online quantities
   // h_triggerJetPt = ibooker_.book1D("triggerJetPt", "Trigger Jet Pt; 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);
   h_triggerMetPt = ibooker_.book1D("triggerMetPt", "Trigger Met Pt; GeV", 20, 0.0, 500.0);
   h_triggerMetPhi = ibooker_.book1D("triggerMetPhi", "Trigger Met Phi", 20, -3.5, 3.5);
   h_triggerHT = ibooker_.book1D("triggerHT", "Trigger HT; GeV", 30, 0.0, 1500.0);
 
   // num and den hists to be divided in harvesting step to make turn on curves
   h_pfMetTurnOn_num = ibooker_.book1D("pfMetTurnOn_num", "PF MET Turn On Numerator", 20, 0.0, 500.0);
   h_pfMetTurnOn_den = ibooker_.book1D("pfMetTurnOn_den", "PF MET Turn OnDenominator", 20, 0.0, 500.0);
   h_pfHTTurnOn_num = ibooker_.book1D("pfHTTurnOn_num", "PF HT Turn On Numerator", 30, 0.0, 1500.0);
   h_pfHTTurnOn_den = ibooker_.book1D("pfHTTurnOn_den", "PF HT Turn On Denominator", 30, 0.0, 1500.0);
 
   ibooker_.cd();
 }

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

overrideprotected

Definition at line 27 of file SUSY_HLT_InclusiveHT.cc.

                                                                             {
   edm::LogInfo("SUSY_HLT_InclusiveHT") << "SUSY_HLT_InclusiveHT::beginRun" << std::endl;
 }

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

overrideprotected

Definition at line 172 of file SUSY_HLT_InclusiveHT.cc.

                                                                                 {
   edm::LogInfo("SUSY_HLT_InclusiveHT") << "SUSY_HLT_InclusiveHT::endRun" << std::endl;
 }