EgammaObjects.h

Go to the documentation of this file.

#ifndef EgammaObjects_h
#define EgammaObjects_h

#include <string>

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

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

#include "DataFormats/EgammaCandidates/interface/Photon.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"

#include "TFile.h"
#include "TH1D.h"
#include "TH2D.h"

namespace edm {
  class HepMCProduct;
}

class EgammaObjects : public edm::EDAnalyzer {
public:
  explicit EgammaObjects(const edm::ParameterSet&);
  ~EgammaObjects() override;

  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void beginJob() override;
  void endJob() override;

private:
  TFile* rootFile_;

  edm::EDGetTokenT<edm::HepMCProduct> MCTruthCollectionT_;
  edm::EDGetTokenT<reco::GsfElectronCollection> RecoCollectionT_;

  int particleID;
  std::string particleString;

  int EtCut;

  double hist_min_Et_;
  double hist_max_Et_;
  int hist_bins_Et_;

  double hist_min_E_;
  double hist_max_E_;
  int hist_bins_E_;

  double hist_min_Eta_;
  double hist_max_Eta_;
  int hist_bins_Eta_;

  double hist_min_Phi_;
  double hist_max_Phi_;
  int hist_bins_Phi_;

  double hist_min_EtOverTruth_;
  double hist_max_EtOverTruth_;
  int hist_bins_EtOverTruth_;

  double hist_min_EOverTruth_;
  double hist_max_EOverTruth_;
  int hist_bins_EOverTruth_;

  double hist_min_EtaOverTruth_;
  double hist_max_EtaOverTruth_;
  int hist_bins_EtaOverTruth_;

  double hist_min_PhiOverTruth_;
  double hist_max_PhiOverTruth_;
  int hist_bins_PhiOverTruth_;

  double hist_min_deltaEta_;
  double hist_max_deltaEta_;
  int hist_bins_deltaEta_;

  double hist_min_deltaPhi_;
  double hist_max_deltaPhi_;
  int hist_bins_deltaPhi_;

  double hist_min_recoMass_;
  double hist_max_recoMass_;
  int hist_bins_recoMass_;

  TH1D* hist_Et_;
  TH1D* hist_EtOverTruth_;
  TH1D* hist_EtEfficiency_;
  TH1D* hist_EtNumRecoOverNumTrue_;
  TH1D* hist_EtOverTruthVsEt_;
  TH1D* hist_EtOverTruthVsE_;
  TH1D* hist_EtOverTruthVsEta_;
  TH1D* hist_EtOverTruthVsPhi_;
  TH1D* hist_resolutionEtVsEt_;
  TH1D* hist_resolutionEtVsE_;
  TH1D* hist_resolutionEtVsEta_;
  TH1D* hist_resolutionEtVsPhi_;

  TH1D* hist_E_;
  TH1D* hist_EOverTruth_;
  TH1D* hist_EEfficiency_;
  TH1D* hist_ENumRecoOverNumTrue_;
  TH1D* hist_EOverTruthVsEt_;
  TH1D* hist_EOverTruthVsE_;
  TH1D* hist_EOverTruthVsEta_;
  TH1D* hist_EOverTruthVsPhi_;
  TH1D* hist_resolutionEVsEt_;
  TH1D* hist_resolutionEVsE_;
  TH1D* hist_resolutionEVsEta_;
  TH1D* hist_resolutionEVsPhi_;

  TH1D* hist_Eta_;
  TH1D* hist_EtaOverTruth_;
  TH1D* hist_EtaEfficiency_;
  TH1D* hist_EtaNumRecoOverNumTrue_;
  TH1D* hist_deltaEtaVsEt_;
  TH1D* hist_deltaEtaVsE_;
  TH1D* hist_deltaEtaVsEta_;
  TH1D* hist_deltaEtaVsPhi_;
  TH1D* hist_resolutionEtaVsEt_;
  TH1D* hist_resolutionEtaVsE_;
  TH1D* hist_resolutionEtaVsEta_;
  TH1D* hist_resolutionEtaVsPhi_;

  TH1D* hist_Phi_;
  TH1D* hist_PhiOverTruth_;
  TH1D* hist_PhiEfficiency_;
  TH1D* hist_PhiNumRecoOverNumTrue_;
  TH1D* hist_deltaPhiVsEt_;
  TH1D* hist_deltaPhiVsE_;
  TH1D* hist_deltaPhiVsEta_;
  TH1D* hist_deltaPhiVsPhi_;
  TH1D* hist_resolutionPhiVsEt_;
  TH1D* hist_resolutionPhiVsE_;
  TH1D* hist_resolutionPhiVsEta_;
  TH1D* hist_resolutionPhiVsPhi_;

  TH1D* hist_All_recoMass_;
  TH1D* hist_BarrelOnly_recoMass_;
  TH1D* hist_EndcapOnly_recoMass_;
  TH1D* hist_Mixed_recoMass_;

  TH1D* hist_recoMass_withBackgroud_NoEtCut_;
  TH1D* hist_recoMass_withBackgroud_5EtCut_;
  TH1D* hist_recoMass_withBackgroud_10EtCut_;
  TH1D* hist_recoMass_withBackgroud_20EtCut_;

  TH2D* _TEMP_scatterPlot_EtOverTruthVsEt_;
  TH2D* _TEMP_scatterPlot_EtOverTruthVsE_;
  TH2D* _TEMP_scatterPlot_EtOverTruthVsEta_;
  TH2D* _TEMP_scatterPlot_EtOverTruthVsPhi_;

  TH2D* _TEMP_scatterPlot_EOverTruthVsEt_;
  TH2D* _TEMP_scatterPlot_EOverTruthVsE_;
  TH2D* _TEMP_scatterPlot_EOverTruthVsEta_;
  TH2D* _TEMP_scatterPlot_EOverTruthVsPhi_;

  TH2D* _TEMP_scatterPlot_deltaEtaVsEt_;
  TH2D* _TEMP_scatterPlot_deltaEtaVsE_;
  TH2D* _TEMP_scatterPlot_deltaEtaVsEta_;
  TH2D* _TEMP_scatterPlot_deltaEtaVsPhi_;

  TH2D* _TEMP_scatterPlot_deltaPhiVsEt_;
  TH2D* _TEMP_scatterPlot_deltaPhiVsE_;
  TH2D* _TEMP_scatterPlot_deltaPhiVsEta_;
  TH2D* _TEMP_scatterPlot_deltaPhiVsPhi_;

  void loadCMSSWObjects(const edm::ParameterSet& ps);
  void loadHistoParameters(const edm::ParameterSet& ps);

  void createBookedHistoObjects();
  void createTempHistoObjects();

  void analyzePhotons(const edm::Event&, const edm::EventSetup&);
  void analyzeElectrons(const edm::Event&, const edm::EventSetup&);

  void getEfficiencyHistosViaDividing();
  void getDeltaResHistosViaSlicing();
  void fitHistos();

  void applyLabels();
  void setDrawOptions();
  void saveHistos();

  double findRecoMass(const reco::Photon& pOne, const reco::Photon& pTwo);
  double findRecoMass(const reco::GsfElectron& eOne, const reco::GsfElectron& eTwo);

  float ecalEta(float EtaParticle, float Zvertex, float plane_Radius);
};
#endif