#include <PFCandidateMonitor.h>

Inheritance diagram for PFCandidateMonitor:

Public Member Functions
template<class T , class C >
void	fill (const T &candidateCollection, const C &matchedCandCollection, float &minVal, float &maxVal, const edm::ParameterSet &parameterSet)
	fill histograms with all particle More...

template<class T , class C , class M >
void	fill (const T &candidateCollection, const C &matchedCandCollection, float &minVal, float &maxVal, const edm::ParameterSet &parameterSet, const M &muonMatchedCandCollection)

void	fillOne (const reco::Candidate &cand)

	PFCandidateMonitor (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 refHistoFlag)
	set the parameters locally More...

void	setParameters (const edm::ParameterSet &parameterSet)
	set the parameters accessing them from ParameterSet More...

void	setup (DQMStore::IBooker &b)
	book histograms More...

void	setup (DQMStore::IBooker &b, const edm::ParameterSet &parameterSet)

virtual	~PFCandidateMonitor ()

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	createEfficiencyHistos_

bool	createReferenceHistos_

TH1F *	deltaR_

float	dRMax_

TH1F *	eta_gen_

TH1F *	eta_ref_

bool	histogramBooked_

MatchCandidateBenchmark	matchCandBench_

bool	matchCharge_

bool	matching_done_

TH1F *	phi_gen_

TH1F *	phi_ref_

TH1F *	pt_gen_

TH1F *	pt_ref_

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 }

Protected Member Functions inherited from Benchmark
TH1F *	book1D (DQMStore::IBooker &b, const char histname, const char title, int nbins, float xmin, float xmax)
	book a 1D histogram, either with DQM or plain root depending if DQM_ has been initialized in a child analyzer or not. More...

TH2F *	book2D (DQMStore::IBooker &b, 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 depending if DQM_ has been initialized in a child analyzer or not. More...

TH2F *	book2D (DQMStore::IBooker &b, 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 depending if DQM_ has been initialized in a child analyzer or not. More...

TProfile *	bookProfile (DQMStore::IBooker &b, const char histname, const char title, int nbinsx, float xmin, float xmax, float ymin, float ymax, const char *option)
	book a TProfile histogram, either with DQM or plain root depending if DQM_ has been initialized in a child analyzer or not. More...

TProfile *	bookProfile (DQMStore::IBooker &b, const char histname, const char title, int nbinsx, float xbins, float ymin, float ymax, const char option)
	book a TProfile histogram, either with DQM or plain root depending if DQM_ has been initialized in a child analyzer or not. More...

Detailed Description

Definition at line 16 of file PFCandidateMonitor.h.

Constructor & Destructor Documentation

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

Definition at line 16 of file PFCandidateMonitor.cc.

References createReferenceHistos_, deltaR_, eta_gen_, eta_ref_, histogramBooked_, phi_gen_, phi_ref_, pt_gen_, pt_ref_, 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);
   
   pt_gen_   = 0;
   eta_gen_  = 0;
   phi_gen_  = 0;
 
   pt_ref_   = 0;
   eta_ref_  = 0;
   phi_ref_  = 0;
 
   deltaR_   = 0;
 
   createReferenceHistos_ = false;
   histogramBooked_ = false;
 }  

PFCandidateMonitor::~PFCandidateMonitor ( )

virtual

Definition at line 43 of file PFCandidateMonitor.cc.

43 {}

Member Function Documentation

template<class T , class C >

void PFCandidateMonitor::fill	(	const T &	candidateCollection,
		const C &	matchedCandCollection,
		float &	minVal,
		float &	maxVal,
		const edm::ParameterSet &	parameterSet
	)

fill histograms with all particle

Definition at line 78 of file PFCandidateMonitor.h.

