PFJetMonitor Class Reference

#include <PFJetMonitor.h>

Inheritance diagram for PFJetMonitor:

Public Member Functions
template<class T , class C >
void	fill (const T &jetCollection, const C &matchedJetCollection, float &minVal, float &maxVal)
	fill histograms with all particle More...
void	fillOne (const reco::Jet &jet, const reco::Jet &matchedJet)
	PFJetMonitor (float dRMax=0.3, bool matchCharge=true, Benchmark::Mode mode=Benchmark::DEFAULT)
void	setDirectory (TDirectory *dir)
	set directory (to use in ROOT) More...
void	setParameters (const edm::ParameterSet &parameterSet)
	set the parameters More...
void	setup ()
	book histograms More...
void	setup (const edm::ParameterSet &parameterSet)
	book histograms More...
virtual	~PFJetMonitor ()
Public Member Functions inherited from Benchmark
	Benchmark (Mode mode=DEFAULT)
bool	isInRange (float pt, float eta, float phi) const
void	setParameters (Mode mode)
void	setRange (float ptMin, float ptMax, float etaMin, float etaMax, float phiMin, float phiMax)
void	write ()
	write to the TFile, in plain ROOT mode. No need to call this function in DQM mode More...
virtual	~Benchmark ()

Protected Attributes
CandidateBenchmark	candBench_
bool	createPFractionHistos_
TH2F *	delta_frac_VS_frac_charged_hadron_
TH2F *	delta_frac_VS_frac_electron_
TH2F *	delta_frac_VS_frac_muon_
TH2F *	delta_frac_VS_frac_neutral_hadron_
TH2F *	delta_frac_VS_frac_photon_
float	dRMax_
MatchCandidateBenchmark	matchCandBench_
bool	matchCharge_
Protected Attributes inherited from Benchmark
TDirectory *	dir_
float	etaMax_
float	etaMin_
Mode	mode_
float	phiMax_
float	phiMin_
float	ptMax_
float	ptMin_

Additional Inherited Members
Public Types inherited from Benchmark
enum	Mode { DEFAULT, DQMOFFLINE, VALIDATION }
Static Public Attributes inherited from Benchmark
static DQMStore *	DQM_ = 0
Protected Member Functions inherited from Benchmark
TH1F *	book1D (const char *histname, const char *title, int nbins, float xmin, float xmax)
	book a 1D histogram, either with DQM or plain root. More...
TH2F *	book2D (const char *histname, const char *title, int nbinsx, float xmin, float xmax, int nbinsy, float ymin, float ymax)
	book a 2D histogram, either with DQM or plain root. More...
TH2F *	book2D (const char *histname, const char *title, int nbinsx, float *xbins, int nbinsy, float ymin, float ymax)
	book a 2D histogram, either with DQM or plain root. More...

Detailed Description

Definition at line 12 of file PFJetMonitor.h.

Constructor & Destructor Documentation

PFJetMonitor::PFJetMonitor ( float  dRMax = 0.3, bool  matchCharge = true, Benchmark::Mode  mode = Benchmark::DEFAULT  )

inline

Definition at line 16 of file PFJetMonitor.h.

    : 
    Benchmark(mode), 
    candBench_(mode), matchCandBench_(mode), 
    dRMax_(dRMax), matchCharge_(matchCharge) {}

PFJetMonitor::~PFJetMonitor ( )

virtual

Definition at line 13 of file PFJetMonitor.cc.

{}

Member Function Documentation

template<class T , class C >

void PFJetMonitor::fill	(	const T &	jetCollection,
		const C &	matchedJetCollection,
		float &	minVal,
		float &	maxVal
	)

fill histograms with all particle

Definition at line 65 of file PFJetMonitor.h.

