#include <JetMETHLTOfflineSource.h>

Inheritance diagram for JetMETHLTOfflineSource:

Classes
class	PathInfo

class	PathInfoCollection

Public Member Functions
	JetMETHLTOfflineSource (const edm::ParameterSet &)

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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

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

virtual void	fillMEforEffAllTrigger (const edm::Event &iEvent, const edm::EventSetup &)

virtual void	fillMEforMonAllTrigger (const edm::Event &iEvent, const edm::EventSetup &)

virtual void	fillMEforMonTriggerSummary (const edm::Event &iEvent, const edm::EventSetup &)

virtual void	fillMEforTriggerNTfired ()

const std::string	getL1ConditionModuleName (const std::string &pathname)

virtual bool	isBarrel (double eta)

virtual bool	isEndCap (double eta)

virtual bool	isForward (double eta)

virtual bool	isHLTPathAccepted (std::string pathName)

virtual bool	isTriggerObjectFound (std::string objectName)

virtual double	TriggerPosition (std::string trigName)

virtual bool	validPathHLT (std::string path)

Private Attributes
double	_fEMF

double	_feta

double	_fHPD

double	_max_CEMF

double	_max_CHEF

double	_max_NEMF

double	_max_NHEF

double	_min_CEMF

double	_min_CHEF

double	_min_NEMF

double	_min_NHEF

double	_n90Hits

double	_pfMHT

reco::CaloJetCollection	calojet

edm::Handle< reco::CaloJetCollection >	calojetColl_

edm::EDGetTokenT< reco::JetCorrector >	CaloJetCorToken_

double	CaloJetEMF [2]

double	CaloJetEta [2]

double	CaloJetfHPD [2]

double	CaloJetn90 [2]

double	CaloJetPhi [2]

double	CaloJetPt [2]

double	CaloJetPx [2]

double	CaloJetPy [2]

edm::EDGetTokenT< reco::CaloJetCollection >	caloJetsToken

edm::Handle< reco::CaloMETCollection >	calometColl_

edm::EDGetTokenT< reco::CaloMETCollection >	caloMetToken

MonitorElement *	correlation_All

MonitorElement *	correlation_AllWrtMB

MonitorElement *	correlation_AllWrtMu

std::vector< std::pair< std::string, std::string > >	custompathnamepairs_

std::string	dirname_

HLTConfigProvider	hltConfig_

PathInfoCollection	hltPathsAll_

PathInfoCollection	hltPathsAllTriggerSummary_

PathInfoCollection	hltPathsEff_

bool	isSetup_

reco::helper::JetIDHelper *	jetID

std::vector< std::string >	MBTrigPaths_

std::vector< std::string >	MuonTrigPaths_

bool	nameForEff_

MonitorElement *	NVertices

std::vector< std::string >	pathFilter_

std::vector< std::string >	pathRejectKeyword_

reco::PFJetCollection	pfjet

double	PFJetCEMF [2]

double	PFJetCHEF [2]

edm::Handle< reco::PFJetCollection >	pfjetColl_

edm::EDGetTokenT< reco::JetCorrector >	PFJetCorToken_

double	PFJetEta [2]

double	PFJetNEMF [2]

double	PFJetNHEF [2]

double	PFJetPhi [2]

double	PFJetPt [2]

double	PFJetPx [2]

double	PFJetPy [2]

edm::EDGetTokenT< reco::PFJetCollection >	pfJetsToken

edm::Handle< reco::PFMETCollection >	pfmetColl_

edm::EDGetTokenT< reco::PFMETCollection >	pfMetToken

double	pfMHTx_All

double	pfMHTy_All

bool	plotAll_

bool	plotEff_

std::vector< int >	prescUsed_

std::string	processname_

MonitorElement *	PVZ

MonitorElement *	rate_All

MonitorElement *	rate_AllWrtMB

MonitorElement *	rate_AllWrtMu

bool	runStandalone_

edm::TriggerNames	triggerNames_

edm::Handle< trigger::TriggerEvent >	triggerObj_

edm::Handle< edm::TriggerResults >	triggerResults_

edm::EDGetTokenT< edm::TriggerResults >	triggerResultsFUToken

edm::InputTag	triggerResultsLabel_

edm::EDGetTokenT< edm::TriggerResults >	triggerResultsToken

edm::EDGetTokenT< trigger::TriggerEvent >	triggerSummaryFUToken

edm::InputTag	triggerSummaryLabel_

edm::EDGetTokenT< trigger::TriggerEvent >	triggerSummaryToken

bool	verbose_

edm::EDGetTokenT< reco::VertexCollection >	vertexToken

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

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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 Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 67 of file JetMETHLTOfflineSource.h.

Constructor & Destructor Documentation

◆ JetMETHLTOfflineSource()

JetMETHLTOfflineSource::JetMETHLTOfflineSource ( const edm::ParameterSet & iConfig )

explicit

Definition at line 38 of file JetMETHLTOfflineSource.cc.

References _fEMF, _feta, _fHPD, _max_CEMF, _max_CHEF, _max_NEMF, _max_NHEF, _min_CEMF, _min_CHEF, _min_NEMF, _min_NHEF, _n90Hits, CaloJetCorToken_, caloJetsToken, caloMetToken, custompathnamepairs_, dirname_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ProducerED_cfi::InputTag, edm::InputTag::instance(), jetID, edm::InputTag::label(), LogDebug, MBTrigPaths_, MuonTrigPaths_, nameForEff_, castor_dqm_sourceclient_file_cfg::path, pathFilter_, pathRejectKeyword_, Skims_PA_cff::paths, PFJetCorToken_, pfJetsToken, pfMetToken, plotAll_, plotEff_, processname_, runStandalone_, AlCaHLTBitMon_QueryRunRegistry::string, triggerResultsFUToken, triggerResultsLabel_, triggerResultsToken, triggerSummaryFUToken, triggerSummaryLabel_, triggerSummaryToken, verbose_, and vertexToken.

                                                                              : isSetup_(false) {
   LogDebug("JetMETHLTOfflineSource") << "constructor....";
 
   //
   dirname_ = iConfig.getUntrackedParameter("dirname", std::string("HLT/JetMET/"));
   processname_ = iConfig.getParameter<std::string>("processname");
   triggerSummaryLabel_ = iConfig.getParameter<edm::InputTag>("triggerSummaryLabel");
   triggerResultsLabel_ = iConfig.getParameter<edm::InputTag>("triggerResultsLabel");
   triggerSummaryToken = consumes<trigger::TriggerEvent>(triggerSummaryLabel_);
   triggerResultsToken = consumes<edm::TriggerResults>(triggerResultsLabel_);
   triggerSummaryFUToken = consumes<trigger::TriggerEvent>(
       edm::InputTag(triggerSummaryLabel_.label(), triggerSummaryLabel_.instance(), std::string("FU")));
   triggerResultsFUToken = consumes<edm::TriggerResults>(
       edm::InputTag(triggerResultsLabel_.label(), triggerResultsLabel_.instance(), std::string("FU")));
   //
   verbose_ = iConfig.getUntrackedParameter<bool>("verbose", false);
   runStandalone_ = iConfig.getUntrackedParameter<bool>("runStandalone", false);
   //
   plotAll_ = iConfig.getUntrackedParameter<bool>("plotAll", true);
   plotEff_ = iConfig.getUntrackedParameter<bool>("plotEff", true);
   nameForEff_ = iConfig.getUntrackedParameter<bool>("nameForEff", true);
   MuonTrigPaths_ = iConfig.getUntrackedParameter<vector<std::string> >("pathnameMuon");
   MBTrigPaths_ = iConfig.getUntrackedParameter<vector<std::string> >("pathnameMB");
   //CaloJet, CaloMET
   caloJetsToken = consumes<reco::CaloJetCollection>(iConfig.getParameter<edm::InputTag>("CaloJetCollectionLabel"));
   caloMetToken = consumes<reco::CaloMETCollection>(iConfig.getParameter<edm::InputTag>("CaloMETCollectionLabel"));
   //PFJet, PFMET
   pfJetsToken = consumes<reco::PFJetCollection>(iConfig.getParameter<edm::InputTag>("PFJetCollectionLabel"));
   pfMetToken = consumes<reco::PFMETCollection>(iConfig.getParameter<edm::InputTag>("PFMETCollectionLabel"));
   //pfmhtTag_       = iConfig.getParameter<edm::InputTag>("PFMHTCollectionLabel");
   //Vertex info
   vertexToken = consumes<reco::VertexCollection>(std::string("offlinePrimaryVertices"));
   //
   CaloJetCorToken_ = consumes<reco::JetCorrector>(iConfig.getParameter<edm::InputTag>("CaloJetCorLabel"));
   PFJetCorToken_ = consumes<reco::JetCorrector>(iConfig.getParameter<edm::InputTag>("PFJetCorLabel"));
   //JetID
   jetID = new reco::helper::JetIDHelper(iConfig.getParameter<ParameterSet>("JetIDParams"), consumesCollector());
   _fEMF = iConfig.getUntrackedParameter<double>("fEMF", 0.01);
   _feta = iConfig.getUntrackedParameter<double>("feta", 2.60);
   _fHPD = iConfig.getUntrackedParameter<double>("fHPD", 0.98);
   _n90Hits = iConfig.getUntrackedParameter<double>("n90Hits", 1.0);
   _min_NHEF = iConfig.getUntrackedParameter<double>("minNHEF", 0.);
   _max_NHEF = iConfig.getUntrackedParameter<double>("maxNHEF", 0.99);
   _min_CHEF = iConfig.getUntrackedParameter<double>("minCHEF", 0.);
   _max_CHEF = iConfig.getUntrackedParameter<double>("maxCHEF", 1.);
   _min_NEMF = iConfig.getUntrackedParameter<double>("minNEMF", 0.);
   _max_NEMF = iConfig.getUntrackedParameter<double>("maxNEMF", 0.99);
   _min_CEMF = iConfig.getUntrackedParameter<double>("minCEMF", 0.);
   _max_CEMF = iConfig.getUntrackedParameter<double>("maxCEMF", 0.99);
   //Paths
   pathFilter_ = iConfig.getUntrackedParameter<vector<std::string> >("pathFilter");
   pathRejectKeyword_ = iConfig.getUntrackedParameter<vector<std::string> >("pathRejectKeyword");
   std::vector<edm::ParameterSet> paths = iConfig.getParameter<std::vector<edm::ParameterSet> >("pathPairs");
   for (auto& path : paths) {
     custompathnamepairs_.push_back(
         make_pair(path.getParameter<std::string>("pathname"), path.getParameter<std::string>("denompathname")));
   }
 }

◆ ~JetMETHLTOfflineSource()

JetMETHLTOfflineSource::~JetMETHLTOfflineSource ( )

override

Definition at line 98 of file JetMETHLTOfflineSource.cc.

References jetID.

                                                 {
   //
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
   delete jetID;
 }

Member Function Documentation

◆ analyze()

void JetMETHLTOfflineSource::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 106 of file JetMETHLTOfflineSource.cc.

References reco::helper::JetIDHelper::calculate(), calojet, calojetColl_, CaloJetCorToken_, CaloJetEMF, CaloJetEta, CaloJetfHPD, CaloJetn90, CaloJetPhi, CaloJetPt, CaloJetPx, CaloJetPy, caloJetsToken, calometColl_, caloMetToken, reco::JetCorrector::correction(), gather_cfg::cout, reco::helper::JetIDHelper::fHPD(), fillMEforEffAllTrigger(), fillMEforMonAllTrigger(), fillMEforMonTriggerSummary(), fillMEforTriggerNTfired(), mps_fire::i, iEvent, edm::HandleBase::isValid(), jetID, reco::helper::JetIDHelper::n90Hits(), pfjet, PFJetCEMF, PFJetCHEF, pfjetColl_, PFJetCorToken_, PFJetEta, PFJetNEMF, PFJetNHEF, PFJetPhi, PFJetPt, PFJetPx, PFJetPy, pfJetsToken, pfmetColl_, pfMetToken, pfMHTx_All, pfMHTy_All, plotAll_, plotEff_, runStandalone_, pfClustersFromCombinedCaloHF_cfi::scale, triggerNames_, triggerObj_, triggerResults_, triggerResultsFUToken, triggerResultsToken, triggerSummaryFUToken, triggerSummaryToken, and verbose_.

                                                                                         {
   if (verbose_) {
     cout << endl;
     cout << "============================================================" << endl;
     cout << " New event" << endl << endl;
   }
 
   //---------- triggerResults ----------
   iEvent.getByToken(triggerResultsToken, triggerResults_);
   if (!triggerResults_.isValid()) {
     iEvent.getByToken(triggerResultsFUToken, triggerResults_);
     if (!triggerResults_.isValid()) {
       if (verbose_)
         cout << " triggerResults not valid" << endl;
       edm::LogInfo("JetMETHLTOfflineSource") << "TriggerResults not found, "
                                                 "skipping event";
       return;
     }
   }
   if (verbose_)
     cout << " done triggerResults" << endl;
 
   //---------- triggerResults ----------
   triggerNames_ = iEvent.triggerNames(*triggerResults_);
 
   //---------- triggerSummary ----------
   iEvent.getByToken(triggerSummaryToken, triggerObj_);
   if (!triggerObj_.isValid()) {
     iEvent.getByToken(triggerSummaryFUToken, triggerObj_);
     if (!triggerObj_.isValid()) {
       edm::LogInfo("JetMETHLTOfflineSource") << "TriggerEvent not found, "
                                                 "skipping event";
       return;
     }
   }
   if (verbose_)
     cout << " done triggerSummary" << endl;
 
   if (verbose_) {
     cout << endl;
     cout << "============================================================" << endl;
     cout << " Reading in offline objects" << endl << endl;
   }
 
   //------------ Offline Objects -------
   iEvent.getByToken(caloJetsToken, calojetColl_);
   if (!calojetColl_.isValid())
     return;
   calojet = *calojetColl_;
   //std::stable_sort( calojet.begin(), calojet.end(), PtSorter() );
 
   if (verbose_)
     cout << " done calo" << endl;
 
   iEvent.getByToken(pfJetsToken, pfjetColl_);
   if (!pfjetColl_.isValid())
     return;
   pfjet = *pfjetColl_;
   //std::stable_sort( pfjet.begin(), pfjet.end(), PtSorter() );
 
   if (verbose_)
     cout << " done pf" << endl;
 
   iEvent.getByToken(caloMetToken, calometColl_);
   if (!calometColl_.isValid())
     return;
 
   iEvent.getByToken(pfMetToken, pfmetColl_);
   if (!pfmetColl_.isValid())
     return;
 
   if (verbose_) {
     cout << endl;
     cout << "============================================================" << endl;
     cout << " Read in offline objects" << endl << endl;
   }
 
   //---------- Event counting (DEBUG) ----------
   if (verbose_ && iEvent.id().event() % 10000 == 0)
     cout << "Run = " << iEvent.id().run() << ", LS = " << iEvent.luminosityBlock()
          << ", Event = " << iEvent.id().event() << endl;
 
   //Define on-the-fly correction Jet
   for (int i = 0; i < 2; i++) {
     CaloJetPx[i] = 0.;
     CaloJetPy[i] = 0.;
     CaloJetPt[i] = 0.;
     CaloJetEta[i] = 0.;
     CaloJetPhi[i] = 0.;
     CaloJetEMF[i] = 0.;
     CaloJetfHPD[i] = 0.;
     CaloJetn90[i] = 0.;
     PFJetPx[i] = 0.;
     PFJetPy[i] = 0.;
     PFJetPt[i] = 0.;
     PFJetEta[i] = 0.;
     PFJetPhi[i] = 0.;
     PFJetNHEF[i] = 0.;
     PFJetCHEF[i] = 0.;
     PFJetNEMF[i] = 0.;
     PFJetCEMF[i] = 0.;
   }
 
   //---------- CaloJet Correction (on-the-fly) ----------
   edm::Handle<reco::JetCorrector> calocorrector;
   iEvent.getByToken(CaloJetCorToken_, calocorrector);
   auto calojet_ = calojet.begin();
   for (; calojet_ != calojet.end(); ++calojet_) {
     double scale = calocorrector->correction(*calojet_);
     jetID->calculate(iEvent, iSetup, *calojet_);
 
     if (scale * calojet_->pt() > CaloJetPt[0]) {
       CaloJetPt[1] = CaloJetPt[0];
       CaloJetPx[1] = CaloJetPx[0];
       CaloJetPy[1] = CaloJetPy[0];
       CaloJetEta[1] = CaloJetEta[0];
       CaloJetPhi[1] = CaloJetPhi[0];
       CaloJetEMF[1] = CaloJetEMF[0];
       CaloJetfHPD[1] = CaloJetfHPD[0];
       CaloJetn90[1] = CaloJetn90[0];
       //
       CaloJetPt[0] = scale * calojet_->pt();
       CaloJetPx[0] = scale * calojet_->px();
       CaloJetPy[0] = scale * calojet_->py();
       CaloJetEta[0] = calojet_->eta();
       CaloJetPhi[0] = calojet_->phi();
       CaloJetEMF[0] = calojet_->emEnergyFraction();
       CaloJetfHPD[0] = jetID->fHPD();
       CaloJetn90[0] = jetID->n90Hits();
     } else if (scale * calojet_->pt() < CaloJetPt[0] && scale * calojet_->pt() > CaloJetPt[1]) {
       CaloJetPt[1] = scale * calojet_->pt();
       CaloJetPx[1] = scale * calojet_->px();
       CaloJetPy[1] = scale * calojet_->py();
       CaloJetEta[1] = calojet_->eta();
       CaloJetPhi[1] = calojet_->phi();
       CaloJetEMF[1] = calojet_->emEnergyFraction();
       CaloJetfHPD[1] = jetID->fHPD();
       CaloJetn90[1] = jetID->n90Hits();
     } else {
     }
   }
 
   //---------- PFJet Correction (on-the-fly) ----------
   pfMHTx_All = 0.;
   pfMHTy_All = 0.;
   edm::Handle<reco::JetCorrector> pfcorrector;
   iEvent.getByToken(PFJetCorToken_, pfcorrector);
   auto pfjet_ = pfjet.begin();
   for (; pfjet_ != pfjet.end(); ++pfjet_) {
     double scale = pfcorrector->correction(*pfjet_);
     pfMHTx_All = pfMHTx_All + scale * pfjet_->px();
     pfMHTy_All = pfMHTy_All + scale * pfjet_->py();
     if (scale * pfjet_->pt() > PFJetPt[0]) {
       PFJetPt[1] = PFJetPt[0];
       PFJetPx[1] = PFJetPx[0];
       PFJetPy[1] = PFJetPy[0];
       PFJetEta[1] = PFJetEta[0];
       PFJetPhi[1] = PFJetPhi[0];
       PFJetNHEF[1] = PFJetNHEF[0];
       PFJetCHEF[1] = PFJetCHEF[0];
       PFJetNEMF[1] = PFJetNEMF[0];
       PFJetCEMF[1] = PFJetCEMF[0];
       //
       PFJetPt[0] = scale * pfjet_->pt();
       PFJetPx[0] = scale * pfjet_->px();
       PFJetPy[0] = scale * pfjet_->py();
       PFJetEta[0] = pfjet_->eta();
       PFJetPhi[0] = pfjet_->phi();
       PFJetNHEF[0] = pfjet_->neutralHadronEnergyFraction();
       PFJetCHEF[0] = pfjet_->chargedHadronEnergyFraction();
       PFJetNEMF[0] = pfjet_->neutralEmEnergyFraction();
       PFJetCEMF[0] = pfjet_->chargedEmEnergyFraction();
     } else if (scale * pfjet_->pt() < PFJetPt[0] && scale * pfjet_->pt() > PFJetPt[1]) {
       PFJetPt[1] = scale * pfjet_->pt();
       PFJetPx[1] = scale * pfjet_->px();
       PFJetPy[1] = scale * pfjet_->py();
       PFJetEta[1] = pfjet_->eta();
       PFJetPhi[1] = pfjet_->phi();
       PFJetNHEF[1] = pfjet_->neutralHadronEnergyFraction();
       PFJetCHEF[1] = pfjet_->chargedHadronEnergyFraction();
       PFJetNEMF[1] = pfjet_->neutralEmEnergyFraction();
       PFJetCEMF[1] = pfjet_->chargedEmEnergyFraction();
     } else {
     }
   }
 
   if (verbose_) {
     for (int i = 0; i < 2; i++) {
       cout << "CaloJet-0: " << CaloJetPt[i] << ", Eta = " << CaloJetEta[i] << ", Phi = " << CaloJetPhi[i] << endl;
       cout << "fHPD = " << CaloJetfHPD[0] << ", n90 = " << CaloJetn90[0] << endl;
     }
     for (int i = 0; i < 2; i++) {
       cout << "PFJet-0: " << PFJetPt[i] << ", Eta = " << PFJetEta[i] << ", Phi = " << PFJetPhi[i] << endl;
     }
   }
 
   //---------- RUN ----------
   fillMEforMonTriggerSummary(iEvent, iSetup);
   if (plotAll_)
     fillMEforMonAllTrigger(iEvent, iSetup);
   if (plotEff_)
     fillMEforEffAllTrigger(iEvent, iSetup);
   if (runStandalone_)
     fillMEforTriggerNTfired();
 }

◆ bookHistograms()

