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

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () 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 22 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

References HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::hltLeptonTag, HLT_2018_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 ( )

override

Definition at line 41 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 47 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), defaultModuleLabel(), and HLT_2018_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
	)

overrideprivate

Definition at line 64 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

References PbPb_ZMuSkimMuonDPG_cff::deltaR, spr::find(), edm::Event::getByToken(), trigger::TriggerRefsCollections::getObjects(), metsig::jet, singleTopDQM_cfi::jets, HLT_2018_cff::jetType, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theJetToken, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theLeptonToken, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::minDeltaR_, eostools::move(), edm::Event::put(), 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;
 }