To plot Candidate quantities. More...

#include <MatchCandidateBenchmark.h>

Inheritance diagram for MatchCandidateBenchmark:

Public Member Functions
void	fillOne (const reco::Candidate &candidate, const reco::Candidate &matchedCandidate)
	fill histograms with a given particle
	MatchCandidateBenchmark (Mode mode)
void	setup (const edm::ParameterSet &parameterSet)
	book histograms
void	setup ()
	book histograms
virtual	~MatchCandidateBenchmark ()
Protected Attributes
TH2F *	delta_et_Over_et_VS_et_
TH2F *	delta_et_VS_et_
TH2F *	delta_eta_VS_et_
TH2F *	delta_phi_VS_et_
bool	histogramBooked_

Detailed Description

To plot Candidate quantities.

Definition at line 13 of file MatchCandidateBenchmark.h.

Constructor & Destructor Documentation

MatchCandidateBenchmark::MatchCandidateBenchmark ( Mode mode )

Definition at line 17 of file MatchCandidateBenchmark.cc.

References delta_et_Over_et_VS_et_, delta_et_VS_et_, delta_eta_VS_et_, delta_phi_VS_et_, and histogramBooked_.

                                                           : Benchmark(mode) {
  delta_et_Over_et_VS_et_ = 0; 
  delta_et_VS_et_         = 0; 
  delta_eta_VS_et_        = 0; 
  delta_phi_VS_et_        = 0;

  histogramBooked_ = false;
}

MatchCandidateBenchmark::~MatchCandidateBenchmark ( ) [virtual]

Definition at line 26 of file MatchCandidateBenchmark.cc.

{}

Member Function Documentation

void MatchCandidateBenchmark::fillOne	(	const reco::Candidate &	candidate,
		const reco::Candidate &	matchedCandidate
	)

fill histograms with a given particle

Definition at line 131 of file MatchCandidateBenchmark.cc.

References delta_et_Over_et_VS_et_, delta_et_VS_et_, delta_eta_VS_et_, delta_phi_VS_et_, reco::Candidate::eta(), histogramBooked_, Benchmark::isInRange(), reco::Candidate::phi(), and reco::Candidate::pt().

Referenced by PFCandidateManager::fill(), PFJetMonitor::fill(), and PFMETMonitor::fillOne().

                                                                        {
  
  if( !isInRange(cand.pt(), cand.eta(), cand.phi() ) ) return;

  if (histogramBooked_) {
    if (delta_et_Over_et_VS_et_) delta_et_Over_et_VS_et_->Fill( matchedCand.pt(), (cand.pt() - matchedCand.pt())/matchedCand.pt() );
    if (delta_et_VS_et_) delta_et_VS_et_->Fill( matchedCand.pt(), cand.pt() - matchedCand.pt() );
    if (delta_eta_VS_et_) delta_eta_VS_et_->Fill( matchedCand.pt(), cand.eta() - matchedCand.eta() );
    if (delta_phi_VS_et_) delta_phi_VS_et_->Fill( matchedCand.pt(), cand.phi() - matchedCand.phi() );
  }  

}

void MatchCandidateBenchmark::setup ( void )

book histograms

Definition at line 29 of file MatchCandidateBenchmark.cc.

References Benchmark::book2D(), delta_et_Over_et_VS_et_, delta_et_VS_et_, delta_eta_VS_et_, delta_phi_VS_et_, Benchmark::DQMOFFLINE, histogramBooked_, Benchmark::PhaseSpace::m, Benchmark::PhaseSpace::M, Benchmark::mode_, Benchmark::PhaseSpace::n, jptDQMConfig_cff::ptBins, and Benchmark::VALIDATION.

