#include <HLTJetL1MatchProducer.h>

Inheritance diagram for HLTJetL1MatchProducer:

Public Member Functions
virtual void	beginJob ()

	HLTJetL1MatchProducer (const edm::ParameterSet &)

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

	~HLTJetL1MatchProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Private Attributes
double	DeltaR_

edm::InputTag	jetsInput_

edm::InputTag	L1CenJets_

edm::InputTag	L1ForJets_

edm::InputTag	L1TauJets_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 9 of file HLTJetL1MatchProducer.h.

Constructor & Destructor Documentation

HLTJetL1MatchProducer::HLTJetL1MatchProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 9 of file HLTJetL1MatchProducer.cc.

References DeltaR_, edm::ParameterSet::getParameter(), jetsInput_, L1CenJets_, L1ForJets_, and L1TauJets_.

 {
   jetsInput_ = iConfig.getParameter<edm::InputTag>("jetsInput");
   L1TauJets_ = iConfig.getParameter<edm::InputTag>("L1TauJets");
   L1CenJets_ = iConfig.getParameter<edm::InputTag>("L1CenJets");
   L1ForJets_ = iConfig.getParameter<edm::InputTag>("L1ForJets");
   DeltaR_ = iConfig.getParameter<double>("DeltaR");
 
   produces< reco::CaloJetCollection > ();
 
 }

HLTJetL1MatchProducer::~HLTJetL1MatchProducer ( )

Definition at line 26 of file HLTJetL1MatchProducer.cc.

 {
 
 }

Member Function Documentation

void HLTJetL1MatchProducer::beginJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 21 of file HLTJetL1MatchProducer.cc.

 {
 
 }

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

virtual

Implements edm::EDProducer.

Definition at line 31 of file HLTJetL1MatchProducer.cc.

References Geom::deltaPhi(), DeltaR_, edm::Event::getByLabel(), jetsInput_, L1CenJets_, L1ForJets_, L1TauJets_, edm::Event::put(), query::result, and mathSSE::sqrt().

 {
 
   edm::Handle<reco::CaloJetCollection> calojets;
   iEvent.getByLabel(jetsInput_, calojets);
 
   std::auto_ptr<reco::CaloJetCollection> result (new reco::CaloJetCollection);
 
 
   edm::Handle<l1extra::L1JetParticleCollection> l1TauJets;
   iEvent.getByLabel(L1TauJets_,l1TauJets);
 
   edm::Handle<l1extra::L1JetParticleCollection> l1CenJets;
   iEvent.getByLabel(L1CenJets_,l1CenJets);
 
   edm::Handle<l1extra::L1JetParticleCollection> l1ForJets;
   iEvent.getByLabel(L1ForJets_,l1ForJets);
 
 
   for (reco::CaloJetCollection::const_iterator calojetc = calojets->begin(); 
        calojetc != calojets->end(); ++calojetc) {
 
     bool isMatched=false;
 
     //std::cout << "FL: l1TauJets.size  = " << l1TauJets->size() << std::endl;
     for (unsigned int jetc=0;jetc<l1TauJets->size();++jetc)
     {
       const double deltaeta=calojetc->eta()-(*l1TauJets)[jetc].eta();
       const double deltaphi=deltaPhi(calojetc->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=calojetc->eta()-(*l1CenJets)[jetc].eta();
       const double deltaphi=deltaPhi(calojetc->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=calojetc->eta()-(*l1ForJets)[jetc].eta();
       const double deltaphi=deltaPhi(calojetc->phi(),(*l1ForJets)[jetc].phi());
       if (sqrt(deltaeta*deltaeta+deltaphi*deltaphi) < DeltaR_) isMatched=true;
     }
 
 
     if (isMatched==true) result->push_back(*calojetc);
 
   } // calojetc
 
   iEvent.put( result);
 
 }