References candBench_, createPFractionHistos_, dRMax_, reco::Candidate::eta(), reco::LeafCandidate::eta(), MatchCandidateBenchmark::fillOne(), CandidateBenchmark::fillOne(), fillOne(), i, Benchmark::isInRange(), metsig::jet, PFB::match(), matchCandBench_, matchCharge_, reco::Candidate::phi(), reco::LeafCandidate::phi(), reco::Candidate::pt(), reco::LeafCandidate::pt(), and findQualityFiles::size.

Referenced by PFJetDQMAnalyzer::analyze().

                                                                         {
  

  std::vector<int> matchIndices;
  PFB::match( jetCollection, matchedJetCollection, matchIndices, 
          matchCharge_, dRMax_ );

  for (unsigned int i = 0; i < (jetCollection).size(); i++) {
    const reco::Jet& jet = jetCollection[i];

    if( !isInRange(jet.pt(), jet.eta(), jet.phi() ) ) continue;
    
    int iMatch = matchIndices[i];
    assert(iMatch< static_cast<int>(matchedJetCollection.size()));

    if( iMatch!=-1 ) {
      const reco::Candidate& matchedJet = matchedJetCollection[ iMatch ];
      if(!isInRange(matchedJet.pt(),matchedJet.eta(),matchedJet.phi() ) ) continue;
      float ptRes = (jet.pt() - matchedJet.pt())/matchedJet.pt();
      
      if (ptRes > maxVal) maxVal = ptRes;
      if (ptRes < minVal) minVal = ptRes;
 
      candBench_.fillOne(jet);
      matchCandBench_.fillOne(jet, matchedJetCollection[ iMatch ]);
      if (createPFractionHistos_) fillOne(jet, matchedJetCollection[ iMatch ]);
    }
  }
}

void PFJetMonitor::fillOne	(	const reco::Jet &	jet,
		const reco::Jet &	matchedJet
	)

Definition at line 78 of file PFJetMonitor.cc.

References reco::PFJet::chargedHadronMultiplicity(), createPFractionHistos_, delta_frac_VS_frac_charged_hadron_, delta_frac_VS_frac_electron_, delta_frac_VS_frac_muon_, delta_frac_VS_frac_neutral_hadron_, delta_frac_VS_frac_photon_, reco::PFJet::electronMultiplicity(), metsig::jet, reco::PFJet::muonMultiplicity(), reco::PFJet::neutralHadronMultiplicity(), and reco::PFJet::photonMultiplicity().

Referenced by fill().

                                       {

  const reco::PFJet* pfJet = dynamic_cast<const reco::PFJet*>(&jet);
  const reco::PFJet* pfMatchedJet = dynamic_cast<const reco::PFJet*>(&matchedJet);
  if (pfJet && pfMatchedJet && createPFractionHistos_) {
    float frac_muon = -99.9;
    float frac_elec = -99.9; 
    float frac_phot = -99.9;
    float frac_ch_had = -99.9;
    float frac_neu_had = -99.9;
 
    if (pfMatchedJet->muonMultiplicity() > 0) frac_muon = (pfJet->muonMultiplicity() - pfMatchedJet->muonMultiplicity())*1.0/pfMatchedJet->muonMultiplicity(); 
    if (pfMatchedJet->chargedHadronMultiplicity() > 0) frac_ch_had = (pfJet->chargedHadronMultiplicity() - pfMatchedJet->chargedHadronMultiplicity())*1.0/pfMatchedJet->chargedHadronMultiplicity(); 
    if (pfMatchedJet->neutralHadronMultiplicity() > 0) frac_neu_had = (pfJet->neutralHadronMultiplicity() - pfMatchedJet->neutralHadronMultiplicity())*1.0/pfMatchedJet->neutralHadronMultiplicity(); 
    if (pfMatchedJet->photonMultiplicity() > 0) frac_phot = (pfJet->photonMultiplicity() - pfMatchedJet->photonMultiplicity())*1.0/pfMatchedJet->photonMultiplicity(); 
    if (pfMatchedJet->electronMultiplicity() > 0) frac_elec = (pfJet->electronMultiplicity() - pfMatchedJet->electronMultiplicity())*1.0/pfMatchedJet->electronMultiplicity(); 

    delta_frac_VS_frac_muon_->Fill(frac_muon);
    delta_frac_VS_frac_electron_->Fill(frac_elec);
    delta_frac_VS_frac_photon_->Fill(frac_phot);
    delta_frac_VS_frac_charged_hadron_->Fill(frac_ch_had);
    delta_frac_VS_frac_neutral_hadron_->Fill(frac_neu_had);
  }
}

