HLTJetCollectionsForBoostedLeptonPlusJets< jetType > Class Template Reference

#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	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () 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

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

References PbPb_ZMuSkimMuonDPG_cff::deltaR, spr::find(), trigger::TriggerRefsCollections::getObjects(), iEvent, metsig::jet, singleTopDQM_cfi::jets, HLT_FULL_cff::jetType, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theJetToken, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theLeptonToken, HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::minDeltaR_, eostools::move(), trigger::TriggerCluster, trigger::TriggerElectron, trigger::TriggerMuon, trigger::TriggerPhoton, and trackerHitRTTI::vector.

Member Data Documentation

hltLeptonTag

template<typename jetType >

edm::InputTag HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::hltLeptonTag

private

Definition at line 51 of file HLTJetCollectionsForBoostedLeptonPlusJets.h.

Referenced by HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::HLTJetCollectionsForBoostedLeptonPlusJets().

m_theJetToken

template<typename jetType >

edm::EDGetTokenT<std::vector<jetType> > HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theJetToken

private

Definition at line 50 of file HLTJetCollectionsForBoostedLeptonPlusJets.h.

Referenced by HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::HLTJetCollectionsForBoostedLeptonPlusJets(), and HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::produce().

m_theLeptonToken

template<typename jetType >

edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::m_theLeptonToken

private

Definition at line 49 of file HLTJetCollectionsForBoostedLeptonPlusJets.h.

Referenced by HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::HLTJetCollectionsForBoostedLeptonPlusJets(), and HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::produce().

minDeltaR_

template<typename jetType >

double HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::minDeltaR_

private

Definition at line 54 of file HLTJetCollectionsForBoostedLeptonPlusJets.h.

Referenced by HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::produce().

sourceJetTag

template<typename jetType >

edm::InputTag HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::sourceJetTag

private

Definition at line 52 of file HLTJetCollectionsForBoostedLeptonPlusJets.h.

Referenced by HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::HLTJetCollectionsForBoostedLeptonPlusJets().