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

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 29 of file HLTJetCollectionsForBoostedLeptonPlusJets.h.

Constructor & Destructor Documentation

◆ HLTJetCollectionsForBoostedLeptonPlusJets()

template<typename jetType >

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

explicit

Definition at line 23 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

References HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::hltLeptonTag, HLT_2024v14_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")),
       // minimum delta-R^2 threshold with sign
       minDeltaR2_(iConfig.getParameter<double>("minDeltaR") * std::abs(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>();
 }

◆ ~HLTJetCollectionsForBoostedLeptonPlusJets()

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

◆ fillDescriptions()

template<typename jetType >

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

static

Definition at line 49 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

References edm::ConfigurationDescriptions::add(), defaultModuleLabel(), submitPVResolutionJobs::desc, and ProducerED_cfi::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);
 }

◆ produce()

template<typename jetType >

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

overrideprivate

Definition at line 66 of file HLTJetCollectionsForBoostedLeptonPlusJets.cc.

References reco::deltaR2(), MillePedeFileConverter_cfg::e, spr::find(), trigger::TriggerRefsCollections::getObjects(), iEvent, metsig::jet, PDWG_EXODelayedJetMET_cff::jets, HLT_2024v14_cff::jetType, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theJetToken, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theLeptonToken, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::minDeltaR2_, eostools::move(), jetUpdater_cfi::sort, trigger::TriggerCluster, trigger::TriggerElectron, trigger::TriggerMuon, trigger::TriggerPhoton, and trackerHitRTTI::vector.

                                                                                                               {
   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 (reco::deltaR2((*muonCands[candNr]), cleanedJet) <= minDeltaR2_) {
           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 (reco::deltaR2((*muonCands[candNr]), (**i_candidate)) < 1e-6) {
               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 (reco::deltaR2((*eleCands[candNr]), cleanedJet) <= minDeltaR2_) {
           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 (reco::deltaR2((*eleCands[candNr]), (**i_candidate)) < 1e-6) {
               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 (reco::deltaR2((*photonCands[candNr]), cleanedJet) <= minDeltaR2_) {
           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 (reco::deltaR2((*photonCands[candNr]), (**i_candidate)) < 1e-6) {
               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 (reco::deltaR2((*clusCands[candNr]), cleanedJet) <= minDeltaR2_) {
           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 (reco::deltaR2((*clusCands[candNr]), (**i_candidate)) < 1e-6) {
               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;
   cleanedJetRefs.reserve(jets.size());
   for (unsigned iJet = 0; iJet < jets.size(); ++iJet) {
     cleanedJetRefs.push_back(JetRef(cleanedJetHandle, iJet));
   }
 
   product->emplace_back(cleanedJetRefs);
   iEvent.put(std::move(product));
 
   return;
 }