reco::CorrectedJetProducer< T > Class Template Reference

#include <CorrectedJetProducer.h>

Inheritance diagram for reco::CorrectedJetProducer< T >:

Public Types
typedef std::vector< T >	JetCollection
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
	CorrectedJetProducer (const edm::ParameterSet &fParameters)
virtual void	produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override
virtual	~CorrectedJetProducer ()
Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default
Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducerBase ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Attributes
const std::vector < edm::EDGetTokenT < reco::JetCorrector > >	mCorrectorTokens
const edm::EDGetTokenT < JetCollection >	mInput
const bool	mVerbose

Additional Inherited Members
Static Public Member Functions inherited from edm::global::EDProducerBase
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)
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

template<class T>
class reco::CorrectedJetProducer< T >

Definition at line 28 of file CorrectedJetProducer.h.

Member Typedef Documentation

template<class T >

typedef std::vector<T> reco::CorrectedJetProducer< T >::JetCollection

Definition at line 31 of file CorrectedJetProducer.h.

Constructor & Destructor Documentation

template<class T >

reco::CorrectedJetProducer< T >::CorrectedJetProducer ( const edm::ParameterSet &  fParameters )

explicit

Definition at line 47 of file CorrectedJetProducer.h.

References GlobalPosition_Frontier_DevDB_cff::tag.

    : mInput(consumes<JetCollection>(fConfig.getParameter <edm::InputTag> ("src")))
    , mCorrectorTokens(edm::vector_transform(fConfig.getParameter<std::vector<edm::InputTag> >("correctors"), [this](edm::InputTag const & tag){return consumes<reco::JetCorrector>(tag);}))

template<class T >

virtual reco::CorrectedJetProducer< T >::~CorrectedJetProducer ( )

inlinevirtual

Definition at line 33 of file CorrectedJetProducer.h.

{}

Member Function Documentation

template<class T >

void reco::CorrectedJetProducer< T >::produce ( edm::StreamID  , edm::Event &  fEvent, const edm::EventSetup &  fSetup  ) const

overridevirtual

Implements edm::global::EDProducerBase.

Definition at line 60 of file CorrectedJetProducer.h.

References HLT_25ns14e33_v1_cff::correctors, gather_cfg::cout, edm::Event::getByToken(), patZpeak::handle, i, cmsHarvester::index, metsig::jet, fwrapper::jets, reco::mVerbose(), edm::Handle< T >::product(), edm::Event::put(), query::result, pileupReCalc_HLTpaths::scale, and python.multivaluedict::sort().

  {
    // FIXME - use something more efficient instead of an std::vector
    std::vector<reco::JetCorrector const *> correctors(mCorrectorTokens.size(), nullptr);

    // look for correctors
    for (unsigned i = 0; i < mCorrectorTokens.size(); i++)
      {
        edm::Handle <reco::JetCorrector> handle;
        fEvent.getByToken (mCorrectorTokens [i], handle);
        correctors[i] = handle.product();
      }
    edm::Handle<JetCollection> jets;                         //Define Inputs
    fEvent.getByToken (mInput, jets);                        //Get Inputs
    std::auto_ptr<JetCollection> result (new JetCollection); //Corrected jets
    typename JetCollection::const_iterator jet;
    for (jet = jets->begin(); jet != jets->end(); jet++)
      {
    const T* referenceJet = &*jet;
    int index = jet-jets->begin();
    edm::RefToBase<reco::Jet> jetRef(edm::Ref<JetCollection>(jets,index));
    T correctedJet = *jet; //copy original jet
    if (mVerbose)
      std::cout<<"CorrectedJetProducer::produce-> original jet: "
           <<jet->print()<<std::endl;
    for (unsigned i = 0; i < mCorrectorTokens.size(); ++i)
      {
        if ( !(correctors[i]->vectorialCorrection()) ) {
          // Scalar correction
              double scale = 1.;
              if (!(correctors[i]->refRequired()))
            scale = correctors[i]->correction (*referenceJet);
              else
                scale = correctors[i]->correction (*referenceJet,jetRef);
          if (mVerbose)
        std::cout<<"CorrectedJetProducer::produce-> Corrector # "
             <<i<<", correction factor: "<<scale<<std::endl;
          correctedJet.scaleEnergy (scale); // apply scalar correction
          referenceJet = &correctedJet;
        } else {
          // Vectorial correction
          reco::JetCorrector::LorentzVector corrected;
          double scale = correctors[i]->correction (*referenceJet, jetRef, corrected);
          if (mVerbose)
        std::cout<<"CorrectedJetProducer::produce-> Corrector # "
             <<i<<", correction factor: "<<scale<<std::endl;
          correctedJet.setP4( corrected ); // apply vectorial correction
          referenceJet = &correctedJet;
        }
      }
    if (mVerbose)
      std::cout<<"CorrectedJetProducer::produce-> corrected jet: "
           <<correctedJet.print ()<<std::endl;
    result->push_back (correctedJet);
      }
    NumericSafeGreaterByPt<T> compJets;
    // reorder corrected jets
    std::sort (result->begin (), result->end (), compJets);
    // put corrected jet collection into event
    fEvent.put(result);
  }

Member Data Documentation

template<class T >

const std::vector<edm::EDGetTokenT<reco::JetCorrector> > reco::CorrectedJetProducer< T >::mCorrectorTokens

private

Definition at line 37 of file CorrectedJetProducer.h.

template<class T >

const edm::EDGetTokenT<JetCollection> reco::CorrectedJetProducer< T >::mInput

private

Definition at line 36 of file CorrectedJetProducer.h.

template<class T >

const bool reco::CorrectedJetProducer< T >::mVerbose

private

Definition at line 38 of file CorrectedJetProducer.h.