Checks what a re-reco changes in PFCandidates. More...

Inheritance diagram for PFCandidateChecker:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

	PFCandidateChecker (const edm::ParameterSet &)

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

	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

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

	EDAnalyzerBase ()

	EDAnalyzerBase (const EDAnalyzerBase &)=delete

ModuleDescription const &	moduleDescription () const

const EDAnalyzerBase &	operator= (const EDAnalyzerBase &)=delete

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

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	printElementsInBlocks (const reco::PFCandidate &cand, std::ostream &out=std::cout) const

void	printJets (const reco::PFJetCollection &pfJetsReco, const reco::PFJetCollection &pfJetsReReco) const

void	printMet (const reco::PFCandidateCollection &pfReco, const reco::PFCandidateCollection &pfReReco) const

Static Private Member Functions
static bool	greaterPt (const reco::PFCandidate &a, const reco::PFCandidate &b)

Private Attributes
double	deltaEMax_
	Cuts for comparison. More...

double	deltaEtaMax_

double	deltaPhiMax_

unsigned	entry_
	Counter. More...

edm::InputTag	inputTagPFCandidatesReco_
	PFCandidates in which we'll look for pile up particles. More...

edm::InputTag	inputTagPFCandidatesReReco_

edm::InputTag	inputTagPFJetsReco_

edm::InputTag	inputTagPFJetsReReco_

bool	printBlocks_
	print the blocks associated to a given candidate ? More...

bool	rankByPt_
	rank the candidates by Pt More...

bool	verbose_
	verbose ? More...

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer<>
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

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Checks what a re-reco changes in PFCandidates.

Author: Patrick Janot

Date: August 2011

Definition at line 23 of file PFCandidateChecker.cc.

Constructor & Destructor Documentation

◆ PFCandidateChecker()

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

explicit

Definition at line 68 of file PFCandidateChecker.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and LogDebug.

                                                                      {
   inputTagPFCandidatesReco_ = iConfig.getParameter<InputTag>("pfCandidatesReco");
 
   inputTagPFCandidatesReReco_ = iConfig.getParameter<InputTag>("pfCandidatesReReco");
 
   inputTagPFJetsReco_ = iConfig.getParameter<InputTag>("pfJetsReco");
 
   inputTagPFJetsReReco_ = iConfig.getParameter<InputTag>("pfJetsReReco");
 
   deltaEMax_ = iConfig.getParameter<double>("deltaEMax");
 
   deltaEtaMax_ = iConfig.getParameter<double>("deltaEtaMax");
 
   deltaPhiMax_ = iConfig.getParameter<double>("deltaPhiMax");
 
   verbose_ = iConfig.getUntrackedParameter<bool>("verbose", false);
 
   printBlocks_ = iConfig.getUntrackedParameter<bool>("printBlocks", false);
 
   rankByPt_ = iConfig.getUntrackedParameter<bool>("rankByPt", false);
 
   entry_ = 0;
 
   LogDebug("PFCandidateChecker") << " input collections : " << inputTagPFCandidatesReco_ << " "
                                  << inputTagPFCandidatesReReco_;
 }

Member Function Documentation

◆ analyze()

void PFCandidateChecker::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 95 of file PFCandidateChecker.cc.

