DTNoiseCalibration.h

Go to the documentation of this file.

#ifndef DTNoiseCalibration_H
#define DTNoiseCalibration_H

/*
 * \file DTNoiseCalibration.h
 *
 * \author G. Mila - INFN Torino
 *         A. Vilela Pereira - INFN Torino
 *
*/

#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "CondFormats/DataRecord/interface/DTTtrigRcd.h"

#include <string>
#include <vector>
#include <map>
#include <ctime>

class DTGeometry;
class DTChamberId;
class DTSuperLayerId;
class DTLayerId;
class DTWireId;
class DTTtrig;
class TFile;
class TH2F;
class TH1F;

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

  void beginJob() override;
  void beginRun(const edm::Run& run, const edm::EventSetup& setup) override;
  void analyze(const edm::Event& e, const edm::EventSetup& c) override;
  void endRun(const edm::Run& run, const edm::EventSetup& setup) override{};
  void endJob() override;

private:
  std::string getChannelName(const DTWireId&) const;
  // Get the name of the layer
  std::string getLayerName(const DTLayerId&) const;
  // Get the name of the superLayer
  std::string getSuperLayerName(const DTSuperLayerId&) const;
  // Get the name of the chamber
  std::string getChamberName(const DTChamberId&) const;

  edm::InputTag digiLabel_;
  bool useTimeWindow_;
  double triggerWidth_;
  int timeWindowOffset_;
  double maximumNoiseRate_;
  bool useAbsoluteRate_;

  /*bool fastAnalysis;
  int wh;
  int sect;*/

  bool readDB_;
  int defaultTtrig_;

  std::vector<DTWireId> wireIdWithHisto_;
  unsigned int lumiMax_;

  int nevents_;
  //int counter;
  time_t runBeginTime_;
  time_t runEndTime_;

  // Get the DT Geometry
  edm::ESHandle<DTGeometry> dtGeom_;
  const edm::ESGetToken<DTGeometry, MuonGeometryRecord> dtGeomToken_;

  // tTrig map
  edm::ESHandle<DTTtrig> tTrigMap_;
  const edm::ESGetToken<DTTtrig, DTTtrigRcd> ttrigToken_;

  TFile* rootFile_;
  // TDC digi distribution
  TH1F* hTDCTriggerWidth_;
  // Map of the occupancy histograms by layer
  std::map<DTLayerId, TH1F*> theHistoOccupancyMap_;
  // Map of occupancy by lumi by wire
  std::map<DTWireId, TH1F*> theHistoOccupancyVsLumiMap_;
  // Map of occupancy by lumi by chamber
  std::map<DTChamberId, TH1F*> chamberOccupancyVsLumiMap_;
  // Map of occupancy by time by chamber
  std::map<DTChamberId, TH1F*> chamberOccupancyVsTimeMap_;
  // Map of the histograms with the number of events per evt per wire
  //std::map<DTLayerId, TH2F*> theHistoEvtPerWireMap_;
  // Map of skipped histograms
  //std::map<DTLayerId, int> skippedPlot;
};
#endif