Referenced by PFMETMonitor::setup(), and PFJetMonitor::setup().

                                    {
  if (!histogramBooked_) {
    PhaseSpace ptPS;
    PhaseSpace dptOvptPS;
    PhaseSpace dptPS;
    PhaseSpace detaPS;
    PhaseSpace dphiPS;
    switch(mode_) {
    case VALIDATION:
      ptPS = PhaseSpace(100,0,1000);
      dptOvptPS = PhaseSpace( 200, -1, 1);
      dphiPS = PhaseSpace( 200, -1, 1);
      detaPS = PhaseSpace( 200, -1, 1);
      dptPS = PhaseSpace( 100, -100, 100);
      break;
    case DQMOFFLINE:
    default:
      ptPS = PhaseSpace(50,0,100);
      dptOvptPS = PhaseSpace( 50, -1, 1);
      dphiPS = PhaseSpace( 50, -1, 1);
      detaPS = PhaseSpace( 50, -1, 1);
      dptPS = PhaseSpace( 50, -50, 50);
      break;
    }
    float ptBins[11] = {0, 1, 2, 5, 10, 20, 50, 100, 200, 400, 1000};
    
    delta_et_Over_et_VS_et_ = book2D("delta_et_Over_et_VS_et_", 
                                     ";E_{T, true} (GeV);#DeltaE_{T}/E_{T}",
                                     10, ptBins, 
                                     dptOvptPS.n, dptOvptPS.m, dptOvptPS.M );
    
    
    delta_et_VS_et_ = book2D("delta_et_VS_et_", 
                             ";E_{T, true} (GeV);#DeltaE_{T}",
                             10, ptBins,
                             dptPS.n, dptPS.m, dptPS.M );
    
    delta_eta_VS_et_ = book2D("delta_eta_VS_et_", 
                              ";#E_{T, true} (GeV);#Delta#eta",
                              10, ptBins,
                              detaPS.n, detaPS.m, detaPS.M );
    
    delta_phi_VS_et_ = book2D("delta_phi_VS_et_", 
                              ";E_{T, true} (GeV);#Delta#phi",
                              10, ptBins,
                              dphiPS.n, dphiPS.m, dphiPS.M );
    histogramBooked_ = true;
  } 
}

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

book histograms

Definition at line 78 of file MatchCandidateBenchmark.cc.

References Benchmark::book2D(), delta_et_Over_et_VS_et_, delta_et_VS_et_, delta_eta_VS_et_, delta_phi_VS_et_, edm::ParameterSet::getParameter(), histogramBooked_, i, and jptDQMConfig_cff::ptBins.

                                                                       {

  if (histogramBooked_) return;

  edm::ParameterSet dptPS = parameterSet.getParameter<edm::ParameterSet>("DeltaPtHistoParameter");
  edm::ParameterSet dptOvptPS = parameterSet.getParameter<edm::ParameterSet>("DeltaPtOvPtHistoParameter");
  edm::ParameterSet detaPS = parameterSet.getParameter<edm::ParameterSet>("DeltaEtaHistoParameter");
  edm::ParameterSet dphiPS = parameterSet.getParameter<edm::ParameterSet>("DeltaPhiHistoParameter");

  std::vector<double> ptBinsPS = parameterSet.getParameter< std::vector<double> >( "VariablePtBins" );
  float* ptBins = new float[ptBinsPS.size()];
  for (size_t i = 0; i < ptBinsPS.size(); i++) {
    ptBins[i] = ptBinsPS[i];
  }
  
  if (dptOvptPS.getParameter<bool>("switchOn")) {
    delta_et_Over_et_VS_et_ = book2D("delta_et_Over_et_VS_et_", 
                                     ";E_{T, true} (GeV);#DeltaE_{T}/E_{T}",
                                     ptBinsPS.size()-1, ptBins, 
                                     dptOvptPS.getParameter<int32_t>("nBin"), 
                                     dptOvptPS.getParameter<double>("xMin"), 
                                     dptOvptPS.getParameter<double>("xMax"));
  }
        
  if (dptPS.getParameter<bool>("switchOn")) {
    delta_et_VS_et_ = book2D("delta_et_VS_et_", 
                             ";E_{T, true} (GeV);#DeltaE_{T}",
                             ptBinsPS.size()-1, ptBins,
                             dptPS.getParameter<int32_t>("nBin"), 
                             dptPS.getParameter<double>("xMin"), 
                             dptPS.getParameter<double>("xMax"));
  }
    
  if (detaPS.getParameter<bool>("switchOn")) {
    delta_eta_VS_et_ = book2D("delta_eta_VS_et_", 
                              ";#E_{T, true} (GeV);#Delta#eta",
                              ptBinsPS.size()-1, ptBins,
                              detaPS.getParameter<int32_t>("nBin"), 
                              detaPS.getParameter<double>("xMin"), 
                              detaPS.getParameter<double>("xMax"));
  }
   
  if (dphiPS.getParameter<bool>("switchOn")) {
    delta_phi_VS_et_ = book2D("delta_phi_VS_et_", 
                              ";E_{T, true} (GeV);#Delta#phi",
                              ptBinsPS.size()-1, ptBins,
                              dphiPS.getParameter<int32_t>("nBin"), 
                              dphiPS.getParameter<double>("xMin"),
                              dphiPS.getParameter<double>("xMax"));
  }
  histogramBooked_ = true;
}