#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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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<>
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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

Member Function Documentation

template<typename T >

void HLTJetL1MatchProducer< T >::beginJob ( void )

virtual

Definition at line 29 of file HLTJetL1MatchProducer.cc.

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(), submitPVResolutionJobs::desc, and HLT_FULL_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 srCondWrite_cfg::deltaPhi, edm::Event::getByToken(), trackerHitRTTI::isMatched(), fwrapper::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));
 }