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

#include <PFCandidateChecker.h>

Inheritance diagram for PFCandidateChecker:

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

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

	PFCandidateChecker (const edm::ParameterSet &)

	~PFCandidateChecker ()

Public Member Functions inherited from edm::stream::EDAnalyzer<>
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

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

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

Checks what a re-reco changes in PFCandidates.

Author: Patrick Janot

Date: August 2011

Definition at line 31 of file PFCandidateChecker.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 17 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_;
    
 }

PFCandidateChecker::~PFCandidateChecker ( )

Definition at line 61 of file PFCandidateChecker.cc.

61 { }

Member Function Documentation

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

virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 71 of file PFCandidateChecker.cc.

References gather_cfg::cout, spr::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), edm::EventID::event(), edm::Event::getByLabel(), i, edm::EventBase::id(), LogDebug, reco::LeafCandidate::phi(), and edm::EventID::run().

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

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

virtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 66 of file PFCandidateChecker.cc.

67 { }

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

inlinestaticprivate

Definition at line 78 of file PFCandidateChecker.h.

References reco::LeafCandidate::pt().

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

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

private

Definition at line 214 of file PFCandidateChecker.cc.

References assert(), ecal_dqm_sourceclient-live_cfg::cerr, reco::PFCandidate::elementsInBlocks(), i, and edm::Ref< C, T, F >::isNull().

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

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

private

Definition at line 177 of file PFCandidateChecker.cc.

References gather_cfg::cout, spr::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, reco::LeafCandidate::et(), reco::LeafCandidate::eta(), 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;
     }
   }
 
 }

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

private

Definition at line 154 of file PFCandidateChecker.cc.

References gather_cfg::cout, i, objects.METAnalyzer::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

double PFCandidateChecker::deltaEMax_

private

Cuts for comparison.

Definition at line 62 of file PFCandidateChecker.h.

double PFCandidateChecker::deltaEtaMax_

private

Definition at line 63 of file PFCandidateChecker.h.

double PFCandidateChecker::deltaPhiMax_

private

Definition at line 64 of file PFCandidateChecker.h.

unsigned PFCandidateChecker::entry_

private

Counter.

Definition at line 76 of file PFCandidateChecker.h.

edm::InputTag PFCandidateChecker::inputTagPFCandidatesReco_

private

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

Definition at line 56 of file PFCandidateChecker.h.

edm::InputTag PFCandidateChecker::inputTagPFCandidatesReReco_

private

Definition at line 57 of file PFCandidateChecker.h.

edm::InputTag PFCandidateChecker::inputTagPFJetsReco_

private

Definition at line 58 of file PFCandidateChecker.h.

edm::InputTag PFCandidateChecker::inputTagPFJetsReReco_

private

Definition at line 59 of file PFCandidateChecker.h.

bool PFCandidateChecker::printBlocks_

private

print the blocks associated to a given candidate ?

Definition at line 70 of file PFCandidateChecker.h.

bool PFCandidateChecker::rankByPt_

private

rank the candidates by Pt

Definition at line 73 of file PFCandidateChecker.h.

bool PFCandidateChecker::verbose_

private

verbose ?

Definition at line 67 of file PFCandidateChecker.h.

Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation