#include <HLTJetCollectionsForBoostedLeptonPlusJets.h>

Inheritance diagram for HLTJetCollectionsForBoostedLeptonPlusJets< jetType >:

Public Member Functions
	HLTJetCollectionsForBoostedLeptonPlusJets (const edm::ParameterSet &)

	~HLTJetCollectionsForBoostedLeptonPlusJets ()

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

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Private Attributes
edm::InputTag	hltLeptonTag

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

edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	m_theLeptonToken

double	minDeltaR_

edm::InputTag	sourceJetTag

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

template<typename jetType>
class HLTJetCollectionsForBoostedLeptonPlusJets< jetType >

This class is an EDProducer implementing an HLT trigger for lepton and jet objects, cutting on variables relating to the jet 4-momentum representation. The producer checks for overlaps between leptons and jets and if a combination of one lepton + jets cleaned against this leptons satisfy the cuts. These jets are then added to a cleaned jet collection which is put into the event.

Author: Lukasz Kreczko

Definition at line 40 of file HLTJetCollectionsForBoostedLeptonPlusJets.h.

Constructor & Destructor Documentation

template<typename jetType >

HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::HLTJetCollectionsForBoostedLeptonPlusJets ( const edm::ParameterSet & iConfig )

explicit

Definition at line 25 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

References HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::hltLeptonTag, HLT_25ns14e33_v1_cff::jetType, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theJetToken, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theLeptonToken, and HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::sourceJetTag.

                                                                                                                            :
   hltLeptonTag(iConfig.getParameter< edm::InputTag > ("HltLeptonTag")),
   sourceJetTag(iConfig.getParameter< edm::InputTag > ("SourceJetTag")),
   minDeltaR_(iConfig.getParameter< double > ("minDeltaR"))
 {
   using namespace edm;
   using namespace std;
 
   typedef vector<RefVector<vector<jetType>,jetType,refhelper::FindUsingAdvance<vector<jetType>,jetType> > > JetCollectionVector;
   typedef vector<jetType> JetCollection;
 
   m_theLeptonToken = consumes<trigger::TriggerFilterObjectWithRefs>(hltLeptonTag);
   m_theJetToken = consumes<std::vector<jetType>>(sourceJetTag);
   produces<JetCollectionVector>();
   produces<JetCollection>();
 }

template<typename jetType >

HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::~HLTJetCollectionsForBoostedLeptonPlusJets ( )

