#include <SUSY_HLT_Muon_Hadronic.h>

Inheritance diagram for SUSY_HLT_Muon_Hadronic:

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

	~SUSY_HLT_Muon_Hadronic () 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	etaMuonOffline_

double	etaThrJet_

HLTConfigProvider	fHltConfig

MonitorElement *	h_MuTurnOn_den

MonitorElement *	h_MuTurnOn_num

MonitorElement *	h_pfHTTurnOn_den

MonitorElement *	h_pfHTTurnOn_num

MonitorElement *	h_pfMetTurnOn_den

MonitorElement *	h_pfMetTurnOn_num

MonitorElement *	h_triggerMuEta

MonitorElement *	h_triggerMuPhi

MonitorElement *	h_triggerMuPt

std::string	HLTProcess_

double	HTOffline_

double	METOffline_

double	ptMuonOffline_

double	ptThrJet_

edm::EDGetTokenT< reco::CaloJetCollection >	theCaloJetCollection_

edm::EDGetTokenT< reco::MuonCollection >	theMuonCollection_

edm::EDGetTokenT< reco::PFJetCollection >	thePfJetCollection_

edm::EDGetTokenT< reco::PFMETCollection >	thePfMETCollection_

edm::EDGetTokenT< trigger::TriggerEvent >	theTrigSummary_

edm::InputTag	triggerFilter_

std::string	triggerPath_

std::string	triggerPathAuxiliaryForHadronic_

std::string	triggerPathAuxiliaryForMuon_

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_Muon_Hadronic.h.

Constructor & Destructor Documentation

◆ SUSY_HLT_Muon_Hadronic()

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

Definition at line 10 of file SUSY_HLT_Muon_Hadronic.cc.

                                                                         {
   edm::LogInfo("SUSY_HLT_Muon_Hadronic") << "Constructor SUSY_HLT_Muon_Hadronic::SUSY_HLT_Muon_Hadronic " << std::endl;
   // Get parameters from configuration file
   theTrigSummary_ = consumes<trigger::TriggerEvent>(ps.getParameter<edm::InputTag>("trigSummary"));
   theMuonCollection_ = consumes<reco::MuonCollection>(ps.getParameter<edm::InputTag>("MuonCollection"));
   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"));
   HLTProcess_ = ps.getParameter<std::string>("HLTProcess");
   triggerPath_ = ps.getParameter<std::string>("TriggerPath");
   triggerPathAuxiliaryForMuon_ = ps.getParameter<std::string>("TriggerPathAuxiliaryForMuon");
   triggerPathAuxiliaryForHadronic_ = ps.getParameter<std::string>("TriggerPathAuxiliaryForHadronic");
   triggerFilter_ = ps.getParameter<edm::InputTag>("TriggerFilter");
   ptMuonOffline_ = ps.getUntrackedParameter<double>("ptMuonOffline");
   etaMuonOffline_ = ps.getUntrackedParameter<double>("etaMuonOffline");
   HTOffline_ = ps.getUntrackedParameter<double>("HTOffline");
   METOffline_ = ps.getUntrackedParameter<double>("METOffline");
   ptThrJet_ = ps.getUntrackedParameter<double>("PtThrJet");
   etaThrJet_ = ps.getUntrackedParameter<double>("EtaThrJet");
 }

References etaMuonOffline_, etaThrJet_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLTProcess_, HTOffline_, METOffline_, ptMuonOffline_, ptThrJet_, AlCaHLTBitMon_QueryRunRegistry::string, theCaloJetCollection_, theMuonCollection_, thePfJetCollection_, thePfMETCollection_, theTrigSummary_, triggerFilter_, triggerPath_, triggerPathAuxiliaryForHadronic_, triggerPathAuxiliaryForMuon_, and triggerResults_.

◆ ~SUSY_HLT_Muon_Hadronic()

SUSY_HLT_Muon_Hadronic::~SUSY_HLT_Muon_Hadronic ( )