void PFJetMonitor::setDirectory ( TDirectory *  dir )

virtual

set directory (to use in ROOT)

Reimplemented from Benchmark.

Definition at line 71 of file PFJetMonitor.cc.

References candBench_, matchCandBench_, and Benchmark::setDirectory().

                                               {
  Benchmark::setDirectory(dir);

  candBench_.setDirectory(dir);
  matchCandBench_.setDirectory(dir);
}

void PFJetMonitor::setParameters

(

const edm::ParameterSet &

parameterSet

)

set the parameters

Definition at line 16 of file PFJetMonitor.cc.

References candBench_, createPFractionHistos_, dRMax_, edm::ParameterSet::getParameter(), matchCandBench_, matchCharge_, Benchmark::mode_, Benchmark::setParameters(), and Benchmark::setRange().

Referenced by PFJetDQMAnalyzer::PFJetDQMAnalyzer().

                                                                      {

  dRMax_                 = parameterSet.getParameter<double>( "deltaRMax" );
  matchCharge_           = parameterSet.getParameter<bool>( "matchCharge" );
  mode_                  = (Benchmark::Mode) parameterSet.getParameter<int>( "mode" );
  createPFractionHistos_ = parameterSet.getParameter<bool>( "CreatePFractionHistos" );

  
  setRange( parameterSet.getParameter<double>("ptMin"),
        parameterSet.getParameter<double>("ptMax"),
        parameterSet.getParameter<double>("etaMin"),
        parameterSet.getParameter<double>("etaMax"),
        parameterSet.getParameter<double>("phiMin"),
        parameterSet.getParameter<double>("phiMax") );

  candBench_.setParameters(mode_);
  matchCandBench_.setParameters(mode_);
}

void PFJetMonitor::setup

(

void

)

book histograms

Definition at line 53 of file PFJetMonitor.cc.

References Benchmark::book2D(), candBench_, createPFractionHistos_, delta_frac_VS_frac_charged_hadron_, delta_frac_VS_frac_electron_, delta_frac_VS_frac_muon_, delta_frac_VS_frac_neutral_hadron_, delta_frac_VS_frac_photon_, matchCandBench_, CandidateBenchmark::setup(), and MatchCandidateBenchmark::setup().

Referenced by PFJetDQMAnalyzer::beginJob().

                         {
  candBench_.setup();
  matchCandBench_.setup();

  if (createPFractionHistos_) {
    delta_frac_VS_frac_muon_ = book2D("delta_frac_VS_frac_muon_", 
 "#DeltaFraction_Vs_Fraction(muon)", 100, 0.0, 1.0, 100, -1.0, 1.0);
    delta_frac_VS_frac_photon_ = book2D("delta_frac_VS_frac_photon_", 
 "#DeltaFraction_Vs_Fraction(photon)", 100, 0.0, 1.0, 100, -1.0, 1.0);
    delta_frac_VS_frac_electron_ = book2D("delta_frac_VS_frac_electron_", 
 "#DeltaFraction_Vs_Fraction(electron)", 100, 0.0, 1.0, 100, -1.0, 1.0);
    delta_frac_VS_frac_charged_hadron_ = book2D("delta_frac_VS_frac_charged_hadron_", 
 "#DeltaFraction_Vs_Fraction(charged hadron)", 100, 0.0, 1.0, 100, -1.0, 1.0);
    delta_frac_VS_frac_neutral_hadron_ = book2D("delta_frac_VS_frac_neutral_hadron_", 
 "#DeltaFraction_Vs_Fraction(neutral hadron)", 100, 0.0, 1.0, 100, -1.0, 1.0);
  }

}

