MCPhotonAnalyzer.h

Go to the documentation of this file.

#ifndef MCPhotonAnalyzer_H
#define MCPhotonAnalyzer_H
#include "RecoEgamma/EgammaMCTools/interface/PhotonMCTruthFinder.h"
 
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
 
#include <map>
#include <vector>
 
// forward declarations
class TFile;
class TH1F;
class TH2F;
class TProfile;
class TTree;
class SimVertex;
class SimTrack;
 
class MCPhotonAnalyzer : public edm::one::EDAnalyzer<> {
public:
  //
  explicit MCPhotonAnalyzer(const edm::ParameterSet&);
  ~MCPhotonAnalyzer() override;
 
  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void beginJob() override;
  void endJob() override;
 
private:
  float etaTransformation(float a, float b);
  float phiNormalization(float& a);
 
  //
  PhotonMCTruthFinder* thePhotonMCTruthFinder_;
 
  const TrackerGeometry* trackerGeom;
 
  std::string fOutputFileName_;
  TFile* fOutputFile_;
 
  int nEvt_;
  int nMatched_;
 
 
  double mcPhi_;
  double mcEta_;
 
  std::string HepMCLabel;
  std::string SimTkLabel;
  std::string SimVtxLabel;
  std::string SimHitLabel;
 
  // all photons
  TH1F* h_MCPhoE_;
  TH1F* h_MCPhoEta_;
  TH1F* h_MCPhoEta1_;
  TH1F* h_MCPhoEta2_;
  TH1F* h_MCPhoEta3_;
  TH1F* h_MCPhoEta4_;
  TH1F* h_MCPhoPhi_;
  // Conversion
  TH1F* h_MCConvPhoE_;
  TH1F* h_MCConvPhoEta_;
  TH1F* h_MCConvPhoPhi_;
  TH1F* h_MCConvPhoR_;
  TH1F* h_MCConvPhoREta1_;
  TH1F* h_MCConvPhoREta2_;
  TH1F* h_MCConvPhoREta3_;
  TH1F* h_MCConvPhoREta4_;
  TH1F* h_convFracEta1_;
  TH1F* h_convFracEta2_;
  TH1F* h_convFracEta3_;
  TH1F* h_convFracEta4_;
 
  TH1F* h_MCConvPhoTwoTracksE_;
  TH1F* h_MCConvPhoTwoTracksEta_;
  TH1F* h_MCConvPhoTwoTracksPhi_;
  TH1F* h_MCConvPhoTwoTracksR_;
  TH1F* h_MCConvPhoOneTrackE_;
  TH1F* h_MCConvPhoOneTrackEta_;
  TH1F* h_MCConvPhoOneTrackPhi_;
  TH1F* h_MCConvPhoOneTrackR_;
 
  TH1F* h_MCEleE_;
  TH1F* h_MCEleEta_;
  TH1F* h_MCElePhi_;
  TH1F* h_BremFrac_;
  TH1F* h_BremEnergy_;
  TH2F* h_EleEvsPhoE_;
  TH2F* h_bremEvsEleE_;
 
  TProfile* p_BremVsR_;
  TProfile* p_BremVsEta_;
 
  TProfile* p_BremVsConvR_;
  TProfile* p_BremVsConvEta_;
 
  TH2F* h_bremFracVsConvR_;
};
 
#endif