References gather_cfg::cout, EgHLTOffHistBins_cfi::deltaE, spr::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), mps_fire::i, iEvent, LogDebug, reco::LeafCandidate::phi(), and jetUpdater_cfi::sort.

                                                                               {
   LogDebug("PFCandidateChecker") << "START event: " << iEvent.id().event() << " in run " << iEvent.id().run() << endl;
 
   // get PFCandidates
 
   Handle<PFCandidateCollection> pfCandidatesReco;
   iEvent.getByLabel(inputTagPFCandidatesReco_, pfCandidatesReco);
 
   Handle<PFCandidateCollection> pfCandidatesReReco;
   iEvent.getByLabel(inputTagPFCandidatesReReco_, pfCandidatesReReco);
 
   Handle<PFJetCollection> pfJetsReco;
   iEvent.getByLabel(inputTagPFJetsReco_, pfJetsReco);
 
   Handle<PFJetCollection> pfJetsReReco;
   iEvent.getByLabel(inputTagPFJetsReReco_, pfJetsReReco);
 
   reco::PFCandidateCollection pfReco, pfReReco;
 
   // to sort, one needs to copy
   if (rankByPt_) {
     pfReco = *pfCandidatesReco;
     pfReReco = *pfCandidatesReReco;
     sort(pfReco.begin(), pfReco.end(), greaterPt);
     sort(pfReReco.begin(), pfReReco.end(), greaterPt);
   }
 
   unsigned minSize = pfReco.size() < pfReReco.size() ? pfReco.size() : pfReReco.size();
   bool differentCand = false;
   bool differentSize = pfReco.size() != pfReReco.size();
   if (differentSize)
     std::cout << "+++WARNING+++ PFCandidate size changed for entry " << entry_ << " !" << endl
               << " - RECO    size : " << pfReco.size() << endl
               << " - Re-RECO size : " << pfReReco.size() << endl;
 
   unsigned npr = 0;
   for (unsigned i = 0; i < minSize; i++) {
     const reco::PFCandidate& candReco = (rankByPt_) ? pfReco[i] : (*pfCandidatesReco)[i];
     const reco::PFCandidate& candReReco = (rankByPt_) ? pfReReco[i] : (*pfCandidatesReReco)[i];
 
     double deltaE = (candReReco.energy() - candReco.energy()) / (candReReco.energy() + candReco.energy());
     double deltaEta = candReReco.eta() - candReco.eta();
     double deltaPhi = candReReco.phi() - candReco.phi();
     if (fabs(deltaE) > deltaEMax_ || fabs(deltaEta) > deltaEtaMax_ || fabs(deltaPhi) > deltaPhiMax_) {
       differentCand = true;
       std::cout << "+++WARNING+++ PFCandidate " << i << " changed  for entry " << entry_ << " ! " << std::endl
                 << " - RECO     : " << candReco << std::endl
                 << " - Re-RECO  : " << candReReco << std::endl
                 << " DeltaE   = : " << deltaE << std::endl
                 << " DeltaEta = : " << deltaEta << std::endl
                 << " DeltaPhi = : " << deltaPhi << std::endl
                 << std::endl;
       if (printBlocks_) {
         std::cout << "Elements in Block for RECO: " << std::endl;
         printElementsInBlocks(candReco);
         std::cout << "Elements in Block for Re-RECO: " << std::endl;
         printElementsInBlocks(candReReco);
       }
       if (++npr == 5)
         break;
     }
   }
 
   if (differentSize || differentCand) {
     printJets(*pfJetsReco, *pfJetsReReco);
     printMet(pfReco, pfReReco);
   }
 
   ++entry_;
   LogDebug("PFCandidateChecker") << "STOP event: " << iEvent.id().event() << " in run " << iEvent.id().run()
                                  << std::endl;
 }

◆ greaterPt()

static bool PFCandidateChecker::greaterPt	(	const reco::PFCandidate &	a,
		const reco::PFCandidate &	b
	)

inlinestaticprivate

Definition at line 59 of file PFCandidateChecker.cc.

References a, and b.

59 { return (a.pt() > b.pt()); }

b

double b

Definition: hdecay.h:120

a

double a

Definition: hdecay.h:121

◆ printElementsInBlocks()

void PFCandidateChecker::printElementsInBlocks	(	const reco::PFCandidate &	cand,
		std::ostream &	out = `std::cout`
	)		const

private

Definition at line 221 of file PFCandidateChecker.cc.

References cms::cuda::assert(), DMR_cfg::cerr, mps_fire::i, edm::Ref< C, T, F >::isNull(), and MillePedeFileConverter_cfg::out.

                                                                                           {
   if (!out)
     return;
 
   PFBlockRef firstRef;
 
   assert(!cand.elementsInBlocks().empty());
   for (unsigned i = 0; i < cand.elementsInBlocks().size(); i++) {
     PFBlockRef blockRef = cand.elementsInBlocks()[i].first;
 
     if (blockRef.isNull()) {
       cerr << "ERROR! no block ref!";
       continue;
     }
 
     if (!i) {
       out << (*blockRef);
       firstRef = blockRef;
     } else if (blockRef != firstRef) {
       cerr << "WARNING! This PFCandidate is not made from a single block" << endl;
     }
 
     out << "\t" << cand.elementsInBlocks()[i].second << endl;
   }
 }

◆ printJets()

void PFCandidateChecker::printJets	(	const reco::PFJetCollection &	pfJetsReco,
		const reco::PFJetCollection &	pfJetsReReco
	)		const

private

Definition at line 188 of file PFCandidateChecker.cc.