override

Definition at line 32 of file SUSY_HLT_Muon_Hadronic.cc.

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

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 72 of file SUSY_HLT_Muon_Hadronic.cc.

                                                                                    {
   edm::LogInfo("SUSY_HLT_Muon_Hadronic") << "SUSY_HLT_Muon_Hadronic::analyze" << std::endl;
  
   //-------------------------------
   //--- MET
   //-------------------------------
   edm::Handle<reco::PFMETCollection> pfMETCollection;
   e.getByToken(thePfMETCollection_, pfMETCollection);
   if (!pfMETCollection.isValid()) {
     edm::LogError("SUSY_HLT_Muon_Hadronic") << "invalid collection: PFMET"
                                             << "\n";
     return;
   }
   //-------------------------------
   //--- Jets
   //-------------------------------
   edm::Handle<reco::PFJetCollection> pfJetCollection;
   e.getByToken(thePfJetCollection_, pfJetCollection);
   if (!pfJetCollection.isValid()) {
     edm::LogError("SUSY_HLT_Muon_Hadronic") << "invalid collection: PFJets"
                                             << "\n";
     return;
   }
   edm::Handle<reco::CaloJetCollection> caloJetCollection;
   e.getByToken(theCaloJetCollection_, caloJetCollection);
   if (!caloJetCollection.isValid()) {
     edm::LogError("SUSY_HLT_Muon_Hadronic") << "invalid collection: CaloJets"
                                             << "\n";
     return;
   }
  
   //-------------------------------
   //--- Muon
   //-------------------------------
   edm::Handle<reco::MuonCollection> MuonCollection;
   e.getByToken(theMuonCollection_, MuonCollection);
   if (!MuonCollection.isValid()) {
     edm::LogError("SUSY_HLT_Mu_Hadronic") << "invalid collection: Muons "
                                           << "\n";
     return;
   }
  
   //-------------------------------
   //--- Trigger
   //-------------------------------
   edm::Handle<edm::TriggerResults> hltresults;
   e.getByToken(triggerResults_, hltresults);
   if (!hltresults.isValid()) {
     edm::LogError("SUSY_HLT_Muon_Hadronic") << "invalid collection: TriggerResults"
                                             << "\n";
     return;
   }
   edm::Handle<trigger::TriggerEvent> triggerSummary;
   e.getByToken(theTrigSummary_, triggerSummary);
   if (!triggerSummary.isValid()) {
     edm::LogError("SUSY_HLT_Muon_Hadronic") << "invalid collection: TriggerSummary"
                                             << "\n";
     return;
   }
  
   // get online objects
   std::vector<Lepton> onlineMuons;
   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 (fabs(foundObject.id()) == 13) {  // It's a muon
         h_triggerMuPt->Fill(foundObject.pt());
         h_triggerMuEta->Fill(foundObject.eta());
         h_triggerMuPhi->Fill(foundObject.phi());
         Lepton theMuon;
         theMuon.pt = foundObject.pt();
         theMuon.phi = foundObject.phi();
         theMuon.eta = foundObject.eta();
         onlineMuons.push_back(theMuon);
       }
     }
   }
  
   bool hasFired = false;
   bool hasFiredAuxiliaryForMuonLeg = 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(triggerPathAuxiliaryForMuon_) != std::string::npos &&
         hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex))
       hasFiredAuxiliaryForMuonLeg = true;
     if (trigNames.triggerName(hltIndex).find(triggerPathAuxiliaryForHadronic_) != std::string::npos &&
         hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex))
       hasFiredAuxiliaryForHadronicLeg = true;
   }
  
   if (hasFiredAuxiliaryForMuonLeg || hasFiredAuxiliaryForHadronicLeg || !e.isRealData()) {
     std::vector<Lepton> offlineMuons;
     for (reco::MuonCollection::const_iterator muon = MuonCollection->begin(); muon != MuonCollection->end(); ++muon) {
       if (fabs(muon->eta()) > etaMuonOffline_)
         continue;
       Lepton theMuon;
       theMuon.pt = muon->pt();
       theMuon.phi = muon->phi();
       theMuon.eta = muon->eta();
       offlineMuons.push_back(theMuon);
     }
  
     std::sort(onlineMuons.begin(), onlineMuons.end(), comparePt);
     std::sort(offlineMuons.begin(), offlineMuons.end(), comparePt);
     /*
     std::cout << "OFFLINE MUON" << std::endl;
     for(std::vector<Lepton>::iterator offline_muon = offlineMuons.begin();
     offline_muon != offlineMuons.end(); offline_muon++) { std::cout <<
     offline_muon->pt << std::endl;
     }
     std::cout << "ONLINE MUON" << std::endl;
     for(std::vector<Lepton>::iterator online_muon = onlineMuons.begin();
     online_muon != onlineMuons.end(); online_muon++) { std::cout <<
     online_muon->pt << std::endl;
     }
     */
     // Matching the muon
     int indexOfMatchedMuon = -1;
     int offlineCounter = 0;
     for (std::vector<Lepton>::iterator offline_muon = offlineMuons.begin();
          (offline_muon != offlineMuons.end() && indexOfMatchedMuon == -1);
          offline_muon++) {
       for (std::vector<Lepton>::iterator online_muon = onlineMuons.begin(); online_muon != onlineMuons.end();
            online_muon++) {
         if (sqrt((offline_muon->phi - online_muon->phi) * (offline_muon->phi - online_muon->phi) +
                  (offline_muon->eta - online_muon->eta) * (offline_muon->eta - online_muon->eta)) < 0.5) {
           indexOfMatchedMuon = offlineCounter;
           // std::cout << "Offline " << offline_muon->pt << " " << "Online " <<
           // online_muon->pt << std::endl;
           break;
         }
       }
       offlineCounter++;
     }
     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();
     }
     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;
       caloHT += i_calojet->pt();
     }
     // Need to apply the MET offline cut to be in the MET plateau
     if ((hasFiredAuxiliaryForMuonLeg || !e.isRealData()) && !offlineMuons.empty() &&
         pfMETCollection->begin()->et() > METOffline_ && pfHT > HTOffline_) {
       if (hasFired && indexOfMatchedMuon >= 0) {
         h_MuTurnOn_num->Fill(offlineMuons.at(indexOfMatchedMuon).pt);
         h_MuTurnOn_den->Fill(offlineMuons.at(indexOfMatchedMuon).pt);
       } else {
         h_MuTurnOn_den->Fill(offlineMuons.at(0).pt);
       }
     }
     // Need to apply the pt offline cut to be in the muon pt plateau
     if ((hasFiredAuxiliaryForHadronicLeg || !e.isRealData()) && indexOfMatchedMuon >= 0 &&
         offlineMuons.at(indexOfMatchedMuon).pt > ptMuonOffline_ && pfHT > HTOffline_) {
       if (hasFired) {
         h_pfMetTurnOn_num->Fill(pfMETCollection->begin()->et());
       }
       h_pfMetTurnOn_den->Fill(pfMETCollection->begin()->et());
     }
     if ((hasFiredAuxiliaryForHadronicLeg || !e.isRealData()) && indexOfMatchedMuon >= 0 &&
         offlineMuons.at(indexOfMatchedMuon).pt > ptMuonOffline_ && pfMETCollection->begin()->et() > METOffline_) {
       if (hasFired) {
         h_pfHTTurnOn_num->Fill(pfHT);
       }
       h_pfHTTurnOn_den->Fill(pfHT);
     }
   }
 }

