EcalEBTrigPrimAnalyzer.h

Go to the documentation of this file.

#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Utilities/interface/ESGetToken.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "Geometry/CaloTopology/interface/EcalTrigTowerConstituentsMap.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"
#include "Geometry/Records/interface/EcalBarrelGeometryRecord.h"
#include "CalibCalorimetry/EcalTPGTools/interface/EcalTPGScale.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include <vector>
#include <string>
#include <TH1I.h>
#include <TFile.h>
#include <TTree.h>
#include <TH2F.h>
#include <TH1F.h>

//
// class declaration
//

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

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

private:
  int nEvents_;
  void init(const edm::EventSetup &);

  // for histos of nr of hits
  std::vector<std::string> ecal_parts_;
  TH1I *ecal_et_[2];
  TH1I *ecal_tt_[2];
  TH1I *ecal_fgvb_[2];
  TH1I *histEndc, *histBar;
  TFile *histfile_;
  TH1F *hAllTPperEvt_;
  TH1F *hTPperEvt_;
  TH2F *hTPvsRechit_;
  TH1F *hTPoverRechit_;
  TH1F *hDeltaEt_;
  TH1F *hAllRechitEt_;
  TH1F *hRechitEt_;
  TH1F *hTPEt_;
  TH1F *hRatioEt_;
  TH2F *hTPvsTow_ieta_;
  TH2F *hTPvsTow_eta_;
  TH2F *hTP_iphiVsieta_;
  TH2F *hRH_iphiVsieta_;
  TH2F *hTP_iphiVsieta_fullrange_;
  TH2F *hRH_iphiVsieta_fullrange_;

  TTree *tree_;

  int tpIphi_, tpIeta_, tpgADC_, ttf_, fg_;
  int rhIeta_, rhIphi_;
  float eRec_, tpgGeV_;

  //edm::InputTag label_;
  edm::EDGetTokenT<EcalEBTrigPrimDigiCollection> primToken_;
  edm::EDGetTokenT<EcalRecHitCollection> rechits_labelEB_;
  edm::EDGetTokenT<EBDigiCollection> tokenEBdigi_;
  bool recHits_;
  bool debug_;
  edm::ESHandle<EcalTrigTowerConstituentsMap> eTTmap_;
  edm::ESGetToken<EcalTrigTowerConstituentsMap, IdealGeometryRecord> eTTmapToken_;
  edm::ESGetToken<CaloSubdetectorGeometry, EcalBarrelGeometryRecord> barrelGeomToken_;
  EcalTPGScale::Tokens tokens_;

  int getIndex(const EBDigiCollection *, EcalTrigTowerDetId &id) { return id.hashedIndex(); }
};