#include <JetCorrectionProducer.h>

Inheritance diagram for cms::JetCorrectionProducer< T >:

Public Types
typedef std::vector< T >	JetCollection

Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Member Functions
	JetCorrectionProducer (const edm::ParameterSet &fParameters)

void	produce (edm::Event &, const edm::EventSetup &) override

	~JetCorrectionProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Attributes
unsigned long long	mCacheId

std::vector< std::string >	mCorrectorNames

std::vector< const JetCorrector * >	mCorrectors

edm::EDGetTokenT< JetCollection >	mInput

bool	mVerbose

Detailed Description

template<class T>
class cms::JetCorrectionProducer< T >

Definition at line 31 of file JetCorrectionProducer.h.

Member Typedef Documentation

template<class T >

typedef std::vector<T> cms::JetCorrectionProducer< T >::JetCollection

Definition at line 34 of file JetCorrectionProducer.h.

Constructor & Destructor Documentation

template<class T >

cms::JetCorrectionProducer< T >::JetCorrectionProducer ( const edm::ParameterSet & fParameters )

explicit

Definition at line 53 of file JetCorrectionProducer.h.

References egammaCTFFinalFitWithMaterial_cff::alias, edm::ParameterSet::getUntrackedParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

     : mInput(consumes<JetCollection>(fConfig.getParameter <edm::InputTag> ("src")))
     , mCorrectorNames(fConfig.getParameter<std::vector<std::string> >("correctors"))
     , mCorrectors(mCorrectorNames.size(), nullptr)
     , mCacheId (0)
     , mVerbose (fConfig.getUntrackedParameter <bool> ("verbose", false))
   {
     std::string alias = fConfig.getUntrackedParameter <std::string> ("alias", "");
     if (alias.empty ())
       produces <JetCollection>();
     else
       produces <JetCollection>().setBranchAlias(alias);
   }

template<class T >

cms::JetCorrectionProducer< T >::~JetCorrectionProducer ( )

inlineoverride

Definition at line 36 of file JetCorrectionProducer.h.

36 {}

Member Function Documentation

template<class T >

void cms::JetCorrectionProducer< T >::produce	(	edm::Event &	fEvent,
		const edm::EventSetup &	fSetup
	)

override

Definition at line 68 of file JetCorrectionProducer.h.

References edm::eventsetup::EventSetupRecord::cacheIdentifier(), gather_cfg::cout, edm::EventSetup::get(), edm::eventsetup::EventSetupRecordImplementation< T >::get(), edm::Event::getByToken(), cmsBatch::handle, mps_fire::i, metsig::jet, fwrapper::jets, cms::JetCorrectionProducer< T >::mCacheId, cms::JetCorrectionProducer< T >::mCorrectorNames, cms::JetCorrectionProducer< T >::mCorrectors, cms::JetCorrectionProducer< T >::mInput, eostools::move(), cms::JetCorrectionProducer< T >::mVerbose, edm::Event::put(), record, mps_fire::result, Scenarios_cff::scale, and jetUpdater_cfi::sort.

   {
     // look for correctors
     const JetCorrectionsRecord& record = fSetup.get <JetCorrectionsRecord> ();
     if (record.cacheIdentifier() != mCacheId)
       { // need to renew cache
     for (unsigned i = 0; i < mCorrectorNames.size(); i++)
       {
         edm::ESHandle <JetCorrector> handle;
         record.get (mCorrectorNames [i], handle);
         mCorrectors [i] = &*handle;
       }
     mCacheId = record.cacheIdentifier();
       }
     edm::Handle<JetCollection> jets;                         //Define Inputs
     fEvent.getByToken (mInput, jets);                        //Get Inputs
     std::unique_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<<"JetCorrectionProducer::produce-> original jet: "
            <<jet->print()<<std::endl;
     for (unsigned i = 0; i < mCorrectors.size(); ++i)
       {
         if ( !(mCorrectors[i]->vectorialCorrection()) ) {
           // Scalar correction
               double scale = 1.;
               if (!(mCorrectors[i]->refRequired()))
             scale = mCorrectors[i]->correction (*referenceJet,fEvent,fSetup);
               else
                 scale = mCorrectors[i]->correction (*referenceJet,jetRef,fEvent,fSetup);
           if (mVerbose)
         std::cout<<"JetCorrectionProducer::produce-> Corrector # "
              <<i<<", correction factor: "<<scale<<std::endl;
           correctedJet.scaleEnergy (scale); // apply scalar correction
           referenceJet = &correctedJet;
         } else {
           // Vectorial correction
           JetCorrector::LorentzVector corrected;
           double scale = mCorrectors[i]->correction (*referenceJet, jetRef,
                              fEvent, fSetup, corrected);
           if (mVerbose)
         std::cout<<"JetCorrectionProducer::produce-> Corrector # "
              <<i<<", correction factor: "<<scale<<std::endl;
           correctedJet.setP4( corrected ); // apply vectorial correction
           referenceJet = &correctedJet;
         }
       }
     if (mVerbose)
       std::cout<<"JetCorrectionProducer::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(std::move(result));
   }