References pdwgLeptonRecoSkim_cfi::caloJetCollection, comparePt(), MillePedeFileConverter_cfg::e, trigger::TriggerObject::eta(), Lepton::eta, etaMuonOffline_, etaThrJet_, dqm::impl::MonitorElement::Fill(), h_MuTurnOn_den, h_MuTurnOn_num, h_pfHTTurnOn_den, h_pfHTTurnOn_num, h_pfMetTurnOn_den, h_pfMetTurnOn_num, h_triggerMuEta, h_triggerMuPhi, h_triggerMuPt, HLTBitAnalyser_cfi::hltresults, HTOffline_, trigger::TriggerObject::id(), dqmiolumiharvest::j, relativeConstraints::keys, METOffline_, HLT_2018_cff::muon, ExoticaDQM_cfi::pfJetCollection, B2GDQM_cfi::pfMETCollection, trigger::TriggerObject::phi(), Lepton::phi, trigger::TriggerObject::pt(), Lepton::pt, ptMuonOffline_, ptThrJet_, mathSSE::sqrt(), theCaloJetCollection_, theMuonCollection_, thePfJetCollection_, thePfMETCollection_, theTrigSummary_, triggerFilter_, triggerMatchMonitor_cfi::triggerObjects, triggerPath_, triggerPathAuxiliaryForHadronic_, triggerPathAuxiliaryForMuon_, triggerResults_, HLTMuonOfflineAnalyzer_cfi::triggerSummary, and trigNames.