References gather_cfg::cout, EgHLTOffHistBins_cfi::deltaE, spr::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, reco::LeafCandidate::et(), reco::LeafCandidate::eta(), mps_fire::i, and reco::LeafCandidate::phi().

                                                                                                                {
   bool differentSize = pfJetsReco.size() != pfJetsReReco.size();
   if (differentSize)
     std::cout << "+++WARNING+++ PFJet size changed for entry " << entry_ << " !" << endl
               << " - RECO    size : " << pfJetsReco.size() << endl
               << " - Re-RECO size : " << pfJetsReReco.size() << endl;
   unsigned minSize = pfJetsReco.size() < pfJetsReReco.size() ? pfJetsReco.size() : pfJetsReReco.size();
   unsigned npr = 0;
   for (unsigned i = 0; i < minSize; ++i) {
     const reco::PFJet& candReco = pfJetsReco[i];
     const reco::PFJet& candReReco = pfJetsReReco[i];
     if (candReco.et() < 20. && candReReco.et() < 20.)
       break;
     double deltaE = (candReReco.et() - candReco.et()) / (candReReco.et() + candReco.et());
     double deltaEta = candReReco.eta() - candReco.eta();
     double deltaPhi = candReReco.phi() - candReco.phi();
     if (fabs(deltaE) > deltaEMax_ || fabs(deltaEta) > deltaEtaMax_ || fabs(deltaPhi) > deltaPhiMax_) {
       std::cout << "+++WARNING+++ PFJet " << i << " changed  for entry " << entry_ << " ! " << std::endl
                 << " - RECO     : " << candReco.et() << " " << candReco.eta() << " " << candReco.phi() << std::endl
                 << " - Re-RECO  : " << candReReco.et() << " " << candReReco.eta() << " " << candReReco.phi()
                 << std::endl
                 << " DeltaE   = : " << deltaE << std::endl
                 << " DeltaEta = : " << deltaEta << std::endl
                 << " DeltaPhi = : " << deltaPhi << std::endl
                 << std::endl;
       if (++npr == 5)
         break;
     } else {
       std::cout << "Jet " << i << " " << candReco.et() << std::endl;
     }
   }
 }

◆ printMet()

void PFCandidateChecker::printMet	(	const reco::PFCandidateCollection &	pfReco,
		const reco::PFCandidateCollection &	pfReReco
	)		const

private

Definition at line 168 of file PFCandidateChecker.cc.

References gather_cfg::cout, mps_fire::i, BTaggingMonitor_cfi::met, and mathSSE::sqrt().

                                                                                                                   {
   double metX = 0.;
   double metY = 0.;
   for (unsigned i = 0; i < pfReco.size(); i++) {
     metX += pfReco[i].px();
     metY += pfReco[i].py();
   }
   double met = std::sqrt(metX * metX + metY * metY);
   std::cout << "MET RECO    = " << metX << " " << metY << " " << met << std::endl;
 
   metX = 0.;
   metY = 0.;
   for (unsigned i = 0; i < pfReReco.size(); i++) {
     metX += pfReReco[i].px();
     metY += pfReReco[i].py();
   }
   met = std::sqrt(metX * metX + metY * metY);
   std::cout << "MET Re-RECO = " << metX << " " << metY << " " << met << std::endl;
 }

Member Data Documentation

◆ deltaEMax_

double PFCandidateChecker::deltaEMax_

private

Cuts for comparison.

Definition at line 43 of file PFCandidateChecker.cc.

◆ deltaEtaMax_

double PFCandidateChecker::deltaEtaMax_

private

Definition at line 44 of file PFCandidateChecker.cc.

◆ deltaPhiMax_

double PFCandidateChecker::deltaPhiMax_

private

Definition at line 45 of file PFCandidateChecker.cc.

◆ entry_

unsigned PFCandidateChecker::entry_

private

Counter.

Definition at line 57 of file PFCandidateChecker.cc.

◆ inputTagPFCandidatesReco_

edm::InputTag PFCandidateChecker::inputTagPFCandidatesReco_

private

PFCandidates in which we'll look for pile up particles.

Definition at line 37 of file PFCandidateChecker.cc.

◆ inputTagPFCandidatesReReco_

edm::InputTag PFCandidateChecker::inputTagPFCandidatesReReco_

private

Definition at line 38 of file PFCandidateChecker.cc.

◆ inputTagPFJetsReco_

edm::InputTag PFCandidateChecker::inputTagPFJetsReco_

private

Definition at line 39 of file PFCandidateChecker.cc.

◆ inputTagPFJetsReReco_

edm::InputTag PFCandidateChecker::inputTagPFJetsReReco_

private

Definition at line 40 of file PFCandidateChecker.cc.

◆ printBlocks_

bool PFCandidateChecker::printBlocks_

private

print the blocks associated to a given candidate ?

Definition at line 51 of file PFCandidateChecker.cc.

◆ rankByPt_

bool PFCandidateChecker::rankByPt_

private

rank the candidates by Pt

Definition at line 54 of file PFCandidateChecker.cc.

◆ verbose_

bool PFCandidateChecker::verbose_

private

verbose ?

Definition at line 48 of file PFCandidateChecker.cc.

Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PFCandidateChecker()

Member Function Documentation

◆ analyze()

◆ greaterPt()

◆ printElementsInBlocks()

◆ printJets()

◆ printMet()

Member Data Documentation

◆ deltaEMax_

◆ deltaEtaMax_

◆ deltaPhiMax_

◆ entry_

◆ inputTagPFCandidatesReco_

◆ inputTagPFCandidatesReReco_

◆ inputTagPFJetsReco_

◆ inputTagPFJetsReReco_

◆ printBlocks_

◆ rankByPt_

◆ verbose_