void JetMETHLTOfflineSource::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	run,
		edm::EventSetup const &	c
	)

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 1622 of file JetMETHLTOfflineSource.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookFloat(), correlation_All, correlation_AllWrtMB, correlation_AllWrtMu, TrackerOfflineValidation_Dqm_cff::dirName, dirname_, dqm::impl::MonitorElement::getTH1(), hltPathsAll_, hltPathsAllTriggerSummary_, hltPathsEff_, mps_fire::i, isSetup_, N, LaserClient_cfi::nbins, NVertices, VtxSmearedParameters_cfi::Phi, plotAll_, plotEff_, PVZ, rate_All, rate_AllWrtMB, rate_AllWrtMu, HLTConfigProvider::removeVersion(), runStandalone_, dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, runGCPTkAlMap::title, trigger::TriggerJet, trigger::TriggerMET, HLTObjectsMonitor_cfi::TriggerSummary, trigger::TriggerTET, and findQualityFiles::v.

                                                                                                                {
   if (!isSetup_) {
     iBooker.setCurrentFolder(dirname_);
 
     //-----------------------------------------------------------------
     //---book trigger summary histos
     if (!isSetup_) {
       std::string foldernm = "/TriggerSummary/";
       iBooker.setCurrentFolder(dirname_ + foldernm);
 
       int TrigBins_ = hltPathsAllTriggerSummary_.size();
       double TrigMin_ = -0.5;
       double TrigMax_ = hltPathsAllTriggerSummary_.size() - 0.5;
 
       std::string histonm = "JetMET_TriggerRate";
       std::string histot = "JetMET TriggerRate Summary";
       rate_All = iBooker.book1D(histonm.c_str(), histot.c_str(), TrigBins_, TrigMin_, TrigMax_);
 
       histonm = "JetMET_TriggerRate_Correlation";
       histot = "JetMET TriggerRate Correlation Summary;y&&!x;x&&y";
       correlation_All =
           iBooker.book2D(histonm.c_str(), histot.c_str(), TrigBins_, TrigMin_, TrigMax_, TrigBins_, TrigMin_, TrigMax_);
 
       histonm = "JetMET_NVertices";
       histot = "No. of vertices";
       NVertices = iBooker.book1D(histonm.c_str(), histot.c_str(), 100, 0, 50);
 
       histonm = "JetMET_PVZ";
       histot = "Primary Vertex Z pos";
       PVZ = iBooker.book1D(histonm.c_str(), histot.c_str(), 100, -50., 50.);
 
       if (runStandalone_) {
         histonm = "JetMET_TriggerRate_WrtMuTrigger";
         histot = "JetMET TriggerRate Summary Wrt Muon Trigger ";
         rate_AllWrtMu = iBooker.book1D(histonm.c_str(), histot.c_str(), TrigBins_, TrigMin_, TrigMax_);
 
         histonm = "JetMET_TriggerRate_Correlation_WrtMuTrigger";
         histot = "JetMET TriggerRate Correlation Summary Wrt Muon Trigger;y&&!x;x&&y";
         correlation_AllWrtMu = iBooker.book2D(
             histonm.c_str(), histot.c_str(), TrigBins_, TrigMin_, TrigMax_, TrigBins_, TrigMin_, TrigMax_);
 
         histonm = "JetMET_TriggerRate_WrtMBTrigger";
         histot = "JetMET TriggerRate Summary Wrt MB Trigger";
         rate_AllWrtMB = iBooker.book1D(histonm.c_str(), histot.c_str(), TrigBins_, TrigMin_, TrigMax_);
 
         histonm = "JetMET_TriggerRate_Correlation_WrtMBTrigger";
         histot = "JetMET TriggerRate Correlation Wrt MB Trigger;y&&!x;x&&y";
         correlation_AllWrtMB = iBooker.book2D(
             histonm.c_str(), histot.c_str(), TrigBins_, TrigMin_, TrigMax_, TrigBins_, TrigMin_, TrigMax_);
       }
       isSetup_ = true;
     }
     //---Set bin label
 
     for (auto& v : hltPathsAllTriggerSummary_) {
       std::string labelnm("dummy");
       labelnm = v.getPath();
       int nbins = rate_All->getTH1()->GetNbinsX();
       for (int ibin = 1; ibin < nbins + 1; ibin++) {
         const char* binLabel = rate_All->getTH1()->GetXaxis()->GetBinLabel(ibin);
         std::string binLabel_str = string(binLabel);
         if (binLabel_str == labelnm)
           break;
         if (binLabel[0] == '\0') {
           rate_All->setBinLabel(ibin, labelnm);
           correlation_All->setBinLabel(ibin, labelnm, 1);
           correlation_All->setBinLabel(ibin, labelnm, 2);
           if (runStandalone_) {
             rate_AllWrtMu->setBinLabel(ibin, labelnm);
             rate_AllWrtMB->setBinLabel(ibin, labelnm);
             correlation_AllWrtMu->setBinLabel(ibin, labelnm, 1);
             correlation_AllWrtMB->setBinLabel(ibin, labelnm, 1);
             correlation_AllWrtMu->setBinLabel(ibin, labelnm, 2);
             correlation_AllWrtMB->setBinLabel(ibin, labelnm, 2);
           }
           break;
         }
       }
     }
 
     // Now define histos for All triggers
     if (plotAll_) {
       //
       int Nbins_ = 10;
       double Nmin_ = -0.5;
       double Nmax_ = 9.5;
       //
       int Ptbins_ = 100;
       if (runStandalone_)
         Ptbins_ = 1000;
       double PtMin_ = 0.;
       double PtMax_ = 1000.;
       //
       int Etabins_ = 50;
       if (runStandalone_)
         Etabins_ = 100;
       double EtaMin_ = -5.;
       double EtaMax_ = 5.;
       //
       int Phibins_ = 35;
       double PhiMin_ = -3.5;
       double PhiMax_ = 3.5;
 
       int Resbins_ = 30;
       double ResMin_ = -1.5;
       double ResMax_ = 1.5;
       //
       std::string dirName = dirname_ + "/MonitorAllTriggers/";
       for (auto& v : hltPathsAll_) {
         //
         std::string trgPathName = HLTConfigProvider::removeVersion(v.getPath());
         std::string subdirName = dirName + trgPathName;
         std::string trigPath = "(" + trgPathName + ")";
         iBooker.setCurrentFolder(subdirName);
 
         std::string labelname("ME");
         std::string histoname(labelname + "");
         std::string title(labelname + "");
 
         MonitorElement* dummy;
         dummy = iBooker.bookFloat("dummy");
 
         if (v.getObjectType() == trigger::TriggerJet && v.getTriggerType() == "SingleJet_Trigger") {
           histoname = labelname + "_recObjN";
           title = labelname + "_recObjN;Reco multiplicity()" + trigPath;
           MonitorElement* N = iBooker.book1D(histoname.c_str(), title.c_str(), Nbins_, Nmin_, Nmax_);
           N->getTH1();
 
           histoname = labelname + "_recObjPt";
           title = labelname + "_recObjPt; Reco Pt[GeV/c]" + trigPath;
           MonitorElement* Pt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           Pt->getTH1();
 
           histoname = labelname + "_recObjPtBarrel";
           title = labelname + "_recObjPtBarrel;Reco Pt[GeV/c]" + trigPath;
           MonitorElement* PtBarrel = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           PtBarrel->getTH1();
 
           histoname = labelname + "_recObjPtEndcap";
           title = labelname + "_recObjPtEndcap;Reco Pt[GeV/c]" + trigPath;
           MonitorElement* PtEndcap = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           PtEndcap->getTH1();
 
           histoname = labelname + "_recObjPtForward";
           title = labelname + "_recObjPtForward;Reco Pt[GeV/c]" + trigPath;
           MonitorElement* PtForward = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           PtForward->getTH1();
 
           histoname = labelname + "_recObjEta";
           title = labelname + "_recObjEta;Reco #eta" + trigPath;
           MonitorElement* Eta = iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           Eta->getTH1();
 
           histoname = labelname + "_recObjPhi";
           title = labelname + "_recObjPhi;Reco #Phi" + trigPath;
           MonitorElement* Phi = iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           Phi->getTH1();
 
           histoname = labelname + "_recObjEtaPhi";
           title = labelname + "_recObjEtaPhi;Reco #eta;Reco #Phi" + trigPath;
           MonitorElement* EtaPhi =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           EtaPhi->getTH1();
 
           histoname = labelname + "_l1ObjPt";
           title = labelname + "_l1ObjPt;L1 Pt[GeV/c]" + trigPath;
           MonitorElement* Pt_L1 = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           Pt_L1->getTH1();
 
           histoname = labelname + "_l1ObjEta";
           title = labelname + "_l1ObjEta;L1 #eta" + trigPath;
           MonitorElement* Eta_L1 = iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           Eta_L1->getTH1();
 
           histoname = labelname + "_l1ObjPhi";
           title = labelname + "_l1ObjPhi;L1 #Phi" + trigPath;
           MonitorElement* Phi_L1 = iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           Phi_L1->getTH1();
 
           histoname = labelname + "_l1ObjEtaPhi";
           title = labelname + "_l1ObjEtaPhi;L1 #eta;L1 #Phi" + trigPath;
           MonitorElement* EtaPhi_L1 =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           EtaPhi_L1->getTH1();
 
           histoname = labelname + "_l1ObjN";
           title = labelname + "_l1ObjN;L1 multiplicity" + trigPath;
           MonitorElement* N_L1 = iBooker.book1D(histoname.c_str(), title.c_str(), Nbins_, Nmin_, Nmax_);
           N_L1->getTH1();
 
           histoname = labelname + "_l1ObjPtBarrel";
           title = labelname + "_l1ObjPtBarrel;L1 Pt[GeV/c]" + trigPath;
           MonitorElement* PtBarrel_L1 = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           PtBarrel_L1->getTH1();
 
           histoname = labelname + "_l1ObjPtEndcap";
           title = labelname + "_l1ObjPtEndcap;L1 Pt[GeV/c]" + trigPath;
           MonitorElement* PtEndcap_L1 = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           PtEndcap_L1->getTH1();
 
           histoname = labelname + "_l1ObjPtForward";
           title = labelname + "_l1ObjPtForward;L1 Pt[GeV/c]" + trigPath;
           MonitorElement* PtForward_L1 = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           PtForward_L1->getTH1();
 
           histoname = labelname + "_hltObjN";
           title = labelname + "_hltObjN;HLT multiplicity" + trigPath;
           MonitorElement* N_HLT = iBooker.book1D(histoname.c_str(), title.c_str(), Nbins_, Nmin_, Nmax_);
           N_HLT->getTH1();
 
           histoname = labelname + "_hltObjPtBarrel";
           title = labelname + "_hltObjPtBarrel;HLT Pt[GeV/c]" + trigPath;
           MonitorElement* PtBarrel_HLT = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           PtBarrel_HLT->getTH1();
 
           histoname = labelname + "_hltObjPtEndcap";
           title = labelname + "_hltObjPtEndcap;HLT Pt[GeV/c]" + trigPath;
           MonitorElement* PtEndcap_HLT = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           PtEndcap_HLT->getTH1();
 
           histoname = labelname + "_hltObjPtForward";
           title = labelname + "_hltObjPtForward;HLT Pt[GeV/c]" + trigPath;
           MonitorElement* PtForward_HLT = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           PtForward_HLT->getTH1();
 
           histoname = labelname + "_hltObjPt";
           title = labelname + "_hltObjPt;HLT Pt[GeV/c]" + trigPath;
           MonitorElement* Pt_HLT = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           Pt_HLT->getTH1();
 
           histoname = labelname + "_hltObjEta";
           title = labelname + "_hltObjEta;HLT #eta" + trigPath;
           MonitorElement* Eta_HLT = iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           Eta_HLT->getTH1();
 
           histoname = labelname + "_hltObjPhi";
           title = labelname + "_hltObjPhi;HLT #Phi" + trigPath;
           MonitorElement* Phi_HLT = iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           Phi_HLT->getTH1();
 
           histoname = labelname + "_hltObjEtaPhi";
           title = labelname + "_hltObjEtaPhi;HLT #eta;HLT #Phi" + trigPath;
           MonitorElement* EtaPhi_HLT =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           EtaPhi_HLT->getTH1();
 
           histoname = labelname + "_l1HLTPtResolution";
           title = labelname + "_l1HLTPtResolution;(Pt(L1)-Pt(HLT))/Pt(L1)" + trigPath;
           MonitorElement* PtResolution_L1HLT =
               iBooker.book1D(histoname.c_str(), title.c_str(), Resbins_, ResMin_, ResMax_);
           PtResolution_L1HLT->getTH1();
 
           histoname = labelname + "_l1HLTEtaResolution";
           title = labelname + "_l1HLTEtaResolution;#eta(L1)-#eta(HLT)" + trigPath;
           MonitorElement* EtaResolution_L1HLT =
               iBooker.book1D(histoname.c_str(), title.c_str(), Resbins_, ResMin_, ResMax_);
           EtaResolution_L1HLT->getTH1();
 
           histoname = labelname + "_l1HLTPhiResolution";
           title = labelname + "_l1HLTPhiResolution;#Phi(L1)-#Phi(HLT)" + trigPath;
           MonitorElement* PhiResolution_L1HLT =
               iBooker.book1D(histoname.c_str(), title.c_str(), Resbins_, ResMin_, ResMax_);
           PhiResolution_L1HLT->getTH1();
 
           histoname = labelname + "_l1HLTPtCorrelation";
           title = labelname + "_l1HLTPtCorrelation;Pt(L1)[GeV/c];Pt(HLT)[GeV/c]" + trigPath;
           MonitorElement* PtCorrelation_L1HLT =
               iBooker.book2D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_, Ptbins_, PtMin_, PtMax_);
           PtCorrelation_L1HLT->getTH1();
 
           histoname = labelname + "_l1HLTEtaCorrelation";
           title = labelname + "_l1HLTEtaCorrelation;#eta(L1);#eta(HLT)" + trigPath;
           MonitorElement* EtaCorrelation_L1HLT =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Etabins_, EtaMin_, EtaMax_);
           EtaCorrelation_L1HLT->getTH1();
 
           histoname = labelname + "_l1HLTPhiCorrelation";
           title = labelname + "_l1HLTPhiCorrelation;#Phi(L1);#Phi(HLT)" + trigPath;
           MonitorElement* PhiCorrelation_L1HLT =
               iBooker.book2D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_, Phibins_, PhiMin_, PhiMax_);
           PhiCorrelation_L1HLT->getTH1();
 
           histoname = labelname + "_hltRecObjPtResolution";
           title = labelname + "_hltRecObjPtResolution;(Pt(HLT)-Pt(Reco))/Pt(HLT)" + trigPath;
           MonitorElement* PtResolution_HLTRecObj =
               iBooker.book1D(histoname.c_str(), title.c_str(), Resbins_, ResMin_, ResMax_);
           PtResolution_HLTRecObj->getTH1();
 
           histoname = labelname + "_hltRecObjEtaResolution";
           title = labelname + "_hltRecObjEtaResolution;#eta(HLT)-#eta(Reco)" + trigPath;
           MonitorElement* EtaResolution_HLTRecObj =
               iBooker.book1D(histoname.c_str(), title.c_str(), Resbins_, ResMin_, ResMax_);
           EtaResolution_HLTRecObj->getTH1();
 
           histoname = labelname + "_hltRecObjPhiResolution";
           title = labelname + "_hltRecObjPhiResolution;#Phi(HLT)-#Phi(Reco)" + trigPath;
           MonitorElement* PhiResolution_HLTRecObj =
               iBooker.book1D(histoname.c_str(), title.c_str(), Resbins_, ResMin_, ResMax_);
           PhiResolution_HLTRecObj->getTH1();
 
           histoname = labelname + "_hltRecObjPtCorrelation";
           title = labelname + "_hltRecObjPtCorrelation;Pt(HLT)[GeV/c];Pt(Reco)[GeV/c]" + trigPath;
           MonitorElement* PtCorrelation_HLTRecObj =
               iBooker.book2D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_, Ptbins_, PtMin_, PtMax_);
           PtCorrelation_HLTRecObj->getTH1();
 
           histoname = labelname + "_hltRecObjEtaCorrelation";
           title = labelname + "_hltRecObjEtaCorrelation;#eta(HLT);#eta(Reco)" + trigPath;
           MonitorElement* EtaCorrelation_HLTRecObj =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Etabins_, EtaMin_, EtaMax_);
           EtaCorrelation_HLTRecObj->getTH1();
 
           histoname = labelname + "_hltRecObjPhiCorrelation";
           title = labelname + "_hltRecObjPhiCorrelation;#Phi(HLT);#Phi(Reco)" + trigPath;
           MonitorElement* PhiCorrelation_HLTRecObj =
               iBooker.book2D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_, Phibins_, PhiMin_, PhiMax_);
           PhiCorrelation_HLTRecObj->getTH1();
 
           v.setHistos(N,
                       Pt,
                       PtBarrel,
                       PtEndcap,
                       PtForward,
                       Eta,
                       Phi,
                       EtaPhi,
                       N_L1,
                       Pt_L1,
                       PtBarrel_L1,
                       PtEndcap_L1,
                       PtForward_L1,
                       Eta_L1,
                       Phi_L1,
                       EtaPhi_L1,
                       N_HLT,
                       Pt_HLT,
                       PtBarrel_HLT,
                       PtEndcap_HLT,
                       PtForward_HLT,
                       Eta_HLT,
                       Phi_HLT,
                       EtaPhi_HLT,
                       PtResolution_L1HLT,
                       EtaResolution_L1HLT,
                       PhiResolution_L1HLT,
                       PtResolution_HLTRecObj,
                       EtaResolution_HLTRecObj,
                       PhiResolution_HLTRecObj,
                       PtCorrelation_L1HLT,
                       EtaCorrelation_L1HLT,
                       PhiCorrelation_L1HLT,
                       PtCorrelation_HLTRecObj,
                       EtaCorrelation_HLTRecObj,
                       PhiCorrelation_HLTRecObj,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy);
         }  // histos for SingleJet Triggers
 
         if (v.getObjectType() == trigger::TriggerJet && v.getTriggerType() == "DiJet_Trigger") {
           histoname = labelname + "_RecObjAveragePt";
           title = labelname + "_RecObjAveragePt;Reco Average Pt[GeV/c]" + trigPath;
           MonitorElement* jetAveragePt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           jetAveragePt->getTH1();
 
           histoname = labelname + "_RecObjAverageEta";
           title = labelname + "_RecObjAverageEta;Reco Average #eta" + trigPath;
           MonitorElement* jetAverageEta = iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           jetAverageEta->getTH1();
 
           histoname = labelname + "_RecObjPhiDifference";
           title = labelname + "_RecObjPhiDifference;Reco #Delta#Phi" + trigPath;
           MonitorElement* jetPhiDifference =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           jetPhiDifference->getTH1();
 
           histoname = labelname + "_hltObjAveragePt";
           title = labelname + "_hltObjAveragePt;HLT Average Pt[GeV/c]" + trigPath;
           MonitorElement* hltAveragePt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           hltAveragePt->getTH1();
 
           histoname = labelname + "_hltObjAverageEta";
           title = labelname + "_hltObjAverageEta;HLT Average #eta" + trigPath;
           MonitorElement* hltAverageEta = iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           hltAverageEta->getTH1();
 
           histoname = labelname + "_hltObjPhiDifference";
           title = labelname + "_hltObjPhiDifference;Reco #Delta#Phi" + trigPath;
           MonitorElement* hltPhiDifference =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           hltPhiDifference->getTH1();
 
           v.setHistos(dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       jetAveragePt,
                       jetAverageEta,
                       jetPhiDifference,
                       hltAveragePt,
                       hltAverageEta,
                       hltPhiDifference,
                       dummy,
                       dummy,
                       dummy);
         }  // histos for DiJet Triggers
 
         if (v.getObjectType() == trigger::TriggerMET || (v.getObjectType() == trigger::TriggerTET)) {
           histoname = labelname + "_recObjPt";
           title = labelname + "_recObjPt;Reco Pt[GeV/c]" + trigPath;
           MonitorElement* Pt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           Pt->getTH1();
 
           histoname = labelname + "_recObjPhi";
           title = labelname + "_recObjPhi;Reco #Phi" + trigPath;
           MonitorElement* Phi = iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           Phi->getTH1();
 
           histoname = labelname + "_l1ObjPt";
           title = labelname + "_l1ObjPt;L1 Pt[GeV/c]" + trigPath;
           MonitorElement* Pt_L1 = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           Pt_L1->getTH1();
 
           histoname = labelname + "_l1ObjPhi";
           title = labelname + "_l1ObjPhi;L1 #Phi" + trigPath;
           MonitorElement* Phi_L1 = iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           Phi_L1->getTH1();
 
           histoname = labelname + "_hltObjPt";
           title = labelname + "_hltObjPt;HLT Pt[GeV/c]" + trigPath;
           MonitorElement* Pt_HLT = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           Pt_HLT->getTH1();
 
           histoname = labelname + "_hltObjPhi";
           title = labelname + "_hltObjPhi;HLT #Phi" + trigPath;
           MonitorElement* Phi_HLT = iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           Phi_HLT->getTH1();
 
           histoname = labelname + "_l1HLTPtResolution";
           title = labelname + "_l1HLTPtResolution;(Pt(L1)-Pt(HLT))/Pt(L1)" + trigPath;
           MonitorElement* PtResolution_L1HLT =
               iBooker.book1D(histoname.c_str(), title.c_str(), Resbins_, ResMin_, ResMax_);
           PtResolution_L1HLT->getTH1();
 
           histoname = labelname + "_l1HLTPhiResolution";
           title = labelname + "_l1HLTPhiResolution;#Phi(L1)-#Phi(HLT)" + trigPath;
           MonitorElement* PhiResolution_L1HLT =
               iBooker.book1D(histoname.c_str(), title.c_str(), Resbins_, ResMin_, ResMax_);
           PhiResolution_L1HLT->getTH1();
 
           histoname = labelname + "_l1HLTPtCorrelation";
           title = labelname + "_l1HLTPtCorrelation;Pt(L1)[GeV/c];Pt(HLT)[GeV/c]" + trigPath;
           MonitorElement* PtCorrelation_L1HLT =
               iBooker.book2D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_, Ptbins_, PtMin_, PtMax_);
           PtCorrelation_L1HLT->getTH1();
 
           histoname = labelname + "_l1HLTPhiCorrelation";
           title = labelname + "_l1HLTPhiCorrelation;#Phi(L1);#Phi(HLT)" + trigPath;
           MonitorElement* PhiCorrelation_L1HLT =
               iBooker.book2D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_, Phibins_, PhiMin_, PhiMax_);
           PhiCorrelation_L1HLT->getTH1();
 
           histoname = labelname + "_hltRecObjPtResolution";
           title = labelname + "_hltRecObjPtResolution;(Pt(HLT)-Pt(Reco))/Pt(HLT)" + trigPath;
           MonitorElement* PtResolution_HLTRecObj =
               iBooker.book1D(histoname.c_str(), title.c_str(), Resbins_, ResMin_, ResMax_);
           PtResolution_HLTRecObj->getTH1();
 
           histoname = labelname + "_hltRecObjPhiResolution";
           title = labelname + "_hltRecObjPhiResolution;#Phi(HLT)-#Phi(Reco)" + trigPath;
           MonitorElement* PhiResolution_HLTRecObj =
               iBooker.book1D(histoname.c_str(), title.c_str(), Resbins_, ResMin_, ResMax_);
           PhiResolution_HLTRecObj->getTH1();
 
           histoname = labelname + "_hltRecObjPtCorrelation";
           title = labelname + "_hltRecObjPtCorrelation;Pt(HLT)[GeV/c];Pt(Reco)[GeV/c]" + trigPath;
           MonitorElement* PtCorrelation_HLTRecObj =
               iBooker.book2D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_, Ptbins_, PtMin_, PtMax_);
           PtCorrelation_HLTRecObj->getTH1();
 
           histoname = labelname + "_hltRecObjPhiCorrelation";
           title = labelname + "_hltRecObjPhiCorrelation;#Phi(HLT);#Phi(Reco)" + trigPath;
           MonitorElement* PhiCorrelation_HLTRecObj =
               iBooker.book2D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_, Phibins_, PhiMin_, PhiMax_);
           PhiCorrelation_HLTRecObj->getTH1();
 
           v.setHistos(dummy,
                       Pt,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       Phi,
                       dummy,
                       dummy,
                       Pt_L1,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       Phi_L1,
                       dummy,
                       dummy,
                       Pt_HLT,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       Phi_HLT,
                       dummy,
                       PtResolution_L1HLT,
                       dummy,
                       PhiResolution_L1HLT,
                       PtResolution_HLTRecObj,
                       dummy,
                       PhiResolution_HLTRecObj,
                       PtCorrelation_L1HLT,
                       dummy,
                       PhiCorrelation_L1HLT,
                       PtCorrelation_HLTRecObj,
                       dummy,
                       PhiCorrelation_HLTRecObj,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy,
                       dummy);
         }  // histos for MET Triggers
       }
     }  //plotAll_
 
     //-------Now Efficiency histos--------
     if (plotEff_) {
       int Ptbins_ = 100;
       if (runStandalone_)
         Ptbins_ = 1000;
       double PtMin_ = 0.;
       double PtMax_ = 1000.;
       //
       int Etabins_ = 50;
       double EtaMin_ = -5.;
       double EtaMax_ = 5.;
       //
       int Phibins_ = 35;
       double PhiMin_ = -3.5;
       double PhiMax_ = 3.5;
       // Now define histos wrt lower threshold trigger
       std::string dirName1 = dirname_ + "/RelativeTriggerEff/";
       for (auto& v : hltPathsEff_) {
         //
         std::string trgPathName = HLTConfigProvider::removeVersion(v.getPath());
         std::string trgPathNameD = HLTConfigProvider::removeVersion(v.getDenomPath());
         //
         std::string labelname("ME");
         std::string subdirName = dirName1 + trgPathName + "_wrt_" + trgPathNameD;
         iBooker.setCurrentFolder(subdirName);
         //
         std::string histoname(labelname + "");
         std::string title(labelname + "");
 
         MonitorElement* dummy;
         dummy = iBooker.bookFloat("dummy");
 
         if ((v.getObjectType() == trigger::TriggerJet) && (v.getTriggerType() == "SingleJet_Trigger")) {
           histoname = labelname + "_NumeratorPt";
           title = labelname + "NumeratorPt;Calo Pt[GeV/c]";
           MonitorElement* NumeratorPt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           NumeratorPt->getTH1();
 
           histoname = labelname + "_NumeratorPtBarrel";
           title = labelname + "NumeratorPtBarrel;Calo Pt[GeV/c] ";
           MonitorElement* NumeratorPtBarrel = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           NumeratorPtBarrel->getTH1();
 
           histoname = labelname + "_NumeratorPtEndcap";
           title = labelname + "NumeratorPtEndcap;Calo Pt[GeV/c]";
           MonitorElement* NumeratorPtEndcap = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           NumeratorPtEndcap->getTH1();
 
           histoname = labelname + "_NumeratorPtForward";
           title = labelname + "NumeratorPtForward;Calo Pt[GeV/c]";
           MonitorElement* NumeratorPtForward =
               iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           NumeratorPtForward->getTH1();
 
           histoname = labelname + "_NumeratorEta";
           title = labelname + "NumeratorEta;Calo #eta ";
           MonitorElement* NumeratorEta = iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorEta->getTH1();
 
           histoname = labelname + "_NumeratorPhi";
           title = labelname + "NumeratorPhi;Calo #Phi";
           MonitorElement* NumeratorPhi = iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPhi->getTH1();
 
           histoname = labelname + "_NumeratorEtaPhi";
           title = labelname + "NumeratorEtaPhi;Calo #eta;Calo #Phi";
           MonitorElement* NumeratorEtaPhi =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           NumeratorEtaPhi->getTH1();
 
           histoname = labelname + "_NumeratorEtaBarrel";
           title = labelname + "NumeratorEtaBarrel;Calo #eta ";
           MonitorElement* NumeratorEtaBarrel =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorEtaBarrel->getTH1();
 
           histoname = labelname + "_NumeratorPhiBarrel";
           title = labelname + "NumeratorPhiBarrel;Calo #Phi";
           MonitorElement* NumeratorPhiBarrel =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPhiBarrel->getTH1();
 
           histoname = labelname + "_NumeratorEtaEndcap";
           title = labelname + "NumeratorEtaEndcap;Calo #eta ";
           MonitorElement* NumeratorEtaEndcap =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorEtaEndcap->getTH1();
 
           histoname = labelname + "_NumeratorPhiEndcap";
           title = labelname + "NumeratorPhiEndcap;Calo #Phi";
           MonitorElement* NumeratorPhiEndcap =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPhiEndcap->getTH1();
 
           histoname = labelname + "_NumeratorEtaForward";
           title = labelname + "NumeratorEtaForward;Calo #eta ";
           MonitorElement* NumeratorEtaForward =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorEtaForward->getTH1();
 
           histoname = labelname + "_NumeratorPhiForward";
           title = labelname + "NumeratorPhiForward;Calo #Phi";
           MonitorElement* NumeratorPhiForward =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPhiForward->getTH1();
 
           histoname = labelname + "_NumeratorEta_LowpTcut";
           title = labelname + "NumeratorEta_LowpTcut;Calo #eta ";
           MonitorElement* NumeratorEta_LowpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorEta_LowpTcut->getTH1();
 
           histoname = labelname + "_NumeratorPhi_LowpTcut";
           title = labelname + "NumeratorPhi_LowpTcut;Calo #Phi";
           MonitorElement* NumeratorPhi_LowpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPhi_LowpTcut->getTH1();
 
           histoname = labelname + "_NumeratorEtaPhi_LowpTcut";
           title = labelname + "NumeratorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
           MonitorElement* NumeratorEtaPhi_LowpTcut =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           NumeratorEtaPhi_LowpTcut->getTH1();
 
           histoname = labelname + "_NumeratorEta_MedpTcut";
           title = labelname + "NumeratorEta_MedpTcut;Calo #eta ";
           MonitorElement* NumeratorEta_MedpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorEta_MedpTcut->getTH1();
 
           histoname = labelname + "_NumeratorPhi_MedpTcut";
           title = labelname + "NumeratorPhi_MedpTcut;Calo #Phi";
           MonitorElement* NumeratorPhi_MedpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPhi_MedpTcut->getTH1();
 
           histoname = labelname + "_NumeratorEtaPhi_MedpTcut";
           title = labelname + "NumeratorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
           MonitorElement* NumeratorEtaPhi_MedpTcut =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           NumeratorEtaPhi_MedpTcut->getTH1();
 
           histoname = labelname + "_NumeratorEta_HighpTcut";
           title = labelname + "NumeratorEta_HighpTcut;Calo #eta ";
           MonitorElement* NumeratorEta_HighpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorEta_HighpTcut->getTH1();
 
           histoname = labelname + "_NumeratorPhi_HighpTcut";
           title = labelname + "NumeratorPhi_HighpTcut;Calo #Phi";
           MonitorElement* NumeratorPhi_HighpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPhi_HighpTcut->getTH1();
 
           histoname = labelname + "_NumeratorEtaPhi_HighpTcut";
           title = labelname + "NumeratorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
           MonitorElement* NumeratorEtaPhi_HighpTcut =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           NumeratorEtaPhi_HighpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPt";
           title = labelname + "DenominatorPt;Calo Pt[GeV/c]";
           MonitorElement* DenominatorPt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           DenominatorPt->getTH1();
 
           histoname = labelname + "_DenominatorPtBarrel";
           title = labelname + "DenominatorPtBarrel;Calo Pt[GeV/c]";
           MonitorElement* DenominatorPtBarrel =
               iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           DenominatorPtBarrel->getTH1();
 
           histoname = labelname + "_DenominatorPtEndcap";
           title = labelname + "DenominatorPtEndcap;Calo Pt[GeV/c]";
           MonitorElement* DenominatorPtEndcap =
               iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           DenominatorPtEndcap->getTH1();
 
           histoname = labelname + "_DenominatorPtForward";
           title = labelname + "DenominatorPtForward;Calo Pt[GeV/c] ";
           MonitorElement* DenominatorPtForward =
               iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           DenominatorPtForward->getTH1();
 
           histoname = labelname + "_DenominatorEta";
           title = labelname + "DenominatorEta;Calo #eta ";
           MonitorElement* DenominatorEta = iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorEta->getTH1();
 
           histoname = labelname + "_DenominatorPhi";
           title = labelname + "DenominatorPhi;Calo #Phi";
           MonitorElement* DenominatorPhi = iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPhi->getTH1();
 
           histoname = labelname + "_DenominatorEtaPhi";
           title = labelname + "DenominatorEtaPhi;Calo #eta; Calo #Phi";
           MonitorElement* DenominatorEtaPhi =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           DenominatorEtaPhi->getTH1();
 
           histoname = labelname + "_DenominatorEtaBarrel";
           title = labelname + "DenominatorEtaBarrel;Calo #eta ";
           MonitorElement* DenominatorEtaBarrel =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorEtaBarrel->getTH1();
 
           histoname = labelname + "_DenominatorPhiBarrel";
           title = labelname + "DenominatorPhiBarrel;Calo #Phi";
           MonitorElement* DenominatorPhiBarrel =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPhiBarrel->getTH1();
 
           histoname = labelname + "_DenominatorEtaEndcap";
           title = labelname + "DenominatorEtaEndcap;Calo #eta ";
           MonitorElement* DenominatorEtaEndcap =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorEtaEndcap->getTH1();
 
           histoname = labelname + "_DenominatorPhiEndcap";
           title = labelname + "DenominatorPhiEndcap;Calo #Phi";
           MonitorElement* DenominatorPhiEndcap =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPhiEndcap->getTH1();
 
           histoname = labelname + "_DenominatorEtaForward";
           title = labelname + "DenominatorEtaForward;Calo #eta ";
           MonitorElement* DenominatorEtaForward =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorEtaForward->getTH1();
 
           histoname = labelname + "_DenominatorPhiForward";
           title = labelname + "DenominatorPhiForward;Calo #Phi";
           MonitorElement* DenominatorPhiForward =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPhiForward->getTH1();
 
           histoname = labelname + "_DenominatorEta_LowpTcut";
           title = labelname + "DenominatorEta_LowpTcut;Calo #eta ";
           MonitorElement* DenominatorEta_LowpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorEta_LowpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPhi_LowpTcut";
           title = labelname + "DenominatorPhi_LowpTcut;Calo #Phi";
           MonitorElement* DenominatorPhi_LowpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPhi_LowpTcut->getTH1();
 
           histoname = labelname + "_DenominatorEtaPhi_LowpTcut";
           title = labelname + "DenominatorEtaPhi_LowpTcut;Calo #eta;Calo #Phi";
           MonitorElement* DenominatorEtaPhi_LowpTcut =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           DenominatorEtaPhi_LowpTcut->getTH1();
 
           histoname = labelname + "_DenominatorEta_MedpTcut";
           title = labelname + "DenominatorEta_MedpTcut;Calo #eta ";
           MonitorElement* DenominatorEta_MedpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorEta_MedpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPhi_MedpTcut";
           title = labelname + "DenominatorPhi_MedpTcut;Calo #Phi";
           MonitorElement* DenominatorPhi_MedpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPhi_MedpTcut->getTH1();
 
           histoname = labelname + "_DenominatorEtaPhi_MedpTcut";
           title = labelname + "DenominatorEtaPhi_MedpTcut;Calo #eta;Calo #Phi";
           MonitorElement* DenominatorEtaPhi_MedpTcut =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           DenominatorEtaPhi_MedpTcut->getTH1();
 
           histoname = labelname + "_DenominatorEta_HighpTcut";
           title = labelname + "DenominatorEta_HighpTcut;Calo #eta ";
           MonitorElement* DenominatorEta_HighpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorEta_HighpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPhi_HighpTcut";
           title = labelname + "DenominatorPhi_HighpTcut;Calo #Phi";
           MonitorElement* DenominatorPhi_HighpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPhi_HighpTcut->getTH1();
 
           histoname = labelname + "_DenominatorEtaPhi_HighpTcut";
           title = labelname + "DenominatorEtaPhi_HighpTcut;Calo #eta;Calo #Phi";
           MonitorElement* DenominatorEtaPhi_HighpTcut =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           DenominatorEtaPhi_HighpTcut->getTH1();
 
           histoname = labelname + "_DeltaR";
           title = labelname + "DeltaR;";
           MonitorElement* DeltaR = iBooker.book1D(histoname.c_str(), title.c_str(), 100, 0., 1.5);
           DeltaR->getTH1();
 
           histoname = labelname + "_DeltaPhi";
           title = labelname + "DeltaPhi;";
           MonitorElement* DeltaPhi = iBooker.book1D(histoname.c_str(), title.c_str(), 500, -5.0, 5.0);
           DeltaPhi->getTH1();
 
           histoname = labelname + "_NumeratorPFMHT";
           title = labelname + "NumeratorPFMHT;PFMHT[GeV/c]";
           MonitorElement* NumeratorPFMHT = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           NumeratorPFMHT->getTH1();
 
           histoname = labelname + "_NumeratorPFPt";
           title = labelname + "NumeratorPFPt;PF Pt[GeV/c]";
           MonitorElement* NumeratorPFPt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           NumeratorPFPt->getTH1();
 
           histoname = labelname + "_NumeratorPFPtBarrel";
           title = labelname + "NumeratorPFPtBarrel;PF Pt[GeV/c] ";
           MonitorElement* NumeratorPFPtBarrel =
               iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           NumeratorPFPtBarrel->getTH1();
 
           histoname = labelname + "_NumeratorPFPtEndcap";
           title = labelname + "NumeratorPFPtEndcap;PF Pt[GeV/c]";
           MonitorElement* NumeratorPFPtEndcap =
               iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           NumeratorPFPtEndcap->getTH1();
 
           histoname = labelname + "_NumeratorPFPtForward";
           title = labelname + "NumeratorPFPtForward;PF Pt[GeV/c]";
           MonitorElement* NumeratorPFPtForward =
               iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           NumeratorPFPtForward->getTH1();
 
           histoname = labelname + "_NumeratorPFEta";
           title = labelname + "NumeratorPFEta;PF #eta ";
           MonitorElement* NumeratorPFEta = iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorPFEta->getTH1();
 
           histoname = labelname + "_NumeratorPFPhi";
           title = labelname + "NumeratorPFPhi;Calo #Phi";
           MonitorElement* NumeratorPFPhi = iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPFPhi->getTH1();
 
           histoname = labelname + "_NumeratorPFEtaPhi";
           title = labelname + "NumeratorPFEtaPhi;PF #eta;Calo #Phi";
           MonitorElement* NumeratorPFEtaPhi =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           NumeratorPFEtaPhi->getTH1();
 
           histoname = labelname + "_NumeratorPFEtaBarrel";
           title = labelname + "NumeratorPFEtaBarrel;PF #eta ";
           MonitorElement* NumeratorPFEtaBarrel =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorPFEtaBarrel->getTH1();
 
           histoname = labelname + "_NumeratorPFPhiBarrel";
           title = labelname + "NumeratorPFPhiBarrel;PF #Phi";
           MonitorElement* NumeratorPFPhiBarrel =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPFPhiBarrel->getTH1();
 
           histoname = labelname + "_NumeratorPFEtaEndcap";
           title = labelname + "NumeratorPFEtaEndcap;Calo #eta ";
           MonitorElement* NumeratorPFEtaEndcap =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorPFEtaEndcap->getTH1();
 
           histoname = labelname + "_NumeratorPFPhiEndcap";
           title = labelname + "NumeratorPFPhiEndcap;PF #Phi";
           MonitorElement* NumeratorPFPhiEndcap =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPFPhiEndcap->getTH1();
 
           histoname = labelname + "_NumeratorPFEtaForward";
           title = labelname + "NumeratorPFEtaForward;Calo #eta ";
           MonitorElement* NumeratorPFEtaForward =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorPFEtaForward->getTH1();
 
           histoname = labelname + "_NumeratorPFPhiForward";
           title = labelname + "NumeratorPFPhiForward;PF #Phi";
           MonitorElement* NumeratorPFPhiForward =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPFPhiForward->getTH1();
 
           histoname = labelname + "_NumeratorPFEta_LowpTcut";
           title = labelname + "NumeratorPFEta_LowpTcut;PF #eta ";
           MonitorElement* NumeratorPFEta_LowpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorPFEta_LowpTcut->getTH1();
 
           histoname = labelname + "_NumeratorPFPhi_LowpTcut";
           title = labelname + "NumeratorPFPhi_LowpTcut;PF #Phi";
           MonitorElement* NumeratorPFPhi_LowpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPFPhi_LowpTcut->getTH1();
 
           histoname = labelname + "_NumeratorPFEtaPhi_LowpTcut";
           title = labelname + "NumeratorPFEtaPhi_LowpTcut;PF #eta;Calo #Phi";
           MonitorElement* NumeratorPFEtaPhi_LowpTcut =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           NumeratorPFEtaPhi_LowpTcut->getTH1();
 
           histoname = labelname + "_NumeratorPFEta_MedpTcut";
           title = labelname + "NumeratorPFEta_MedpTcut;PF #eta ";
           MonitorElement* NumeratorPFEta_MedpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorPFEta_MedpTcut->getTH1();
 
           histoname = labelname + "_NumeratorPFPhi_MedpTcut";
           title = labelname + "NumeratorPFPhi_MedpTcut;PF #Phi";
           MonitorElement* NumeratorPFPhi_MedpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPFPhi_MedpTcut->getTH1();
 
           histoname = labelname + "_NumeratorPFEtaPhi_MedpTcut";
           title = labelname + "NumeratorPFEtaPhi_MedpTcut;PF #eta;PF #Phi";
           MonitorElement* NumeratorPFEtaPhi_MedpTcut =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           NumeratorPFEtaPhi_MedpTcut->getTH1();
 
           histoname = labelname + "_NumeratorPFEta_HighpTcut";
           title = labelname + "NumeratorPFEta_HighpTcut;Calo #eta ";
           MonitorElement* NumeratorPFEta_HighpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorPFEta_HighpTcut->getTH1();
 
           histoname = labelname + "_NumeratorPFPhi_HighpTcut";
           title = labelname + "NumeratorPFPhi_HighpTcut;PF #Phi";
           MonitorElement* NumeratorPFPhi_HighpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPFPhi_HighpTcut->getTH1();
 
           histoname = labelname + "_NumeratorPFEtaPhi_HighpTcut";
           title = labelname + "NumeratorPFEtaPhi_HighpTcut;PF #eta;PF #Phi";
           MonitorElement* NumeratorPFEtaPhi_HighpTcut =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           NumeratorPFEtaPhi_HighpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPFMHT";
           title = labelname + "DenominatorPFMHT;PF Pt[GeV/c]";
           MonitorElement* DenominatorPFMHT = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           DenominatorPFMHT->getTH1();
 
           histoname = labelname + "_DenominatorPFPt";
           title = labelname + "DenominatorPFPt;PF Pt[GeV/c]";
           MonitorElement* DenominatorPFPt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           DenominatorPFPt->getTH1();
 
           histoname = labelname + "_DenominatorPFPtBarrel";
           title = labelname + "DenominatorPFPtBarrel;Calo Pt[GeV/c]";
           MonitorElement* DenominatorPFPtBarrel =
               iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           DenominatorPFPtBarrel->getTH1();
 
           histoname = labelname + "_DenominatorPFPtEndcap";
           title = labelname + "DenominatorPFPtEndcap;PF Pt[GeV/c]";
           MonitorElement* DenominatorPFPtEndcap =
               iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           DenominatorPFPtEndcap->getTH1();
 
           histoname = labelname + "_DenominatorPFPtForward";
           title = labelname + "DenominatorPFPtForward;PF Pt[GeV/c] ";
           MonitorElement* DenominatorPFPtForward =
               iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           DenominatorPFPtForward->getTH1();
 
           histoname = labelname + "_DenominatorPFEta";
           title = labelname + "DenominatorPFEta;PF #eta ";
           MonitorElement* DenominatorPFEta =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorPFEta->getTH1();
 
           histoname = labelname + "_DenominatorPFPhi";
           title = labelname + "DenominatorPFPhi;PF #Phi";
           MonitorElement* DenominatorPFPhi =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPFPhi->getTH1();
 
           histoname = labelname + "_DenominatorPFEtaPhi";
           title = labelname + "DenominatorPFEtaPhi;PF #eta; Calo #Phi";
           MonitorElement* DenominatorPFEtaPhi =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           DenominatorPFEtaPhi->getTH1();
 
           histoname = labelname + "_DenominatorPFEtaBarrel";
           title = labelname + "DenominatorPFEtaBarrel;Calo #eta ";
           MonitorElement* DenominatorPFEtaBarrel =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorPFEtaBarrel->getTH1();
 
           histoname = labelname + "_DenominatorPFPhiBarrel";
           title = labelname + "DenominatorPFPhiBarrel;PF #Phi";
           MonitorElement* DenominatorPFPhiBarrel =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPFPhiBarrel->getTH1();
 
           histoname = labelname + "_DenominatorPFEtaEndcap";
           title = labelname + "DenominatorPFEtaEndcap;PF #eta ";
           MonitorElement* DenominatorPFEtaEndcap =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorPFEtaEndcap->getTH1();
 
           histoname = labelname + "_DenominatorPFPhiEndcap";
           title = labelname + "DenominatorPFPhiEndcap;Calo #Phi";
           MonitorElement* DenominatorPFPhiEndcap =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPFPhiEndcap->getTH1();
 
           histoname = labelname + "_DenominatorPFEtaForward";
           title = labelname + "DenominatorPFEtaForward;PF #eta ";
           MonitorElement* DenominatorPFEtaForward =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorPFEtaForward->getTH1();
 
           histoname = labelname + "_DenominatorPFPhiForward";
           title = labelname + "DenominatorPFPhiForward;PF #Phi";
           MonitorElement* DenominatorPFPhiForward =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPFPhiForward->getTH1();
 
           histoname = labelname + "_DenominatorPFEta_LowpTcut";
           title = labelname + "DenominatorPFEta_LowpTcut;PF #eta ";
           MonitorElement* DenominatorPFEta_LowpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorPFEta_LowpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPFPhi_LowpTcut";
           title = labelname + "DenominatorPFPhi_LowpTcut;PF #Phi";
           MonitorElement* DenominatorPFPhi_LowpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPFPhi_LowpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPFEtaPhi_LowpTcut";
           title = labelname + "DenominatorPFEtaPhi_LowpTcut;PF #eta;Calo #Phi";
           MonitorElement* DenominatorPFEtaPhi_LowpTcut =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           DenominatorPFEtaPhi_LowpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPFEta_MedpTcut";
           title = labelname + "DenominatorPFEta_MedpTcut;PF #eta ";
           MonitorElement* DenominatorPFEta_MedpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorPFEta_MedpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPFPhi_MedpTcut";
           title = labelname + "DenominatorPFPhi_MedpTcut;PF #Phi";
           MonitorElement* DenominatorPFPhi_MedpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPFPhi_MedpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPFEtaPhi_MedpTcut";
           title = labelname + "DenominatorPFEtaPhi_MedpTcut;PF #eta;Calo #Phi";
           MonitorElement* DenominatorPFEtaPhi_MedpTcut =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           DenominatorPFEtaPhi_MedpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPFEta_HighpTcut";
           title = labelname + "DenominatorPFEta_HighpTcut;PF #eta ";
           MonitorElement* DenominatorPFEta_HighpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorPFEta_HighpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPFPhi_HighpTcut";
           title = labelname + "DenominatorPFPhi_HighpTcut;PF #Phi";
           MonitorElement* DenominatorPFPhi_HighpTcut =
               iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPFPhi_HighpTcut->getTH1();
 
           histoname = labelname + "_DenominatorPFEtaPhi_HighpTcut";
           title = labelname + "DenominatorPFEtaPhi_HighpTcut;PF #eta;Calo #Phi";
           MonitorElement* DenominatorPFEtaPhi_HighpTcut =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Phibins_, PhiMin_, PhiMax_);
           DenominatorPFEtaPhi_HighpTcut->getTH1();
 
           histoname = labelname + "_PFDeltaR";
           title = labelname + "PFDeltaR;";
           MonitorElement* PFDeltaR = iBooker.book1D(histoname.c_str(), title.c_str(), 100, 0., 1.5);
           PFDeltaR->getTH1();
 
           histoname = labelname + "_PFDeltaPhi";
           title = labelname + "PFDeltaPhi;";
           MonitorElement* PFDeltaPhi = iBooker.book1D(histoname.c_str(), title.c_str(), 500, -5.0, 5.0);
           PFDeltaPhi->getTH1();
 
           v.setEffHistos(NumeratorPt,
                          NumeratorPtBarrel,
                          NumeratorPtEndcap,
                          NumeratorPtForward,
                          NumeratorEta,
                          NumeratorPhi,
                          NumeratorEtaPhi,
                          //
                          NumeratorEtaBarrel,
                          NumeratorPhiBarrel,
                          NumeratorEtaEndcap,
                          NumeratorPhiEndcap,
                          NumeratorEtaForward,
                          NumeratorPhiForward,
                          NumeratorEta_LowpTcut,
                          NumeratorPhi_LowpTcut,
                          NumeratorEtaPhi_LowpTcut,
                          NumeratorEta_MedpTcut,
                          NumeratorPhi_MedpTcut,
                          NumeratorEtaPhi_MedpTcut,
                          NumeratorEta_HighpTcut,
                          NumeratorPhi_HighpTcut,
                          NumeratorEtaPhi_HighpTcut,
                          //
                          DenominatorPt,
                          DenominatorPtBarrel,
                          DenominatorPtEndcap,
                          DenominatorPtForward,
                          DenominatorEta,
                          DenominatorPhi,
                          DenominatorEtaPhi,
                          //
                          DenominatorEtaBarrel,
                          DenominatorPhiBarrel,
                          DenominatorEtaEndcap,
                          DenominatorPhiEndcap,
                          DenominatorEtaForward,
                          DenominatorPhiForward,
                          DenominatorEta_LowpTcut,
                          DenominatorPhi_LowpTcut,
                          DenominatorEtaPhi_LowpTcut,
                          DenominatorEta_MedpTcut,
                          DenominatorPhi_MedpTcut,
                          DenominatorEtaPhi_MedpTcut,
                          DenominatorEta_HighpTcut,
                          DenominatorPhi_HighpTcut,
                          DenominatorEtaPhi_HighpTcut,
                          DeltaR,
                          DeltaPhi,
                          //
                          NumeratorPFPt,
                          NumeratorPFMHT,
                          NumeratorPFPtBarrel,
                          NumeratorPFPtEndcap,
                          NumeratorPFPtForward,
                          NumeratorPFEta,
                          NumeratorPFPhi,
                          NumeratorPFEtaPhi,
                          NumeratorPFEtaBarrel,
                          NumeratorPFPhiBarrel,
                          NumeratorPFEtaEndcap,
                          NumeratorPFPhiEndcap,
                          NumeratorPFEtaForward,
                          NumeratorPFPhiForward,
                          NumeratorPFEta_LowpTcut,
                          NumeratorPFPhi_LowpTcut,
                          NumeratorPFEtaPhi_LowpTcut,
                          NumeratorPFEta_MedpTcut,
                          NumeratorPFPhi_MedpTcut,
                          NumeratorPFEtaPhi_MedpTcut,
                          NumeratorPFEta_HighpTcut,
                          NumeratorPFPhi_HighpTcut,
                          NumeratorPFEtaPhi_HighpTcut,
                          DenominatorPFPt,
                          DenominatorPFMHT,
                          DenominatorPFPtBarrel,
                          DenominatorPFPtEndcap,
                          DenominatorPFPtForward,
                          DenominatorPFEta,
                          DenominatorPFPhi,
                          DenominatorPFEtaPhi,
                          DenominatorPFEtaBarrel,
                          DenominatorPFPhiBarrel,
                          DenominatorPFEtaEndcap,
                          DenominatorPFPhiEndcap,
                          DenominatorPFEtaForward,
                          DenominatorPFPhiForward,
                          DenominatorPFEta_LowpTcut,
                          DenominatorPFPhi_LowpTcut,
                          DenominatorPFEtaPhi_LowpTcut,
                          DenominatorPFEta_MedpTcut,
                          DenominatorPFPhi_MedpTcut,
                          DenominatorPFEtaPhi_MedpTcut,
                          DenominatorPFEta_HighpTcut,
                          DenominatorPFPhi_HighpTcut,
                          DenominatorPFEtaPhi_HighpTcut,
                          PFDeltaR,
                          PFDeltaPhi);
 
         }  // Loop over Jet Trigger
 
         if ((v.getObjectType() == trigger::TriggerJet) && (v.getTriggerType() == "DiJet_Trigger")) {
           histoname = labelname + "_NumeratorAvrgPt";
           title = labelname + "NumeratorAvrgPt;Calo Pt[GeV/c]";
           MonitorElement* NumeratorPt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           NumeratorPt->getTH1();
 
           histoname = labelname + "_NumeratorAvrgEta";
           title = labelname + "NumeratorAvrgEta;Calo #eta";
           MonitorElement* NumeratorEta = iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorEta->getTH1();
 
           histoname = labelname + "_DenominatorAvrgPt";
           title = labelname + "DenominatorAvrgPt;Calo Pt[GeV/c] ";
           MonitorElement* DenominatorPt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           DenominatorPt->getTH1();
 
           histoname = labelname + "_DenominatorAvrgEta";
           title = labelname + "DenominatorAvrgEta;Calo #eta";
           MonitorElement* DenominatorEta = iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorEta->getTH1();
 
           histoname = labelname + "_DeltaR";
           title = labelname + "DeltaR;";
           MonitorElement* DeltaR = iBooker.book1D(histoname.c_str(), title.c_str(), 100, 0., 1.5);
           DeltaR->getTH1();
 
           histoname = labelname + "_DeltaPhi";
           title = labelname + "DeltaPhi;";
           MonitorElement* DeltaPhi = iBooker.book1D(histoname.c_str(), title.c_str(), 500, -5., 5.);
           DeltaPhi->getTH1();
 
           //add PF histo: SJ
           histoname = labelname + "_NumeratorAvrgPFPt";
           title = labelname + "NumeratorAvrgPFPt;PF Pt[GeV/c]";
           MonitorElement* NumeratorPFPt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           NumeratorPFPt->getTH1();
 
           histoname = labelname + "_NumeratorAvrgPFEta";
           title = labelname + "NumeratorAvrgPFEta;PF #eta";
           MonitorElement* NumeratorPFEta = iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           NumeratorPFEta->getTH1();
 
           histoname = labelname + "_DenominatorAvrgPFPt";
           title = labelname + "DenominatorAvrgPFPt;PF Pt[GeV/c] ";
           MonitorElement* DenominatorPFPt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           DenominatorPFPt->getTH1();
 
           histoname = labelname + "_DenominatorAvrgPFEta";
           title = labelname + "DenominatorAvrgPFEta;PF #eta";
           MonitorElement* DenominatorPFEta =
               iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           DenominatorPFEta->getTH1();
 
           histoname = labelname + "_PFDeltaR";
           title = labelname + "PFDeltaR;";
           MonitorElement* PFDeltaR = iBooker.book1D(histoname.c_str(), title.c_str(), 100, 0., 1.5);
           PFDeltaR->getTH1();
 
           histoname = labelname + "_PFDeltaPhi";
           title = labelname + "PFDeltaPhi;";
           MonitorElement* PFDeltaPhi = iBooker.book1D(histoname.c_str(), title.c_str(), 500, -5., 5.);
           PFDeltaPhi->getTH1();
 
           v.setEffHistos(dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy);
         }
 
         if (v.getObjectType() == trigger::TriggerMET || (v.getObjectType() == trigger::TriggerTET)) {
           histoname = labelname + "_NumeratorPt";
           if (v.getPath().find("HLT_PFMET") == std::string::npos)
             title = labelname + "NumeratorPt; CaloMET[GeV/c]";
           else
             title = labelname + "NumeratorPt; PFMET[GeV/c]";
           MonitorElement* NumeratorPt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           NumeratorPt->getTH1();
 
           histoname = labelname + "_NumeratorPhi";
           title = labelname + "NumeratorPhi; #Phi";
           MonitorElement* NumeratorPhi = iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           NumeratorPhi->getTH1();
 
           histoname = labelname + "_DenominatorPt";
           if (v.getPath().find("HLT_PFMET") == std::string::npos)
             title = labelname + "DenominatorPt; CaloMET[GeV/c]";
           else
             title = labelname + "DenominatorPt; PFMET[GeV/c]";
           MonitorElement* DenominatorPt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           DenominatorPt->getTH1();
 
           histoname = labelname + "_DenominatorPhi";
           title = labelname + "DenominatorPhi; #Phi";
           MonitorElement* DenominatorPhi = iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           DenominatorPhi->getTH1();
 
           v.setEffHistos(NumeratorPt,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          NumeratorPhi,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          DenominatorPt,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          DenominatorPhi,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy,
                          dummy);
         }  // Loop over MET trigger
       }
     }  //plotEff_
 
     if (runStandalone_) {  //runStandalone
       //--------Histos to see WHY trigger is NOT fired----------
       int Nbins_ = 10;
       int Nmin_ = 0;
       int Nmax_ = 10;
       int Ptbins_ = 1000;
       int Etabins_ = 40;
       int Phibins_ = 35;
       double PtMin_ = 0.;
       double PtMax_ = 1000.;
       double EtaMin_ = -5.;
       double EtaMax_ = 5.;
       double PhiMin_ = -3.14159;
       double PhiMax_ = 3.14159;
 
       std::string dirName4_ = dirname_ + "/TriggerNotFired/";
       //      iBooker.setCurrentFolder(dirName4);
 
       for (auto& v : hltPathsAll_) {
         MonitorElement* dummy;
         dummy = iBooker.bookFloat("dummy");
 
         std::string labelname("ME");
         std::string histoname(labelname + "");
         std::string title(labelname + "");
         iBooker.setCurrentFolder(dirName4_ + v.getPath());
 
         histoname = labelname + "_TriggerSummary";
         title = labelname + "Summary of trigger levels";
         MonitorElement* TriggerSummary = iBooker.book1D(histoname.c_str(), title.c_str(), 7, -0.5, 6.5);
 
         std::vector<std::string> trigger;
         trigger.emplace_back("Nevt");
         trigger.emplace_back("L1 failed");
         trigger.emplace_back("L1 & HLT failed");
         trigger.emplace_back("L1 failed but not HLT");
         trigger.emplace_back("L1 passed");
         trigger.emplace_back("L1 & HLT passed");
         trigger.emplace_back("L1 passed but not HLT");
 
         for (unsigned int i = 0; i < trigger.size(); i++)
           TriggerSummary->setBinLabel(i + 1, trigger[i]);
 
         if ((v.getTriggerType() == "SingleJet_Trigger")) {
           histoname = labelname + "_JetPt";
           title = labelname + "Leading jet pT;Pt[GeV/c]";
           MonitorElement* JetPt = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           JetPt->getTH1();
 
           histoname = labelname + "_JetEtaVsPt";
           title = labelname + "Leading jet #eta vs pT;#eta;Pt[GeV/c]";
           MonitorElement* JetEtaVsPt =
               iBooker.book2D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_, Ptbins_, PtMin_, PtMax_);
           JetEtaVsPt->getTH1();
 
           histoname = labelname + "_JetPhiVsPt";
           title = labelname + "Leading jet #Phi vs pT;#Phi;Pt[GeV/c]";
           MonitorElement* JetPhiVsPt =
               iBooker.book2D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_, Ptbins_, PtMin_, PtMax_);
           JetPhiVsPt->getTH1();
 
           v.setDgnsHistos(
               TriggerSummary, dummy, JetPt, JetEtaVsPt, JetPhiVsPt, dummy, dummy, dummy, dummy, dummy, dummy);
         }  // single Jet trigger
 
         if ((v.getTriggerType() == "DiJet_Trigger")) {
           histoname = labelname + "_JetSize";
           title = labelname + "Jet Size;multiplicity";
           MonitorElement* JetSize = iBooker.book1D(histoname.c_str(), title.c_str(), Nbins_, Nmin_, Nmax_);
           JetSize->getTH1();
 
           histoname = labelname + "_AvergPt";
           title = labelname + "Average Pt;Pt[GeV/c]";
           MonitorElement* Pt12 = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           Pt12->getTH1();
 
           histoname = labelname + "_AvergEta";
           title = labelname + "Average Eta;#eta";
           MonitorElement* Eta12 = iBooker.book1D(histoname.c_str(), title.c_str(), Etabins_, EtaMin_, EtaMax_);
           Eta12->getTH1();
 
           histoname = labelname + "_PhiDifference";
           title = labelname + "#Delta#Phi;#Delta#Phi";
           MonitorElement* Phi12 = iBooker.book1D(histoname.c_str(), title.c_str(), Phibins_, PhiMin_, PhiMax_);
           Phi12->getTH1();
 
           histoname = labelname + "_Pt3Jet";
           title = labelname + "Pt of 3rd Jet;Pt[GeV/c]";
           MonitorElement* Pt3 = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           Pt3->getTH1();
 
           histoname = labelname + "_Pt12VsPt3Jet";
           title = labelname + "Pt of 3rd Jet vs Average Pt of leading jets;Avergage Pt[GeV/c]; Pt of 3rd Jet [GeV/c]";
           MonitorElement* Pt12Pt3 =
               iBooker.book2D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_, Ptbins_, PtMin_, PtMax_);
           Pt12Pt3->getTH1();
 
           histoname = labelname + "_Pt12VsPhi12";
           title = labelname +
                   "Average Pt of leading jets vs #Delta#Phi between leading jets;Avergage Pt[GeV/c]; #Delta#Phi";
           MonitorElement* Pt12Phi12 =
               iBooker.book2D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_, Phibins_, PhiMin_, PhiMax_);
           Pt12Phi12->getTH1();
 
           v.setDgnsHistos(TriggerSummary, JetSize, dummy, dummy, dummy, Pt12, Eta12, Phi12, Pt3, Pt12Pt3, Pt12Phi12);
         }  // Dijet Jet trigger
 
         if ((v.getTriggerType() == "MET_Trigger")) {
           histoname = labelname + "_MET";
           title = labelname + "MET;Pt[GeV/c]";
           MonitorElement* MET = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           MET->getTH1();
 
           v.setDgnsHistos(TriggerSummary, dummy, MET, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy);
         }  // MET trigger
 
         if ((v.getTriggerType() == "TET_Trigger")) {
           histoname = labelname + "_TET";
           title = labelname + "TET;Pt[GeV/c]";
           MonitorElement* TET = iBooker.book1D(histoname.c_str(), title.c_str(), Ptbins_, PtMin_, PtMax_);
           TET->getTH1();
 
           v.setDgnsHistos(TriggerSummary, dummy, TET, dummy, dummy, dummy, dummy, dummy, dummy, dummy, dummy);
         }  // TET trigger
       }
     }  //runStandalone
   }
 }

