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 final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const 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

bool	wantsStreamRuns () const

	~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

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

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices 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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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_