MatchCandidateBenchmark.h

Go to the documentation of this file.

#ifndef RecoParticleFlow_Benchmark_MatchCandidateBenchmark_h
#define RecoParticleFlow_Benchmark_MatchCandidateBenchmark_h

#include "DQMOffline/PFTau/interface/Benchmark.h"

#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include <vector>

class MatchCandidateBenchmark : public Benchmark {
public:
  MatchCandidateBenchmark(Mode mode);

  ~MatchCandidateBenchmark() override;

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

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

  void fillOne(const reco::Candidate &candidate,
               const reco::Candidate &matchedCandidate,
               const edm::ParameterSet &parameterSet);

protected:
  TH2F *delta_et_Over_et_VS_et_;
  TH2F *delta_et_VS_et_;
  TH2F *delta_eta_VS_et_;
  TH2F *delta_phi_VS_et_;

  TH2F *BRdelta_et_Over_et_VS_et_;
  TH2F *ERdelta_et_Over_et_VS_et_;
  std::vector<TH1F *> pTRes_;
  std::vector<TH1F *> BRpTRes_;
  std::vector<TH1F *> ERpTRes_;
  std::vector<float> ptBins_;

  bool histogramBooked_;
  double eta_min_barrel_;
  double eta_max_barrel_;
  double eta_min_endcap_;
  double eta_max_endcap_;

private:
  void computePtBins(const edm::ParameterSet &, const edm::ParameterSet &);
  bool inEtaRange(double, bool);
  inline bool inBarrelRange(double value) { return inEtaRange(value, true); }
  inline bool inEndcapRange(double value) { return inEtaRange(value, false); }
};

#endif