◆ dqmBeginRun()

void JetMETHLTOfflineSource::dqmBeginRun	(	edm::Run const &	run,
		edm::EventSetup const &	c
	)

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 1323 of file JetMETHLTOfflineSource.cc.

References HltBtagPostValidation_cff::c, config::checkPath, gather_cfg::cout, custompathnamepairs_, JetMETHLTOfflineSource_cfi::denompathname, getL1ConditionModuleName(), hltConfig_, hltPathsAll_, hltPathsAllTriggerSummary_, hltPathsEff_, mps_fire::i, HLTConfigProvider::init(), isSetup_, LogDebug, BTaggingMonitor_cfi::met, HLTConfigProvider::moduleLabels(), HLTConfigProvider::moduleType(), MuonTrigPaths_, dqmiodumpmetadata::n, nameForEff_, pathFilter_, BTVHLTOfflineSource_cfi::pathname, pathRejectKeyword_, processname_, writedatasetfile::run, HLTConfigProvider::size(), AlCaHLTBitMon_QueryRunRegistry::string, trigger::TriggerJet, trigger::TriggerMET, HLTConfigProvider::triggerName(), trigger::TriggerTET, HLT_2023v12_cff::triggerType, and verbose_.

                                                                                   {
   if (!isSetup_) {
     //--- htlConfig_
     bool changed(true);
     if (!hltConfig_.init(run, c, processname_, changed)) {
       LogDebug("HLTJetMETDQMSource") << "HLTConfigProvider failed to initialize.";
     }
 
     /*
       Here we select the Single Jet, DiJet, MET trigger. SingleJet and DiJet trigger are saved under same object type "TriggerJet". 
       We can easily separate out single and di jet trigger later. For the first trigger in the list, denominator trigger is dummy 
       (empty) whereas for other triggers denom is previous trigger of same type. e.g. SingleJet50 has singleJet30 as denominator.
     */
 
     const unsigned int n(hltConfig_.size());
     int singleJet = 0;
     int diJet = 0;
     int met = 0;
     int tet = 0;
     for (unsigned int i = 0; i != n; ++i) {
       bool denomFound = false;
       bool numFound = false;
       bool checkPath = false;
 
       //Look for paths if "path name fraction" is found in the pathname
       std::string pathname = hltConfig_.triggerName(i);
       //Filter only paths JetMET triggers are interested in
       auto controlPathname = pathFilter_.begin();
       for (; controlPathname != pathFilter_.end(); ++controlPathname) {
         if (pathname.find((*controlPathname)) != std::string::npos) {
           checkPath = true;
           break;
         }
       }
       if (checkPath == false)
         continue;
 
       //Reject if keyword(s) is found in the pathname
       auto rejectPathname = pathRejectKeyword_.begin();
       for (; rejectPathname != pathRejectKeyword_.end(); ++rejectPathname) {
         if (pathname.find((*rejectPathname)) != std::string::npos) {
           checkPath = false;
           break;
         }
       }
       if (checkPath == false)
         continue;
 
       //
       if (verbose_)
         cout << "==pathname==" << pathname << endl;
       std::string dpathname = MuonTrigPaths_[0];
       std::string l1pathname = "dummy";
       std::string denompathname = "";
       unsigned int usedPrescale = 1;
       unsigned int objectType = 0;
       std::string triggerType = "";
       std::string filtername("dummy");
       std::string Denomfiltername("denomdummy");
 
       if (pathname.find("Jet") != std::string::npos && !(pathname.find("DoubleJet") != std::string::npos) &&
           !(pathname.find("DiJet") != std::string::npos) && !(pathname.find("DiPFJet") != std::string::npos) &&
           !(pathname.find("BTag") != std::string::npos) && !(pathname.find("Mu") != std::string::npos) &&
           !(pathname.find("Fwd") != std::string::npos)) {
         triggerType = "SingleJet_Trigger";
         objectType = trigger::TriggerJet;
       }
       if (pathname.find("DiJet") != std::string::npos || pathname.find("DiPFJet") != std::string::npos ||
           pathname.find("DoubleJet") != std::string::npos) {
         triggerType = "DiJet_Trigger";
         objectType = trigger::TriggerJet;
       }
       if (pathname.find("MET") != std::string::npos) {
         triggerType = "MET_Trigger";
         objectType = trigger::TriggerMET;
       }
       if (pathname.find("HT") != std::string::npos) {
         triggerType = "TET_Trigger";
         objectType = trigger::TriggerTET;
       }
 
       //
       if (objectType == trigger::TriggerJet && !(pathname.find("DiJet") != std::string::npos) &&
           !(pathname.find("DiPFJet") != std::string::npos) && !(pathname.find("DoubleJet") != std::string::npos)) {
         singleJet++;
         if (singleJet > 1)
           dpathname = dpathname = hltConfig_.triggerName(i - 1);
         if (singleJet == 1)
           dpathname = MuonTrigPaths_[0];
       }
       if (objectType == trigger::TriggerJet &&
           ((pathname.find("DiJet") != std::string::npos) || (pathname.find("DiPFJet") != std::string::npos))) {
         diJet++;
         if (diJet > 1)
           dpathname = dpathname = hltConfig_.triggerName(i - 1);
         if (diJet == 1)
           dpathname = MuonTrigPaths_[0];
       }
       if (objectType == trigger::TriggerMET) {
         met++;
         if (met > 1)
           dpathname = dpathname = hltConfig_.triggerName(i - 1);
         if (met == 1)
           dpathname = MuonTrigPaths_[0];
       }
       if (objectType == trigger::TriggerTET) {
         tet++;
         if (tet > 1)
           dpathname = dpathname = hltConfig_.triggerName(i - 1);
         if (tet == 1)
           dpathname = MuonTrigPaths_[0];
       }
 
       // find L1 condition for numpath with numpath objecttype
       // find PSet for L1 global seed for numpath,sss
       // list module labels for numpath
 
       // Checking if the trigger exist in HLT table or not
       for (unsigned int i = 0; i != n; ++i) {
         std::string HLTname = hltConfig_.triggerName(i);
         if (HLTname == pathname)
           numFound = true;
         if (HLTname == dpathname)
           denomFound = true;
       }
 
       if (numFound) {  //make trigger exist in the menu
         //ml needs change l1pathname
         l1pathname = getL1ConditionModuleName(pathname);  //ml added L1conditionmodulename
         //ml added
         std::vector<std::string> numpathmodules = hltConfig_.moduleLabels(pathname);
         for (auto& numpathmodule : numpathmodules) {
           edm::InputTag testTag(numpathmodule, "", processname_);
           if ((hltConfig_.moduleType(numpathmodule) == "HLT1CaloJet") ||
               (hltConfig_.moduleType(numpathmodule) == "HLT1PFJet") ||
               (hltConfig_.moduleType(numpathmodule) == "HLTDiJetAveFilter") ||
               (hltConfig_.moduleType(numpathmodule) == "HLTDiPFJetAveFilter") ||
               (hltConfig_.moduleType(numpathmodule) == "HLT1CaloMET") ||
               (hltConfig_.moduleType(numpathmodule) == "HLTMhtFilter") ||
               (hltConfig_.moduleType(numpathmodule) == "HLTPrescaler"))
             filtername = numpathmodule;
         }
       }
 
       if (objectType != 0 && denomFound) {
         std::vector<std::string> numpathmodules = hltConfig_.moduleLabels(dpathname);
         for (auto& numpathmodule : numpathmodules) {
           edm::InputTag testTag(numpathmodule, "", processname_);
           if ((hltConfig_.moduleType(numpathmodule) == "HLT1CaloJet") ||
               (hltConfig_.moduleType(numpathmodule) == "HLT1PFJet") ||
               (hltConfig_.moduleType(numpathmodule) == "HLTDiJetAveFilter") ||
               (hltConfig_.moduleType(numpathmodule) == "HLTDiPFJetAveFilter") ||
               (hltConfig_.moduleType(numpathmodule) == "HLT1CaloMET") ||
               (hltConfig_.moduleType(numpathmodule) == "HLTMhtFilter") ||
               (hltConfig_.moduleType(numpathmodule) == "HLTPrescaler"))
             Denomfiltername = numpathmodule;
         }
       }
 
       if (objectType != 0 && numFound) {
         if (verbose_)
           cout << "=Pathname= " << pathname << " | =Denompath= " << dpathname << " | =Filtername= " << filtername
                << " | =Denomfiltername= " << Denomfiltername << " | =L1pathname= " << l1pathname
                << " | =ObjectType= " << objectType << endl;
         if (!((pathname.find("HT") != std::string::npos) || (pathname.find("Quad") != std::string::npos))) {
           hltPathsAll_.push_back(PathInfo(usedPrescale,
                                           dpathname,
                                           pathname,
                                           l1pathname,
                                           filtername,
                                           Denomfiltername,
                                           processname_,
                                           objectType,
                                           triggerType));
           if (!nameForEff_ && denomFound)
             hltPathsEff_.push_back(PathInfo(usedPrescale,
                                             dpathname,
                                             pathname,
                                             l1pathname,
                                             filtername,
                                             Denomfiltername,
                                             processname_,
                                             objectType,
                                             triggerType));
         }
         hltPathsAllTriggerSummary_.push_back(PathInfo(usedPrescale,
                                                       dpathname,
                                                       pathname,
                                                       l1pathname,
                                                       filtername,
                                                       Denomfiltername,
                                                       processname_,
                                                       objectType,
                                                       triggerType));
       }
     }  //Loop over paths
 
     if (verbose_)
       cout << "get names for efficicncy------------------" << endl;
     //---------bool to pick trigger names pair from config file-------------
     if (nameForEff_) {
       std::string l1pathname = "dummy";
       std::string denompathname = "";
       unsigned int usedPrescale = 1;
       unsigned int objectType = 0;
       std::string triggerType = "";
       std::string filtername("dummy");
       std::string Denomfiltername("denomdummy");
       for (auto& custompathnamepair : custompathnamepairs_) {
         std::string pathname = custompathnamepair.first;
         std::string dpathname = custompathnamepair.second;
         bool numFound = false;
         bool denomFound = false;
         // Checking if the trigger exist in HLT table or not
         for (unsigned int i = 0; i != n; ++i) {
           std::string HLTname = hltConfig_.triggerName(i);
           if (HLTname.find(pathname) != std::string::npos) {
             numFound = true;
             pathname = HLTname;
           }  //changed to get versions
           if (HLTname.find(dpathname) != std::string::npos) {
             denomFound = true;
             dpathname = HLTname;
           }
         }
         if (numFound && denomFound) {
           if (pathname.find("Jet") != std::string::npos && !(pathname.find("DiJet") != std::string::npos) &&
               !(pathname.find("DiPFJet") != std::string::npos) && !(pathname.find("DoubleJet") != std::string::npos) &&
               !(pathname.find("BTag") != std::string::npos) && !(pathname.find("Mu") != std::string::npos) &&
               !(pathname.find("Fwd") != std::string::npos)) {
             triggerType = "SingleJet_Trigger";
             objectType = trigger::TriggerJet;
           }
           if (pathname.find("DiJet") != std::string::npos || pathname.find("DiPFJet") != std::string::npos ||
               pathname.find("DoubleJet") != std::string::npos) {
             triggerType = "DiJet_Trigger";
             objectType = trigger::TriggerJet;
           }
           if (pathname.find("MET") != std::string::npos) {
             triggerType = "MET_Trigger";
             objectType = trigger::TriggerMET;
           }
           if (pathname.find("TET") != std::string::npos) {
             triggerType = "TET_Trigger";
             objectType = trigger::TriggerTET;
           }
 
           l1pathname = getL1ConditionModuleName(pathname);  //ml added L1conditionmodulename
           std::vector<std::string> numpathmodules = hltConfig_.moduleLabels(pathname);
           for (auto& numpathmodule : numpathmodules) {
             edm::InputTag testTag(numpathmodule, "", processname_);
             if ((hltConfig_.moduleType(numpathmodule) == "HLT1CaloJet") ||
                 (hltConfig_.moduleType(numpathmodule) == "HLT1PFJet") ||
                 (hltConfig_.moduleType(numpathmodule) == "HLTDiJetAveFilter") ||
                 (hltConfig_.moduleType(numpathmodule) == "HLTDiPFJetAveFilter") ||
                 (hltConfig_.moduleType(numpathmodule) == "HLT1CaloMET") ||
                 (hltConfig_.moduleType(numpathmodule) == "HLTMhtFilter") ||
                 (hltConfig_.moduleType(numpathmodule) == "HLTPrescaler"))
               filtername = numpathmodule;
           }
 
           if (objectType != 0) {
             std::vector<std::string> numpathmodules = hltConfig_.moduleLabels(dpathname);
             for (auto& numpathmodule : numpathmodules) {
               edm::InputTag testTag(numpathmodule, "", processname_);
               if ((hltConfig_.moduleType(numpathmodule) == "HLT1CaloJet") ||
                   (hltConfig_.moduleType(numpathmodule) == "HLT1PFJet") ||
                   (hltConfig_.moduleType(numpathmodule) == "HLTDiJetAveFilter") ||
                   (hltConfig_.moduleType(numpathmodule) == "HLTDiPFJetAveFilter") ||
                   (hltConfig_.moduleType(numpathmodule) == "HLT1CaloMET") ||
                   (hltConfig_.moduleType(numpathmodule) == "HLTMhtFilter") ||
                   (hltConfig_.moduleType(numpathmodule) == "HLTPrescaler"))
                 Denomfiltername = numpathmodule;
             }
 
             if (verbose_)
               cout << "==pathname==" << pathname << "==denompath==" << dpathname << "==filtername==" << filtername
                    << "==denomfiltername==" << Denomfiltername << "==l1pathname==" << l1pathname
                    << "==objectType==" << objectType << endl;
             hltPathsEff_.push_back(PathInfo(usedPrescale,
                                             dpathname,
                                             pathname,
                                             l1pathname,
                                             filtername,
                                             Denomfiltername,
                                             processname_,
                                             objectType,
                                             triggerType));
           }
         }
       }
     }
 
     if (verbose_)
       cout << "== end hltPathsEff_.push_back ======" << endl;
   }
 }

