PatPhotonSimpleAnalyzer.h

Go to the documentation of this file.

#ifndef RecoEgamma_PhotonIdentification_PatPhotonSimpleAnalyzer_H
#define RecoEgamma_PhotonIdentification_PatPhotonSimpleAnalyzer_H

//

//

// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include <string>
#include "TH1.h"
#include "TTree.h"

class TFile;

//
// class declaration
//
class PatPhotonSimpleAnalyzer : public edm::one::EDAnalyzer<> {
public:
  explicit PatPhotonSimpleAnalyzer(const edm::ParameterSet&);
  ~PatPhotonSimpleAnalyzer() override;

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

private:
  std::string outputFile_;  // output file
  double minPhotonEt_;      // minimum photon Et
  double minPhotonAbsEta_;  // min and
  double maxPhotonAbsEta_;  // max abs(eta)
  double minPhotonR9_;      // minimum R9 = E(3x3)/E(SuperCluster)
  double maxPhotonHoverE_;  // maximum HCAL / ECAL
  bool createPhotonTTree_;  // Create a TTree of photon variables

  // Will be used for creating TTree of photons.
  // These names did not have to match those from a phtn->...
  // but do match for clarity.
  struct struct_recPhoton {
    float isolationEcalRecHit;
    float isolationHcalRecHit;
    float isolationSolidTrkCone;
    float isolationHollowTrkCone;
    float nTrkSolidCone;
    float nTrkHollowCone;
    float isEBGap;
    float isEEGap;
    float isEBEEGap;
    float r9;
    float et;
    float eta;
    float phi;
    float hadronicOverEm;
    float ecalIso;
    float hcalIso;
    float trackIso;
  };
  struct_recPhoton recPhoton;

  // root file to store histograms
  TFile* rootFile_;

  // data members for histograms to be filled

  // PhotonID Histograms
  TH1F* h_isoEcalRecHit_;
  TH1F* h_isoHcalRecHit_;
  TH1F* h_trk_pt_solid_;
  TH1F* h_trk_pt_hollow_;
  TH1F* h_ntrk_solid_;
  TH1F* h_ntrk_hollow_;
  TH1F* h_ebgap_;
  TH1F* h_eeGap_;
  TH1F* h_ebeeGap_;
  TH1F* h_r9_;

  // Photon Histograms
  TH1F* h_photonEt_;
  TH1F* h_photonEta_;
  TH1F* h_photonPhi_;
  TH1F* h_hadoverem_;

  // Photon's SuperCluster Histograms
  TH1F* h_photonScEt_;
  TH1F* h_photonScEta_;
  TH1F* h_photonScPhi_;
  TH1F* h_photonScEtaWidth_;

  // Composite or Other Histograms
  TH1F* h_photonInAnyGap_;
  TH1F* h_nPassingPho_;
  TH1F* h_nPho_;

  // TTree
  TTree* tree_PhotonAll_;
};
#endif