References assert(), candBench_, createEfficiencyHistos_, createReferenceHistos_, reco::deltaR(), deltaR_, dRMax_, eta, reco::Candidate::eta(), CandidateBenchmark::fillOne(), MatchCandidateBenchmark::fillOne(), fillOne(), i, Benchmark::isInRange(), j, PFB::match(), matchCandBench_, matchCharge_, matching_done_, phi, reco::Candidate::phi(), EnergyCorrector::pt, reco::Candidate::pt(), and findQualityFiles::size.

Referenced by PFCandidateDQMAnalyzer::analyze(), and PFMuonDQMAnalyzer::analyze().

                                                       {
   
   matching_done_ = false;
   if ( createEfficiencyHistos_ ) {
     for( unsigned i=0; i<candCollection.size(); ++i) {
       if( !isInRange(candCollection[i].pt(), candCollection[i].eta(), candCollection[i].phi() ) ) continue;
       fillOne(candCollection[i]); // fill pt_gen, eta_gen and phi_gen histos for UNMATCHED generated candidate
 
       for( unsigned j=0; j<matchedCandCollection.size(); ++j)  // for DeltaR spectrum
     if (deltaR_) deltaR_->Fill( reco::deltaR( candCollection[i], matchedCandCollection[j] ) ) ;
     }
   }
 
   std::vector<int> matchIndices;
   PFB::match( candCollection, matchedCandCollection, matchIndices, matchCharge_, dRMax_ );
   //PFB::match( candCollection, matchedCandCollection, matchIndices, parameterSet, matchCharge_, dRMax_ );
   // now matchIndices[i] stores the j-th closest matched jet 
   matching_done_ = true;
 
   for (unsigned int i = 0; i < (candCollection).size(); i++) {
      const reco::Candidate& cand = candCollection[i];
 
     if( !isInRange(cand.pt(), cand.eta(), cand.phi() ) ) continue;
     
     int iMatch = matchIndices[i];
     assert(iMatch< static_cast<int>(matchedCandCollection.size()));
 
     if( iMatch!=-1 ) {
       const reco::Candidate& matchedCand = matchedCandCollection[ iMatch ];
       if(!isInRange(matchedCand.pt(),matchedCand.eta(),matchedCand.phi() ) ) continue;
       //std::cout <<"PFJet pT " <<cand.pt() <<" eta " <<cand.eta() <<" phi " <<cand.phi() ;
       //std::cout <<"\nmatched genJet pT " <<matchedCand.pt() <<" eta " <<matchedCand.eta() <<" phi " <<matchedCand.phi() <<"\n" <<std::endl ;
       float ptRes = (cand.pt() - matchedCand.pt()) / matchedCand.pt();
       
       if (ptRes > maxVal) maxVal = ptRes;
       if (ptRes < minVal) minVal = ptRes;
  
       if ( !createEfficiencyHistos_ ) {
     candBench_.fillOne(cand);   // fill pt, eta phi and charge histos for MATCHED candidate
     //matchCandBench_.fillOne(cand, matchedCand);  // fill delta_x_VS_y histos for matched couple
     matchCandBench_.fillOne(cand, matchedCand, parameterSet);  // fill delta_x_VS_y histos for matched couple
     if (createReferenceHistos_) fillOne(matchedCand); // fill pt_ref, eta_ref and phi_ref histos for MATCHED reference candidate
       }
       else {
     candBench_.fillOne(matchedCand);   // fill pt, eta phi and charge histos for MATCHED candidate
     //matchCandBench_.fillOne(matchedCand, cand);  // fill delta_x_VS_y histos for matched couple
     matchCandBench_.fillOne(cand, matchedCand, parameterSet);  // fill delta_x_VS_y histos for matched couple
     if (createReferenceHistos_) fillOne(cand); // fill pt_ref, eta_ref and phi_ref histos for MATCHED reference candidate
       }
 
     }
   }
 }

template<class T , class C , class M >

void PFCandidateMonitor::fill	(	const T &	candidateCollection,
		const C &	matchedCandCollection,
		float &	minVal,
		float &	maxVal,
		const edm::ParameterSet &	parameterSet,
		const M &	muonMatchedCandCollection
	)