◆ fillMEforEffAllTrigger()

void JetMETHLTOfflineSource::fillMEforEffAllTrigger	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Definition at line 816 of file JetMETHLTOfflineSource.cc.

References _fEMF, _feta, _fHPD, _max_CEMF, _max_CHEF, _min_CEMF, _min_CHEF, _min_NEMF, _min_NHEF, _n90Hits, _pfMHT, funct::abs(), edm::HLTGlobalStatus::accept(), calojet, calojetColl_, CaloJetEMF, CaloJetEta, CaloJetfHPD, CaloJetn90, CaloJetPhi, CaloJetPt, calometColl_, gather_cfg::cout, SiPixelRawToDigiRegional_cfi::deltaPhi, PbPb_ZMuSkimMuonDPG_cff::deltaR, HLT_2023v12_cff::dPhi, HGC3DClusterGenMatchSelector_cfi::dR, HLT_2023v12_cff::dRmin, PVValHelper::eta, trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::getObjects(), hltPathsEff_, EgHLTOfflineClient_cfi::hltTag, createfilelist::int, isBarrel(), isEndCap(), isForward(), edm::HandleBase::isValid(), BTaggingMonitor_cfi::met, pfjet, PFJetCEMF, PFJetCHEF, pfjetColl_, PFJetEta, PFJetNEMF, PFJetNHEF, PFJetPhi, PFJetPt, pfmetColl_, pfMHTx_All, pfMHTy_All, phi, processname_, edm::Handle< T >::product(), edm::TriggerNames::size(), trigger::TriggerEvent::sizeFilters(), mathSSE::sqrt(), ClusterTask_cfi::trig, edm::TriggerNames::triggerIndex(), trigger::TriggerJet, trigger::TriggerMET, triggerNames_, triggerObj_, triggerResults_, trigger::TriggerTET, findQualityFiles::v, and verbose_.