void PFJetMonitor::setup

(

const edm::ParameterSet &

parameterSet

)

book histograms

Definition at line 35 of file PFJetMonitor.cc.

References Benchmark::book2D(), candBench_, createPFractionHistos_, delta_frac_VS_frac_charged_hadron_, delta_frac_VS_frac_electron_, delta_frac_VS_frac_muon_, delta_frac_VS_frac_neutral_hadron_, delta_frac_VS_frac_photon_, matchCandBench_, CandidateBenchmark::setup(), and MatchCandidateBenchmark::setup().

                                                             {
  candBench_.setup(parameterSet);
  matchCandBench_.setup(parameterSet);

  if (createPFractionHistos_) {
    delta_frac_VS_frac_muon_ = book2D("delta_frac_VS_frac_muon_", 
 "#DeltaFraction_Vs_Fraction(muon)", 100, 0.0, 1.0, 100, -1.0, 1.0);
    delta_frac_VS_frac_photon_ = book2D("delta_frac_VS_frac_photon_", 
 "#DeltaFraction_Vs_Fraction(photon)", 100, 0.0, 1.0, 100, -1.0, 1.0);
    delta_frac_VS_frac_electron_ = book2D("delta_frac_VS_frac_electron_", 
 "#DeltaFraction_Vs_Fraction(electron)", 100, 0.0, 1.0, 100, -1.0, 1.0);
    delta_frac_VS_frac_charged_hadron_ = book2D("delta_frac_VS_frac_charged_hadron_", 
 "#DeltaFraction_Vs_Fraction(charged hadron)", 100, 0.0, 1.0, 100, -1.0, 1.0);
    delta_frac_VS_frac_neutral_hadron_ = book2D("delta_frac_VS_frac_neutral_hadron_", 
 "#DeltaFraction_Vs_Fraction(neutral hadron)", 100, 0.0, 1.0, 100, -1.0, 1.0);
  }

}

Member Data Documentation

CandidateBenchmark PFJetMonitor::candBench_

protected

Definition at line 48 of file PFJetMonitor.h.

Referenced by fill(), setDirectory(), setParameters(), and setup().

bool PFJetMonitor::createPFractionHistos_

protected

Definition at line 59 of file PFJetMonitor.h.

Referenced by fill(), fillOne(), setParameters(), and setup().

TH2F* PFJetMonitor::delta_frac_VS_frac_charged_hadron_

protected

Definition at line 54 of file PFJetMonitor.h.

Referenced by fillOne(), and setup().

TH2F* PFJetMonitor::delta_frac_VS_frac_electron_

protected

Definition at line 53 of file PFJetMonitor.h.

Referenced by fillOne(), and setup().

TH2F* PFJetMonitor::delta_frac_VS_frac_muon_

protected

Definition at line 51 of file PFJetMonitor.h.

Referenced by fillOne(), and setup().

TH2F* PFJetMonitor::delta_frac_VS_frac_neutral_hadron_

protected

Definition at line 55 of file PFJetMonitor.h.

Referenced by fillOne(), and setup().

TH2F* PFJetMonitor::delta_frac_VS_frac_photon_

protected

Definition at line 52 of file PFJetMonitor.h.

Referenced by fillOne(), and setup().

float PFJetMonitor::dRMax_

protected

Definition at line 57 of file PFJetMonitor.h.

Referenced by fill(), and setParameters().

MatchCandidateBenchmark PFJetMonitor::matchCandBench_

protected

Definition at line 49 of file PFJetMonitor.h.

Referenced by fill(), setDirectory(), setParameters(), and setup().

bool PFJetMonitor::matchCharge_

protected

Definition at line 58 of file PFJetMonitor.h.

Referenced by fill(), and setParameters().