#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

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

◆ HLTJetCollectionsForBoostedLeptonPlusJets()

template<typename jetType >

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

explicit

Definition at line 22 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

     : 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>();
 }

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

◆ ~HLTJetCollectionsForBoostedLeptonPlusJets()

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

◆ fillDescriptions()

template<typename jetType >

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

static

Definition at line 47 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

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

References edm::ConfigurationDescriptions::add(), defaultModuleLabel(), submitPVResolutionJobs::desc, and HLT_FULL_cff::InputTag.

◆ produce()

template<typename jetType >

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

overrideprivate

Definition at line 64 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

                                                                                                               {
   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;
 }