#include <JVFJetIdProducer.h>

Inheritance diagram for JVFJetIdProducer:

Public Member Functions
	JVFJetIdProducer (const edm::ParameterSet &)

	~JVFJetIdProducer () 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

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
double	dZcut_

double	JVFcut_

double	minTrackPt_

int	neutralJetOption_

edm::EDGetTokenT < reco::VertexCollection >	srcHardScatterVertex_

edm::EDGetTokenT < reco::PFJetCollection >	srcJets_

edm::EDGetTokenT < reco::PFCandidateCollection >	srcPFCandidates_

edm::EDGetTokenT < PFCandToVertexAssMap >	srcPFCandToVertexAssociations_

int	verbosity_

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

Discriminate jets originating from the hard-scatter event from pile-up jets, based on the fraction of tracks within the jet that are associated to the hard-scatter vertex. Jets outside the tracking acceptance are considered to originate from the hard-scatter event per default. Optionally, they can be classified as pile-up.

Authors: Christian Veelken, LLR

Definition at line 25 of file JVFJetIdProducer.h.

Constructor & Destructor Documentation

JVFJetIdProducer::JVFJetIdProducer ( const edm::ParameterSet & cfg )

Definition at line 21 of file JVFJetIdProducer.cc.

References dZcut_, Exception, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), JVFcut_, kNeutralJetNoPU, kNeutralJetPU, minTrackPt_, neutralJetOption_, srcHardScatterVertex_, srcJets_, srcPFCandidates_, srcPFCandToVertexAssociations_, AlCaHLTBitMon_QueryRunRegistry::string, and verbosity_.

                                                              {
   srcJets_ = consumes<reco::PFJetCollection>(cfg.getParameter<edm::InputTag>("srcJets"));
 
   srcPFCandidates_ = consumes<reco::PFCandidateCollection>(cfg.getParameter<edm::InputTag>("srcPFCandidates"));
   srcPFCandToVertexAssociations_ =
       consumes<PFCandToVertexAssMap>(cfg.getParameter<edm::InputTag>("srcPFCandToVertexAssociations"));
   srcHardScatterVertex_ = consumes<reco::VertexCollection>(cfg.getParameter<edm::InputTag>("srcHardScatterVertex"));
   minTrackPt_ = cfg.getParameter<double>("minTrackPt");
   dZcut_ = cfg.getParameter<double>("dZcut");
 
   JVFcut_ = cfg.getParameter<double>("JVFcut");
 
   std::string neutralJetOption_string = cfg.getParameter<std::string>("neutralJetOption");
   if (neutralJetOption_string == "PU")
     neutralJetOption_ = kNeutralJetPU;
   else if (neutralJetOption_string == "noPU")
     neutralJetOption_ = kNeutralJetNoPU;
   else
     throw cms::Exception("JVFJetIdProducer")
         << "Invalid Configuration Parameter 'neutralJetOption' = " << neutralJetOption_string << " !!\n";
 
   verbosity_ = (cfg.exists("verbosity")) ? cfg.getParameter<int>("verbosity") : 0;
 
   produces<edm::ValueMap<double>>("Discriminant");
   produces<edm::ValueMap<int>>("Id");
 }

JVFJetIdProducer::~JVFJetIdProducer ( )

override

Definition at line 48 of file JVFJetIdProducer.cc.

                                     {
   // nothing to be done yet...
 }

Member Function Documentation

void JVFJetIdProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

overrideprivate

Definition at line 120 of file JVFJetIdProducer.cc.

References dZcut_, edm::Event::getByToken(), metsig::jet, fwrapper::jets, JVFcut_, kNeutralJetNoPU, minTrackPt_, eostools::move(), neutralJetOption_, slimmedMuons_cfi::pfCandidates, edm::Event::put(), srcHardScatterVertex_, srcJets_, srcPFCandidates_, srcPFCandToVertexAssociations_, and verbosity_.

                                                                      {
   // get jets
   edm::Handle<reco::PFJetCollection> jets;
   evt.getByToken(srcJets_, jets);
 
   // get PFCandidates
   edm::Handle<reco::PFCandidateCollection> pfCandidates;
   evt.getByToken(srcPFCandidates_, pfCandidates);
 
   // get PFCandidate-to-vertex associations and "the" hard-scatter vertex
   edm::Handle<PFCandToVertexAssMap> pfCandToVertexAssociations;
   evt.getByToken(srcPFCandToVertexAssociations_, pfCandToVertexAssociations);
 
   edm::Handle<reco::VertexCollection> hardScatterVertex;
   evt.getByToken(srcHardScatterVertex_, hardScatterVertex);
 
   std::vector<double> jetIdDiscriminants;
   std::vector<int> jetIdFlags;
 
   size_t numJets = jets->size();
   for (size_t iJet = 0; iJet < numJets; ++iJet) {
     reco::PFJetRef jet(jets, iJet);
 
     double jetJVF = computeJVF(
         *jet, *pfCandToVertexAssociations, *hardScatterVertex, dZcut_, minTrackPt_, verbosity_ && jet->pt() > 20.);
     jetIdDiscriminants.push_back(jetJVF);
 
     int jetIdFlag = 0;
     if (jetJVF > JVFcut_)
       jetIdFlag = 255;
     else if (jetJVF < -0.5 && neutralJetOption_ == kNeutralJetNoPU)
       jetIdFlag = 255;
     jetIdFlags.push_back(jetIdFlag);
   }
 
   auto jetIdDiscriminants_ptr = std::make_unique<edm::ValueMap<double>>();
   edm::ValueMap<double>::Filler jetIdDiscriminantFiller(*jetIdDiscriminants_ptr);
   jetIdDiscriminantFiller.insert(jets, jetIdDiscriminants.begin(), jetIdDiscriminants.end());
   jetIdDiscriminantFiller.fill();
 
   auto jetIdFlags_ptr = std::make_unique<edm::ValueMap<int>>();
   edm::ValueMap<int>::Filler jetIdFlagFiller(*jetIdFlags_ptr);
   jetIdFlagFiller.insert(jets, jetIdFlags.begin(), jetIdFlags.end());
   jetIdFlagFiller.fill();
 
   evt.put(std::move(jetIdDiscriminants_ptr), "Discriminant");
   evt.put(std::move(jetIdFlags_ptr), "Id");
 }