PFTauElecRejectionBenchmark.h

Go to the documentation of this file.

#ifndef RecoParticleFlow_Benchmark_PFTauElecRejectionBenchmark_h
#define RecoParticleFlow_Benchmark_PFTauElecRejectionBenchmark_h

#include "DataFormats/TauReco/interface/PFTau.h"
#include "DataFormats/TauReco/interface/PFTauDiscriminator.h"
#include "DataFormats/TauReco/interface/PFTauDiscriminatorByIsolation.h"

#include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidateFwd.h"

#include "Math/GenVector/VectorUtil.h"
#include "Math/GenVector/PxPyPzE4D.h"
#include "TLorentzVector.h"

#include "CommonTools/UtilAlgos/interface/TFileService.h"
#include "TH1F.h"
#include "TH2F.h"
#include <string>
#include <TFile.h>

#include "DQMServices/Core/interface/DQMStore.h"

#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"

class PFTauElecRejectionBenchmark;

class PFTauElecRejectionBenchmark {
public:
  typedef dqm::legacy::DQMStore DQMStore;
  typedef dqm::legacy::MonitorElement MonitorElement;

  PFTauElecRejectionBenchmark();
  virtual ~PFTauElecRejectionBenchmark();

  void setup(std::string Filename,
             std::string benchmarkLabel,
             double maxDeltaR,
             double minRecoPt,
             double maxRecoAbsEta,
             double minMCPt,
             double maxMCAbsEta,
             std::string sGenMatchObjectLabel,
             bool applyEcalCrackCut,
             DQMStore *db_store);
  void process(edm::Handle<edm::HepMCProduct> mcevt,
               edm::Handle<reco::PFTauCollection> pfTaus,
               edm::Handle<reco::PFTauDiscriminator> pfTauIsoDiscr,
               edm::Handle<reco::PFTauDiscriminator> pfTauElecDiscr);
  void write();

private:
  bool isInEcalCrack(double eta) const;

  TFile *file_;
  std::string outputFile_;
  std::string benchmarkLabel_;
  double maxDeltaR_;
  double minMCPt_;
  double maxMCAbsEta_;
  double minRecoPt_;
  double maxRecoAbsEta_;
  std::string sGenMatchObjectLabel_;
  bool applyEcalCrackCut_;

  // histograms
  TH1F *hEoverP;
  TH1F *hHoverP;
  TH1F *hEmfrac;

  TH1F *hEoverP_barrel;
  TH1F *hHoverP_barrel;
  TH1F *hEmfrac_barrel;

  TH1F *hEoverP_endcap;
  TH1F *hHoverP_endcap;
  TH1F *hEmfrac_endcap;

  TH1F *hEoverP_preid0;
  TH1F *hHoverP_preid0;
  TH1F *hEmfrac_preid0;

  TH1F *hEoverP_preid1;
  TH1F *hHoverP_preid1;
  TH1F *hEmfrac_preid1;

  TH1F *hElecPreID;
  TH1F *hElecMVA;
  TH1F *hTauElecDiscriminant;

  TH2F *hHoPvsEoP;
  TH2F *hHoPvsEoP_preid0;
  TH2F *hHoPvsEoP_preid1;

  TH2F *hEmfracvsEoP;
  TH2F *hEmfracvsEoP_preid0;
  TH2F *hEmfracvsEoP_preid1;

  TH1F *hpfcand_deltaEta;
  TH1F *hpfcand_deltaEta_weightE;
  TH1F *hpfcand_deltaPhiOverQ;
  TH1F *hpfcand_deltaPhiOverQ_weightE;

  TH1F *hleadTk_pt;
  TH1F *hleadTk_eta;
  TH1F *hleadTk_phi;

  // to be filled yet!
  TH1F *hleadGsfTk_pt;
  TH1F *hleadGsfTk_eta;
  TH1F *hleadGsfTk_phi;

  std::vector<TLorentzVector> _GenObjects;

protected:
  DQMStore *db_;
};

#endif  // RecoParticleFlow_Benchmark_PFTauElecRejectionBenchmark_h