Definition at line 137 of file PFCandidateMonitor.h.

References assert(), candBench_, createEfficiencyHistos_, createReferenceHistos_, reco::deltaR(), deltaR_, dRMax_, eta, reco::Candidate::eta(), CandidateBenchmark::fillOne(), MatchCandidateBenchmark::fillOne(), fillOne(), i, Benchmark::isInRange(), j, PFB::match(), matchCandBench_, matchCharge_, matching_done_, phi, reco::Candidate::phi(), EnergyCorrector::pt, reco::Candidate::pt(), and findQualityFiles::size.

                                                                                           {
 
   
   matching_done_ = false;
   if ( createEfficiencyHistos_ ) {
     for( unsigned i=0; i<candCollection.size(); ++i) {
       if( !isInRange(candCollection[i].pt(), candCollection[i].eta(), candCollection[i].phi() ) ) continue;
       fillOne(candCollection[i]); // fill pt_gen, eta_gen and phi_gen histos for UNMATCHED generated candidate
 
       for( unsigned j=0; j<matchedCandCollection.size(); ++j)  // for DeltaR spectrum
     if (deltaR_) deltaR_->Fill( reco::deltaR( candCollection[i], matchedCandCollection[j] ) ) ;
     }
   }
 
   std::vector<int> matchIndices;
   //PFB::match( candCollection, matchedCandCollection, matchIndices, matchCharge_, dRMax_ );
   //PFB::match( candCollection, matchedCandCollection, matchIndices, parameterSet, matchCharge_, dRMax_ );
   PFB::match( candCollection, matchedCandCollection, matchIndices, parameterSet, muonMatchedCandCollection, matchCharge_, dRMax_ );
   // now matchIndices[i] stores the j-th closest matched jet 
   matching_done_ = true;
 
   for (unsigned int i = 0; i < (candCollection).size(); i++) {
      const reco::Candidate& cand = candCollection[i];
 
     if( !isInRange(cand.pt(), cand.eta(), cand.phi() ) ) continue;
     
     int iMatch = matchIndices[i];
     assert(iMatch< static_cast<int>(matchedCandCollection.size()));
 
     if( iMatch!=-1 ) {
       const reco::Candidate& matchedCand = matchedCandCollection[ iMatch ];
       if(!isInRange(matchedCand.pt(),matchedCand.eta(),matchedCand.phi() ) ) continue;
       //std::cout <<"PFJet pT " <<cand.pt() <<" eta " <<cand.eta() <<" phi " <<cand.phi() ;
       //std::cout <<"\nmatched genJet pT " <<matchedCand.pt() <<" eta " <<matchedCand.eta() <<" phi " <<matchedCand.phi() <<"\n" <<std::endl ;
       float ptRes = (cand.pt() - matchedCand.pt()) / matchedCand.pt();
       
       if (ptRes > maxVal) maxVal = ptRes;
       if (ptRes < minVal) minVal = ptRes;
  
       if ( !createEfficiencyHistos_ ) {
     candBench_.fillOne(cand);   // fill pt, eta phi and charge histos for MATCHED candidate
     matchCandBench_.fillOne(cand, matchedCand, parameterSet);  // fill delta_x_VS_y histos for matched couple
     if (createReferenceHistos_) fillOne(matchedCand); // fill pt_ref, eta_ref and phi_ref histos for MATCHED reference candidate
       }
       else {
     candBench_.fillOne(matchedCand);   // fill pt, eta phi and charge histos for MATCHED candidate
     matchCandBench_.fillOne(cand, matchedCand, parameterSet);  // fill delta_x_VS_y histos for matched couple
     if (createReferenceHistos_) fillOne(cand); // fill pt_ref, eta_ref and phi_ref histos for MATCHED reference candidate
       }
 
     }
   }
 }

void PFCandidateMonitor::fillOne ( const reco::Candidate & cand )

Definition at line 192 of file PFCandidateMonitor.cc.

