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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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

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_, m_thePFJetToken, and ~HLTPFJetIDProducer().

    : 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

Referenced by HLTPFJetIDProducer().

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(), and HLT_2018_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_, 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().