#include <HLTJetCollForElePlusJets.h>

Inheritance diagram for HLTJetCollForElePlusJets< T >:

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

	~HLTJetCollForElePlusJets () 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	hltElectronTag

edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs >	m_theElectronToken

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

double	maxAbsJetEta_

double	minDeltaEta_

double	minDeltaR_

double	minJetPt_

unsigned int	minNJets_

double	minSoftJetPt_

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 T>
class HLTJetCollForElePlusJets< T >

This class is an EDProducer implementing an HLT trigger for electron and jet objects, cutting on variables relating to the jet 4-momentum representation. The producer checks for overlaps between electrons and jets and if a combination of one electron + jets cleaned against this electrons 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 HLTJetCollForElePlusJets.h.

Constructor & Destructor Documentation

◆ HLTJetCollForElePlusJets()

template<typename T >

HLTJetCollForElePlusJets< T >::HLTJetCollForElePlusJets ( const edm::ParameterSet & iConfig )

explicit

Definition at line 15 of file HLTJetCollForElePlusJets.cc.

References HLTJetCollForElePlusJets< T >::hltElectronTag, HLTJetCollForElePlusJets< T >::m_theElectronToken, HLTJetCollForElePlusJets< T >::m_theJetToken, and HLTJetCollForElePlusJets< T >::sourceJetTag.

     : hltElectronTag(iConfig.getParameter<edm::InputTag>("HltElectronTag")),
       sourceJetTag(iConfig.getParameter<edm::InputTag>("SourceJetTag")),
       minJetPt_(iConfig.getParameter<double>("MinJetPt")),
       maxAbsJetEta_(iConfig.getParameter<double>("MaxAbsJetEta")),
       minNJets_(iConfig.getParameter<unsigned int>("MinNJets")),
       minDeltaR_(iConfig.getParameter<double>("minDeltaR")),
       //Only for VBF
       minSoftJetPt_(iConfig.getParameter<double>("MinSoftJetPt")),
       minDeltaEta_(iConfig.getParameter<double>("MinDeltaEta")) {
   typedef std::vector<T> TCollection;
   m_theElectronToken = consumes<trigger::TriggerFilterObjectWithRefs>(hltElectronTag);
   m_theJetToken = consumes<TCollection>(sourceJetTag);
   produces<TCollection>();
 }

◆ ~HLTJetCollForElePlusJets()

template<typename T >

HLTJetCollForElePlusJets< T >::~HLTJetCollForElePlusJets ( )

override

Definition at line 32 of file HLTJetCollForElePlusJets.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 T >

void HLTJetCollForElePlusJets< T >::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 38 of file HLTJetCollForElePlusJets.cc.

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

                                                                                              {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("HltElectronTag", edm::InputTag("triggerFilterObjectWithRefs"));
   desc.add<edm::InputTag>("SourceJetTag", edm::InputTag("jetCollection"));
   desc.add<double>("MinJetPt", 30.);
   desc.add<double>("MaxAbsJetEta", 2.6);
   desc.add<unsigned int>("MinNJets", 1);
   desc.add<double>("minDeltaR", 0.5);
   //Only for VBF
   desc.add<double>("MinSoftJetPt", 25.);
   desc.add<double>("MinDeltaEta", -1.);
   descriptions.add(defaultModuleLabel<HLTJetCollForElePlusJets<T>>(), desc);
 }

◆ produce()

template<typename T >

void HLTJetCollForElePlusJets< T >::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 58 of file HLTJetCollForElePlusJets.cc.

