Inheritance diagram for PFJetAnalyzerDQM:

Classes
class	Plot1DInBin

class	Plot1DInBinVariable

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

	PFJetAnalyzerDQM (const edm::ParameterSet &)

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

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

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	fillJetResponse (edm::View< pat::Jet > &recoJetCollection, edm::View< reco::Jet > &genJetCollection)

void	prepareGenJetPlots (const std::vector< edm::ParameterSet > &genjet_plots_pset)

void	prepareJetResponsePlots (const std::vector< edm::ParameterSet > &genjet_plots_pset)

Private Attributes
std::vector< Plot1DInBinVariable >	genJetPlots

edm::InputTag	genJetsLabel

bool	genJetsOn

edm::EDGetTokenT< edm::View< reco::Jet > >	genJetsToken

bool	isMC

std::string	jetCollectionName

float	jetDeltaR

std::vector< Plot1DInBin >	jetResponsePlots

std::vector< Plot1DInBin >	jetResponsePlots_noJEC

edm::InputTag	recoJetsLabel

edm::EDGetTokenT< edm::View< pat::Jet > >	recoJetsToken

edm::EDGetTokenT< reco::CandViewMatchMap >	srcRefToJetMap

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 21 of file PFJetAnalyzerDQM.cc.

Constructor & Destructor Documentation

◆ PFJetAnalyzerDQM()

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

Definition at line 191 of file PFJetAnalyzerDQM.cc.

                                                                  {
   recoJetsLabel = iConfig.getParameter<edm::InputTag>("recoJetCollection");
   genJetsLabel = iConfig.getParameter<edm::InputTag>("genJetCollection");
  
   //label for making new folder
   jetCollectionName = recoJetsLabel.label();
  
   //DeltaR for reco to gen jet matching
   jetDeltaR = iConfig.getParameter<double>("jetDeltaR");
  
   //for turn genJet on/off
   genJetsOn = iConfig.getParameter<bool>("genJetsOn");
  
   //Create all jet response plots in bins of genjet pt and eta
   const auto& response_plots = iConfig.getParameter<std::vector<edm::ParameterSet>>("responsePlots");
   prepareJetResponsePlots(response_plots);
  
   const auto& genjet_plots = iConfig.getParameter<std::vector<edm::ParameterSet>>("genJetPlots");
   prepareGenJetPlots(genjet_plots);
  
   recoJetsToken = consumes<edm::View<pat::Jet>>(recoJetsLabel);
   genJetsToken = consumes<edm::View<reco::Jet>>(genJetsLabel);
 }

References genJetsLabel, genJetsOn, genJetsToken, edm::ParameterSet::getParameter(), jetCollectionName, jetDeltaR, edm::InputTag::label(), prepareGenJetPlots(), prepareJetResponsePlots(), recoJetsLabel, and recoJetsToken.

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 281 of file PFJetAnalyzerDQM.cc.

                                                                            {
   edm::Handle<edm::View<pat::Jet>> recoJetCollectionHandle;
   iEvent.getByToken(recoJetsToken, recoJetCollectionHandle);
   auto recoJetCollection = *recoJetCollectionHandle;
  
   isMC = !iEvent.isRealData();
  
   if (isMC) {
     edm::Handle<edm::View<reco::Jet>> genJetCollectionHandle;
     iEvent.getByToken(genJetsToken, genJetCollectionHandle);
     auto genJetCollection = *genJetCollectionHandle;
  
     fillJetResponse(recoJetCollection, genJetCollection);
   }
 }

References fillJetResponse(), particleFlowDQM_cff::genJetCollection, genJetsToken, iEvent, isMC, particleFlowDQM_cff::recoJetCollection, and recoJetsToken.

◆ bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 261 of file PFJetAnalyzerDQM.cc.

                                                                                                   {
   booker.setCurrentFolder("ParticleFlow/JetResponse/" + jetCollectionName + "/JEC/");
   for (auto& plot : jetResponsePlots) {
     plot.book(booker);
   }
   //Book plots for noJEC
   booker.setCurrentFolder("ParticleFlow/JetResponse/" + jetCollectionName + "/noJEC/");
   for (auto& plot : jetResponsePlots_noJEC) {
     plot.book(booker);
   }
  
   //Book plots for gen-jet pt spectra
   if (genJetsOn) {
     booker.setCurrentFolder("ParticleFlow/GenJets/");
     for (auto& plot : genJetPlots) {
       plot.book(booker);
     }
   }
 }

References genJetPlots, genJetsOn, jetCollectionName, jetResponsePlots, jetResponsePlots_noJEC, plotFactory::plot, and dqm::implementation::NavigatorBase::setCurrentFolder().

◆ fillJetResponse()

void PFJetAnalyzerDQM::fillJetResponse	(	edm::View< pat::Jet > &	recoJetCollection,
		edm::View< reco::Jet > &	genJetCollection
	)

private

