Module to analyze pat::Jets in the context of a more complex exercise. More...

Inheritance diagram for PatJetAnalyzer:

Public Member Functions
	PatJetAnalyzer (const edm::ParameterSet &cfg)
	default contructor More...

	~PatJetAnalyzer () override
	default destructor More...

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup) override
	everything that needs to be done during the even loop More...

bool	booked (const std::string histName) const
	check if histogram was booked More...

void	fill (const std::string histName, double value) const
	fill histogram if it had been booked before More...

void	print (edm::View< pat::Jet >::const_iterator &jet, unsigned int idx)

Private Attributes
std::string	corrLevel_
	correction level for pat jet More...

std::map< std::string, TH1F * >	hists_
	management of 1d histograms More...

edm::EDGetTokenT< edm::View< pat::Jet > >	jetsToken_
	pat jets More...

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Module to analyze pat::Jets in the context of a more complex exercise.

Basic quantities of jets like the transverse momentum, eta and phi as well as the invariant dijet mass are plotted. Basic histograms for a jet energy response plot as a function of the pt of the reference object are filled. As reference matched partons are chosen. Input parameters are:

src –> input for the pat jet collection (edm::InputTag).
corrLevel –> string for the pat jet correction level.

Definition at line 42 of file PatJetAnalyzer.cc.

Constructor & Destructor Documentation

◆ PatJetAnalyzer()

PatJetAnalyzer::PatJetAnalyzer ( const edm::ParameterSet & cfg )

explicit

default contructor

Definition at line 72 of file PatJetAnalyzer.cc.

References BINS, edmScanValgrind::buffer, compareTotals::fs, hists_, heavyIonCSV_trainingSettings::idx, TFileService::kSharedResource, MAXBIN, and runGCPTkAlMap::title.

     : corrLevel_(cfg.getParameter<std::string>("corrLevel")),
       jetsToken_(consumes<edm::View<pat::Jet> >(cfg.getParameter<edm::InputTag>("src"))) {
   usesResource(TFileService::kSharedResource);
   // register TFileService
   edm::Service<TFileService> fs;
 
   // jet multiplicity
   hists_["mult"] = fs->make<TH1F>("mult", "N_{Jet}", 15, 0., 15.);
   // jet pt (for all jets)
   hists_["pt"] = fs->make<TH1F>("pt", "p_{T}(Jet) [GeV]", 60, 0., 300.);
   // jet eta (for all jets)
   hists_["eta"] = fs->make<TH1F>("eta", "#eta (Jet)", 60, -3., 3.);
   // jet phi (for all jets)
   hists_["phi"] = fs->make<TH1F>("phi", "#phi (Jet)", 60, 3.2, 3.2);
   // dijet mass (if available)
   hists_["mass"] = fs->make<TH1F>("mass", "M_{jj} [GeV]", 50, 0., 500.);
   // basic histograms for jet energy response
   for (unsigned int idx = 0; idx < MAXBIN; ++idx) {
     char buffer[10];
     sprintf(buffer, "jes_%i", idx);
     char title[50];
     sprintf(title, "p_{T}^{rec}/p_{T}^{gen} [%i GeV - %i GeV]", (int)BINS[idx], (int)BINS[idx + 1]);
     hists_[buffer] = fs->make<TH1F>(buffer, title, 100, 0., 2.);
   }
 }

◆ ~PatJetAnalyzer()

PatJetAnalyzer::~PatJetAnalyzer ( )

inlineoverride

default destructor

Definition at line 47 of file PatJetAnalyzer.cc.

47 {};

Member Function Documentation

◆ analyze()

void PatJetAnalyzer::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overrideprivatevirtual

everything that needs to be done during the even loop

Implements edm::one::EDAnalyzerBase.

Definition at line 99 of file PatJetAnalyzer.cc.

