#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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

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::JVFJetIdProducer ( const edm::ParameterSet & cfg )

Definition at line 21 of file JVFJetIdProducer.cc.

                                                              {
   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");
 }

References looper::cfg, dZcut_, Exception, JVFcut_, kNeutralJetNoPU, kNeutralJetPU, minTrackPt_, neutralJetOption_, srcHardScatterVertex_, srcJets_, srcPFCandidates_, srcPFCandToVertexAssociations_, AlCaHLTBitMon_QueryRunRegistry::string, and verbosity_.

◆ ~JVFJetIdProducer()

JVFJetIdProducer::~JVFJetIdProducer ( )

override

Definition at line 48 of file JVFJetIdProducer.cc.

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

Member Function Documentation

◆ produce()

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

overrideprivate

Definition at line 120 of file JVFJetIdProducer.cc.

                                                                      {
   // 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");
 }