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

	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

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 >
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 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 ( const edm::ParameterSet & iConfig )

Definition at line 191 of file PFJetAnalyzerDQM.cc.

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

                                                                  {
   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);
 }

Member Function Documentation

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 281 of file PFJetAnalyzerDQM.cc.

References fillJetResponse(), genJetsToken, edm::Event::getByToken(), isMC, edm::EventBase::isRealData(), and recoJetsToken.

                                                                            {
   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);
   }
 }

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 261 of file PFJetAnalyzerDQM.cc.

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

                                                                                                   {
   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);
     }
   }
 }

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

private

Definition at line 215 of file PFJetAnalyzerDQM.cc.

References funct::abs(), edm::View< T >::at(), genJetPlots, genJetsOn, mps_fire::i, jetDeltaR, jetResponsePlots, jetResponsePlots_noJEC, PFB::match(), bigModule::plot(), uploadConditions::response, and edm::View< T >::size().

Referenced by analyze().

                                                                                                              {
   //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);
         }
       }
     }
   }
 }

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

private

Definition at line 168 of file PFJetAnalyzerDQM.cc.

References genJetPlots, mergeVDriftHistosByStation::name, TrackValidation_cff::pset, defaults_cfi::ptbins, AlCaHLTBitMon_QueryRunRegistry::string, and runGCPTkAlMap::title.

Referenced by PFJetAnalyzerDQM().

                                                                                              {
   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));
   }
 }

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

private

Definition at line 138 of file PFJetAnalyzerDQM.cc.

References jetResponsePlots, jetResponsePlots_noJEC, mergeVDriftHistosByStation::name, TrackValidation_cff::pset, AlCaHLTBitMon_QueryRunRegistry::string, and runGCPTkAlMap::title.

Referenced by PFJetAnalyzerDQM().

                                                                                                {
   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.");
   }
 }