PhotonPostprocessing.h

Go to the documentation of this file.

#ifndef PhotonPostprocessing_H
#define PhotonPostprocessing_H

#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CommonTools/UtilAlgos/interface/TFileService.h"
//
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Exception.h"

#include "TFile.h"
#include "TH1.h"
#include "TH2.h"
#include "TTree.h"
#include "TVector3.h"
#include "TProfile.h"
//

#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"

//DQM services
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"

//

#include <vector>

// forward declarations
class TFile;
class TH1F;
class TH2F;
class TProfile;
class TTree;

class PhotonPostprocessing : public edm::one::EDAnalyzer<edm::one::WatchRuns> {
public:
  typedef dqm::legacy::DQMStore DQMStore;
  typedef dqm::legacy::MonitorElement MonitorElement;

  //
  explicit PhotonPostprocessing(const edm::ParameterSet& pset);
  ~PhotonPostprocessing() override;

  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void beginJob() override;
  void endJob() override;
  void beginRun(const edm::Run&, const edm::EventSetup&) override {}
  void endRun(const edm::Run&, const edm::EventSetup&) override;

private:
  //

  void dividePlots(MonitorElement* dividend, MonitorElement* numerator, MonitorElement* denominator, std::string type);
  void dividePlots(MonitorElement* dividend, MonitorElement* numerator, double denominator);
  virtual void runPostprocessing();

  DQMStore* dbe_;
  int verbosity_;

  edm::ParameterSet parameters_;

  bool standAlone_;
  bool batch_;
  bool isRunCentrally_;
  bool fastSim_;
  std::string analyzerName_;
  std::string outputFileName_;
  std::string inputFileName_;

  std::stringstream currentFolder_;

  double etMin;
  double etMax;
  int etBin;
  double etaMin;
  double etaMax;
  int etaBin;
  int etaBin2;
  double phiMin;
  double phiMax;
  int phiBin;
  double rMin;
  double rMax;
  int rBin;
  double zMin;
  double zMax;
  int zBin;

  MonitorElement* phoRecoEffEta_;
  MonitorElement* phoRecoEffPhi_;
  MonitorElement* phoRecoEffEt_;

  MonitorElement* phoDeadChEta_;
  MonitorElement* phoDeadChPhi_;
  MonitorElement* phoDeadChEt_;

  MonitorElement* convEffEtaTwoTracks_;
  MonitorElement* convEffPhiTwoTracks_;
  MonitorElement* convEffRTwoTracks_;
  MonitorElement* convEffZTwoTracks_;
  MonitorElement* convEffEtTwoTracks_;

  MonitorElement* convEffEtaTwoTracksAndVtxProbGT0_;
  MonitorElement* convEffEtaTwoTracksAndVtxProbGT0005_;
  MonitorElement* convEffRTwoTracksAndVtxProbGT0_;
  MonitorElement* convEffRTwoTracksAndVtxProbGT0005_;

  MonitorElement* convEffEtaOneTrack_;
  MonitorElement* convEffPhiOneTrack_;
  MonitorElement* convEffROneTrack_;
  MonitorElement* convEffZOneTrack_;
  MonitorElement* convEffEtOneTrack_;

  MonitorElement* h_simConvEtaMTotal_;
  MonitorElement* h_simConvPhiMTotal_;
  MonitorElement* h_simConvRMTotal_;
  MonitorElement* h_simConvZMTotal_;
  MonitorElement* h_simConvEtMTotal_;

  MonitorElement* convEffEtaTotal_;
  MonitorElement* convEffPhiTotal_;
  MonitorElement* convEffRTotal_;
  MonitorElement* convEffZTotal_;
  MonitorElement* convEffEtTotal_;

  MonitorElement* convFakeRateEtaTwoTracks_;
  MonitorElement* convFakeRatePhiTwoTracks_;
  MonitorElement* convFakeRateRTwoTracks_;
  MonitorElement* convFakeRateZTwoTracks_;
  MonitorElement* convFakeRateEtTwoTracks_;

  MonitorElement* bkgRecoEffEta_;
  MonitorElement* bkgRecoEffPhi_;
  MonitorElement* bkgRecoEffEt_;

  MonitorElement* bkgDeadChEta_;
  MonitorElement* bkgDeadChPhi_;
  MonitorElement* bkgDeadChEt_;

  MonitorElement* convVsEt_[2];
};

#endif