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/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::EDAnalyzer {
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 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