Definition at line 43 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

 {
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

template<typename jetType >

void HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 52 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

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

                                                                                                                {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag> ("HltLeptonTag", edm::InputTag("triggerFilterObjectWithRefs"));
   //(2)
   desc.add<edm::InputTag> ("SourceJetTag", edm::InputTag("jetCollection"));
   //(2)
   desc.add<double> ("minDeltaR", 0.5);
   descriptions.add(defaultModuleLabel<HLTJetCollectionsForBoostedLeptonPlusJets<jetType>>(), desc);
 }

template<typename jetType >

void HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::stream::EDProducerBase.

Definition at line 71 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

References deltaR(), spr::find(), edm::Event::getByToken(), metsig::jet, fwrapper::jets, HLT_25ns14e33_v1_cff::jetType, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theJetToken, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theLeptonToken, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::minDeltaR_, edm::Event::put(), dt_dqm_sourceclient_common_cff::reco, python.multivaluedict::sort(), trigger::TriggerCluster, trigger::TriggerElectron, trigger::TriggerMuon, and trigger::TriggerPhoton.

 {
   using namespace edm;
   using namespace std;
   //(3)
   using namespace reco;
   //(3)
   typedef vector<RefVector<vector<jetType>,jetType,refhelper::FindUsingAdvance<vector<jetType>,jetType> > > JetCollectionVector;
   typedef vector<jetType> JetCollection;
   typedef edm::Ref<JetCollection> JetRef;
   typedef edm::RefVector<JetCollection> JetRefVector;
 
   Handle<trigger::TriggerFilterObjectWithRefs> PrevFilterOutput;
   iEvent.getByToken(m_theLeptonToken,PrevFilterOutput);
  
   //its easier on the if statement flow if I try everything at once, shouldnt add to timing
   // Electrons can be stored as objects of types TriggerCluster, TriggerElectron, or TriggerPhoton
   vector<reco::RecoChargedCandidateRef> muonCands;
   PrevFilterOutput->getObjects(trigger::TriggerMuon,muonCands);
 
   vector<Ref<reco::ElectronCollection> > eleCands;
   PrevFilterOutput->getObjects(trigger::TriggerElectron,eleCands);
   
   trigger::VRphoton photonCands;
   PrevFilterOutput->getObjects(trigger::TriggerPhoton, photonCands);
   
   vector<Ref<reco::RecoEcalCandidateCollection> > clusCands;
   PrevFilterOutput->getObjects(trigger::TriggerCluster,clusCands);
 
   Handle<JetCollection> theJetCollectionHandle;
   iEvent.getByToken(m_theJetToken, theJetCollectionHandle);
   
   typename JetCollection::const_iterator jet;
   
   auto_ptr<JetCollection> allSelections(new JetCollection);
   auto_ptr<JetCollectionVector> product(new JetCollectionVector);
 
   std::vector<size_t> usedCands;
 
   if(!muonCands.empty()){ // muons
     for (jet = theJetCollectionHandle->begin(); jet != theJetCollectionHandle->end(); jet++) {
       //const jetType* referenceJet = &*jet;
       jetType cleanedJet = *jet; //copy original jet
       for(size_t candNr=0;candNr<muonCands.size();candNr++){
         if (std::find(usedCands.begin(),usedCands.end(),candNr)!=usedCands.end()) continue;
         if (deltaR((*muonCands[candNr]),cleanedJet) <= minDeltaR_) {
           std::vector<edm::Ptr<reco::PFCandidate> > pfConstituents = cleanedJet.getPFConstituents();
           for(std::vector<edm::Ptr<reco::PFCandidate> >::const_iterator i_candidate = pfConstituents.begin(); i_candidate != pfConstituents.end(); ++i_candidate){
             if (deltaR((*muonCands[candNr]),(**i_candidate))<0.001) {
           cleanedJet.setP4( cleanedJet.p4() - muonCands[candNr]->p4());
               usedCands.push_back(candNr);
               break;
             }//if constituent matched
           }//for constituents
         }//if dR<min
       }//for cands
       allSelections->push_back(cleanedJet);
     }//for jets
   }//if cands
 
   if(!eleCands.empty()){ // electrons
     for (jet = theJetCollectionHandle->begin(); jet != theJetCollectionHandle->end(); jet++) {
       //const jetType* referenceJet = &*jet;
       jetType cleanedJet = *jet; //copy original jet
       for(size_t candNr=0;candNr<eleCands.size();candNr++){
         if (std::find(usedCands.begin(),usedCands.end(),candNr)!=usedCands.end()) continue;
         if (deltaR((*eleCands[candNr]),cleanedJet) <= minDeltaR_) {
           std::vector<edm::Ptr<reco::PFCandidate> > pfConstituents = cleanedJet.getPFConstituents();
           for(std::vector<edm::Ptr<reco::PFCandidate> >::const_iterator i_candidate = pfConstituents.begin(); i_candidate != pfConstituents.end(); ++i_candidate){
             if (deltaR((*eleCands[candNr]),(**i_candidate))<0.001) {
           cleanedJet.setP4( cleanedJet.p4() - eleCands[candNr]->p4());
               usedCands.push_back(candNr);
               break;
             }//if constituent matched
           }//for constituents
         }//if dR<min
       }//for cands
       allSelections->push_back(cleanedJet);
     }//for jets
   }//if cands
 
   if(!photonCands.empty()){ // photons
     for (jet = theJetCollectionHandle->begin(); jet != theJetCollectionHandle->end(); jet++) {
       //const jetType* referenceJet = &*jet;
       jetType cleanedJet = *jet; //copy original jet
       for(size_t candNr=0;candNr<photonCands.size();candNr++){
         if (std::find(usedCands.begin(),usedCands.end(),candNr)!=usedCands.end()) continue;
         if (deltaR((*photonCands[candNr]),cleanedJet) <= minDeltaR_) {
           std::vector<edm::Ptr<reco::PFCandidate> > pfConstituents = cleanedJet.getPFConstituents();
           for(std::vector<edm::Ptr<reco::PFCandidate> >::const_iterator i_candidate = pfConstituents.begin(); i_candidate != pfConstituents.end(); ++i_candidate){
             if (deltaR((*photonCands[candNr]),(**i_candidate))<0.001) {
           cleanedJet.setP4( cleanedJet.p4() - photonCands[candNr]->p4());
               usedCands.push_back(candNr);
               break;
             }//if constituent matched
           }//for constituents
         }//if dR<min
       }//for cands
       allSelections->push_back(cleanedJet);
     }//for jets
   }//if cands
 
   if(!clusCands.empty()){ // trigger clusters
     for (jet = theJetCollectionHandle->begin(); jet != theJetCollectionHandle->end(); jet++) {
       //const jetType* referenceJet = &*jet;
       jetType cleanedJet = *jet; //copy original jet
       for(size_t candNr=0;candNr<clusCands.size();candNr++){
         if (std::find(usedCands.begin(),usedCands.end(),candNr)!=usedCands.end()) continue;
         if (deltaR((*clusCands[candNr]),cleanedJet) <= minDeltaR_) {
           std::vector<edm::Ptr<reco::PFCandidate> > pfConstituents = cleanedJet.getPFConstituents();
           for(std::vector<edm::Ptr<reco::PFCandidate> >::const_iterator i_candidate = pfConstituents.begin(); i_candidate != pfConstituents.end(); ++i_candidate){
             if (deltaR((*clusCands[candNr]),(**i_candidate))<0.001) {
           cleanedJet.setP4( cleanedJet.p4() - clusCands[candNr]->p4());
               usedCands.push_back(candNr);
               break;
             }//if constituent matched
           }//for constituents
         }//if dR<min
       }//for cands
       allSelections->push_back(cleanedJet);
     }//for jets
   }//if cands
 
   NumericSafeGreaterByPt<jetType> compJets;  
   // reorder cleaned jets
   std::sort (allSelections->begin(), allSelections->end(), compJets);
   edm::OrphanHandle<JetCollection> cleanedJetHandle = iEvent.put(allSelections);
 
   JetCollection const & jets = *cleanedJetHandle;
 
   JetRefVector cleanedJetRefs;
 
   for (unsigned iJet = 0; iJet < jets.size(); ++iJet) {
     cleanedJetRefs.push_back(JetRef(cleanedJetHandle, iJet));
   }
 
   product->emplace_back(cleanedJetRefs);
   iEvent.put(product);
 
   return;
   
 }