Referenced by analyze().

                                                                                                     {
   if (!triggerResults_.isValid())
     return;
 
   bool denompassed = false;
   bool numpassed = false;
   const trigger::TriggerObjectCollection& toc(triggerObj_->getObjects());
 
   for (auto& v : hltPathsEff_) {
     denompassed = false;
     numpassed = false;
 
     unsigned indexNum = triggerNames_.triggerIndex(v.getPath());
     unsigned indexDenom = triggerNames_.triggerIndex(v.getDenomPath());
 
     if (indexNum < triggerNames_.size() && triggerResults_->accept(indexNum))
       numpassed = true;
     if (indexDenom < triggerNames_.size() && triggerResults_->accept(indexDenom))
       denompassed = true;
 
     if (denompassed == false)
       continue;
 
     //if(numpassed==true){
     edm::InputTag hltTag(v.getLabel(), "", processname_);
     const int hltIndex = triggerObj_->filterIndex(hltTag);
     edm::InputTag l1Tag(v.getl1Path(), "", processname_);
     const int l1Index = triggerObj_->filterIndex(l1Tag);
     //}
 
     //----------------------------------------------------------------------
     //double pTcut = 0;
     double trigLowpTcut = 0;
     double trigMedpTcut = 0;
     double trigHighpTcut = 0;
     double trigLowpTcutFwd = 0;
     double trigMedpTcutFwd = 0;
     double trigHighpTcutFwd = 0;
     //
     //double pTPFcut = 0 ;
     double trigLowpTPFcut = 0;
     double trigMedpTPFcut = 0;
     double trigHighpTPFcut = 0;
     double trigLowpTPFcutFwd = 0;
     double trigMedpTPFcutFwd = 0;
     double trigHighpTPFcutFwd = 0;
     //
     //cout<<"pre-path" << v->getPath()<<endl;
     size_t jetstrfound = v.getPath().find("Jet");
     //size_t censtrfound = v->getPath().find("Central"); //shoouldn't be needed?
     string tpath = v.getPath();
     string jetTrigVal;
     float jetVal = 0.;
     //
     if (jetstrfound != string::npos) {  // && ustrfound != string::npos ){
       //cout<<v->getPath()<<endl;
       for (int trig = int(jetstrfound) + 3; trig < int(jetstrfound) + 7; trig++) {  // int(ustrfound); trig++){
         if (!isdigit(tpath[trig]))
           break;
         jetTrigVal += tpath[trig];
       }
       jetVal = atof(jetTrigVal.c_str());
       //
       if (jetVal > 0.) {
         if (jetVal < 50.) {
           //pTcut = jetVal / 2.;
           trigMedpTcut = jetVal + 5.;
           trigHighpTcut = jetVal + 10.;
           //
           trigLowpTcutFwd = jetVal + 9.;
           trigMedpTcutFwd = jetVal + 15.;
           trigHighpTcutFwd = jetVal + 21.;
         } else {
           //pTcut = jetVal - 20. ;
           trigMedpTcut = jetVal + 2.;
           trigHighpTcut = jetVal + 60.;
           //
           trigLowpTcutFwd = jetVal + 22.;
           trigMedpTcutFwd = jetVal + 25.;
           trigHighpTcutFwd = jetVal + 110.;
         }
         trigLowpTcut = jetVal;
       }
       //
       if (jetVal > 0.) {
         if (jetVal < 50.) {
           //pTPFcut = jetVal ;
           trigMedpTPFcut = jetVal + 20.;
           trigHighpTPFcut = jetVal + 40.;
           //
           trigLowpTPFcutFwd = jetVal + 60.;
           trigMedpTPFcutFwd = jetVal + 80.;
           trigHighpTPFcutFwd = jetVal + 100.;
         } else {
           //pTPFcut = jetVal  ;
           trigMedpTPFcut = jetVal + 40.;
           trigHighpTPFcut = jetVal + 140.;
           //
           trigLowpTPFcutFwd = jetVal + 110.;
           trigMedpTPFcutFwd = jetVal + 130.;
           trigHighpTPFcutFwd = jetVal + 190.;
         }
         trigLowpTPFcut = jetVal;
       }
     }
     //----------------------------------------------------------------------
 
     //CaloJet paths
     if (verbose_)
       std::cout << "fillMEforEffAllTrigger: CaloJet -------------------" << std::endl;
     if (calojetColl_.isValid() && (v.getObjectType() == trigger::TriggerJet)) {
       //cout<<"   - CaloJet "<<endl;
       //&& (v->getPath().find("HLT_PFJet")==std::string::npos)
       //&& (v->getPath().find("HLT_DiPFJet")==std::string::npos)){
       bool jetIDbool = false;
       double leadjpt = CaloJetPt[0];
       double leadjeta = CaloJetEta[0];
       double leadjphi = CaloJetPhi[0];
       //double ljemf    = CaloJetEMF[0];
       double ljfhpd = CaloJetfHPD[0];
       double ljn90 = CaloJetn90[0];
       if ((v.getTriggerType() == "SingleJet_Trigger") && !calojet.empty()) {  //this line stops the central jets
         if ((ljfhpd < _fHPD) && (ljn90 > _n90Hits)) {
           if (verbose_)
             cout << "Passed CaloJet ID -------------------" << endl;
           jetIDbool = true;
           //Denominator fill
           v.getMEhisto_DenominatorPt()->Fill(leadjpt);
           if (isBarrel(leadjeta))
             v.getMEhisto_DenominatorPtBarrel()->Fill(leadjpt);
           if (isEndCap(leadjeta))
             v.getMEhisto_DenominatorPtEndcap()->Fill(leadjpt);
           if (isForward(leadjeta))
             v.getMEhisto_DenominatorPtForward()->Fill(leadjpt);
           v.getMEhisto_DenominatorEta()->Fill(leadjeta);
           v.getMEhisto_DenominatorPhi()->Fill(leadjphi);
           v.getMEhisto_DenominatorEtaPhi()->Fill(leadjeta, leadjphi);
           if (isBarrel(leadjeta)) {
             v.getMEhisto_DenominatorEtaBarrel()->Fill(leadjeta);
             v.getMEhisto_DenominatorPhiBarrel()->Fill(leadjphi);
           }
           if (isEndCap(leadjeta)) {
             v.getMEhisto_DenominatorEtaEndcap()->Fill(leadjeta);
             v.getMEhisto_DenominatorPhiEndcap()->Fill(leadjphi);
           }
           if (isForward(leadjeta)) {
             v.getMEhisto_DenominatorEtaForward()->Fill(leadjeta);
             v.getMEhisto_DenominatorPhiForward()->Fill(leadjphi);
           }
           if ((leadjpt > trigLowpTcut && !isForward(leadjeta)) || (leadjpt > trigLowpTcutFwd && isForward(leadjeta))) {
             v.getMEhisto_DenominatorEta_LowpTcut()->Fill(leadjeta);
             v.getMEhisto_DenominatorPhi_LowpTcut()->Fill(leadjphi);
             v.getMEhisto_DenominatorEtaPhi_LowpTcut()->Fill(leadjeta, leadjphi);
           }
           if ((leadjpt > trigMedpTcut && !isForward(leadjeta)) || (leadjpt > trigMedpTcutFwd && isForward(leadjeta))) {
             v.getMEhisto_DenominatorEta_MedpTcut()->Fill(leadjeta);
             v.getMEhisto_DenominatorPhi_MedpTcut()->Fill(leadjphi);
             v.getMEhisto_DenominatorEtaPhi_MedpTcut()->Fill(leadjeta, leadjphi);
           }
           if ((leadjpt > trigHighpTcut && !isForward(leadjeta)) ||
               (leadjpt > trigHighpTcutFwd && isForward(leadjeta))) {
             v.getMEhisto_DenominatorEta_HighpTcut()->Fill(leadjeta);
             v.getMEhisto_DenominatorPhi_HighpTcut()->Fill(leadjphi);
             v.getMEhisto_DenominatorEtaPhi_HighpTcut()->Fill(leadjeta, leadjphi);
           }
 
           //Numerator fill
           if (numpassed) {
             //
             double dRmin = 99999.;
             double dPhimin = 9999.;
             if (v.getPath().find("L1") != std::string::npos) {
               if (l1Index >= triggerObj_->sizeFilters()) {
                 edm::LogInfo("JetMETHLTOfflineSource") << "no index hlt" << hltIndex << " of that name ";
               } else {
                 const trigger::Keys& kl1 = triggerObj_->filterKeys(l1Index);
                 for (unsigned short ki : kl1) {
                   double dR = deltaR(toc[ki].eta(), toc[ki].phi(), leadjeta, leadjphi);
                   if (dR < dRmin) {
                     dRmin = dR;
                   }
                 }
               }
             } else {
               if (hltIndex >= triggerObj_->sizeFilters()) {
                 edm::LogInfo("JetMETHLTOfflineSource") << "no index hlt" << hltIndex << " of that name ";
               } else {
                 const trigger::Keys& khlt = triggerObj_->filterKeys(hltIndex);
                 auto kj = khlt.begin();
                 for (; kj != khlt.end(); ++kj) {
                   double dR = deltaR(toc[*kj].eta(), toc[*kj].phi(), leadjeta, leadjphi);
                   if (dR < dRmin) {
                     dRmin = dR;
                   }
                   double dPhi = deltaPhi(toc[*kj].phi(), leadjphi);
                   if (dPhi < dPhimin) {
                     dPhimin = dPhi;
                   }
                 }
                 //v->getMEhisto_DeltaPhi()->Fill(dPhimin);
                 v.getMEhisto_DeltaPhi()->Fill(dPhimin);
                 v.getMEhisto_DeltaR()->Fill(dRmin);
               }
             }
             if (dRmin < 0.1 || (v.getPath().find("L1") != std::string::npos && dRmin < 0.4)) {
               v.getMEhisto_NumeratorPt()->Fill(leadjpt);
               if (isBarrel(leadjeta))
                 v.getMEhisto_NumeratorPtBarrel()->Fill(leadjpt);
               if (isEndCap(leadjeta))
                 v.getMEhisto_NumeratorPtEndcap()->Fill(leadjpt);
               if (isForward(leadjeta))
                 v.getMEhisto_NumeratorPtForward()->Fill(leadjpt);
               v.getMEhisto_NumeratorEta()->Fill(leadjeta);
               v.getMEhisto_NumeratorPhi()->Fill(leadjphi);
               v.getMEhisto_NumeratorEtaPhi()->Fill(leadjeta, leadjphi);
               if (isBarrel(leadjeta)) {
                 v.getMEhisto_NumeratorEtaBarrel()->Fill(leadjeta);
                 v.getMEhisto_NumeratorPhiBarrel()->Fill(leadjphi);
               }
               if (isEndCap(leadjeta)) {
                 v.getMEhisto_NumeratorEtaEndcap()->Fill(leadjeta);
                 v.getMEhisto_NumeratorPhiEndcap()->Fill(leadjphi);
               }
               if (isForward(leadjeta)) {
                 v.getMEhisto_NumeratorEtaForward()->Fill(leadjeta);
                 v.getMEhisto_NumeratorPhiForward()->Fill(leadjphi);
               }
               if ((leadjpt > trigLowpTcut && !isForward(leadjeta)) ||
                   (leadjpt > trigLowpTcutFwd && isForward(leadjeta))) {
                 v.getMEhisto_NumeratorEta_LowpTcut()->Fill(leadjeta);
                 v.getMEhisto_NumeratorPhi_LowpTcut()->Fill(leadjphi);
                 v.getMEhisto_NumeratorEtaPhi_LowpTcut()->Fill(leadjeta, leadjphi);
               }
               if ((leadjpt > trigMedpTcut && !isForward(leadjeta)) ||
                   (leadjpt > trigMedpTcutFwd && isForward(leadjeta))) {
                 v.getMEhisto_NumeratorEta_MedpTcut()->Fill(leadjeta);
                 v.getMEhisto_NumeratorPhi_MedpTcut()->Fill(leadjphi);
                 v.getMEhisto_NumeratorEtaPhi_MedpTcut()->Fill(leadjeta, leadjphi);
               }
               if ((leadjpt > trigHighpTcut && !isForward(leadjeta)) ||
                   (leadjpt > trigHighpTcutFwd && isForward(leadjeta))) {
                 v.getMEhisto_NumeratorEta_HighpTcut()->Fill(leadjeta);
                 v.getMEhisto_NumeratorPhi_HighpTcut()->Fill(leadjphi);
                 v.getMEhisto_NumeratorEtaPhi_HighpTcut()->Fill(leadjeta, leadjphi);
               }
             }
           }  //numpassed
         }    //CalojetID filter
       }
 
       if (jetIDbool == true && (v.getTriggerType() == "DiJet_Trigger") && calojet.size() > 1) {
         if (((CaloJetEMF[1] > _fEMF || std::abs(CaloJetEta[1]) > _feta) && CaloJetfHPD[0] < _fHPD &&
              CaloJetn90[0] > _n90Hits)) {
           v.getMEhisto_DenominatorPt()->Fill((CaloJetPt[0] + CaloJetPt[1]) / 2.);
           v.getMEhisto_DenominatorEta()->Fill((CaloJetEta[0] + CaloJetEta[1]) / 2.);
           if (numpassed == true) {
             v.getMEhisto_NumeratorPt()->Fill((CaloJetPt[0] + CaloJetPt[1]) / 2.);
             v.getMEhisto_NumeratorEta()->Fill((CaloJetEta[0] + CaloJetEta[1]) / 2.);
           }
         }
       }
     }  // Jet trigger and valid jet collection
 
     //PFJet paths
     if (verbose_)
       std::cout << "fillMEforEffAllTrigger: PFJet -------------------" << std::endl;
     if (pfjetColl_.isValid() && (v.getObjectType() == trigger::TriggerJet)) {
       //cout<<"   - PFJet "<<endl;
       //&& (v->getPath().find("HLT_PFJet")!=std::string::npos)
       //&& (v->getPath().find("HLT_DiPFJet")!=std::string::npos)){
       bool jetIDbool = false;
       double leadjpt = PFJetPt[0];
       double leadjeta = PFJetEta[0];
       double leadjphi = PFJetPhi[0];
       double ljNHEF = PFJetNHEF[0];
       double ljCHEF = PFJetCHEF[0];
       double ljNEMF = PFJetNEMF[0];
       double ljCEMF = PFJetCEMF[0];
       //double sleadjpt  = PFJetPt[1];
       //double sleadjeta = PFJetEta[1];
       //double sleadjphi = PFJetPhi[1];
       double sljNHEF = PFJetNHEF[1];
       double sljCHEF = PFJetCHEF[1];
       double sljNEMF = PFJetNEMF[1];
       double sljCEMF = PFJetCEMF[1];
       //
       double pfMHTx = pfMHTx_All;
       double pfMHTy = pfMHTy_All;
       //
       if ((v.getTriggerType() == "SingleJet_Trigger") && !pfjet.empty()) {  //this line stops the central jets
 
         //======get pfmht
         _pfMHT = sqrt(pfMHTx * pfMHTx + pfMHTy * pfMHTy);
         v.getMEhisto_DenominatorPFMHT()->Fill(_pfMHT);
 
         if (ljNHEF >= _min_NHEF && ljNHEF <= _max_NHEF && ljCHEF >= _min_CHEF && ljCHEF <= _max_CHEF &&
             ljNEMF >= _min_NEMF && ljNEMF <= _max_NEMF && ljCEMF >= _min_CEMF && ljCEMF <= _max_CEMF) {
           if (verbose_)
             cout << "Passed PFJet ID -------------------" << endl;
           jetIDbool = true;
           v.getMEhisto_DenominatorPFPt()->Fill(leadjpt);
           if (isBarrel(leadjeta))
             v.getMEhisto_DenominatorPFPtBarrel()->Fill(leadjpt);
           if (isEndCap(leadjeta))
             v.getMEhisto_DenominatorPFPtEndcap()->Fill(leadjpt);
           if (isForward(leadjeta))
             v.getMEhisto_DenominatorPFPtForward()->Fill(leadjpt);
           v.getMEhisto_DenominatorPFEta()->Fill(leadjeta);
           v.getMEhisto_DenominatorPFPhi()->Fill(leadjphi);
           v.getMEhisto_DenominatorPFEtaPhi()->Fill(leadjeta, leadjphi);
           if (isBarrel(leadjeta)) {
             v.getMEhisto_DenominatorPFEtaBarrel()->Fill(leadjeta);
             v.getMEhisto_DenominatorPFPhiBarrel()->Fill(leadjphi);
           }
           if (isEndCap(leadjeta)) {
             v.getMEhisto_DenominatorPFEtaEndcap()->Fill(leadjeta);
             v.getMEhisto_DenominatorPFPhiEndcap()->Fill(leadjphi);
           }
           if (isForward(leadjeta)) {
             v.getMEhisto_DenominatorPFEtaForward()->Fill(leadjeta);
             v.getMEhisto_DenominatorPFPhiForward()->Fill(leadjphi);
           }
           if ((leadjpt > trigLowpTPFcut && !isForward(leadjeta)) ||
               (leadjpt > trigLowpTPFcutFwd && isForward(leadjeta))) {
             v.getMEhisto_DenominatorPFEta_LowpTcut()->Fill(leadjeta);
             v.getMEhisto_DenominatorPFPhi_LowpTcut()->Fill(leadjphi);
             v.getMEhisto_DenominatorPFEtaPhi_LowpTcut()->Fill(leadjeta, leadjphi);
           }
           if ((leadjpt > trigMedpTPFcut && !isForward(leadjeta)) ||
               (leadjpt > trigMedpTPFcutFwd && isForward(leadjeta))) {
             v.getMEhisto_DenominatorPFEta_MedpTcut()->Fill(leadjeta);
             v.getMEhisto_DenominatorPFPhi_MedpTcut()->Fill(leadjphi);
             v.getMEhisto_DenominatorPFEtaPhi_MedpTcut()->Fill(leadjeta, leadjphi);
           }
           if ((leadjpt > trigHighpTPFcut && !isForward(leadjeta)) ||
               (leadjpt > trigHighpTPFcutFwd && isForward(leadjeta))) {
             v.getMEhisto_DenominatorPFEta_HighpTcut()->Fill(leadjeta);
             v.getMEhisto_DenominatorPFPhi_HighpTcut()->Fill(leadjphi);
             v.getMEhisto_DenominatorPFEtaPhi_HighpTcut()->Fill(leadjeta, leadjphi);
           }
 
           //Numerator fill
           if (numpassed) {
             double dRmin = 99999.;
             double dPhimin = 9999.;
             if (v.getPath().find("L1") != std::string::npos) {
               if (l1Index >= triggerObj_->sizeFilters()) {
                 edm::LogInfo("JetMETHLTOfflineSource") << "no index hlt" << hltIndex << " of that name ";
               } else {
                 const trigger::Keys& kl1 = triggerObj_->filterKeys(l1Index);
                 for (unsigned short ki : kl1) {
                   double dR = deltaR(toc[ki].eta(), toc[ki].phi(), leadjeta, leadjphi);
                   if (dR < dRmin) {
                     dRmin = dR;
                   }
                 }
               }
             } else {
               if (hltIndex >= triggerObj_->sizeFilters()) {
                 edm::LogInfo("JetMETHLTOfflineSource") << "no index hlt" << hltIndex << " of that name ";
               } else {
                 const trigger::Keys& khlt = triggerObj_->filterKeys(hltIndex);
                 for (unsigned short kj : khlt) {
                   double dR = deltaR(toc[kj].eta(), toc[kj].phi(), leadjeta, leadjphi);
                   if (dR < dRmin) {
                     dRmin = dR;
                   }
                   double dPhi = deltaPhi(toc[kj].phi(), leadjphi);
                   if (dPhi < dPhimin) {
                     dPhimin = dPhi;
                   }
                 }
                 v.getMEhisto_PFDeltaPhi()->Fill(dPhimin);
                 v.getMEhisto_PFDeltaR()->Fill(dRmin);
               }
             }
             if (dRmin < 0.1 || (v.getPath().find("L1") != std::string::npos && dRmin < 0.4)) {
               v.getMEhisto_NumeratorPFPt()->Fill(leadjpt);
               if (isBarrel(leadjeta))
                 v.getMEhisto_NumeratorPFPtBarrel()->Fill(leadjpt);
               if (isEndCap(leadjeta))
                 v.getMEhisto_NumeratorPFPtEndcap()->Fill(leadjpt);
               if (isForward(leadjeta))
                 v.getMEhisto_NumeratorPFPtForward()->Fill(leadjpt);
               v.getMEhisto_NumeratorPFEta()->Fill(leadjeta);
               v.getMEhisto_NumeratorPFPhi()->Fill(leadjphi);
               v.getMEhisto_NumeratorPFEtaPhi()->Fill(leadjeta, leadjphi);
               if (isBarrel(leadjeta)) {
                 v.getMEhisto_NumeratorPFEtaBarrel()->Fill(leadjeta);
                 v.getMEhisto_NumeratorPFPhiBarrel()->Fill(leadjphi);
               }
               if (isEndCap(leadjeta)) {
                 v.getMEhisto_NumeratorPFEtaEndcap()->Fill(leadjeta);
                 v.getMEhisto_NumeratorPFPhiEndcap()->Fill(leadjphi);
               }
               if (isForward(leadjeta)) {
                 v.getMEhisto_NumeratorPFEtaForward()->Fill(leadjeta);
                 v.getMEhisto_NumeratorPFPhiForward()->Fill(leadjphi);
               }
               if ((leadjpt > trigLowpTPFcut && !isForward(leadjeta)) ||
                   (leadjpt > trigLowpTPFcutFwd && isForward(leadjeta))) {
                 v.getMEhisto_NumeratorPFEta_LowpTcut()->Fill(leadjeta);
                 v.getMEhisto_NumeratorPFPhi_LowpTcut()->Fill(leadjphi);
                 v.getMEhisto_NumeratorPFEtaPhi_LowpTcut()->Fill(leadjeta, leadjphi);
               }
               if ((leadjpt > trigMedpTPFcut && !isForward(leadjeta)) ||
                   (leadjpt > trigMedpTPFcutFwd && isForward(leadjeta))) {
                 v.getMEhisto_NumeratorPFEta_MedpTcut()->Fill(leadjeta);
                 v.getMEhisto_NumeratorPFPhi_MedpTcut()->Fill(leadjphi);
                 v.getMEhisto_NumeratorPFEtaPhi_MedpTcut()->Fill(leadjeta, leadjphi);
               }
               if ((leadjpt > trigHighpTPFcut && !isForward(leadjeta)) ||
                   (leadjpt > trigHighpTPFcutFwd && isForward(leadjeta))) {
                 v.getMEhisto_NumeratorPFEta_HighpTcut()->Fill(leadjeta);
                 v.getMEhisto_NumeratorPFPhi_HighpTcut()->Fill(leadjphi);
                 v.getMEhisto_NumeratorPFEtaPhi_HighpTcut()->Fill(leadjeta, leadjphi);
               }
             }
           }
         }
       }
       if (jetIDbool == true && (v.getTriggerType() == "DiJet_Trigger") && pfjet.size() > 1) {
         if (ljNHEF >= _min_NHEF && ljNHEF <= _max_NHEF && ljCHEF >= _min_CHEF && ljCHEF <= _max_CHEF &&
             ljNEMF >= _min_NEMF && ljNEMF <= _max_NEMF && ljCEMF >= _min_CEMF && ljCEMF <= _max_CEMF &&
             sljNHEF >= _min_NHEF && sljNHEF <= _max_NHEF && sljCHEF >= _min_CHEF && sljCHEF <= _max_CHEF &&
             sljNEMF >= _min_NEMF && sljNEMF <= _max_NEMF && sljCEMF >= _min_CEMF && sljCEMF <= _max_CEMF) {
           v.getMEhisto_DenominatorPFPt()->Fill((PFJetPt[0] + PFJetPt[1]) / 2.);
           v.getMEhisto_DenominatorPFEta()->Fill((PFJetEta[0] + PFJetEta[1]) / 2.);
           if (numpassed) {
             v.getMEhisto_NumeratorPFPt()->Fill((PFJetPt[0] + PFJetPt[1]) / 2.);
             v.getMEhisto_NumeratorPFEta()->Fill((PFJetEta[0] + PFJetEta[1]) / 2.);
           }
         }
       }
     }  // PF Jet trigger and valid jet collection
 
     //CaloMET path
     if (verbose_)
       std::cout << "fillMEforEffAllTrigger: CaloMET -------------------" << std::endl;
     if (calometColl_.isValid() &&
         ((v.getObjectType() == trigger::TriggerMET) || (v.getObjectType() == trigger::TriggerTET)) &&
         (v.getPath().find("HLT_PFMET") == std::string::npos)) {
       const CaloMETCollection* calometcol = calometColl_.product();
       const CaloMET met = calometcol->front();
       v.getMEhisto_DenominatorPt()->Fill(met.et());
       v.getMEhisto_DenominatorPhi()->Fill(met.phi());
       if (numpassed) {
         v.getMEhisto_NumeratorPt()->Fill(met.et());
         v.getMEhisto_NumeratorPhi()->Fill(met.phi());
         if (hltIndex >= triggerObj_->sizeFilters()) {
           edm::LogInfo("JetMETHLTOfflineSource") << "no index hlt" << hltIndex << " of that name ";
         } else {
           double dPhimin = 9999.;  //
           const trigger::Keys& khlt = triggerObj_->filterKeys(hltIndex);
           for (unsigned short kj : khlt) {
             double dPhi = deltaPhi(toc[kj].phi(), met.phi());
             if (dPhi < dPhimin) {
               dPhimin = dPhi;
             }
           }
           v.getMEhisto_DeltaPhi()->Fill(dPhimin);
         }
       }
     }
 
     //PFMET
     if (verbose_)
       std::cout << "fillMEforEffAllTrigger: PFMET -------------------" << std::endl;
     if (pfmetColl_.isValid() &&
         ((v.getObjectType() == trigger::TriggerMET) || (v.getObjectType() == trigger::TriggerTET)) &&
         (v.getPath().find("HLT_PFMET") != std::string::npos)) {
       const PFMETCollection* pfmetcol = pfmetColl_.product();
       const PFMET met = pfmetcol->front();
       v.getMEhisto_DenominatorPt()->Fill(met.et());
       v.getMEhisto_DenominatorPhi()->Fill(met.phi());
       if (numpassed) {
         v.getMEhisto_NumeratorPt()->Fill(met.et());
         v.getMEhisto_NumeratorPhi()->Fill(met.phi());
         if (hltIndex >= triggerObj_->sizeFilters()) {
           edm::LogInfo("JetMETHLTOfflineSource") << "no index hlt" << hltIndex << " of that name ";
         } else {
           double dPhimin = 9999.;  //
           const trigger::Keys& khlt = triggerObj_->filterKeys(hltIndex);
           for (unsigned short kj : khlt) {
             double dPhi = deltaPhi(toc[kj].phi(), met.phi());
             if (dPhi < dPhimin) {
               dPhimin = dPhi;
             }
           }
           v.getMEhisto_DeltaPhi()->Fill(dPhimin);
         }
       }
     }
 
     /*
       if(pfmhtColl_.isValid() && ((v->getObjectType() == trigger::TriggerMET)|| (v->getObjectType() == trigger::TriggerTET))){
       const PFMHTCollection *pfmhtcol = pfmhtColl_.product();
       const PFMHT met = pfmhtcol->front();
       v->getMEhisto_DenominatorPFPt()->Fill(met.pt());
       v->getMEhisto_DenominatorPFPhi()->Fill(met.phi());
       }// PFMHT  trigger and valid MET collection 
     */
   }  // trigger under study
 }

