EcalDisplaysByEvent.h

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:   EcalDisplaysByEvent
// Class:     EcalDisplaysByEvent
//
//
// Original Author:  Seth COOPER
//         Created:  Th Aug 28 5:46:22 CEST 2007
//
//

// system include files
#include <memory>
#include <vector>
#include <map>
#include <set>
#include <string>

// 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/EventSetup.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CommonTools/UtilAlgos/interface/TFileService.h"
#include "Geometry/EcalMapping/interface/EcalMappingRcd.h"

#include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalRawData/interface/EcalRawDataCollections.h"
#include "DataFormats/EcalRawData/interface/EcalDCCHeaderBlock.h"

#include "Geometry/EcalMapping/interface/EcalElectronicsMapping.h"
#include "Geometry/CaloTopology/interface/CaloTopology.h"
#include "Geometry/Records/interface/CaloTopologyRecord.h"

#include "CaloOnlineTools/EcalTools/interface/EcalFedMap.h"

#include "TFile.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TH3F.h"
#include "TGraph.h"
#include "TTree.h"
#include "TCanvas.h"

//
// class declaration
//

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

private:
  void beginRun(edm::Run const&, edm::EventSetup const&) override;
  void endRun(edm::Run const&, edm::EventSetup const&) override;
  void analyze(edm::Event const&, edm::EventSetup const&) override;
  void endJob() override;
  std::string intToString(int num);
  std::string floatToString(float num);
  void initHists(int);
  void initEvtByEvtHists(int naiveEvtNum_, int ievt);
  void deleteEvtByEvtHists();
  void initAllEventHistos();
  enum Ecal2DHistSubDetType { EB_FINE = 0, EB_COARSE = 1, EEM_FINE = 2, EEM_COARSE = 3, EEP_FINE = 4, EEP_COARSE = 5 };
  TH2F* init2DEcalHist(std::string histTypeName, int subDet);
  TH3F* init3DEcalHist(std::string histTypeName, int dubDet);
  TCanvas* init2DEcalCanvas(std::string canvasName);
  void selectHits(edm::Handle<EcalRecHitCollection> hits, int ievt);
  TGraph* selectDigi(DetId det, int ievt);
  int getEEIndex(EcalElectronicsId elecId);
  void makeHistos(edm::Handle<EBDigiCollection> ebDigis);
  void makeHistos(edm::Handle<EEDigiCollection> eeDigis);
  void makeHistos(edm::Handle<EcalRecHitCollection> hits);
  void drawHistos();
  void drawCanvas(TCanvas* canvas, TH1F* hist1, TH1F* hist2, TH1F* hist3);
  void drawCanvas(TCanvas* canvas, TH2F* hist1, TH2F* hist2, TH2F* hist3);
  void drawCanvas(TCanvas* canvas, TH3F* hist1, TH3F* hist2, TH3F* hist3);
  void drawTimingErrors(TProfile2D* profile);
  void drawEELines();

  // ----------member data ---------------------------

  const edm::InputTag EBRecHitCollection_;
  const edm::InputTag EERecHitCollection_;
  const edm::InputTag EBDigis_;
  const edm::InputTag EEDigis_;
  const edm::InputTag headerProducer_;

  edm::Handle<EBDigiCollection> EBdigisHandle;
  edm::Handle<EEDigiCollection> EEdigisHandle;

  const edm::EDGetTokenT<EcalRawDataCollection> rawDataToken_;
  const edm::EDGetTokenT<EcalRecHitCollection> ebRecHitToken_;
  const edm::EDGetTokenT<EcalRecHitCollection> eeRecHitToken_;
  const edm::EDGetTokenT<EBDigiCollection> ebDigiToken_;
  const edm::EDGetTokenT<EEDigiCollection> eeDigiToken_;

  const edm::ESGetToken<EcalElectronicsMapping, EcalMappingRcd> ecalMappingToken_;
  const edm::ESGetToken<CaloTopology, CaloTopologyRecord> topologyToken_;

  int runNum_;
  const int side_;
  const double threshold_;
  const double minTimingAmp_;
  const bool makeDigiGraphs_;
  const bool makeTimingHistos_;
  const bool makeEnergyHistos_;
  const bool makeOccupancyHistos_;
  const double histRangeMin_;
  const double histRangeMax_;
  const double minTimingEnergyEB_;
  const double minTimingEnergyEE_;

  std::set<EBDetId> listEBChannels;
  std::set<EEDetId> listEEChannels;

  int abscissa[10];
  int ordinate[10];

  static float gainRatio[3];
  static edm::Service<TFileService> fileService;

  std::vector<std::string>* names;
  std::vector<std::string>* histoCanvasNamesVector;
  std::vector<int> maskedChannels_;
  std::vector<int> maskedFEDs_;
  std::vector<int> seedCrys_;
  std::vector<std::string> maskedEBs_;
  std::map<int, TH1F*> FEDsAndTimingHists_;
  std::map<int, float> crysAndAmplitudesMap_;
  std::map<int, EcalDCCHeaderBlock> FEDsAndDCCHeaders_;
  std::map<std::string, int> seedFrequencyMap_;

  TH1F* allFedsTimingHist_;
  // For event-by-evet histos
  TH1F* timingEB_;
  TH1F* timingEEM_;
  TH1F* timingEEP_;
  TH1F* energyEB_;
  TH1F* energyEEM_;
  TH1F* energyEEP_;
  TH2F* energyMapEB_;
  TH2F* energyMapEBcoarse_;
  TH2F* energyMapEEM_;
  TH2F* energyMapEEMcoarse_;
  TH2F* energyMapEEP_;
  TH2F* energyMapEEPcoarse_;
  TH2F* recHitOccupancyEB_;
  TH2F* recHitOccupancyEBcoarse_;
  TH2F* recHitOccupancyEEM_;
  TH2F* recHitOccupancyEEMcoarse_;
  TH2F* recHitOccupancyEEP_;
  TH2F* recHitOccupancyEEPcoarse_;
  TH2F* digiOccupancyEB_;
  TH2F* digiOccupancyEBcoarse_;
  TH2F* digiOccupancyEEM_;
  TH2F* digiOccupancyEEMcoarse_;
  TH2F* digiOccupancyEEP_;
  TH2F* digiOccupancyEEPcoarse_;
  TH3F* timingMapEBCoarse_;
  TH3F* timingMapEEMCoarse_;
  TH3F* timingMapEEPCoarse_;
  TH3F* timingMapEB_;
  TH3F* timingMapEEM_;
  TH3F* timingMapEEP_;
  TCanvas* timingCanvas_;
  TCanvas* energyCanvas_;
  TCanvas* energyMapCanvas_;
  TCanvas* energyMapCoarseCanvas_;
  TCanvas* recHitOccupancyCanvas_;
  TCanvas* recHitOccupancyCoarseCanvas_;
  TCanvas* digiOccupancyCanvas_;
  TCanvas* digiOccupancyCoarseCanvas_;
  TCanvas* timingMapCoarseCanvas_;
  TCanvas* timingMapCanvas_;

  // For all-event hitos
  TH1F* timingEBAll_;
  TH1F* timingEEMAll_;
  TH1F* timingEEPAll_;
  TH1F* energyEBAll_;
  TH1F* energyEEMAll_;
  TH1F* energyEEPAll_;
  TH2F* energyMapEBAll_;
  TH2F* energyMapEBcoarseAll_;
  TH2F* energyMapEEMAll_;
  TH2F* energyMapEEMcoarseAll_;
  TH2F* energyMapEEPAll_;
  TH2F* energyMapEEPcoarseAll_;
  TH2F* recHitOccupancyEBAll_;
  TH2F* recHitOccupancyEBcoarseAll_;
  TH2F* recHitOccupancyEEMAll_;
  TH2F* recHitOccupancyEEMcoarseAll_;
  TH2F* recHitOccupancyEEPAll_;
  TH2F* recHitOccupancyEEPcoarseAll_;
  TH2F* digiOccupancyEBAll_;
  TH2F* digiOccupancyEBcoarseAll_;
  TH2F* digiOccupancyEEMAll_;
  TH2F* digiOccupancyEEMcoarseAll_;
  TH2F* digiOccupancyEEPAll_;
  TH2F* digiOccupancyEEPcoarseAll_;
  TH3F* timingMapEBCoarseAll_;
  TH3F* timingMapEEMCoarseAll_;
  TH3F* timingMapEEPCoarseAll_;
  TH3F* timingMapEBAll_;
  TH3F* timingMapEEMAll_;
  TH3F* timingMapEEPAll_;
  TCanvas* timingCanvasAll_;
  TCanvas* energyCanvasAll_;
  TCanvas* energyMapCanvasAll_;
  TCanvas* energyMapCoarseCanvasAll_;
  TCanvas* recHitOccupancyCanvasAll_;
  TCanvas* recHitOccupancyCoarseCanvasAll_;
  TCanvas* digiOccupancyCanvasAll_;
  TCanvas* digiOccupancyCoarseCanvasAll_;
  TCanvas* timingMapCoarseCanvasAll_;
  TCanvas* timingMapCanvasAll_;

  TTree* canvasNames_;
  TTree* histoCanvasNames_;
  EcalFedMap* fedMap_;
  const EcalElectronicsMapping* ecalElectronicsMap_;
  const CaloTopology* caloTopo_;

  int naiveEvtNum_;
};