Inheritance diagram for PartonJetCorrectionExample:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

	PartonJetCorrectionExample (const edm::ParameterSet &fParameters)

	~PartonJetCorrectionExample () override

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

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)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
std::string	m_bJ_CorrectorName

std::string	m_bT_CorrectorName

std::string	m_gJ_CorrectorName

std::string	m_qJ_CorrectorName

edm::InputTag	mInput

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

Definition at line 11 of file PartonJetCorrectionExample.cc.

Constructor & Destructor Documentation

PartonJetCorrectionExample::PartonJetCorrectionExample ( const edm::ParameterSet & fParameters )

explicit

Definition at line 37 of file PartonJetCorrectionExample.cc.

   : mInput (fConfig.getParameter <edm::InputTag> ("src")),
     m_gJ_CorrectorName (fConfig.getParameter <std::string> ("gJetCorrector")),
     m_qJ_CorrectorName (fConfig.getParameter <std::string> ("qJetCorrector")),
     m_bJ_CorrectorName (fConfig.getParameter <std::string> ("bJetCorrector")),
     m_bT_CorrectorName (fConfig.getParameter <std::string> ("bTopCorrector"))
 {}

PartonJetCorrectionExample::~PartonJetCorrectionExample ( )

inlineoverride

Definition at line 14 of file PartonJetCorrectionExample.cc.

References analyze().

14 {}

Member Function Documentation

void PartonJetCorrectionExample::analyze	(	const edm::Event &	fEvent,
		const edm::EventSetup &	fSetup
	)

override

Definition at line 45 of file PartonJetCorrectionExample.cc.

References reco::JetCorrector::correction(), mitigatedMETSequence_cff::corrector, gather_cfg::cout, DEFINE_FWK_MODULE, reco::LeafCandidate::eta(), edm::Event::getByLabel(), JetCorrector::getJetCorrector(), metsig::jet, fwrapper::jets, m_bJ_CorrectorName, m_bT_CorrectorName, m_gJ_CorrectorName, m_qJ_CorrectorName, mInput, reco::LeafCandidate::phi(), and reco::LeafCandidate::pt().

Referenced by ~PartonJetCorrectionExample().

                                                                                             {
   // get all correctors
   const JetCorrector* gJetCorrector = JetCorrector::getJetCorrector (m_gJ_CorrectorName, fSetup);
   const JetCorrector* qJetCorrector = JetCorrector::getJetCorrector (m_qJ_CorrectorName, fSetup);
   const JetCorrector* bJetCorrector = JetCorrector::getJetCorrector (m_bJ_CorrectorName, fSetup);
   const JetCorrector* bTopCorrector = JetCorrector::getJetCorrector (m_bT_CorrectorName, fSetup);
   const JetCorrector* corrector = nullptr;
   
   // get input jets (supposed to be MC corrected already)
   edm::Handle<CaloJetCollection> jets;                    
   fEvent.getByLabel (mInput, jets);
   // loop over jets                      
   for (unsigned ijet = 0; ijet < jets->size(); ++ijet) {
     const CaloJet& jet = (*jets)[ijet];
     std::cout << "PartonJetCorrectionExample::analize-> jet #" << ijet;
     if (ijet%4 == 0) { // assume it is gluon from diJet
       std::cout << ": use gJ corrections" << std::endl;
       corrector = gJetCorrector;
     }
     else if (ijet%4 == 1) { // assume it is light quark from diJet
       std::cout << ": use qJ corrections" << std::endl;
       corrector = qJetCorrector;
     }
     else if (ijet%4 == 2) { // assume it is b quark from diJet
       std::cout << ": use bJ corrections" << std::endl;
       corrector = bJetCorrector;
     }
     else { // assume it is b quark from ttbar
       std::cout << ": use bT corrections" << std::endl;
       corrector = bTopCorrector;
     }
     // get selected correction for the jet
     double correction = corrector->correction (jet);
     // dump it
     std::cout << "  jet pt/eta/phi: " << jet.pt() << '/' <<  jet.eta() << '/' << jet.phi() 
           << " -> correction factor: " << correction 
           << ", corrected pt: " << jet.pt()*correction
           << std::endl;
   }
 }