◆ fillMEforMonAllTrigger()

void JetMETHLTOfflineSource::fillMEforMonAllTrigger	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Definition at line 486 of file JetMETHLTOfflineSource.cc.

References calojetColl_, CaloJetEta, CaloJetPhi, CaloJetPt, CaloJetPx, CaloJetPy, calometColl_, gather_cfg::cout, SiPixelRawToDigiRegional_cfi::deltaPhi, PbPb_ZMuSkimMuonDPG_cff::deltaR, EgHLTOffHistBins_cfi::et, reco::LeafCandidate::et(), PVValHelper::eta, trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::getObjects(), hltPathsAll_, EgHLTOfflineClient_cfi::hltTag, isBarrel(), isEndCap(), isForward(), isHLTPathAccepted(), edm::HandleBase::isValid(), BTaggingMonitor_cfi::met, pfjetColl_, PFJetEta, PFJetPhi, PFJetPt, PFJetPx, PFJetPy, pfmetColl_, phi, reco::LeafCandidate::phi(), processname_, edm::Handle< T >::product(), DiDispStaMuonMonitor_cfi::pt, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, trigger::TriggerEvent::sizeFilters(), trigger::TriggerJet, trigger::TriggerMET, triggerObj_, triggerResults_, trigger::TriggerTET, findQualityFiles::v, and verbose_.