References createEfficiencyHistos_, createReferenceHistos_, reco::Candidate::eta(), eta_gen_, eta_ref_, histogramBooked_, matching_done_, reco::Candidate::phi(), phi_gen_, phi_ref_, reco::Candidate::pt(), pt_gen_, and pt_ref_.

Referenced by fill().

                                                           {
 
   if (matching_done_) {
     if (createReferenceHistos_ && histogramBooked_) {
       if (pt_ref_) pt_ref_->Fill( cand.pt() );
       if (eta_ref_) eta_ref_->Fill( cand.eta() );
       if (phi_ref_) phi_ref_->Fill( cand.phi() );
     }
   } else if (createEfficiencyHistos_ && histogramBooked_) {
     if (pt_gen_) pt_gen_->Fill( cand.pt() );
     if (eta_gen_) eta_gen_->Fill( cand.eta() );
     if (phi_gen_) phi_gen_->Fill( cand.phi() );
   }
 
 }

void PFCandidateMonitor::setDirectory ( TDirectory * dir )

virtual

set directory (to use in ROOT)

Reimplemented from Benchmark.

Definition at line 181 of file PFCandidateMonitor.cc.

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

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

void PFCandidateMonitor::setParameters	(	float	dRMax,
		bool	matchCharge,
		Benchmark::Mode	mode,
		float	ptmin,
		float	ptmax,
		float	etamin,
		float	etamax,
		float	phimin,
		float	phimax,
		bool	refHistoFlag
	)

set the parameters locally

Definition at line 72 of file PFCandidateMonitor.cc.

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

Referenced by PFCandidateDQMAnalyzer::PFCandidateDQMAnalyzer(), and PFMuonDQMAnalyzer::PFMuonDQMAnalyzer().

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

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

set the parameters accessing them from ParameterSet

Definition at line 49 of file PFCandidateMonitor.cc.

References candBench_, createEfficiencyHistos_, createReferenceHistos_, 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" );
   createReferenceHistos_  = parameterSet.getParameter<bool>( "CreateReferenceHistos" );
   createEfficiencyHistos_ = parameterSet.getParameter<bool>( "CreateEfficiencyHistos" );
   
   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 PFCandidateMonitor::setup ( DQMStore::IBooker & b )

book histograms

Definition at line 149 of file PFCandidateMonitor.cc.

References Benchmark::book1D(), candBench_, createEfficiencyHistos_, createReferenceHistos_, eta_gen_, eta_ref_, histogramBooked_, Benchmark::PhaseSpace::m, Benchmark::PhaseSpace::M, matchCandBench_, Benchmark::PhaseSpace::n, phi_gen_, phi_ref_, pt_gen_, pt_ref_, CandidateBenchmark::setup(), and MatchCandidateBenchmark::setup().

Referenced by PFCandidateDQMAnalyzer::bookHistograms(), and PFMuonDQMAnalyzer::bookHistograms().

                                                  {
   candBench_.setup(b);
   matchCandBench_.setup(b);
 
   if (createReferenceHistos_ && !histogramBooked_) {
     PhaseSpace ptPS(100,0,100);
     PhaseSpace phiPS(360, -3.1416, 3.1416);
     PhaseSpace etaPS(100, -5,5);
     
     pt_ref_ = book1D(b, "pt_ref_", "p_{T}_ref;p_{T} (GeV)", ptPS.n, ptPS.m, ptPS.M);
     if (createEfficiencyHistos_) {
       pt_gen_ = book1D(b, "pt_gen_", "p_{T}_gen;p_{T} (GeV)", ptPS.n, ptPS.m, ptPS.M);
     }
 
     eta_ref_ = book1D(b, "eta_ref_", "#eta_ref;#eta", etaPS.n, etaPS.m, etaPS.M);
     if (createEfficiencyHistos_) {
       eta_gen_ = book1D(b, "eta_gen_", "#eta_gen;#eta", etaPS.n, etaPS.m, etaPS.M);
     }
 
     phi_ref_ = book1D(b, "phi_ref_", "#phi_ref;#phi", phiPS.n, phiPS.m, phiPS.M);
     if (createEfficiencyHistos_) {
       phi_gen_ = book1D(b, "phi_gen_", "#phi_gen;#phi", phiPS.n, phiPS.m, phiPS.M);
     }
 
     histogramBooked_ = true;
   }
 }

