#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 (float dRMax, bool matchCharge, Benchmark::Mode mode, float ptmin, float ptmax, float etamin, float etamax, float phimin, float phimax, bool fracHistoFlag=true)
	set the parameters locally More...

void	setParameters (const edm::ParameterSet &parameterSet)
	set the parameters accessing them from ParameterSet 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_

bool	histogramBooked_

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

Definition at line 15 of file PFJetMonitor.cc.

References 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_, histogramBooked_, and Benchmark::setRange().

                                                                              : 
   Benchmark(mode), 
   candBench_(mode), 
   matchCandBench_(mode), 
   dRMax_(dRMax), 
   matchCharge_(matchCharge) {
   
   setRange( 0.0, 10e10, -10.0, 10.0, -3.14, 3.14);
   
   delta_frac_VS_frac_muon_           = 0;
   delta_frac_VS_frac_photon_         = 0;
   delta_frac_VS_frac_electron_       = 0;
   delta_frac_VS_frac_charged_hadron_ = 0;
   delta_frac_VS_frac_neutral_hadron_ = 0;
   
   createPFractionHistos_ = false;
   histogramBooked_ = false;
 }

PFJetMonitor::~PFJetMonitor ( )

virtual

Definition at line 36 of file PFJetMonitor.cc.

36 {}

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 67 of file PFJetMonitor.h.

References candBench_, createPFractionHistos_, dRMax_, reco::Candidate::eta(), reco::LeafCandidate::eta(), MatchCandidateBenchmark::fillOne(), CandidateBenchmark::fillOne(), fillOne(), histogramBooked_, 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(), and PFRootEventManager::processEntry().

                                                                          {
   
 
   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_ && histogramBooked_) fillOne(jet, matchedJetCollection[ iMatch ]);
     }
   }
 }

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

Definition at line 135 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 126 of file PFJetMonitor.cc.

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

Referenced by PFRootEventManager::readOptions().

                                                {
   Benchmark::setDirectory(dir);
 
   candBench_.setDirectory(dir);
   matchCandBench_.setDirectory(dir);
 }

void PFJetMonitor::setParameters	(	float	dRMax,
		bool	matchCharge,
		Benchmark::Mode	mode,
		float	ptmin,
		float	ptmax,
		float	etamin,
		float	etamax,
		float	phimin,
		float	phimax,
		bool	fracHistoFlag = `true`
	)

set the parameters locally

Definition at line 62 of file PFJetMonitor.cc.

References candBench_, createPFractionHistos_, dRMax_, matchCandBench_, matchCharge_, alignBH_cfg::mode, Benchmark::mode_, Benchmark::setParameters(), and Benchmark::setRange().

Referenced by PFJetDQMAnalyzer::PFJetDQMAnalyzer(), and PFRootEventManager::readOptions().

                                                                  {
   dRMax_                 = dRMax;
   matchCharge_           = matchCharge;
   mode_                  = mode;
   createPFractionHistos_ = fracHistoFlag;
   
   setRange( ptmin, ptmax, etamin, etamax, phimin, phimax );
 
   candBench_.setParameters(mode_);
   matchCandBench_.setParameters(mode_);
 
 }

void PFJetMonitor::setParameters ( const edm::ParameterSet & parameterSet )

set the parameters accessing them from ParameterSet

Definition at line 41 of file PFJetMonitor.cc.

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

                                                                       {
 
   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 102 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_, histogramBooked_, matchCandBench_, CandidateBenchmark::setup(), and MatchCandidateBenchmark::setup().

Referenced by PFJetDQMAnalyzer::beginJob(), and PFRootEventManager::readOptions().

                          {
   candBench_.setup();
   matchCandBench_.setup();
 
   if (createPFractionHistos_ && !histogramBooked_) {
     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);
 
     histogramBooked_ = true;
   }
 
 }

void PFJetMonitor::setup ( const edm::ParameterSet & parameterSet )

book histograms

Definition at line 79 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_, histogramBooked_, matchCandBench_, CandidateBenchmark::setup(), and MatchCandidateBenchmark::setup().

                                                              {
   candBench_.setup(parameterSet);
   matchCandBench_.setup(parameterSet);
 
   if (createPFractionHistos_ && !histogramBooked_) {
     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);
   
     histogramBooked_ = true;
   }
 
 }