PhotonOfflineClient.h

Go to the documentation of this file.

#ifndef PhotonOfflineClient_H
#define PhotonOfflineClient_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"
 
// EgammaCoreTools
#include "RecoEcal/EgammaCoreTools/interface/PositionCalc.h"
 
#include "TFile.h"
#include "TH1.h"
#include "TH2.h"
#include "TTree.h"
#include "TVector3.h"
#include "TProfile.h"
//
 
#include "DQMServices/Core/interface/DQMEDHarvester.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 "RooGlobalFunc.h"
#include "RooRealVar.h"
#include "RooDataSet.h"
#include "RooBreitWigner.h"
#include "RooDataHist.h"
#include "RooFitResult.h"
 
#include <vector>
#include <string>
 
// forward declarations
class TFile;
class TH1F;
class TH2F;
class TProfile;
class TTree;
 
class PhotonOfflineClient : public DQMEDHarvester {
public:
  //
  explicit PhotonOfflineClient(const edm::ParameterSet& pset);
  ~PhotonOfflineClient() override;
 
  //  virtual void analyze(const edm::Event&, const edm::EventSetup&  ) ;
  // virtual void beginJob() ;
  //virtual void endJob() ;
  void dqmEndJob(DQMStore::IBooker&, DQMStore::IGetter&) override;
 
  //  virtual void endLuminosityBlock( const edm::LuminosityBlock& , const edm::EventSetup& ) ;
  //virtual void endRun(const edm::Run& , const edm::EventSetup& ) ;
  //virtual void runClient();
 
  virtual void runClient(DQMStore::IBooker& iBooker, DQMStore::IGetter& iGetter);
  MonitorElement* bookHisto(
      DQMStore::IBooker& iBooker, std::string histoName, std::string title, int bin, double min, double max);
 
  void book2DHistoVector(DQMStore::IBooker& iBooker,
                         std::vector<std::vector<MonitorElement*> >& vecOfHist,
                         std::string histoType,
                         std::string histoName,
                         std::string title,
                         int xbin,
                         double xmin,
                         double xmax,
                         int ybin = 1,
                         double ymin = 1,
                         double ymax = 2);
  void book3DHistoVector(DQMStore::IBooker& iBooker,
                         std::vector<std::vector<std::vector<MonitorElement*> > >& vecOfHist,
                         std::string histoType,
                         std::string histoName,
                         std::string title,
                         int xbin,
                         double xmin,
                         double xmax,
                         int ybin = 1,
                         double ymin = 1,
                         double ymax = 2);
 
  MonitorElement* retrieveHisto(DQMStore::IGetter& iGetter, std::string dir, std::string name);
 
private:
  std::string analyzerName_;
  MonitorElement* p_efficiencyVsEtaLoose_;
  MonitorElement* p_efficiencyVsEtLoose_;
  MonitorElement* p_efficiencyVsEtaTight_;
  MonitorElement* p_efficiencyVsEtTight_;
  MonitorElement* p_efficiencyVsEtaHLT_;
  MonitorElement* p_efficiencyVsEtHLT_;
 
  MonitorElement* p_convFractionVsEtaLoose_;
  MonitorElement* p_convFractionVsEtLoose_;
  MonitorElement* p_convFractionVsEtaTight_;
  MonitorElement* p_convFractionVsEtTight_;
 
  std::vector<std::vector<MonitorElement*> > p_convFractionVsEta_;
  std::vector<std::vector<std::vector<MonitorElement*> > > p_convFractionVsPhi_;
  std::vector<std::vector<MonitorElement*> > p_convFractionVsEt_;
 
  std::vector<std::vector<MonitorElement*> > p_badChannelsFractionVsEta_;
  std::vector<std::vector<MonitorElement*> > p_badChannelsFractionVsPhi_;
  std::vector<std::vector<MonitorElement*> > p_badChannelsFractionVsEt_;
 
  MonitorElement* p_vertexReconstructionEfficiencyVsEta_;
 
  void dividePlots(MonitorElement* dividend, MonitorElement* numerator, MonitorElement* denominator);
  void dividePlots(MonitorElement* dividend, MonitorElement* numerator, double denominator);
 
  int verbosity_;
 
  edm::ParameterSet parameters_;
 
  double cutStep_;
  int numberOfSteps_;
  bool minimalSetOfHistos_;
  bool excludeBkgHistos_;
 
  double etMin;
  double etMax;
  int etBin;
  double etaMin;
  double etaMax;
  int etaBin;
  double phiMin;
  double phiMax;
  int phiBin;
 
  bool standAlone_;
  bool batch_;
 
  std::string outputFileName_;
  std::string inputFileName_;
 
  std::stringstream currentFolder_;
 
  int histo_index_photons_;
  int histo_index_conversions_;
  int histo_index_efficiency_;
  int histo_index_invMass_;
 
  std::vector<std::string> types_;
  std::vector<std::string> parts_;
};
 
#endif