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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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

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::vector< ModuleDescription const * > &modules, 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

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::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

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<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)

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 43 of file PatJetAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

default contructor

Definition at line 72 of file PatJetAnalyzer.cc.

References BINS, edmScanValgrind::buffer, hists_, training_settings::idx, TFileService::make(), MAXBIN, and runGCPTkAlMap::title.

                                                         :
   corrLevel_(cfg.getParameter<std::string>("corrLevel")),
   jetsToken_(consumes<edm::View<pat::Jet> >(cfg.getParameter<edm::InputTag>("src")))
 {
   // 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 ( )

inlineoverride

default destructor

Definition at line 49 of file PatJetAnalyzer.cc.

References analyze(), and GeneralSetup::setup().

49 {};

Member Function Documentation

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

overrideprivate

everything that needs to be done during the even loop

Definition at line 98 of file PatJetAnalyzer.cc.

References BINS, edmScanValgrind::buffer, corrLevel_, fill(), training_settings::idx, metsig::jet, fwrapper::jets, jetsToken_, ResonanceBuilder::mass, MAXBIN, and p4.

Referenced by ~PatJetAnalyzer().

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

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

inlineprivate

check if histogram was booked

Definition at line 56 of file PatJetAnalyzer.cc.

References hists_.

Referenced by fill().

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

PatJetAnalyzer::hists_

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

management of 1d histograms

Definition: PatJetAnalyzer.cc:68

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

inlineprivate

fill histogram if it had been booked before

Definition at line 58 of file PatJetAnalyzer.cc.

References booked(), hists_, training_settings::idx, metsig::jet, and print().

Referenced by analyze().

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

PatJetAnalyzer::hists_

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

management of 1d histograms

Definition: PatJetAnalyzer.cc:68

value

Definition: value.py:1

PatJetAnalyzer::booked

bool booked(const std::string histName) const

check if histogram was booked

Definition: PatJetAnalyzer.cc:56

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

private

Definition at line 131 of file PatJetAnalyzer.cc.

References gather_cfg::cout, pat::Jet::currentJECLevel(), DEFINE_FWK_MODULE, training_settings::idx, pat::JetCorrFactors::NONE, reco::LeafCandidate::pt(), edm::View< T >::size(), and pat::JetCorrFactors::UDS.

Referenced by fill().

 {
   //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

std::string PatJetAnalyzer::corrLevel_

private

correction level for pat jet

Definition at line 64 of file PatJetAnalyzer.cc.

Referenced by analyze().

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

private

management of 1d histograms

Definition at line 68 of file PatJetAnalyzer.cc.

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

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

private

pat jets

Definition at line 66 of file PatJetAnalyzer.cc.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation