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	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () 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 28 of file HLTPFJetIDProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 18 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 39 of file HLTPFJetIDProducer.cc.

References edm::ConfigurationDescriptions::add(), and edm::ParameterSetDescription::add().

                                                                                     {
    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 55 of file HLTPFJetIDProducer.cc.

References funct::abs(), CHF_, PVValHelper::eta, edm::Event::getByToken(), m_thePFJetToken, maxEta_, minPt_, eostools::move(), NCH_, NEF_, NTOT_, EnergyCorrector::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 40 of file HLTPFJetIDProducer.h.

double HLTPFJetIDProducer::CHF_

private

charged hadron fraction

Definition at line 38 of file HLTPFJetIDProducer.h.

Referenced by produce().

edm::InputTag HLTPFJetIDProducer::inputTag_

private

input PFJet collection

Definition at line 45 of file HLTPFJetIDProducer.h.

Referenced by HLTPFJetIDProducer().

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

private

Definition at line 47 of file HLTPFJetIDProducer.h.

Referenced by HLTPFJetIDProducer(), and produce().

double HLTPFJetIDProducer::maxCF_

private

total charged energy fraction

Definition at line 42 of file HLTPFJetIDProducer.h.

double HLTPFJetIDProducer::maxEta_

private

Definition at line 37 of file HLTPFJetIDProducer.h.

Referenced by produce().

double HLTPFJetIDProducer::minPt_

private

Definition at line 36 of file HLTPFJetIDProducer.h.

Referenced by produce().

int HLTPFJetIDProducer::NCH_

private

number of charged constituents

Definition at line 43 of file HLTPFJetIDProducer.h.

Referenced by produce().

double HLTPFJetIDProducer::NEF_

private

neutral EM fraction

Definition at line 41 of file HLTPFJetIDProducer.h.

Referenced by produce().

double HLTPFJetIDProducer::NHF_

private

neutral hadron fraction

Definition at line 39 of file HLTPFJetIDProducer.h.

int HLTPFJetIDProducer::NTOT_

private

number of constituents

Definition at line 44 of file HLTPFJetIDProducer.h.

Referenced by produce().