#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	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () 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

◆ HLTJetL1MatchProducer()

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>();
 }

◆ ~HLTJetL1MatchProducer()

template<typename T >

HLTJetL1MatchProducer< T >::~HLTJetL1MatchProducer ( )

overridedefault

Member Function Documentation

◆ beginJob()

template<typename T >

void HLTJetL1MatchProducer< T >::beginJob ( void )

virtual

Definition at line 29 of file HLTJetL1MatchProducer.cc.

29 {}

◆ fillDescriptions()

template<typename T >

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

static

Definition at line 35 of file HLTJetL1MatchProducer.cc.

                                                                                           {
   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);
 }

References edm::ConfigurationDescriptions::add(), defaultModuleLabel(), submitPVResolutionJobs::desc, and HLT_FULL_cff::InputTag.

◆ produce()

template<typename T >

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

override

Definition at line 46 of file HLTJetL1MatchProducer.cc.

                                                                                     {
   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));
 }