#include <HLTJetCollectionsForBoostedLeptonPlusJets.h>

Inheritance diagram for HLTJetCollectionsForBoostedLeptonPlusJets< jetType >:

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

	~HLTJetCollectionsForBoostedLeptonPlusJets () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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

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

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

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, objects.autophobj::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 ( )

override

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(), and defaultModuleLabel().

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

overrideprivate

Definition at line 71 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

References boostedElectronIsolation_cff::deltaR, spr::find(), edm::Event::getByToken(), trigger::TriggerRefsCollections::getObjects(), metsig::jet, fwrapper::jets, objects.autophobj::jetType, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theJetToken, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theLeptonToken, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::minDeltaR_, eostools::move(), edm::Event::put(), jetUpdater_cfi::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;
   
   unique_ptr<JetCollection> allSelections(new JetCollection);
   unique_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(std::move(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(std::move(product));
 
   return;
   
 }