#include <PFJetFilter.h>

Inheritance diagram for PFJetFilter:

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

	~PFJetFilter () override

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDFilter () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	beginJob () override

void	endJob () override

bool	filter (edm::Event &, const edm::EventSetup &) override

double	resolution (double, bool)

double	response (double, bool)

Private Attributes
PFBenchmarkAlgo *	algo_

double	deltaEt_max

double	deltaEt_min

double	deltaR_max

double	deltaR_min

unsigned int	entry

double	eta_max

double	eta_min

double	genPt_cut

edm::EDGetTokenT< edm::View< reco::Candidate > >	inputRecoLabel_

edm::EDGetTokenT< edm::View< reco::Candidate > >	inputTruthLabel_

double	recPt_cut

bool	verbose

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 13 of file PFJetFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 16 of file PFJetFilter.cc.

References mps_splice::entry, MuonTCMETValueMapProducer_cff::eta_max, MuonME0RecHits_cfi::eta_min, and edm::ParameterSet::getParameter().

                                                        {
   inputTruthLabel_ = consumes<edm::View<reco::Candidate>>(iConfig.getParameter<edm::InputTag>("InputTruthLabel"));
   inputRecoLabel_ = consumes<edm::View<reco::Candidate>>(iConfig.getParameter<edm::InputTag>("InputRecoLabel"));
 
   recPt_cut = iConfig.getParameter<double>("recPt");
   genPt_cut = iConfig.getParameter<double>("genPt");
 
   eta_min = iConfig.getParameter<double>("minEta");
   eta_max = iConfig.getParameter<double>("maxEta");
 
   deltaR_min = iConfig.getParameter<double>("deltaRMin");
   deltaR_max = iConfig.getParameter<double>("deltaRMax");
 
   deltaEt_min = iConfig.getParameter<double>("minDeltaEt");
   deltaEt_max = iConfig.getParameter<double>("maxDeltaEt");
 
   verbose = iConfig.getParameter<bool>("verbose");
 
   entry = 0;
 }

PFJetFilter::~PFJetFilter ( )

override

Definition at line 37 of file PFJetFilter.cc.

37 {}

Member Function Documentation

void PFJetFilter::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 39 of file PFJetFilter.cc.

39 {}

void PFJetFilter::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 41 of file PFJetFilter.cc.

41 {}

bool PFJetFilter::filter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iESetup
	)

overrideprivate

Definition at line 43 of file PFJetFilter.cc.

References funct::abs(), gather_cfg::cout, reco::deltaR(), promptTrackCountingComputer_cfi::deltaRmin, mps_splice::entry, reco::Candidate::eta(), MuonTCMETValueMapProducer_cff::eta_max, MuonME0RecHits_cfi::eta_min, edm::Event::getByToken(), mps_fire::i, trackingPlots::other, reco::Candidate::phi(), edm::Handle< T >::product(), EnergyCorrector::pt, reco::Candidate::pt(), ptmin, and fftjetproducer_cfi::resolution.

                                                                        {
   // Typedefs to use views
   typedef edm::View<reco::Candidate> candidateCollection;
   typedef edm::View<reco::Candidate> candidateCollection;
 
   const candidateCollection *truth_candidates;
   const candidateCollection *reco_candidates;
 
   // ==========================================================
   // Retrieve!
   // ==========================================================
 
   // Get Truth Candidates (GenCandidates, GenJets, etc.)
   Handle<candidateCollection> truth_hnd;
   bool isGen = iEvent.getByToken(inputTruthLabel_, truth_hnd);
   if (!isGen) {
     std::cout << "Warning : no Gen jets in input !" << std::endl;
     return false;
   }
 
   truth_candidates = truth_hnd.product();
 
   // Get Reco Candidates (PFlow, CaloJet, etc.)
   Handle<candidateCollection> reco_hnd;
   bool isReco = iEvent.getByToken(inputRecoLabel_, reco_hnd);
   if (!isReco) {
     std::cout << "Warning : no Reco jets in input !" << std::endl;
     return false;
   }
   reco_candidates = reco_hnd.product();
   if (!truth_candidates || !reco_candidates)
     return false;
 
   bool pass = false;
 
   for (unsigned int i = 0; i < reco_candidates->size(); i++) {
     const reco::Candidate *particle = &(*reco_candidates)[i];
 
     // This protection is meant at not being used !
     assert(particle != nullptr);
 
     double rec_pt = particle->pt();
     double rec_eta = particle->eta();
     double rec_phi = particle->phi();
 
     // skip PFjets with pt < recPt_cut GeV
     if (rec_pt < recPt_cut)
       continue;
 
     // skip PFjets with eta > maxEta_cut or eta < minEta_cut
     if (fabs(rec_eta) > eta_max)
       continue;
     if (fabs(rec_eta) < eta_min)
       continue;
 
     bool Barrel = false;
     bool Endcap = false;
     if (std::abs(rec_eta) < 1.4)
       Barrel = true;
     if (std::abs(rec_eta) > 1.4 && std::abs(rec_eta) < 2.6)
       Endcap = true;
 
     // Keep only jets in the barrel or the endcaps, within the tracker
     // acceptance
     if (!Barrel && !Endcap)
       continue;
 
     // Find the closets recJet
     double deltaRmin = 999.;
     double ptmin = 0.;
     for (unsigned int j = 0; j < reco_candidates->size(); j++) {
       if (i == j)
         continue;
       const reco::Candidate *other = &(*reco_candidates)[j];
       double deltaR = algo_->deltaR(particle, other);
       if (deltaR < deltaRmin && other->pt() > 0.25 * particle->pt() && other->pt() > recPt_cut) {
         deltaRmin = deltaR;
         ptmin = other->pt();
       }
       if (deltaRmin < deltaR_min)
         break;
     }
     if (deltaRmin < deltaR_min)
       continue;
 
     // Find the closest genJet.
     const reco::Candidate *gen_particle = algo_->matchByDeltaR(particle, truth_candidates);
 
     // Check there is a genJet associated to the recoJet
     if (gen_particle == nullptr)
       continue;
 
     // check deltaR is small enough
     double deltaR = algo_->deltaR(particle, gen_particle);
     if (deltaR > deltaR_max)
       continue;
 
     // double true_E = gen_particle->p();
     double true_pt = gen_particle->pt();
     double true_eta = gen_particle->eta();
     double true_phi = gen_particle->phi();
 
     // skip PFjets with pt < genPt_cut GeV
     if (true_pt < genPt_cut)
       continue;
 
     // Find the pT residual
     double resPt = (rec_pt - true_pt) / true_pt;
     double sigma = resolution(true_pt, Barrel);
     double avera = response(true_pt, Barrel);
     double nSig = (resPt - avera) / sigma;
 
     if (nSig > deltaEt_max || nSig < deltaEt_min) {
       /* */
       if (verbose)
         std::cout << "Entry " << entry << " resPt = " << resPt << " sigma/avera/nSig = " << sigma << "/" << avera << "/"
                   << nSig << " pT (T/R) " << true_pt << "/" << rec_pt << " Eta (T/R) " << true_eta << "/" << rec_eta
                   << " Phi (T/R) " << true_phi << "/" << rec_phi << " deltaRMin/ptmin " << deltaRmin << "/" << ptmin
                   << std::endl;
       /* */
       pass = true;
     }
 
     if (pass)
       break;
   }
 
   entry++;
   return pass;
 }