◆ bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 66 of file SUSY_HLT_Muon_Hadronic.cc.

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

References bookHistos().

◆ bookHistos()

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

private

Definition at line 261 of file SUSY_HLT_Muon_Hadronic.cc.

                                                                  {
   ibooker_.cd();
   ibooker_.setCurrentFolder("HLT/SUSYBSM/" + triggerPath_);
  
   // offline quantities
  
   // online quantities
   h_triggerMuPt = ibooker_.book1D("triggerMuPt", "Trigger Muon Pt; GeV", 20, 0.0, 500.0);
   h_triggerMuEta = ibooker_.book1D("triggerMuEta", "Trigger Muon Eta", 20, -3.0, 3.0);
   h_triggerMuPhi = ibooker_.book1D("triggerMuPhi", "Trigger Muon 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 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);
   h_MuTurnOn_num = ibooker_.book1D("MuTurnOn_num", "Muon Turn On Numerator", 30, 0.0, 150);
   h_MuTurnOn_den = ibooker_.book1D("MuTurnOn_den", "Muon Turn On Denominator", 30, 0.0, 150.0);
  
   ibooker_.cd();
 }

References dqm::implementation::IBooker::book1D(), dqm::implementation::NavigatorBase::cd(), h_MuTurnOn_den, h_MuTurnOn_num, h_pfHTTurnOn_den, h_pfHTTurnOn_num, h_pfMetTurnOn_den, h_pfMetTurnOn_num, h_triggerMuEta, h_triggerMuPhi, h_triggerMuPt, dqm::implementation::NavigatorBase::setCurrentFolder(), and triggerPath_.

Referenced by bookHistograms().

◆ dqmBeginRun()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 36 of file SUSY_HLT_Muon_Hadronic.cc.

                                                                                   {
   bool changed;
  
   if (!fHltConfig.init(run, e, HLTProcess_, changed)) {
     edm::LogError("SUSY_HLT_Muon_Hadronic") << "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_Muon_Hadronic") << "Path not found"
                                        << "\n";
     return;
   }
   // std::vector<std::string> filtertags = fHltConfig.moduleLabels( triggerPath_
   // ); triggerFilter_ =
   // edm::InputTag(filtertags[filtertags.size()-1],"",fHltConfig.processName());
   // triggerFilter_ = edm::InputTag("hltPFMET120Mu5L3PreFiltered", "",
   // fHltConfig.processName());
  
   edm::LogInfo("SUSY_HLT_Muon_Hadronic") << "SUSY_HLT_Muon_Hadronic::beginRun" << std::endl;
 }