PFCandidateChecker Class Reference

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::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

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::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

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.

{ }

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 71 of file PFCandidateChecker.cc.

References gather_cfg::cout, HLTFastRecoForTau_cff::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), edm::EventID::event(), edm::Event::getByLabel(), i, edm::EventBase::id(), LogDebug, reco::LeafCandidate::phi(), edm::EventID::run(), and python.multivaluedict::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;
}

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 66 of file PFCandidateChecker.cc.

                                            { }

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 dtNoiseDBValidation_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, HLTFastRecoForTau_cff::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, CaloMET_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

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.