PhotonMonitor.h

Go to the documentation of this file.

#ifndef DQMOFFLINE_TRIGGER_PHOTON_H
#define DQMOFFLINE_TRIGGER_PHOTON_H

#include <string>
#include <vector>

#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include <DQMServices/Core/interface/DQMEDAnalyzer.h>

#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/Registry.h"

#include "CommonTools/Utils/interface/StringCutObjectSelector.h"

//DataFormats
#include "DataFormats/METReco/interface/PFMET.h"
#include "DataFormats/METReco/interface/PFMETCollection.h"

#include "DataFormats/JetReco/interface/PFJet.h"
#include "DataFormats/JetReco/interface/PFJetCollection.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/EgammaCandidates/interface/Photon.h"
#include "DataFormats/EgammaCandidates/interface/PhotonFwd.h"

class GenericTriggerEventFlag;

struct MEbinning {
  unsigned int nbins;
  double xmin;
  double xmax;
};

struct PhotonME {
  MonitorElement* numerator;
  MonitorElement* denominator;
};
//
// class declaration
//

class PhotonMonitor : public DQMEDAnalyzer 
{
public:
  PhotonMonitor( const edm::ParameterSet& );
  ~PhotonMonitor() override;
  static void fillDescriptions(edm::ConfigurationDescriptions & descriptions);
  static void fillHistoPSetDescription(edm::ParameterSetDescription & pset);
  static void fillHistoLSPSetDescription(edm::ParameterSetDescription & pset);
protected:

  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
  void bookME(DQMStore::IBooker &, PhotonME& me, const std::string& histname, const std::string& histtitle, unsigned int nbins, double xmin, double xmax);
  void bookME(DQMStore::IBooker &, PhotonME& me, const std::string& histname, const std::string& histtitle, const std::vector<double>& binningX);
  void bookME(DQMStore::IBooker &, PhotonME& me, const std::string& histname, const std::string& histtitle, unsigned int nbinsX, double xmin, double xmax, double ymin, double ymax);
  void bookME(DQMStore::IBooker &, PhotonME& me, const std::string& histname, const std::string& histtitle, unsigned int nbinsX, double xmin, double xmax, unsigned int nbinsY, double ymin, double ymax);
  void bookME(DQMStore::IBooker &, PhotonME& me, const std::string& histname, const std::string& histtitle, const std::vector<double>& binningX, const std::vector<double>& binningY);
  void setTitle(PhotonME& me, const std::string& titleX, const std::string& titleY);

  void analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) override;

private:
  static MEbinning getHistoPSet    (edm::ParameterSet const& pset);
  static MEbinning getHistoLSPSet  (edm::ParameterSet const& pset);

  std::string folderName_;
  std::string histoSuffix_;

  edm::EDGetTokenT<reco::PFMETCollection>       metToken_;
  edm::EDGetTokenT<reco::PFJetCollection>       jetToken_;
  edm::EDGetTokenT<reco::GsfElectronCollection> eleToken_;
  edm::EDGetTokenT<reco::PhotonCollection>     photonToken_;




  std::vector<double> photon_variable_binning_;
  std::vector<double> diphoton_mass_binning_;
  
  MEbinning           photon_binning_;
  MEbinning           ls_binning_;
  
  PhotonME subphotonEtaME_;
  PhotonME subphotonME_;
  PhotonME subphotonPhiME_;
  PhotonME subphotonME_variableBinning_;
  PhotonME subphotonEtaPhiME_;   
  PhotonME subphotonr9ME_;   
  PhotonME subphotonHoverEME_;   
  PhotonME diphotonMassME_;

  PhotonME photonEtaME_;
  PhotonME photonME_;
  PhotonME photonPhiME_;
  PhotonME photonME_variableBinning_;
  PhotonME photonVsLS_;
  PhotonME photonEtaPhiME_;   
  PhotonME photonr9ME_;   
  PhotonME photonHoverEME_;   

  double MAX_PHI1 = 3.2;
  unsigned int N_PHI1 = 64;
  const MEbinning phi_binning_1{
    N_PHI1, -MAX_PHI1, MAX_PHI1
      };


  double MAX_ETA = 1.4442;
  unsigned int N_ETA = 34;
  const MEbinning eta_binning_{
    N_ETA, -MAX_ETA, MAX_ETA
      };



  double MAX_r9 = 1;
  double MIN_r9 = 0;
  unsigned int N_r9 = 50;
  const MEbinning r9_binning_{
    N_r9, MIN_r9, MAX_r9
      };



  double MAX_hoe = 0.02;
  double MIN_hoe= 0;
  const MEbinning hoe_binning_{
    N_r9, MIN_hoe, MAX_hoe
      };





  //  GenericTriggerEventFlag* num_genTriggerEventFlag_;
  //GenericTriggerEventFlag* den_genTriggerEventFlag_;

  std::unique_ptr<GenericTriggerEventFlag> num_genTriggerEventFlag_ ;
  std::unique_ptr<GenericTriggerEventFlag> den_genTriggerEventFlag_ ;

  StringCutObjectSelector<reco::MET,true>         metSelection_;
  StringCutObjectSelector<reco::PFJet,true   >    jetSelection_;
  StringCutObjectSelector<reco::GsfElectron,true> eleSelection_;
  StringCutObjectSelector<reco::Photon,true> photonSelection_;
  unsigned int njets_;
  unsigned int nphotons_;
  unsigned int nelectrons_;

};

#endif // PhotonMonitor_H