Referenced by analyze().

                                                                                                     {
   if (verbose_)
     cout << ">> Inside fillMEforMonAllTrigger " << endl;
   if (!triggerResults_.isValid())
     return;
   if (verbose_)
     cout << "   ... and triggerResults is valid" << endl;
 
   const trigger::TriggerObjectCollection& toc(triggerObj_->getObjects());
   for (auto& v : hltPathsAll_) {
     if (verbose_)
       cout << "   + Checking path " << v.getPath();
     if (isHLTPathAccepted(v.getPath()) == false) {
       if (verbose_)
         cout << " - failed" << endl;
       continue;
     }
     if (verbose_)
       cout << " - PASSED! " << endl;
 
     //New jet collection (after apply JEC)
     std::vector<double> jetPtVec;
     std::vector<double> jetPhiVec;
     std::vector<double> jetEtaVec;
     std::vector<double> jetPxVec;
     std::vector<double> jetPyVec;
     std::vector<double> hltPtVec;
     std::vector<double> hltPhiVec;
     std::vector<double> hltEtaVec;
     std::vector<double> hltPxVec;
     std::vector<double> hltPyVec;
 
     //This will be used to find out punch through trigger
     //bool fillL1HLT = false;
 
     //L1 and HLT indices
     if (verbose_) {
       cout << "     - L1Path = " << v.getl1Path() << endl;
       cout << "     - Label  = " << v.getLabel() << endl;
     }
 
     //edm::InputTag l1Tag(v->getl1Path(),"",processname_);
     edm::InputTag l1Tag(v.getLabel(), "", processname_);
     const int l1Index = triggerObj_->filterIndex(l1Tag);
     if (verbose_)
       cout << "     - l1Index = " << l1Index << " - l1Tag = [" << l1Tag << "]" << endl;
 
     edm::InputTag hltTag(v.getLabel(), "", processname_);
     const int hltIndex = triggerObj_->filterIndex(hltTag);
     if (verbose_)
       cout << "     - hltIndex = " << hltIndex << " - hltTag = [" << hltTag << "]" << endl;
 
     //bool l1TrigBool = false;
     bool hltTrigBool = false;
     bool diJetFire = false;
     int jetsize = 0;
 
     if (l1Index >= triggerObj_->sizeFilters()) {
       edm::LogInfo("JetMETHLTOfflineSource") << "no index " << l1Index << " of that name " << l1Tag;
       if (verbose_)
         cout << "[JetMETHLTOfflineSource::fillMEforMonAllTrigger] - No index l1Index=" << l1Index << " of that name \""
              << l1Tag << "\"" << endl;
     } else {
       //l1TrigBool = true;
       const trigger::Keys& kl1 = triggerObj_->filterKeys(l1Index);
       //
       if (v.getObjectType() == trigger::TriggerJet && v.getTriggerType() == "SingleJet_Trigger")
         v.getMEhisto_N_L1()->Fill(kl1.size());
       //
       auto ki = kl1.begin();
       for (; ki != kl1.end(); ++ki) {
         double l1TrigEta = -100;
         double l1TrigPhi = -100;
         //
         if (v.getObjectType() == trigger::TriggerJet) {
           l1TrigEta = toc[*ki].eta();
           l1TrigPhi = toc[*ki].phi();
           if (v.getTriggerType() == "SingleJet_Trigger") {
             v.getMEhisto_Pt_L1()->Fill(toc[*ki].pt());
             if (isBarrel(toc[*ki].eta()))
               v.getMEhisto_PtBarrel_L1()->Fill(toc[*ki].pt());
             if (isEndCap(toc[*ki].eta()))
               v.getMEhisto_PtEndcap_L1()->Fill(toc[*ki].pt());
             if (isForward(toc[*ki].eta()))
               v.getMEhisto_PtForward_L1()->Fill(toc[*ki].pt());
             v.getMEhisto_Eta_L1()->Fill(toc[*ki].eta());
             v.getMEhisto_Phi_L1()->Fill(toc[*ki].phi());
             v.getMEhisto_EtaPhi_L1()->Fill(toc[*ki].eta(), toc[*ki].phi());
           }
         }
         if (v.getObjectType() == trigger::TriggerMET || v.getObjectType() == trigger::TriggerTET) {
           v.getMEhisto_Pt_L1()->Fill(toc[*ki].pt());
           v.getMEhisto_Phi_L1()->Fill(toc[*ki].phi());
         }
 
         //-----------------------------------------------
         if (hltIndex >= triggerObj_->sizeFilters()) {
           edm::LogInfo("JetMETHLTOfflineSource") << "no index hlt" << hltIndex << " of that name ";
           if (verbose_)
             cout << "[JetMETHLTOfflineSource::fillMEforMonAllTrigger] - No index hltIndex=" << hltIndex
                  << " of that name " << endl;
         } else {
           const trigger::Keys& khlt = triggerObj_->filterKeys(hltIndex);
           if (v.getObjectType() == trigger::TriggerJet && ki == kl1.begin() &&
               v.getTriggerType() == "SingleJet_Trigger")
             v.getMEhisto_N_HLT()->Fill(khlt.size());
           //
           auto kj = khlt.begin();
           //Define hltTrigBool
           for (; kj != khlt.end(); ++kj) {
             if (v.getObjectType() == trigger::TriggerJet) {
               double hltTrigEta = -100;
               double hltTrigPhi = -100;
               hltTrigEta = toc[*kj].eta();
               hltTrigPhi = toc[*kj].phi();
               if ((deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi)) < 0.4 &&
                   (v.getTriggerType() == "DiJet_Trigger"))
                 hltTrigBool = true;
             }
           }
           //
           kj = khlt.begin();
           for (; kj != khlt.end(); ++kj) {
             double hltTrigEta = -100.;
             double hltTrigPhi = -100.;
             //fillL1HLT = true;
             //MET Triggers
             if (verbose_)
               cout << "+ MET Triggers plots" << endl;
             if (v.getObjectType() == trigger::TriggerMET || (v.getObjectType() == trigger::TriggerTET)) {
               v.getMEhisto_Pt_HLT()->Fill(toc[*kj].pt());
               v.getMEhisto_Phi_HLT()->Fill(toc[*kj].phi());
               v.getMEhisto_PtCorrelation_L1HLT()->Fill(toc[*ki].pt(), toc[*kj].pt());
               v.getMEhisto_PhiCorrelation_L1HLT()->Fill(toc[*ki].phi(), toc[*kj].phi());
               v.getMEhisto_PtResolution_L1HLT()->Fill((toc[*ki].pt() - toc[*kj].pt()) / (toc[*ki].pt()));
               v.getMEhisto_PhiResolution_L1HLT()->Fill(toc[*ki].phi() - toc[*kj].phi());
             }
             //Jet Triggers
             if (verbose_)
               cout << "+ Jet Trigger plots" << endl;
             if (v.getObjectType() == trigger::TriggerJet) {
               if (verbose_)
                 cout << "  - Going for those..." << endl;
               hltTrigEta = toc[*kj].eta();
               hltTrigPhi = toc[*kj].phi();
               if ((deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi)) < 0.4) {
                 if (v.getTriggerType() == "SingleJet_Trigger") {
                   v.getMEhisto_PtCorrelation_L1HLT()->Fill(toc[*ki].pt(), toc[*kj].pt());
                   v.getMEhisto_EtaCorrelation_L1HLT()->Fill(toc[*ki].eta(), toc[*kj].eta());
                   v.getMEhisto_PhiCorrelation_L1HLT()->Fill(toc[*ki].phi(), toc[*kj].phi());
                   v.getMEhisto_PtResolution_L1HLT()->Fill((toc[*ki].pt() - toc[*kj].pt()) / (toc[*ki].pt()));
                   v.getMEhisto_EtaResolution_L1HLT()->Fill(toc[*ki].eta() - toc[*kj].eta());
                   v.getMEhisto_PhiResolution_L1HLT()->Fill(toc[*ki].phi() - toc[*kj].phi());
                 }
               }
               if (((deltaR(hltTrigEta, hltTrigPhi, l1TrigEta, l1TrigPhi) < 0.4) ||
                    ((v.getTriggerType() == "DiJet_Trigger") && hltTrigBool)) &&
                   !diJetFire) {
                 if (v.getTriggerType() == "SingleJet_Trigger") {
                   v.getMEhisto_Pt_HLT()->Fill(toc[*kj].pt());
                   if (isBarrel(toc[*kj].eta()))
                     v.getMEhisto_PtBarrel_HLT()->Fill(toc[*kj].pt());
                   if (isEndCap(toc[*kj].eta()))
                     v.getMEhisto_PtEndcap_HLT()->Fill(toc[*kj].pt());
                   if (isForward(toc[*kj].eta()))
                     v.getMEhisto_PtForward_HLT()->Fill(toc[*kj].pt());
                   v.getMEhisto_Eta_HLT()->Fill(toc[*kj].eta());
                   v.getMEhisto_Phi_HLT()->Fill(toc[*kj].phi());
                   v.getMEhisto_EtaPhi_HLT()->Fill(toc[*kj].eta(), toc[*kj].phi());
                 }
 
                 //Calojet
                 if (calojetColl_.isValid() && (v.getObjectType() == trigger::TriggerJet) && (v.getPath() == "PFJet")) {
                   //CaloJetCollection::const_iterator jet = calojet.begin();
                   //for(; jet != calojet.end(); ++jet) {
                   for (int iCalo = 0; iCalo < 2; iCalo++) {
                     if (deltaR(hltTrigEta, hltTrigPhi, CaloJetEta[iCalo], CaloJetPhi[iCalo]) < 0.4) {
                       jetsize++;
                       if (v.getTriggerType() == "SingleJet_Trigger") {
                         v.getMEhisto_Pt()->Fill(CaloJetPt[iCalo]);
                         if (isBarrel(CaloJetEta[iCalo]))
                           v.getMEhisto_PtBarrel()->Fill(CaloJetPt[iCalo]);
                         if (isEndCap(CaloJetEta[iCalo]))
                           v.getMEhisto_PtEndcap()->Fill(CaloJetPt[iCalo]);
                         if (isForward(CaloJetEta[iCalo]))
                           v.getMEhisto_PtForward()->Fill(CaloJetPt[iCalo]);
                         //
                         v.getMEhisto_Eta()->Fill(CaloJetEta[iCalo]);
                         v.getMEhisto_Phi()->Fill(CaloJetPhi[iCalo]);
                         v.getMEhisto_EtaPhi()->Fill(CaloJetEta[iCalo], CaloJetPhi[iCalo]);
                         //
                         v.getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].pt(), CaloJetPt[iCalo]);
                         v.getMEhisto_EtaCorrelation_HLTRecObj()->Fill(toc[*kj].eta(), CaloJetEta[iCalo]);
                         v.getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].phi(), CaloJetPhi[iCalo]);
                         //
                         v.getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].pt() - CaloJetPt[iCalo]) /
                                                                     (toc[*kj].pt()));
                         v.getMEhisto_EtaResolution_HLTRecObj()->Fill(toc[*kj].eta() - CaloJetEta[iCalo]);
                         v.getMEhisto_PhiResolution_HLTRecObj()->Fill(toc[*kj].phi() - CaloJetPhi[iCalo]);
                       }
 
                       //-------------------------------------------------------
                       if ((v.getTriggerType() == "DiJet_Trigger")) {
                         jetPhiVec.push_back(CaloJetPhi[iCalo]);
                         jetPtVec.push_back(CaloJetPt[iCalo]);
                         jetEtaVec.push_back(CaloJetEta[iCalo]);
                         jetPxVec.push_back(CaloJetPx[iCalo]);
                         jetPyVec.push_back(CaloJetPy[iCalo]);
                         //
                         hltPhiVec.push_back(toc[*kj].phi());
                         hltPtVec.push_back(toc[*kj].pt());
                         hltEtaVec.push_back(toc[*kj].eta());
                         hltPxVec.push_back(toc[*kj].px());
                         hltPyVec.push_back(toc[*kj].py());
                       }
                     }  // matching jet
                   }    // Jet Loop
                 }      // valid calojet collection, with calojet trigger
 
                 //PFJet trigger
                 if (pfjetColl_.isValid() && (v.getObjectType() == trigger::TriggerJet) && (v.getPath() != "PFJet")) {
                   //PFJetCollection::const_iterator jet = pfjet.begin();
                   //for(; jet != pfjet.end(); ++jet){
                   for (int iPF = 0; iPF < 2; iPF++) {
                     if (deltaR(hltTrigEta, hltTrigPhi, PFJetEta[iPF], PFJetPhi[iPF]) < 0.4) {
                       jetsize++;
                       if (v.getTriggerType() == "SingleJet_Trigger") {
                         v.getMEhisto_Pt()->Fill(PFJetPt[iPF]);
                         if (isBarrel(PFJetEta[iPF]))
                           v.getMEhisto_PtBarrel()->Fill(PFJetPt[iPF]);
                         if (isEndCap(PFJetEta[iPF]))
                           v.getMEhisto_PtEndcap()->Fill(PFJetPt[iPF]);
                         if (isForward(PFJetEta[iPF]))
                           v.getMEhisto_PtForward()->Fill(PFJetPt[iPF]);
                         //
                         v.getMEhisto_Eta()->Fill(PFJetEta[iPF]);
                         v.getMEhisto_Phi()->Fill(PFJetPhi[iPF]);
                         v.getMEhisto_EtaPhi()->Fill(PFJetEta[iPF], PFJetPhi[iPF]);
                         //
                         v.getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].pt(), PFJetPt[iPF]);
                         v.getMEhisto_EtaCorrelation_HLTRecObj()->Fill(toc[*kj].eta(), PFJetEta[iPF]);
                         v.getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].phi(), PFJetPhi[iPF]);
                         //
                         v.getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].pt() - PFJetPt[iPF]) / (toc[*kj].pt()));
                         v.getMEhisto_EtaResolution_HLTRecObj()->Fill(toc[*kj].eta() - PFJetEta[iPF]);
                         v.getMEhisto_PhiResolution_HLTRecObj()->Fill(toc[*kj].phi() - PFJetPhi[iPF]);
                       }
 
                       //-------------------------------------------------------
                       if ((v.getTriggerType() == "DiJet_Trigger")) {
                         jetPhiVec.push_back(PFJetPhi[iPF]);
                         jetPtVec.push_back(PFJetPt[iPF]);
                         jetEtaVec.push_back(PFJetEta[iPF]);
                         jetPxVec.push_back(PFJetPx[iPF]);
                         jetPyVec.push_back(PFJetPy[iPF]);
                         //
                         hltPhiVec.push_back(toc[*kj].phi());
                         hltPtVec.push_back(toc[*kj].pt());
                         hltEtaVec.push_back(toc[*kj].eta());
                         hltPxVec.push_back(toc[*kj].px());
                         hltPyVec.push_back(toc[*kj].py());
                       }
                     }  // matching jet
                   }    //PFJet loop
                 }      //valid pfjet collection, with pfjet trigger
                        //
               }        // hlt matching with l1
             }          // jet trigger
 
             //------------------------------------------------------
             if (calometColl_.isValid() &&
                 ((v.getObjectType() == trigger::TriggerMET) || (v.getObjectType() == trigger::TriggerTET)) &&
                 (v.getPath().find("HLT_PFMET") == std::string::npos)) {
               const CaloMETCollection* calometcol = calometColl_.product();
               const CaloMET met = calometcol->front();
               //
               v.getMEhisto_Pt()->Fill(met.et());
               v.getMEhisto_Phi()->Fill(met.phi());
               //
               v.getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].et(), met.et());
               v.getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].phi(), met.phi());
               v.getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].et() - met.et()) / (toc[*kj].et()));
               v.getMEhisto_PhiResolution_HLTRecObj()->Fill(toc[*kj].phi() - met.phi());
             }
 
             //--------------------------------------------------------
             if (pfmetColl_.isValid() &&
                 ((v.getObjectType() == trigger::TriggerMET) || (v.getObjectType() == trigger::TriggerTET)) &&
                 (v.getPath().find("HLT_PFMET") != std::string::npos)) {
               const PFMETCollection* pfmetcol = pfmetColl_.product();
               const PFMET pfmet = pfmetcol->front();
               //
               v.getMEhisto_Pt()->Fill(pfmet.et());
               v.getMEhisto_Phi()->Fill(pfmet.phi());
               //
               v.getMEhisto_PtCorrelation_HLTRecObj()->Fill(toc[*kj].et(), pfmet.et());
               v.getMEhisto_PhiCorrelation_HLTRecObj()->Fill(toc[*kj].phi(), pfmet.phi());
               v.getMEhisto_PtResolution_HLTRecObj()->Fill((toc[*kj].et() - pfmet.et()) / (toc[*kj].et()));
               v.getMEhisto_PhiResolution_HLTRecObj()->Fill(toc[*kj].phi() - pfmet.phi());
             }
           }  //Loop over HLT trigger candidates
           if ((v.getTriggerType() == "DiJet_Trigger"))
             diJetFire = true;
         }  // Valid hlt trigger object
       }    // Loop over L1 objects
     }      // Valid L1 trigger object
     v.getMEhisto_N()->Fill(jetsize);
 
     //--------------------------------------------------------
     if ((v.getTriggerType() == "DiJet_Trigger") && jetPtVec.size() > 1) {
       double AveJetPt = (jetPtVec[0] + jetPtVec[1]) / 2;
       double AveJetEta = (jetEtaVec[0] + jetEtaVec[1]) / 2;
       double JetDelPhi = deltaPhi(jetPhiVec[0], jetPhiVec[1]);
       double AveHLTPt = (hltPtVec[0] + hltPtVec[1]) / 2;
       double AveHLTEta = (hltEtaVec[0] + hltEtaVec[1]) / 2;
       double HLTDelPhi = deltaPhi(hltPhiVec[0], hltPhiVec[1]);
       v.getMEhisto_AveragePt_RecObj()->Fill(AveJetPt);
       v.getMEhisto_AverageEta_RecObj()->Fill(AveJetEta);
       v.getMEhisto_DeltaPhi_RecObj()->Fill(JetDelPhi);
       //
       v.getMEhisto_AveragePt_HLTObj()->Fill(AveHLTPt);
       v.getMEhisto_AverageEta_HLTObj()->Fill(AveHLTEta);
       v.getMEhisto_DeltaPhi_HLTObj()->Fill(HLTDelPhi);
     }
   }
   if (verbose_)
     cout << "<< Exiting fillMEforMonAllTrigger " << endl;
 }

◆ fillMEforMonTriggerSummary()

void JetMETHLTOfflineSource::fillMEforMonTriggerSummary	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Definition at line 314 of file JetMETHLTOfflineSource.cc.

References correlation_All, correlation_AllWrtMB, correlation_AllWrtMu, gather_cfg::cout, dqm::impl::MonitorElement::Fill(), hltConfig_, hltPathsAll_, iEvent, isHLTPathAccepted(), dqmiolumiharvest::j, MBTrigPaths_, MuonTrigPaths_, wplusjetsAnalysis_cfi::muTrig, NVertices, BTVHLTOfflineSource_cfi::pathname, PVZ, rate_All, rate_AllWrtMB, rate_AllWrtMu, runStandalone_, HLTConfigProvider::size(), AlCaHLTBitMon_QueryRunRegistry::string, HLTConfigProvider::triggerName(), TriggerPosition(), findQualityFiles::v, verbose_, vertexToken, and w().

Referenced by analyze().

                                                                                                         {
   if (verbose_)
     cout << ">> Inside fillMEforMonTriggerSummary " << endl;
   bool muTrig = false;
 
   for (auto const& MuonTrigPath : MuonTrigPaths_) {
     const unsigned int nPath(hltConfig_.size());
     for (unsigned int j = 0; j != nPath; ++j) {
       std::string pathname = hltConfig_.triggerName(j);
       if (pathname.find(MuonTrigPath) != std::string::npos) {
         if (isHLTPathAccepted(pathname)) {
           muTrig = true;
           if (verbose_)
             cout << "fillMEforMonTriggerSummary: Muon Match" << endl;
         }
       }
       if (muTrig)
         break;
     }
     if (muTrig)
       break;
   }
 
   bool mbTrig = false;
   for (auto const& MBTrigPath : MBTrigPaths_) {
     const unsigned int nPath(hltConfig_.size());
     for (unsigned int j = 0; j != nPath; ++j) {
       std::string pathname = hltConfig_.triggerName(j);
       if (pathname.find(MBTrigPath) != std::string::npos) {
         if (isHLTPathAccepted(pathname)) {
           mbTrig = true;
           if (verbose_)
             cout << "fillMEforMonTriggerSummary: MinBias Match" << endl;
         }
       }
       if (mbTrig)
         break;
     }
     if (mbTrig)
       break;
   }
 
   auto v = hltPathsAll_.begin();
   for (; v != hltPathsAll_.end(); ++v) {
     bool trigFirst = false;
     double binV = TriggerPosition(v->getPath());
     if (isHLTPathAccepted(v->getPath()))
       trigFirst = true;
     if (!trigFirst)
       continue;
     if (trigFirst) {
       rate_All->Fill(binV);
       correlation_All->Fill(binV, binV);
       if (muTrig && runStandalone_) {
         rate_AllWrtMu->Fill(binV);
         correlation_AllWrtMu->Fill(binV, binV);
       }
       if (mbTrig && runStandalone_) {
         rate_AllWrtMB->Fill(binV);
         correlation_AllWrtMB->Fill(binV, binV);
       }
     }
     for (auto w = v + 1; w != hltPathsAll_.end(); ++w) {
       bool trigSec = false;
       double binW = TriggerPosition(w->getPath());
       if (isHLTPathAccepted(w->getPath()))
         trigSec = true;
       if (trigSec && trigFirst) {
         correlation_All->Fill(binV, binW);
         if (muTrig && runStandalone_)
           correlation_AllWrtMu->Fill(binV, binW);
         if (mbTrig && runStandalone_)
           correlation_AllWrtMB->Fill(binV, binW);
       }
       if (!trigSec && trigFirst) {
         correlation_All->Fill(binW, binV);
         if (muTrig && runStandalone_)
           correlation_AllWrtMu->Fill(binW, binV);
         if (mbTrig && runStandalone_)
           correlation_AllWrtMB->Fill(binW, binV);
       }
     }
   }
 
   //Vertex
   edm::Handle<VertexCollection> Vtx;
   iEvent.getByToken(vertexToken, Vtx);
   int vtxcnt = 0;
   for (auto const& itv : *Vtx) {
     //if(vtxcnt>=20) break;
     PVZ->Fill(itv.z());
     //chi2vtx[vtxcnt] = itv->chi2();
     //ndofvtx[vtxcnt] = itv->ndof();
     //ntrkvtx[vtxcnt] = itv->tracksSize();
     vtxcnt++;
   }
   NVertices->Fill(vtxcnt);
 }

◆ fillMEforTriggerNTfired()

void JetMETHLTOfflineSource::fillMEforTriggerNTfired ( )

privatevirtual

Definition at line 414 of file JetMETHLTOfflineSource.cc.

References edm::HLTGlobalStatus::accept(), calojet, calojetColl_, calometColl_, gather_cfg::cout, SiPixelRawToDigiRegional_cfi::deltaPhi, trigger::TriggerEvent::filterIndex(), hltPathsAll_, edm::HandleBase::isValid(), metsig::jet, BTaggingMonitor_cfi::met, processname_, edm::Handle< T >::product(), edm::TriggerNames::size(), edm::TriggerNames::triggerIndex(), triggerNames_, triggerObj_, triggerResults_, findQualityFiles::v, and verbose_.

Referenced by analyze().

                                                      {
   if (verbose_)
     cout << ">> Inside fillMEforTriggerNTfired" << endl;
   if (!triggerResults_.isValid())
     return;
   if (verbose_)
     cout << "   ... and triggerResults is valid" << endl;
 
   //
   for (auto& v : hltPathsAll_) {
     unsigned index = triggerNames_.triggerIndex(v.getPath());
     if (index < triggerNames_.size()) {
       v.getMEhisto_TriggerSummary()->Fill(0.);
       edm::InputTag l1Tag(v.getl1Path(), "", processname_);
       const int l1Index = triggerObj_->filterIndex(l1Tag);
       bool l1found = false;
       if (l1Index < triggerObj_->sizeFilters())
         l1found = true;
       if (!l1found)
         v.getMEhisto_TriggerSummary()->Fill(1.);
       if (!l1found && !(triggerResults_->accept(index)))
         v.getMEhisto_TriggerSummary()->Fill(2.);
       if (!l1found && (triggerResults_->accept(index)))
         v.getMEhisto_TriggerSummary()->Fill(3.);
       if (l1found)
         v.getMEhisto_TriggerSummary()->Fill(4.);
       if (l1found && (triggerResults_->accept(index)))
         v.getMEhisto_TriggerSummary()->Fill(5.);
       if (l1found && !(triggerResults_->accept(index)))
         v.getMEhisto_TriggerSummary()->Fill(6.);
       if (!(triggerResults_->accept(index)) && l1found) {
         //cout<<v->getTriggerType()<<endl;
         if ((v.getTriggerType() == "SingleJet_Trigger") && (calojetColl_.isValid()) && !calojet.empty()) {
           auto jet = calojet.begin();
           v.getMEhisto_JetPt()->Fill(jet->pt());
           v.getMEhisto_EtavsPt()->Fill(jet->eta(), jet->pt());
           v.getMEhisto_PhivsPt()->Fill(jet->phi(), jet->pt());
         }
         // single jet trigger is not fired
 
         if ((v.getTriggerType() == "DiJet_Trigger") && calojetColl_.isValid() && !calojet.empty()) {
           v.getMEhisto_JetSize()->Fill(calojet.size());
           if (calojet.size() >= 2) {
             auto jet = calojet.begin();
             auto jet2 = calojet.begin();
             jet2++;
             double jet3pt = 0.;
             if (calojet.size() > 2) {
               auto jet3 = jet2++;
               jet3pt = jet3->pt();
             }
             v.getMEhisto_Pt12()->Fill((jet->pt() + jet2->pt()) / 2.);
             v.getMEhisto_Eta12()->Fill((jet->eta() + jet2->eta()) / 2.);
             v.getMEhisto_Phi12()->Fill(deltaPhi(jet->phi(), jet2->phi()));
             v.getMEhisto_Pt3()->Fill(jet3pt);
             v.getMEhisto_Pt12Pt3()->Fill((jet->pt() + jet2->pt()) / 2., jet3pt);
             v.getMEhisto_Pt12Phi12()->Fill((jet->pt() + jet2->pt()) / 2., deltaPhi(jet->phi(), jet2->phi()));
           }
         }  // di jet trigger is not fired
 
         if (((v.getTriggerType() == "MET_Trigger") || (v.getTriggerType() == "TET_Trigger")) &&
             calometColl_.isValid()) {
           const CaloMETCollection* calometcol = calometColl_.product();
           const CaloMET met = calometcol->front();
           v.getMEhisto_JetPt()->Fill(met.pt());
         }  //MET trigger is not fired
       }    // L1 is fired
     }      //
   }        // trigger not fired
 }

◆ getL1ConditionModuleName()

const std::string JetMETHLTOfflineSource::getL1ConditionModuleName ( const std::string & pathname )

private

Definition at line 3288 of file JetMETHLTOfflineSource.cc.

References hltConfig_, HLTConfigProvider::moduleLabels(), HLTConfigProvider::moduleType(), BTVHLTOfflineSource_cfi::pathname, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dqmBeginRun().

                                                                                           {
   // find L1 condition for numpath with numpath objecttype
   // find PSet for L1 global seed for numpath,
   // list module labels for numpath
   std::string l1pathname = "dummy";
 
   std::vector<std::string> numpathmodules = hltConfig_.moduleLabels(pathname);
 
   for (auto& numpathmodule : numpathmodules) {
     if (hltConfig_.moduleType(numpathmodule) == "HLTLevel1GTSeed") {
       l1pathname = numpathmodule;
       break;
     }
   }
   return l1pathname;
 }

◆ isBarrel()

bool JetMETHLTOfflineSource::isBarrel ( double eta )

privatevirtual

Definition at line 3306 of file JetMETHLTOfflineSource.cc.

References PVValHelper::eta, and convertSQLitetoXML_cfg::output.

Referenced by fillMEforEffAllTrigger(), and fillMEforMonAllTrigger().

                                                 {
   bool output = false;
   if (fabs(eta) <= 1.3)
     output = true;
   return output;
 }

◆ isEndCap()

bool JetMETHLTOfflineSource::isEndCap ( double eta )

privatevirtual

Definition at line 3314 of file JetMETHLTOfflineSource.cc.

References PVValHelper::eta, and convertSQLitetoXML_cfg::output.