void PFCandidateMonitor::setup	(	DQMStore::IBooker &	b,
		const edm::ParameterSet &	parameterSet
	)

Definition at line 90 of file PFCandidateMonitor.cc.

References Benchmark::book1D(), candBench_, createEfficiencyHistos_, createReferenceHistos_, deltaR_, PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, eta_gen_, eta_ref_, edm::ParameterSet::getParameter(), histogramBooked_, matchCandBench_, phi_gen_, phi_ref_, pt_gen_, pt_ref_, CandidateBenchmark::setup(), and MatchCandidateBenchmark::setup().

                                                                                        {
   candBench_.setup(b, parameterSet);
   matchCandBench_.setup(b, parameterSet);
 
   if (createReferenceHistos_ && !histogramBooked_) {
     edm::ParameterSet ptPS  = parameterSet.getParameter<edm::ParameterSet>("PtHistoParameter");
     edm::ParameterSet etaPS = parameterSet.getParameter<edm::ParameterSet>("EtaHistoParameter");
     edm::ParameterSet phiPS = parameterSet.getParameter<edm::ParameterSet>("PhiHistoParameter");
 
     edm::ParameterSet dR = parameterSet.getParameter<edm::ParameterSet>("DeltaRHistoParameter");
 
     if (ptPS.getParameter<bool>("switchOn")) {
       pt_ref_ = book1D(b, "pt_ref_", "p_{T}_ref;p_{T} (GeV)", ptPS.getParameter<int32_t>("nBin"), 
                ptPS.getParameter<double>("xMin"),
                ptPS.getParameter<double>("xMax"));
       if (createEfficiencyHistos_) {
     pt_gen_ = book1D(b, "pt_gen_", "p_{T}_gen;p_{T} (GeV)", ptPS.getParameter<int32_t>("nBin"), 
              ptPS.getParameter<double>("xMin"),
              ptPS.getParameter<double>("xMax") ) ;
       }
     } 
     
     if (etaPS.getParameter<bool>("switchOn")) {
       eta_ref_ = book1D(b, "eta_ref_", "#eta_ref;#eta", etaPS.getParameter<int32_t>("nBin"), 
             etaPS.getParameter<double>("xMin"),
             etaPS.getParameter<double>("xMax"));
       if (createEfficiencyHistos_) {
     eta_gen_ = book1D(b, "eta_gen_", "#eta_gen;#eta", etaPS.getParameter<int32_t>("nBin"), 
               etaPS.getParameter<double>("xMin"),
               etaPS.getParameter<double>("xMax") ) ;
       }
     }
 
     if (phiPS.getParameter<bool>("switchOn")) {
       phi_ref_ = book1D(b, "phi_ref_", "#phi_ref;#phi", phiPS.getParameter<int32_t>("nBin"), 
             phiPS.getParameter<double>("xMin"),
             phiPS.getParameter<double>("xMax"));
       if (createEfficiencyHistos_) {
     phi_gen_ = book1D(b, "phi_gen_", "#phi_gen;#phi", phiPS.getParameter<int32_t>("nBin"), 
               phiPS.getParameter<double>("xMin"),
               phiPS.getParameter<double>("xMax") ) ;
       }
     }
 
     if ( createEfficiencyHistos_ && dR.getParameter<bool>("switchOn") ) { 
       deltaR_ = book1D(b, "deltaR_", "#DeltaR;#DeltaR",
                dR.getParameter<int32_t>("nBin"), 
                dR.getParameter<double>("xMin"),
                dR.getParameter<double>("xMax"));
     }
 
     histogramBooked_ = true;   
   }
 }