double PFJetFilter::resolution	(	double	pt,
		bool	barrel
	)

private

Definition at line 174 of file PFJetFilter.cc.

References p1, p2, EnergyCorrector::pt, and mathSSE::sqrt().

                                                      {
   double p0 = barrel ? 1.19200e-02 : 8.45341e-03;
   double p1 = barrel ? 1.06138e+00 : 7.96855e-01;
   double p2 = barrel ? -2.05929e+00 : -3.12076e-01;
 
   double resp = p0 + p1 / sqrt(pt) + p2 / pt;
   return resp;
 }

double PFJetFilter::response	(	double	pt,
		bool	barrel
	)

private

Definition at line 183 of file PFJetFilter.cc.

References JetChargeProducer_cfi::exp, p1, and p2.

                                                    {
   double p0 = barrel ? 1.09906E-1 : 6.91939E+1;
   double p1 = barrel ? -1.61443E-1 : -6.92733E+1;
   double p2 = barrel ? 3.45489E+3 : 1.58207E+6;
 
   double resp = p0 + p1 * std::exp(-pt / p2);
   return resp;
 }

Member Data Documentation

PFBenchmarkAlgo* PFJetFilter::algo_

private

Definition at line 40 of file PFJetFilter.h.

double PFJetFilter::deltaEt_max

private

Definition at line 29 of file PFJetFilter.h.

double PFJetFilter::deltaEt_min

private

Definition at line 28 of file PFJetFilter.h.

double PFJetFilter::deltaR_max

private

Definition at line 31 of file PFJetFilter.h.

double PFJetFilter::deltaR_min

private

Definition at line 30 of file PFJetFilter.h.

unsigned int PFJetFilter::entry

private

Definition at line 37 of file PFJetFilter.h.

double PFJetFilter::eta_max

private

Definition at line 33 of file PFJetFilter.h.

double PFJetFilter::eta_min

private

Definition at line 32 of file PFJetFilter.h.

double PFJetFilter::genPt_cut

private

Definition at line 27 of file PFJetFilter.h.

edm::EDGetTokenT<edm::View<reco::Candidate> > PFJetFilter::inputRecoLabel_

private

Definition at line 35 of file PFJetFilter.h.

edm::EDGetTokenT<edm::View<reco::Candidate> > PFJetFilter::inputTruthLabel_

private

Definition at line 34 of file PFJetFilter.h.

double PFJetFilter::recPt_cut

private

Definition at line 26 of file PFJetFilter.h.

bool PFJetFilter::verbose

private

Definition at line 38 of file PFJetFilter.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation