#include <HLTJetL1MatchProducer.h>

Inheritance diagram for HLTJetL1MatchProducer< T >:

Public Member Functions
virtual void	beginJob ()

	HLTJetL1MatchProducer (const edm::ParameterSet &)

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

	~HLTJetL1MatchProducer () override

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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Attributes
double	DeltaR_

edm::InputTag	jetsInput_

edm::InputTag	L1CenJets_

edm::InputTag	L1ForJets_

edm::InputTag	L1TauJets_

edm::EDGetTokenT< std::vector< T > >	m_theJetToken

edm::EDGetTokenT< l1extra::L1JetParticleCollection >	m_theL1CenJetToken

edm::EDGetTokenT< l1extra::L1JetParticleCollection >	m_theL1ForJetToken

edm::EDGetTokenT< l1extra::L1JetParticleCollection >	m_theL1TauJetToken

Additional Inherited Members
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

Detailed Description

template<typename T>
class HLTJetL1MatchProducer< T >

Definition at line 20 of file HLTJetL1MatchProducer.h.

Constructor & Destructor Documentation

template<typename T >

HLTJetL1MatchProducer< T >::HLTJetL1MatchProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 13 of file HLTJetL1MatchProducer.cc.

                                                                               {
   jetsInput_ = iConfig.template getParameter<edm::InputTag>("jetsInput");
   L1TauJets_ = iConfig.template getParameter<edm::InputTag>("L1TauJets");
   L1CenJets_ = iConfig.template getParameter<edm::InputTag>("L1CenJets");
   L1ForJets_ = iConfig.template getParameter<edm::InputTag>("L1ForJets");
   DeltaR_ = iConfig.template getParameter<double>("DeltaR");
 
   typedef std::vector<T> TCollection;
   m_theJetToken = consumes<TCollection>(jetsInput_);
   m_theL1TauJetToken = consumes<l1extra::L1JetParticleCollection>(L1TauJets_);
   m_theL1CenJetToken = consumes<l1extra::L1JetParticleCollection>(L1CenJets_);
   m_theL1ForJetToken = consumes<l1extra::L1JetParticleCollection>(L1ForJets_);
   produces<TCollection>();
 }

template<typename T >

HLTJetL1MatchProducer< T >::~HLTJetL1MatchProducer ( )

overridedefault

Referenced by HLTJetL1MatchProducer< T >::beginJob().

Member Function Documentation

template<typename T >

void HLTJetL1MatchProducer< T >::beginJob ( void )

virtual

Definition at line 29 of file HLTJetL1MatchProducer.cc.

References HLTJetL1MatchProducer< T >::~HLTJetL1MatchProducer().

29 {}

template<typename T >

void HLTJetL1MatchProducer< T >::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 35 of file HLTJetL1MatchProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), defaultModuleLabel(), and HLT_2018_cff::InputTag.

                                                                                           {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("jetsInput", edm::InputTag("hltAntiKT5PFJets"));
   desc.add<edm::InputTag>("L1TauJets", edm::InputTag("hltL1extraParticles", "Tau"));
   desc.add<edm::InputTag>("L1CenJets", edm::InputTag("hltL1extraParticles", "Central"));
   desc.add<edm::InputTag>("L1ForJets", edm::InputTag("hltL1extraParticles", "Forward"));
   desc.add<double>("DeltaR", 0.5);
   descriptions.add(defaultModuleLabel<HLTJetL1MatchProducer<T>>(), desc);
 }

template<typename T >

void HLTJetL1MatchProducer< T >::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 46 of file HLTJetL1MatchProducer.cc.

References SiPixelRawToDigiRegional_cfi::deltaPhi, edm::Event::getByToken(), trackerHitRTTI::isMatched(), singleTopDQM_cfi::jets, eostools::move(), edm::Event::put(), mps_fire::result, and mathSSE::sqrt().

                                                                                     {
   typedef std::vector<T> TCollection;
 
   edm::Handle<TCollection> jets;
   iEvent.getByToken(m_theJetToken, jets);
 
   std::unique_ptr<TCollection> result(new TCollection);
 
   edm::Handle<l1extra::L1JetParticleCollection> l1TauJets;
   iEvent.getByToken(m_theL1TauJetToken, l1TauJets);
 
   edm::Handle<l1extra::L1JetParticleCollection> l1CenJets;
   iEvent.getByToken(m_theL1CenJetToken, l1CenJets);
 
   edm::Handle<l1extra::L1JetParticleCollection> l1ForJets;
   iEvent.getByToken(m_theL1ForJetToken, l1ForJets);
 
   typename TCollection::const_iterator jet_iter;
   for (jet_iter = jets->begin(); jet_iter != jets->end(); ++jet_iter) {
     bool isMatched = false;
 
     //std::cout << "FL: l1TauJets.size  = " << l1TauJets->size() << std::endl;
     for (unsigned int jetc = 0; jetc < l1TauJets->size(); ++jetc) {
       const double deltaeta = jet_iter->eta() - (*l1TauJets)[jetc].eta();
       const double deltaphi = deltaPhi(jet_iter->phi(), (*l1TauJets)[jetc].phi());
       //std::cout << "FL: sqrt(2) = " << sqrt(2) << std::endl;
       if (sqrt(deltaeta * deltaeta + deltaphi * deltaphi) < DeltaR_)
         isMatched = true;
     }
 
     for (unsigned int jetc = 0; jetc < l1CenJets->size(); ++jetc) {
       const double deltaeta = jet_iter->eta() - (*l1CenJets)[jetc].eta();
       const double deltaphi = deltaPhi(jet_iter->phi(), (*l1CenJets)[jetc].phi());
       if (sqrt(deltaeta * deltaeta + deltaphi * deltaphi) < DeltaR_)
         isMatched = true;
     }
 
     for (unsigned int jetc = 0; jetc < l1ForJets->size(); ++jetc) {
       const double deltaeta = jet_iter->eta() - (*l1ForJets)[jetc].eta();
       const double deltaphi = deltaPhi(jet_iter->phi(), (*l1ForJets)[jetc].phi());
       if (sqrt(deltaeta * deltaeta + deltaphi * deltaphi) < DeltaR_)
         isMatched = true;
     }
 
     if (isMatched == true)
       result->push_back(*jet_iter);
 
   }  // jet_iter
 
   iEvent.put(std::move(result));
 }