HLTJetCollForElePlusJets< T > Class Template Reference

#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_2022v12_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, jetUpdater_cfi::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().