HLTPFJetIDProducer Class Reference

This applies PFJet ID and produces a jet collection with jets that pass the ID. More...

#include <HLTPFJetIDProducer.h>

Inheritance diagram for HLTPFJetIDProducer:

Public Member Functions
	HLTPFJetIDProducer (const edm::ParameterSet &iConfig)
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
	~HLTPFJetIDProducer () 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 Attributes
double	CEF_
	charged EM fraction More...
double	CHF_
	charged hadron fraction More...
edm::InputTag	inputTag_
	input PFJet collection More...
edm::EDGetTokenT < reco::PFJetCollection >	m_thePFJetToken
double	maxCF_
	total charged energy fraction More...
double	maxEta_
double	minPt_
int	NCH_
	number of charged constituents More...
double	NEF_
	neutral EM fraction More...
double	NHF_
	neutral hadron fraction More...
int	NTOT_
	number of constituents More...

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

This applies PFJet ID and produces a jet collection with jets that pass the ID.

Author

Michele de Gruttola, Jia Fu Low (Nov 2013)

This receives a PFJet collection, selects jets that pass PFJet ID, and makes an output PFJet collection with only jets that pass.

See header file for documentation

Author

a Jet/MET person

Definition at line 27 of file HLTPFJetIDProducer.h.

Constructor & Destructor Documentation

HLTPFJetIDProducer::HLTPFJetIDProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 17 of file HLTPFJetIDProducer.cc.

References inputTag_, and m_thePFJetToken.

    : minPt_(iConfig.getParameter<double>("minPt")),
      maxEta_(iConfig.getParameter<double>("maxEta")),
      CHF_(iConfig.getParameter<double>("CHF")),
      NHF_(iConfig.getParameter<double>("NHF")),
      CEF_(iConfig.getParameter<double>("CEF")),
      NEF_(iConfig.getParameter<double>("NEF")),
      maxCF_(iConfig.getParameter<double>("maxCF")),
      NCH_(iConfig.getParameter<int>("NCH")),
      NTOT_(iConfig.getParameter<int>("NTOT")),
      inputTag_(iConfig.getParameter<edm::InputTag>("jetsInput")) {
  m_thePFJetToken = consumes<reco::PFJetCollection>(inputTag_);

  // Register the products
  produces<reco::PFJetCollection>();
}

HLTPFJetIDProducer::~HLTPFJetIDProducer ( )

overridedefault

Member Function Documentation

void HLTPFJetIDProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 38 of file HLTPFJetIDProducer.cc.

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

                                                                                    {
  edm::ParameterSetDescription desc;
  desc.add<double>("minPt", 20.);
  desc.add<double>("maxEta", 1e99);
  desc.add<double>("CHF", -99.);
  desc.add<double>("NHF", 99.);
  desc.add<double>("CEF", 99.);
  desc.add<double>("NEF", 99.);
  desc.add<double>("maxCF", 99.);
  desc.add<int>("NCH", 0);
  desc.add<int>("NTOT", 0);
  desc.add<edm::InputTag>("jetsInput", edm::InputTag("hltAntiKT4PFJets"));
  descriptions.add("hltPFJetIDProducer", desc);
}

void HLTPFJetIDProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Definition at line 54 of file HLTPFJetIDProducer.cc.

References funct::abs(), CHF_, PVValHelper::eta, edm::Event::getByToken(), dqmiolumiharvest::j, m_thePFJetToken, maxEta_, minPt_, eostools::move(), NCH_, NEF_, gpuPixelDoublets::ntot, NTOT_, DiDispStaMuonMonitor_cfi::pt, edm::Event::put(), and mps_fire::result.

                                                                              {
  // Create a pointer to the products
  std::unique_ptr<reco::PFJetCollection> result(new reco::PFJetCollection());

  edm::Handle<reco::PFJetCollection> pfjets;
  iEvent.getByToken(m_thePFJetToken, pfjets);

  for (auto const& j : *pfjets) {
    bool pass = false;
    double pt = j.pt();
    double eta = j.eta();
    double abseta = std::abs(eta);

    if (!(pt > 0.))
      continue;  // skip jets with zero or negative pt

    if (pt < minPt_) {
      pass = true;

    } else if (abseta >= maxEta_) {
      pass = true;

    } else {
      double chf = j.chargedHadronEnergyFraction();
      //double nhf  = j->neutralHadronEnergyFraction() + j->HFHadronEnergyFraction();
      double nhf = j.neutralHadronEnergyFraction();
      double cef = j.chargedEmEnergyFraction();
      double nef = j.neutralEmEnergyFraction();
      double cftot = chf + cef + j.chargedMuEnergyFraction();
      int nch = j.chargedMultiplicity();
      int ntot = j.numberOfDaughters();

      pass = true;
      pass = pass && (ntot > NTOT_);
      pass = pass && (nef < NEF_);
      pass = pass && (nhf < NHF_ || abseta >= 2.4);  //NHF-cut does not work in HF anymore with recent PF
      pass = pass && (cef < CEF_ || abseta >= 2.4);
      pass = pass && (chf > CHF_ || abseta >= 2.4);
      pass = pass && (nch > NCH_ || abseta >= 2.4);
      pass = pass && (cftot < maxCF_ || abseta >= 2.4);
    }

    if (pass)
      result->push_back(j);
  }

  // Put the products into the Event
  iEvent.put(std::move(result));
}

Member Data Documentation

double HLTPFJetIDProducer::CEF_

private

charged EM fraction

Definition at line 39 of file HLTPFJetIDProducer.h.

double HLTPFJetIDProducer::CHF_

private

charged hadron fraction

Definition at line 37 of file HLTPFJetIDProducer.h.

Referenced by produce().

edm::InputTag HLTPFJetIDProducer::inputTag_

private

input PFJet collection

Definition at line 44 of file HLTPFJetIDProducer.h.

Referenced by HLTPFJetIDProducer().

edm::EDGetTokenT<reco::PFJetCollection> HLTPFJetIDProducer::m_thePFJetToken

private

Definition at line 46 of file HLTPFJetIDProducer.h.

Referenced by HLTPFJetIDProducer(), and produce().

double HLTPFJetIDProducer::maxCF_

private

total charged energy fraction

Definition at line 41 of file HLTPFJetIDProducer.h.

double HLTPFJetIDProducer::maxEta_

private

Definition at line 36 of file HLTPFJetIDProducer.h.

Referenced by produce().

double HLTPFJetIDProducer::minPt_

private

Definition at line 35 of file HLTPFJetIDProducer.h.

Referenced by produce().

int HLTPFJetIDProducer::NCH_

private

number of charged constituents

Definition at line 42 of file HLTPFJetIDProducer.h.

Referenced by produce().

double HLTPFJetIDProducer::NEF_

private

neutral EM fraction

Definition at line 40 of file HLTPFJetIDProducer.h.

Referenced by produce().

double HLTPFJetIDProducer::NHF_

private

neutral hadron fraction

Definition at line 38 of file HLTPFJetIDProducer.h.

int HLTPFJetIDProducer::NTOT_

private

number of constituents

Definition at line 43 of file HLTPFJetIDProducer.h.

Referenced by produce().