References BINS, edmScanValgrind::buffer, corrLevel_, fill(), heavyIonCSV_trainingSettings::idx, metsig::jet, PDWG_EXODelayedJetMET_cff::jets, jetsToken_, EgHLTOffHistBins_cfi::mass, and MAXBIN.

                                                                               {
   // recieve jet collection label
   edm::Handle<edm::View<pat::Jet> > jets;
   event.getByToken(jetsToken_, jets);
 
   // loop jets
   for (edm::View<pat::Jet>::const_iterator jet = jets->begin(); jet != jets->end(); ++jet) {
     // print jec factors
     // print(jet, jet-jets->begin());
 
     // fill basic kinematics
     fill("pt", jet->correctedJet(corrLevel_).pt());
     fill("eta", jet->eta());
     fill("phi", jet->phi());
     // basic plots for jet responds plot as a function of pt
     if (jet->genJet()) {
       double resp = jet->correctedJet(corrLevel_).pt() / jet->genJet()->pt();
       for (unsigned int idx = 0; idx < MAXBIN; ++idx) {
         if (BINS[idx] <= jet->genJet()->pt() && jet->genJet()->pt() < BINS[idx + 1]) {
           char buffer[10];
           sprintf(buffer, "jes_%i", idx);
           fill(buffer, resp);
         }
       }
     }
   }
   // jet multiplicity
   fill("mult", jets->size());
   // invariant dijet mass for first two leading jets
   if (jets->size() > 1) {
     fill("mass", ((*jets)[0].p4() + (*jets)[1].p4()).mass());
   }
 }

◆ booked()

bool PatJetAnalyzer::booked ( const std::string histName ) const

inlineprivate

check if histogram was booked

Definition at line 54 of file PatJetAnalyzer.cc.

References hists_.

Referenced by fill().

54 { return hists_.find(histName) != hists_.end(); };

PatJetAnalyzer::hists_

std::map< std::string, TH1F * > hists_

management of 1d histograms

Definition: PatJetAnalyzer.cc:69

◆ fill()

void PatJetAnalyzer::fill	(	const std::string	histName,
		double	value
	)		const

inlineprivate

fill histogram if it had been booked before

Definition at line 56 of file PatJetAnalyzer.cc.

References booked(), and hists_.

Referenced by analyze().

                                                           {
     if (booked(histName))
       hists_.find(histName)->second->Fill(value);
   };

◆ print()

void PatJetAnalyzer::print	(	edm::View< pat::Jet >::const_iterator &	jet,
		unsigned int	idx
	)

private

Definition at line 133 of file PatJetAnalyzer.cc.

References gather_cfg::cout, pat::Jet::currentJECLevel(), heavyIonCSV_trainingSettings::idx, metsig::jet, pat::JetCorrFactors::NONE, reco::LeafCandidate::pt(), and pat::JetCorrFactors::UDS.

                                                                                {
   //edm::LogInfo log("JEC");
   std::cout << "[" << idx << "] :: eta=" << std::setw(10) << jet->eta() << " phi=" << std::setw(10) << jet->phi()
             << " size: " << jet->availableJECLevels().size() << std::endl;
   for (unsigned int idx = 0; idx < jet->availableJECLevels().size(); ++idx) {
     pat::Jet correctedJet;
     if (jet->availableJECLevels()[idx].find("L5Flavor") != std::string::npos ||
         jet->availableJECLevels()[idx].find("L7Parton") != std::string::npos) {
       correctedJet = jet->correctedJet(idx, pat::JetCorrFactors::UDS);
     } else {
       correctedJet = jet->correctedJet(idx, pat::JetCorrFactors::NONE);
     }
     std::cout << std::setw(10) << correctedJet.currentJECLevel() << " pt=" << std::setw(10) << correctedJet.pt()
               << std::endl;
   }
 }

Member Data Documentation

◆ corrLevel_

std::string PatJetAnalyzer::corrLevel_

private

correction level for pat jet

Definition at line 65 of file PatJetAnalyzer.cc.

Referenced by analyze().

◆ hists_

std::map<std::string, TH1F*> PatJetAnalyzer::hists_

private

management of 1d histograms

Definition at line 69 of file PatJetAnalyzer.cc.

Referenced by booked(), fill(), and PatJetAnalyzer().

◆ jetsToken_

edm::EDGetTokenT<edm::View<pat::Jet> > PatJetAnalyzer::jetsToken_

private

pat jets

Definition at line 67 of file PatJetAnalyzer.cc.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PatJetAnalyzer()

◆ ~PatJetAnalyzer()

Member Function Documentation

◆ analyze()

◆ booked()

◆ fill()

◆ print()

Member Data Documentation

◆ corrLevel_

◆ hists_

◆ jetsToken_