Definition at line 215 of file PFJetAnalyzerDQM.cc.

                                                                                                              {
   //match gen jets to reco jets, require minimum jetDeltaR, choose closest, do not try to match charge
   std::vector<int> matchIndices;
   PFB::match(genJetCollection, recoJetCollection, matchIndices, false, jetDeltaR);
  
   for (unsigned int i = 0; i < genJetCollection.size(); i++) {
     const auto& genJet = genJetCollection.at(i);
     const auto pt_gen = genJet.pt();
     const auto eta_gen = abs(genJet.eta());
     const int iMatch = matchIndices[i];
  
     //Fill genjet pt if genJetOn
     if (genJetsOn) {
       for (auto& plot : genJetPlots) {
         if (plot.isInEtaBin(eta_gen)) {
           plot.fill(pt_gen);
         }
       }
     }
  
     //If gen jet had a matched reco jet
     if (iMatch != -1) {
       const auto& recoJet = recoJetCollection[iMatch];
       auto pt_reco = recoJet.pt();
  
       const auto response = pt_reco / pt_gen;
       const auto response_raw = pt_reco * recoJet.jecFactor("Uncorrected") / pt_gen;
  
       //Loop linearly through all plots and check if they match the pt and eta bin
       //this is not algorithmically optimal but we don't expect to more than a few hundred plots
       //If this turns out to be a problem, can easily make a 2D-map for indices
       for (auto& plot : jetResponsePlots) {
         if (plot.isInPtEtaBin(pt_gen, eta_gen)) {
           plot.fill(response);
         }
       }
       // this loop should be for NoJEC plots
       for (auto& plot : jetResponsePlots_noJEC) {
         if (plot.isInPtEtaBin(pt_gen, eta_gen)) {
           plot.fill(response_raw);
         }
       }
     }
   }
 }

References funct::abs(), particleFlowDQM_cff::genJetCollection, genJetPlots, genJetsOn, mps_fire::i, jetDeltaR, jetResponsePlots, jetResponsePlots_noJEC, PFB::match(), plotFactory::plot, and particleFlowDQM_cff::recoJetCollection.

Referenced by analyze().

◆ prepareGenJetPlots()

void PFJetAnalyzerDQM::prepareGenJetPlots ( const std::vector< edm::ParameterSet > & genjet_plots_pset )

private

Definition at line 168 of file PFJetAnalyzerDQM.cc.

                                                                                              {
   for (auto& pset : genjet_plots_pset) {
     const auto name = pset.getParameter<std::string>("name");
     const auto title = pset.getParameter<std::string>("title");
  
     //Low and high edges of the eta bins for jets to pass to be filled into this histogram
     const auto ptbins_d = pset.getParameter<std::vector<double>>("ptBins");
     std::vector<float> ptbins(ptbins_d.begin(), ptbins_d.end());
  
     const auto etabin_low = pset.getParameter<double>("etaBinLow");
     const auto etabin_high = pset.getParameter<double>("etaBinHigh");
  
     genJetPlots.push_back(Plot1DInBinVariable(
         name,
         title,
         std::make_unique<TH1F>(name.c_str(), title.c_str(), static_cast<int>(ptbins.size()) - 1, ptbins.data()),
         0.0,
         0.0,
         etabin_low,
         etabin_high));
   }
 }

References genJetPlots, Skims_PA_cff::name, muonDTDigis_cfi::pset, LepHTMonitor_cff::ptbins, AlCaHLTBitMon_QueryRunRegistry::string, and runGCPTkAlMap::title.

Referenced by PFJetAnalyzerDQM().

◆ prepareJetResponsePlots()

void PFJetAnalyzerDQM::prepareJetResponsePlots ( const std::vector< edm::ParameterSet > & genjet_plots_pset )

private

Definition at line 138 of file PFJetAnalyzerDQM.cc.

                                                                                                {
   for (auto& pset : response_plots) {
     //Low and high edges of the pt and eta bins for jets to pass to be filled into this histogram
     const auto ptbin_low = pset.getParameter<double>("ptBinLow");
     const auto ptbin_high = pset.getParameter<double>("ptBinHigh");
     const auto etabin_low = pset.getParameter<double>("etaBinLow");
     const auto etabin_high = pset.getParameter<double>("etaBinHigh");
  
     const auto response_nbins = pset.getParameter<uint32_t>("responseNbins");
     const auto response_low = pset.getParameter<double>("responseLow");
     const auto response_high = pset.getParameter<double>("responseHigh");
  
     const auto name = pset.getParameter<std::string>("name");
     const auto title = pset.getParameter<std::string>("title");
  
     // for title of raw jet response histograms
     auto rawTitle = title;
     rawTitle = rawTitle.replace(rawTitle.begin(), rawTitle.begin(), "Raw ");
  
     jetResponsePlots.push_back(Plot1DInBin(
         name, title, response_nbins, response_low, response_high, ptbin_low, ptbin_high, etabin_low, etabin_high));
  
     jetResponsePlots_noJEC.push_back(Plot1DInBin(
         name, rawTitle, response_nbins, response_low, response_high, ptbin_low, ptbin_high, etabin_low, etabin_high));
   }
   if (jetResponsePlots.size() > 200) {
     throw std::runtime_error("Requested too many jet response plots, aborting as this seems unusual.");
   }
 }