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	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () 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 overlapproblemtsosanalyzer_cfi::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.");
   }
 }