DTTTrigCalibration.h

Go to the documentation of this file.

#ifndef CalibMuon_DTTTrigCalibration_H
#define CalibMuon_DTTTrigCalibration_H

#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "DataFormats/MuonDetId/interface/DTSuperLayerId.h"
#include "DataFormats/MuonDetId/interface/DTLayerId.h"

#include <string>
#include <map>

namespace edm {
  class ParameterSet;
  class Event;
  class EventSetup;
}  // namespace edm

class TFile;
class TH1F;
class DTTimeBoxFitter;
class DTTTrigBaseSync;
class DTTtrig;
class DTStatusFlag;
class DTStatusFlagRcd;

class DTTTrigCalibration : public edm::one::EDAnalyzer<> {
public:
  DTTTrigCalibration(const edm::ParameterSet& pset);

  ~DTTTrigCalibration() override;

  // Operations

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

  void endJob() override;

protected:
private:
  // Generate the time box name
  std::string getTBoxName(const DTSuperLayerId& slId) const;
  // Generate the time box name
  std::string getOccupancyName(const DTLayerId& slId) const;

  // Print computed ttrig
  void dumpTTrigMap(const DTTtrig* tTrig) const;

  // Plot computed ttrig
  void plotTTrig(const DTTtrig* tTrig) const;

  // Debug flag
  bool debug;

  // The label used to retrieve digis from the event
  std::string digiLabel;

  // The TDC time-window
  int maxTDCCounts;
  //The maximum number of digis per layer
  int maxDigiPerLayer;

  // The file which will contain the time boxes
  TFile* theFile;

  // Map of the histograms by SL
  std::map<DTSuperLayerId, TH1F*> theHistoMap;
  std::map<DTLayerId, TH1F*> theOccupancyMap;

  // Switch for t0 subtraction
  bool doSubtractT0;
  // Switch for checking of noisy channels
  bool checkNoisyChannels;
  //card to switch on/off the DB writing
  bool findTMeanAndSigma;
  // the kfactor to be uploaded in the ttrig DB
  double kFactor;

  // The fitter
  std::unique_ptr<DTTimeBoxFitter> theFitter;
  // The module for t0 subtraction
  std::unique_ptr<DTTTrigBaseSync> theSync;  //FIXME: should be const
  edm::ESGetToken<DTStatusFlag, DTStatusFlagRcd> theStatusMapToken;
};
#endif