MCPizeroAnalyzer.h

Go to the documentation of this file.

#ifndef MCPizeroAnalyzer_H
#define MCPizeroAnalyzer_H
#include "RecoEgamma/EgammaMCTools/interface/PizeroMCTruthFinder.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 MCPizeroAnalyzer : public edm::one::EDAnalyzer<> {
public:
  //
  explicit MCPizeroAnalyzer(const edm::ParameterSet&);
  ~MCPizeroAnalyzer() 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);

  //
  PizeroMCTruthFinder* thePizeroMCTruthFinder_;

  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;

  TH1F* h_MCPizE_;
  TH1F* h_MCPizEta_;
  TH1F* h_MCPizUnEta_;
  TH1F* h_MCPiz1ConEta_;
  TH1F* h_MCPiz2ConEta_;
  TH1F* h_MCPizPhi_;
  TH1F* h_MCPizMass1_;
  TH1F* h_MCPizMass2_;

  TH1F* h_MCEleE_;
  TH1F* h_MCEleEta_;
  TH1F* h_MCElePhi_;
  TH1F* h_BremFrac_;
  TH1F* h_BremEnergy_;

  TH2F* h_EleEvsPhoE_;

  TH1F* h_MCPhoE_;
  TH1F* h_MCPhoEta_;
  TH1F* h_MCPhoPhi_;
  TH1F* h_MCConvPhoE_;
  TH1F* h_MCConvPhoEta_;
  TH1F* h_MCConvPhoPhi_;
  TH1F* h_MCConvPhoR_;
};

#endif