Referenced by fillMEforEffAllTrigger(), and fillMEforMonAllTrigger().

                                                 {
   bool output = false;
   if (fabs(eta) <= 3.0 && fabs(eta) > 1.3)
     output = true;
   return output;
 }

◆ isForward()

bool JetMETHLTOfflineSource::isForward ( double eta )

privatevirtual

Definition at line 3322 of file JetMETHLTOfflineSource.cc.

References PVValHelper::eta, and convertSQLitetoXML_cfg::output.

Referenced by fillMEforEffAllTrigger(), and fillMEforMonAllTrigger().

                                                  {
   bool output = false;
   if (fabs(eta) > 3.0)
     output = true;
   return output;
 }

◆ isHLTPathAccepted()

bool JetMETHLTOfflineSource::isHLTPathAccepted ( std::string pathName )

privatevirtual

Definition at line 3341 of file JetMETHLTOfflineSource.cc.

References edm::HLTGlobalStatus::accept(), edm::HandleBase::isValid(), convertSQLitetoXML_cfg::output, hltMonBTagIPClient_cfi::pathName, edm::TriggerNames::size(), edm::TriggerNames::triggerIndex(), triggerNames_, and triggerResults_.

Referenced by fillMEforMonAllTrigger(), and fillMEforMonTriggerSummary().

                                                                  {
   // triggerResults_, triggerNames_ has to be defined first before calling this method
   bool output = false;
   if (triggerResults_.isValid()) {
     unsigned index = triggerNames_.triggerIndex(pathName);
     if (index < triggerNames_.size() && triggerResults_->accept(index))
       output = true;
   }
   return output;
 }

◆ isTriggerObjectFound()

bool JetMETHLTOfflineSource::isTriggerObjectFound ( std::string objectName )

privatevirtual

Definition at line 3375 of file JetMETHLTOfflineSource.cc.

References trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), dqmdumpme::k, summarizeEdmComparisonLogfiles::objectName, convertSQLitetoXML_cfg::output, processname_, trigger::TriggerEvent::sizeFilters(), and triggerObj_.

                                                                       {
   // processname_, triggerObj_ has to be defined before calling this method
   bool output = false;
   edm::InputTag testTag(objectName, "", processname_);
   const int index = triggerObj_->filterIndex(testTag);
   if (index >= triggerObj_->sizeFilters()) {
     edm::LogInfo("JetMETHLTOfflineSource") << "no index " << index << " of that name ";
   } else {
     const trigger::Keys& k = triggerObj_->filterKeys(index);
     if (!k.empty())
       output = true;
   }
   return output;
 }

◆ TriggerPosition()

double JetMETHLTOfflineSource::TriggerPosition ( std::string trigName )

privatevirtual

Definition at line 3354 of file JetMETHLTOfflineSource.cc.

References dqm::impl::MonitorElement::getTH1(), LaserClient_cfi::nbins, rate_All, and cscTnPEfficiencyTask_cfi::trigName.

Referenced by fillMEforMonTriggerSummary().

                                                                  {
   int nbins = rate_All->getTH1()->GetNbinsX();
   double binVal = -100;
   for (int ibin = 1; ibin < nbins + 1; ibin++) {
     const char* binLabel = rate_All->getTH1()->GetXaxis()->GetBinLabel(ibin);
     if (binLabel[0] == '\0')
       continue;
     //       std::string binLabel_str = string(binLabel);
     //       if(binLabel_str.compare(trigName)!=0)continue;
     if (trigName != binLabel)
       continue;
 
     if (trigName == binLabel) {
       binVal = rate_All->getTH1()->GetBinCenter(ibin);
       break;
     }
   }
   return binVal;
 }

◆ validPathHLT()

bool JetMETHLTOfflineSource::validPathHLT ( std::string path )

privatevirtual

Definition at line 3330 of file JetMETHLTOfflineSource.cc.

References hltConfig_, dqmiolumiharvest::j, convertSQLitetoXML_cfg::output, BTVHLTOfflineSource_cfi::pathname, HLTConfigProvider::size(), and HLTConfigProvider::triggerName().

                                                             {
   // hltConfig_ has to be defined first before calling this method
   bool output = false;
   for (unsigned int j = 0; j != hltConfig_.size(); ++j) {
     if (hltConfig_.triggerName(j) == pathname)
       output = true;
   }
   return output;
 }

Member Data Documentation

◆ _fEMF

double JetMETHLTOfflineSource::_fEMF

private

Definition at line 112 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforEffAllTrigger(), and JetMETHLTOfflineSource().

◆ _feta

double JetMETHLTOfflineSource::_feta

private

Definition at line 113 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforEffAllTrigger(), and JetMETHLTOfflineSource().

◆ _fHPD

double JetMETHLTOfflineSource::_fHPD

private

Definition at line 114 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforEffAllTrigger(), and JetMETHLTOfflineSource().

◆ _max_CEMF

double JetMETHLTOfflineSource::_max_CEMF

private

Definition at line 124 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforEffAllTrigger(), and JetMETHLTOfflineSource().

◆ _max_CHEF

double JetMETHLTOfflineSource::_max_CHEF

private

Definition at line 120 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforEffAllTrigger(), and JetMETHLTOfflineSource().

◆ _max_NEMF

double JetMETHLTOfflineSource::_max_NEMF

private

Definition at line 122 of file JetMETHLTOfflineSource.h.

Referenced by JetMETHLTOfflineSource().

◆ _max_NHEF

double JetMETHLTOfflineSource::_max_NHEF

private

Definition at line 118 of file JetMETHLTOfflineSource.h.

Referenced by JetMETHLTOfflineSource().

◆ _min_CEMF

double JetMETHLTOfflineSource::_min_CEMF

private

Definition at line 123 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforEffAllTrigger(), and JetMETHLTOfflineSource().

◆ _min_CHEF

double JetMETHLTOfflineSource::_min_CHEF

private

Definition at line 119 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforEffAllTrigger(), and JetMETHLTOfflineSource().

◆ _min_NEMF

double JetMETHLTOfflineSource::_min_NEMF

private

Definition at line 121 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforEffAllTrigger(), and JetMETHLTOfflineSource().

◆ _min_NHEF

double JetMETHLTOfflineSource::_min_NHEF

private

Definition at line 117 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforEffAllTrigger(), and JetMETHLTOfflineSource().

◆ _n90Hits

double JetMETHLTOfflineSource::_n90Hits

private

Definition at line 115 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforEffAllTrigger(), and JetMETHLTOfflineSource().

◆ _pfMHT

double JetMETHLTOfflineSource::_pfMHT

private

Definition at line 116 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforEffAllTrigger().

◆ calojet

reco::CaloJetCollection JetMETHLTOfflineSource::calojet

private

Definition at line 150 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), and fillMEforTriggerNTfired().

◆ calojetColl_

edm::Handle<reco::CaloJetCollection> JetMETHLTOfflineSource::calojetColl_

private

Definition at line 138 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), fillMEforMonAllTrigger(), and fillMEforTriggerNTfired().

◆ CaloJetCorToken_

edm::EDGetTokenT<reco::JetCorrector> JetMETHLTOfflineSource::CaloJetCorToken_

private

Definition at line 143 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and JetMETHLTOfflineSource().

◆ CaloJetEMF

double JetMETHLTOfflineSource::CaloJetEMF[2]

private

Definition at line 162 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforEffAllTrigger().

◆ CaloJetEta

double JetMETHLTOfflineSource::CaloJetEta[2]

private

Definition at line 160 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), and fillMEforMonAllTrigger().

◆ CaloJetfHPD

double JetMETHLTOfflineSource::CaloJetfHPD[2]

private

Definition at line 163 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforEffAllTrigger().

◆ CaloJetn90

double JetMETHLTOfflineSource::CaloJetn90[2]

private

Definition at line 164 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforEffAllTrigger().

◆ CaloJetPhi

double JetMETHLTOfflineSource::CaloJetPhi[2]

private

Definition at line 161 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), and fillMEforMonAllTrigger().

◆ CaloJetPt

double JetMETHLTOfflineSource::CaloJetPt[2]

private

Definition at line 159 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), and fillMEforMonAllTrigger().

◆ CaloJetPx

double JetMETHLTOfflineSource::CaloJetPx[2]

private

Definition at line 157 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforMonAllTrigger().

◆ CaloJetPy

double JetMETHLTOfflineSource::CaloJetPy[2]

private

Definition at line 158 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforMonAllTrigger().

◆ caloJetsToken

edm::EDGetTokenT<reco::CaloJetCollection> JetMETHLTOfflineSource::caloJetsToken

private

Definition at line 134 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and JetMETHLTOfflineSource().

◆ calometColl_

edm::Handle<reco::CaloMETCollection> JetMETHLTOfflineSource::calometColl_

private

Definition at line 139 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), fillMEforMonAllTrigger(), and fillMEforTriggerNTfired().

◆ caloMetToken

edm::EDGetTokenT<reco::CaloMETCollection> JetMETHLTOfflineSource::caloMetToken

private

Definition at line 135 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and JetMETHLTOfflineSource().

◆ correlation_All

MonitorElement* JetMETHLTOfflineSource::correlation_All

private

Definition at line 948 of file JetMETHLTOfflineSource.h.

Referenced by bookHistograms(), and fillMEforMonTriggerSummary().

◆ correlation_AllWrtMB

MonitorElement* JetMETHLTOfflineSource::correlation_AllWrtMB

private

Definition at line 950 of file JetMETHLTOfflineSource.h.

Referenced by bookHistograms(), and fillMEforMonTriggerSummary().

◆ correlation_AllWrtMu

MonitorElement* JetMETHLTOfflineSource::correlation_AllWrtMu

private

Definition at line 949 of file JetMETHLTOfflineSource.h.

Referenced by bookHistograms(), and fillMEforMonTriggerSummary().

◆ custompathnamepairs_

std::vector<std::pair<std::string, std::string> > JetMETHLTOfflineSource::custompathnamepairs_

private

Definition at line 148 of file JetMETHLTOfflineSource.h.

Referenced by dqmBeginRun(), and JetMETHLTOfflineSource().

◆ dirname_

std::string JetMETHLTOfflineSource::dirname_

private

Definition at line 98 of file JetMETHLTOfflineSource.h.

Referenced by bookHistograms(), and JetMETHLTOfflineSource().

◆ hltConfig_

HLTConfigProvider JetMETHLTOfflineSource::hltConfig_

private

Definition at line 152 of file JetMETHLTOfflineSource.h.

Referenced by dqmBeginRun(), fillMEforMonTriggerSummary(), getL1ConditionModuleName(), and validPathHLT().

◆ hltPathsAll_

PathInfoCollection JetMETHLTOfflineSource::hltPathsAll_

private

Definition at line 941 of file JetMETHLTOfflineSource.h.

Referenced by bookHistograms(), dqmBeginRun(), fillMEforMonAllTrigger(), fillMEforMonTriggerSummary(), and fillMEforTriggerNTfired().

◆ hltPathsAllTriggerSummary_

PathInfoCollection JetMETHLTOfflineSource::hltPathsAllTriggerSummary_

private

Definition at line 940 of file JetMETHLTOfflineSource.h.

Referenced by bookHistograms(), and dqmBeginRun().

◆ hltPathsEff_

PathInfoCollection JetMETHLTOfflineSource::hltPathsEff_

private

Definition at line 942 of file JetMETHLTOfflineSource.h.

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

◆ isSetup_

bool JetMETHLTOfflineSource::isSetup_

private

Definition at line 109 of file JetMETHLTOfflineSource.h.

Referenced by bookHistograms(), and dqmBeginRun().

◆ jetID

reco::helper::JetIDHelper* JetMETHLTOfflineSource::jetID

private

Definition at line 102 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), Jet.Jet::jetID(), JetMETHLTOfflineSource(), Jet.Jet::looseJetId(), and ~JetMETHLTOfflineSource().

◆ MBTrigPaths_

std::vector<std::string> JetMETHLTOfflineSource::MBTrigPaths_

private

Definition at line 95 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforMonTriggerSummary(), and JetMETHLTOfflineSource().

◆ MuonTrigPaths_

std::vector<std::string> JetMETHLTOfflineSource::MuonTrigPaths_

private

Definition at line 94 of file JetMETHLTOfflineSource.h.

Referenced by dqmBeginRun(), fillMEforMonTriggerSummary(), and JetMETHLTOfflineSource().

◆ nameForEff_

bool JetMETHLTOfflineSource::nameForEff_

private

Definition at line 110 of file JetMETHLTOfflineSource.h.

Referenced by dqmBeginRun(), and JetMETHLTOfflineSource().

◆ NVertices

MonitorElement* JetMETHLTOfflineSource::NVertices

private

Definition at line 952 of file JetMETHLTOfflineSource.h.

Referenced by bookHistograms(), and fillMEforMonTriggerSummary().

◆ pathFilter_

std::vector<std::string> JetMETHLTOfflineSource::pathFilter_

private

Definition at line 146 of file JetMETHLTOfflineSource.h.

Referenced by dqmBeginRun(), and JetMETHLTOfflineSource().

◆ pathRejectKeyword_

std::vector<std::string> JetMETHLTOfflineSource::pathRejectKeyword_

private

Definition at line 147 of file JetMETHLTOfflineSource.h.

Referenced by dqmBeginRun(), and JetMETHLTOfflineSource().

◆ pfjet

reco::PFJetCollection JetMETHLTOfflineSource::pfjet

private

Definition at line 151 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforEffAllTrigger().

◆ PFJetCEMF

double JetMETHLTOfflineSource::PFJetCEMF[2]

private

Definition at line 174 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforEffAllTrigger().

◆ PFJetCHEF

double JetMETHLTOfflineSource::PFJetCHEF[2]

private

Definition at line 172 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforEffAllTrigger().

◆ pfjetColl_

edm::Handle<reco::PFJetCollection> JetMETHLTOfflineSource::pfjetColl_

private

Definition at line 140 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), and fillMEforMonAllTrigger().

◆ PFJetCorToken_

edm::EDGetTokenT<reco::JetCorrector> JetMETHLTOfflineSource::PFJetCorToken_

private

Definition at line 144 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and JetMETHLTOfflineSource().

◆ PFJetEta

double JetMETHLTOfflineSource::PFJetEta[2]

private

Definition at line 169 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), and fillMEforMonAllTrigger().

◆ PFJetNEMF

double JetMETHLTOfflineSource::PFJetNEMF[2]

private

Definition at line 173 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforEffAllTrigger().

◆ PFJetNHEF

double JetMETHLTOfflineSource::PFJetNHEF[2]

private

Definition at line 171 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforEffAllTrigger().

◆ PFJetPhi

double JetMETHLTOfflineSource::PFJetPhi[2]

private

Definition at line 170 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), and fillMEforMonAllTrigger().

◆ PFJetPt

double JetMETHLTOfflineSource::PFJetPt[2]

private

Definition at line 168 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), and fillMEforMonAllTrigger().

◆ PFJetPx

double JetMETHLTOfflineSource::PFJetPx[2]

private

Definition at line 166 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforMonAllTrigger().

◆ PFJetPy

double JetMETHLTOfflineSource::PFJetPy[2]

private

Definition at line 167 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforMonAllTrigger().

◆ pfJetsToken

edm::EDGetTokenT<reco::PFJetCollection> JetMETHLTOfflineSource::pfJetsToken

private

Definition at line 132 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and JetMETHLTOfflineSource().

◆ pfmetColl_

edm::Handle<reco::PFMETCollection> JetMETHLTOfflineSource::pfmetColl_

private

Definition at line 141 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), and fillMEforMonAllTrigger().

◆ pfMetToken

edm::EDGetTokenT<reco::PFMETCollection> JetMETHLTOfflineSource::pfMetToken

private

Definition at line 133 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and JetMETHLTOfflineSource().

◆ pfMHTx_All

double JetMETHLTOfflineSource::pfMHTx_All

private

Definition at line 176 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforEffAllTrigger().

◆ pfMHTy_All

double JetMETHLTOfflineSource::pfMHTy_All

private

Definition at line 177 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and fillMEforEffAllTrigger().

◆ plotAll_

bool JetMETHLTOfflineSource::plotAll_

private

Definition at line 106 of file JetMETHLTOfflineSource.h.

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

◆ plotEff_

bool JetMETHLTOfflineSource::plotEff_

private

Definition at line 107 of file JetMETHLTOfflineSource.h.

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

◆ prescUsed_

std::vector<int> JetMETHLTOfflineSource::prescUsed_

private

Definition at line 96 of file JetMETHLTOfflineSource.h.

◆ processname_

std::string JetMETHLTOfflineSource::processname_

private

Definition at line 99 of file JetMETHLTOfflineSource.h.

Referenced by dqmBeginRun(), fillMEforEffAllTrigger(), fillMEforMonAllTrigger(), fillMEforTriggerNTfired(), isTriggerObjectFound(), and JetMETHLTOfflineSource().

◆ PVZ

MonitorElement* JetMETHLTOfflineSource::PVZ

private

Definition at line 951 of file JetMETHLTOfflineSource.h.

Referenced by bookHistograms(), and fillMEforMonTriggerSummary().

◆ rate_All

MonitorElement* JetMETHLTOfflineSource::rate_All

private

Definition at line 944 of file JetMETHLTOfflineSource.h.

Referenced by bookHistograms(), fillMEforMonTriggerSummary(), and TriggerPosition().

◆ rate_AllWrtMB

MonitorElement* JetMETHLTOfflineSource::rate_AllWrtMB

private

Definition at line 946 of file JetMETHLTOfflineSource.h.

Referenced by bookHistograms(), and fillMEforMonTriggerSummary().

◆ rate_AllWrtMu

MonitorElement* JetMETHLTOfflineSource::rate_AllWrtMu

private

Definition at line 945 of file JetMETHLTOfflineSource.h.

Referenced by bookHistograms(), and fillMEforMonTriggerSummary().

◆ runStandalone_

bool JetMETHLTOfflineSource::runStandalone_

private

Definition at line 105 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), bookHistograms(), fillMEforMonTriggerSummary(), and JetMETHLTOfflineSource().

◆ triggerNames_

edm::TriggerNames JetMETHLTOfflineSource::triggerNames_

private

Definition at line 154 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), fillMEforTriggerNTfired(), and isHLTPathAccepted().

◆ triggerObj_

edm::Handle<trigger::TriggerEvent> JetMETHLTOfflineSource::triggerObj_

private

Definition at line 155 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), fillMEforMonAllTrigger(), fillMEforTriggerNTfired(), and isTriggerObjectFound().

◆ triggerResults_

edm::Handle<edm::TriggerResults> JetMETHLTOfflineSource::triggerResults_

private

Definition at line 153 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), fillMEforEffAllTrigger(), fillMEforMonAllTrigger(), fillMEforTriggerNTfired(), and isHLTPathAccepted().

◆ triggerResultsFUToken

edm::EDGetTokenT<edm::TriggerResults> JetMETHLTOfflineSource::triggerResultsFUToken

private

Definition at line 129 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and JetMETHLTOfflineSource().

◆ triggerResultsLabel_

edm::InputTag JetMETHLTOfflineSource::triggerResultsLabel_

private

Definition at line 127 of file JetMETHLTOfflineSource.h.

Referenced by JetMETHLTOfflineSource().

◆ triggerResultsToken

edm::EDGetTokenT<edm::TriggerResults> JetMETHLTOfflineSource::triggerResultsToken

private

Definition at line 128 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and JetMETHLTOfflineSource().

◆ triggerSummaryFUToken

edm::EDGetTokenT<trigger::TriggerEvent> JetMETHLTOfflineSource::triggerSummaryFUToken

private

Definition at line 131 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and JetMETHLTOfflineSource().

◆ triggerSummaryLabel_

edm::InputTag JetMETHLTOfflineSource::triggerSummaryLabel_

private

Definition at line 126 of file JetMETHLTOfflineSource.h.

Referenced by JetMETHLTOfflineSource().

◆ triggerSummaryToken

edm::EDGetTokenT<trigger::TriggerEvent> JetMETHLTOfflineSource::triggerSummaryToken

private

Definition at line 130 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), and JetMETHLTOfflineSource().

◆ verbose_

bool JetMETHLTOfflineSource::verbose_

private

Definition at line 104 of file JetMETHLTOfflineSource.h.

Referenced by analyze(), dqmBeginRun(), fillMEforEffAllTrigger(), fillMEforMonAllTrigger(), fillMEforMonTriggerSummary(), fillMEforTriggerNTfired(), and JetMETHLTOfflineSource().

◆ vertexToken

edm::EDGetTokenT<reco::VertexCollection> JetMETHLTOfflineSource::vertexToken

private

Definition at line 136 of file JetMETHLTOfflineSource.h.

Referenced by fillMEforMonTriggerSummary(), and JetMETHLTOfflineSource().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ JetMETHLTOfflineSource()

◆ ~JetMETHLTOfflineSource()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ dqmBeginRun()

◆ fillMEforEffAllTrigger()

◆ fillMEforMonAllTrigger()

◆ fillMEforMonTriggerSummary()

◆ fillMEforTriggerNTfired()

◆ getL1ConditionModuleName()

◆ isBarrel()

◆ isEndCap()

◆ isForward()

◆ isHLTPathAccepted()

◆ isTriggerObjectFound()

◆ TriggerPosition()

◆ validPathHLT()

Member Data Documentation

◆ _fEMF

◆ _feta

◆ _fHPD

◆ _max_CEMF

◆ _max_CHEF

◆ _max_NEMF

◆ _max_NHEF

◆ _min_CEMF

◆ _min_CHEF

◆ _min_NEMF

◆ _min_NHEF

◆ _n90Hits

◆ _pfMHT

◆ calojet

◆ calojetColl_

◆ CaloJetCorToken_

◆ CaloJetEMF

◆ CaloJetEta

◆ CaloJetfHPD

◆ CaloJetn90

◆ CaloJetPhi

◆ CaloJetPt

◆ CaloJetPx

◆ CaloJetPy

◆ caloJetsToken

◆ calometColl_

◆ caloMetToken

◆ correlation_All

◆ correlation_AllWrtMB

◆ correlation_AllWrtMu

◆ custompathnamepairs_

◆ dirname_

◆ hltConfig_

◆ hltPathsAll_

◆ hltPathsAllTriggerSummary_

◆ hltPathsEff_

◆ isSetup_

◆ jetID

◆ MBTrigPaths_

◆ MuonTrigPaths_

◆ nameForEff_

◆ NVertices

◆ pathFilter_

◆ pathRejectKeyword_

◆ pfjet

◆ PFJetCEMF

◆ PFJetCHEF

◆ pfjetColl_

◆ PFJetCorToken_

◆ PFJetEta

◆ PFJetNEMF

◆ PFJetNHEF

◆ PFJetPhi

◆ PFJetPt