Inheritance diagram for FlavorJetCorrectionExample:

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

	FlavorJetCorrectionExample (const edm::ParameterSet &fParameters)

	~FlavorJetCorrectionExample () 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	mBCorrectorName

std::string	mCCorrectorName

edm::InputTag	mInput

std::string	mUDSCorrectorName

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 FlavorJetCorrectionExample.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 36 of file FlavorJetCorrectionExample.cc.

   : mInput (fConfig.getParameter <edm::InputTag> ("src")),
     mUDSCorrectorName (fConfig.getParameter <std::string> ("UDSQuarksCorrector")),
     mCCorrectorName (fConfig.getParameter <std::string> ("CQuarkCorrector")),
     mBCorrectorName (fConfig.getParameter <std::string> ("BQuarkCorrector"))
 {}

FlavorJetCorrectionExample::~FlavorJetCorrectionExample ( )

inlineoverride

Definition at line 14 of file FlavorJetCorrectionExample.cc.

References analyze().

14 {}

Member Function Documentation

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

override

Definition at line 43 of file FlavorJetCorrectionExample.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, mBCorrectorName, mCCorrectorName, mInput, mUDSCorrectorName, reco::LeafCandidate::phi(), and reco::LeafCandidate::pt().

Referenced by ~FlavorJetCorrectionExample().

                                                                                             {
   // get all correctors
   const JetCorrector* udsJetCorrector = JetCorrector::getJetCorrector (mUDSCorrectorName, fSetup);
   const JetCorrector* cQuarkJetCorrector = JetCorrector::getJetCorrector (mCCorrectorName, fSetup);
   const JetCorrector* bQuarkJetCorrector = JetCorrector::getJetCorrector (mBCorrectorName, 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 << "FlavorJetCorrectionExample::analize-> jet #" << ijet;
     if (ijet%3 == 0) { // assume it is light quark
       std::cout << ": use USD quark corrections" << std::endl;
       corrector = udsJetCorrector;
     }
     else if (ijet%3 == 1) { // assume it is c quark
       std::cout << ": use c quark corrections" << std::endl;
       corrector = cQuarkJetCorrector;
     }
     else { // assume it is b quark
       std::cout << ": use b quark corrections" << std::endl;
       corrector = bQuarkJetCorrector;
     }
     // 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;
   }
 }