PFCandidateChecker Class Reference

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
void	beginRun (const edm::Run &r, const edm::EventSetup &c) override
	PFCandidateChecker (const edm::ParameterSet &)
	~PFCandidateChecker () override
Public Member Functions inherited from edm::stream::EDAnalyzer< T >
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
	~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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
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	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< T >
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
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
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

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 72 of file PFCandidateChecker.cc.

                                                                     {
  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_;
}

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

~PFCandidateChecker()

PFCandidateChecker::~PFCandidateChecker ( )

override

Definition at line 99 of file PFCandidateChecker.cc.

{}

Member Function Documentation

analyze()

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

overridevirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 103 of file PFCandidateChecker.cc.

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

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

beginRun()

void PFCandidateChecker::beginRun ( const edm::Run &  r, const edm::EventSetup &  c  )

overridevirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 101 of file PFCandidateChecker.cc.

{}

greaterPt()

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

inlinestaticprivate

Definition at line 63 of file PFCandidateChecker.cc.

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

References a, and b.

printElementsInBlocks()

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

private

Definition at line 229 of file PFCandidateChecker.cc.

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

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

printJets()

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

private

Definition at line 196 of file PFCandidateChecker.cc.

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

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().

printMet()

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

private

Definition at line 176 of file PFCandidateChecker.cc.

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

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

Member Data Documentation

deltaEMax_

double PFCandidateChecker::deltaEMax_

private

Cuts for comparison.

Definition at line 47 of file PFCandidateChecker.cc.

deltaEtaMax_

double PFCandidateChecker::deltaEtaMax_

private

Definition at line 48 of file PFCandidateChecker.cc.

deltaPhiMax_

double PFCandidateChecker::deltaPhiMax_

private

Definition at line 49 of file PFCandidateChecker.cc.

entry_

unsigned PFCandidateChecker::entry_

private

Counter.

Definition at line 61 of file PFCandidateChecker.cc.

inputTagPFCandidatesReco_

edm::InputTag PFCandidateChecker::inputTagPFCandidatesReco_

private

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

Definition at line 41 of file PFCandidateChecker.cc.

inputTagPFCandidatesReReco_

edm::InputTag PFCandidateChecker::inputTagPFCandidatesReReco_

private

Definition at line 42 of file PFCandidateChecker.cc.

inputTagPFJetsReco_

edm::InputTag PFCandidateChecker::inputTagPFJetsReco_

private

Definition at line 43 of file PFCandidateChecker.cc.

inputTagPFJetsReReco_

edm::InputTag PFCandidateChecker::inputTagPFJetsReReco_

private

Definition at line 44 of file PFCandidateChecker.cc.

printBlocks_

bool PFCandidateChecker::printBlocks_

private

print the blocks associated to a given candidate ?

Definition at line 55 of file PFCandidateChecker.cc.

rankByPt_

bool PFCandidateChecker::rankByPt_

private

rank the candidates by Pt

Definition at line 58 of file PFCandidateChecker.cc.

verbose_

bool PFCandidateChecker::verbose_

private

verbose ?

Definition at line 52 of file PFCandidateChecker.cc.