References funct::abs(), trigger::TriggerRefsCollections::getObjects(), iEvent, dqmiolumiharvest::j, dqmdumpme::k, eostools::move(), multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, jetsAK4_CHS_cff::sort, trigger::TriggerCluster, trigger::TriggerElectron, trigger::TriggerPhoton, and tier0::unique().

                                                                                        {
   using namespace edm;
   using namespace std;
 
   typedef vector<T> TCollection;
   typedef Ref<TCollection> TRef;
   typedef edm::RefVector<TCollection> TRefVector;
   typedef std::vector<edm::RefVector<std::vector<T>, T, edm::refhelper::FindUsingAdvance<std::vector<T>, T>>>
       TCollectionVector;
 
   edm::Handle<trigger::TriggerFilterObjectWithRefs> PrevFilterOutput;
   iEvent.getByToken(m_theElectronToken, 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
   std::vector<edm::Ref<reco::RecoEcalCandidateCollection>> clusCands;
   PrevFilterOutput->getObjects(trigger::TriggerCluster, clusCands);
 
   std::vector<edm::Ref<reco::ElectronCollection>> eleCands;
   PrevFilterOutput->getObjects(trigger::TriggerElectron, eleCands);
 
   trigger::VRphoton photonCands;
   PrevFilterOutput->getObjects(trigger::TriggerPhoton, photonCands);
 
   //prepare the collection of 3-D vector for electron momenta
   std::vector<TVector3> ElePs;
 
   if (!clusCands.empty()) {  //try trigger cluster
     for (auto& clusCand : clusCands) {
       TVector3 positionVector(clusCand->superCluster()->position().x(),
                               clusCand->superCluster()->position().y(),
                               clusCand->superCluster()->position().z());
       ElePs.push_back(positionVector);
     }
   } else if (!eleCands.empty()) {  // try trigger electrons
     for (auto& eleCand : eleCands) {
       TVector3 positionVector(eleCand->superCluster()->position().x(),
                               eleCand->superCluster()->position().y(),
                               eleCand->superCluster()->position().z());
       ElePs.push_back(positionVector);
     }
   } else if (!photonCands.empty()) {  // try trigger photons
     for (auto& photonCand : photonCands) {
       TVector3 positionVector(photonCand->superCluster()->position().x(),
                               photonCand->superCluster()->position().y(),
                               photonCand->superCluster()->position().z());
       ElePs.push_back(positionVector);
     }
   }
 
   edm::Handle<TCollection> theJetCollectionHandle;
   iEvent.getByToken(m_theJetToken, theJetCollectionHandle);
 
   const TCollection& theJetCollection = *theJetCollectionHandle;
 
   std::unique_ptr<TCollection> theFilteredJetCollection(new TCollection);
 
   std::unique_ptr<TCollectionVector> allSelections(new TCollectionVector());
 
   bool foundSolution(false);
 
   for (auto& EleP : ElePs) {
     bool VBFJetPair = false;
     std::vector<int> store_jet;
     TRefVector refVector;
 
     for (unsigned int j = 0; j < theJetCollection.size(); j++) {
       TVector3 JetP(theJetCollection[j].px(), theJetCollection[j].py(), theJetCollection[j].pz());
       double DR = EleP.DeltaR(JetP);
 
       if (JetP.Pt() > minJetPt_ && std::abs(JetP.Eta()) < maxAbsJetEta_ && DR > minDeltaR_) {
         store_jet.push_back(j);
         // The VBF part of the filter
         if (minDeltaEta_ > 0) {
           for (unsigned int k = j + 1; k < theJetCollection.size(); k++) {
             TVector3 SoftJetP(theJetCollection[k].px(), theJetCollection[k].py(), theJetCollection[k].pz());
             double softDR = EleP.DeltaR(SoftJetP);
 
             if (SoftJetP.Pt() > minSoftJetPt_ && std::abs(SoftJetP.Eta()) < maxAbsJetEta_ && softDR > minDeltaR_)
               if (std::abs(SoftJetP.Eta() - JetP.Eta()) > minDeltaEta_) {
                 store_jet.push_back(k);
                 VBFJetPair = true;
               }
           }
         }
       }
     }
 
     // Now remove duplicates from the jet collection to store
     std::sort(store_jet.begin(), store_jet.end());
     store_jet.erase(unique(store_jet.begin(), store_jet.end()), store_jet.end());
 
     // Now save the cleaned jets
     for (int& ijet : store_jet) {
       //store all selections
       refVector.push_back(TRef(theJetCollectionHandle, ijet));
       //store first selection which matches the criteria
       if (!foundSolution)
         theFilteredJetCollection->push_back(theJetCollection[ijet]);
     }
     //store all selections
     allSelections->push_back(refVector);
 
     if (theFilteredJetCollection->size() >= minNJets_ && minDeltaEta_ < 0)
       foundSolution = true;
     else if (VBFJetPair && minDeltaEta_ > 0)
       foundSolution = true;
     else if (!foundSolution)
       theFilteredJetCollection->clear();
   }
 
   iEvent.put(std::move(theFilteredJetCollection));
 
   return;
 }

Member Data Documentation

◆ hltElectronTag

template<typename T >

edm::InputTag HLTJetCollForElePlusJets< T >::hltElectronTag

private

Definition at line 51 of file HLTJetCollForElePlusJets.h.

Referenced by HLTJetCollForElePlusJets< T >::HLTJetCollForElePlusJets().

◆ m_theElectronToken

template<typename T >

edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> HLTJetCollForElePlusJets< T >::m_theElectronToken

private

Definition at line 49 of file HLTJetCollForElePlusJets.h.

Referenced by HLTJetCollForElePlusJets< T >::HLTJetCollForElePlusJets().

◆ m_theJetToken

template<typename T >

edm::EDGetTokenT<std::vector<T> > HLTJetCollForElePlusJets< T >::m_theJetToken

private

Definition at line 50 of file HLTJetCollForElePlusJets.h.

Referenced by HLTJetCollForElePlusJets< T >::HLTJetCollForElePlusJets().

◆ maxAbsJetEta_

template<typename T >

double HLTJetCollForElePlusJets< T >::maxAbsJetEta_

private

Definition at line 55 of file HLTJetCollForElePlusJets.h.

◆ minDeltaEta_

template<typename T >

double HLTJetCollForElePlusJets< T >::minDeltaEta_

private

Definition at line 61 of file HLTJetCollForElePlusJets.h.

◆ minDeltaR_

template<typename T >

double HLTJetCollForElePlusJets< T >::minDeltaR_

private

Definition at line 58 of file HLTJetCollForElePlusJets.h.

◆ minJetPt_

template<typename T >

double HLTJetCollForElePlusJets< T >::minJetPt_

private

Definition at line 54 of file HLTJetCollForElePlusJets.h.

◆ minNJets_

template<typename T >

unsigned int HLTJetCollForElePlusJets< T >::minNJets_

private

Definition at line 56 of file HLTJetCollForElePlusJets.h.

◆ minSoftJetPt_

template<typename T >

double HLTJetCollForElePlusJets< T >::minSoftJetPt_

private

Definition at line 60 of file HLTJetCollForElePlusJets.h.

◆ sourceJetTag

template<typename T >

edm::InputTag HLTJetCollForElePlusJets< T >::sourceJetTag

private

Definition at line 52 of file HLTJetCollForElePlusJets.h.

Referenced by HLTJetCollForElePlusJets< T >::HLTJetCollForElePlusJets().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

template<typename T> class HLTJetCollForElePlusJets< T >

Constructor & Destructor Documentation

◆ HLTJetCollForElePlusJets()

◆ ~HLTJetCollForElePlusJets()

Member Function Documentation

◆ fillDescriptions()

◆ produce()

Member Data Documentation

◆ hltElectronTag

◆ m_theElectronToken

◆ m_theJetToken

◆ maxAbsJetEta_

◆ minDeltaEta_

◆ minDeltaR_

◆ minJetPt_

◆ minNJets_

◆ minSoftJetPt_

◆ sourceJetTag

template<typename T>
class HLTJetCollForElePlusJets< T >