DTT0CalibrationRMS.h

Go to the documentation of this file.

#ifndef CalibMuon_DTT0CalibrationRMS_H
#define CalibMuon_DTT0CalibrationRMS_H

#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "DataFormats/DTDigi/interface/DTDigiCollection.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"

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

class TFile;
class TH1I;
class TH1D;
class DTT0;

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

  ~DTT0CalibrationRMS() override;

  // Operations

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

  void endJob() override;

protected:
private:
  // Generate the histo name
  std::string getHistoName(const DTWireId& wId) const;
  std::string getHistoName(const DTLayerId& lId) const;

  // Debug flag
  bool debug;

  // The token used to retrieve digis from the event
  edm::EDGetTokenT<DTDigiCollection> digiToken;

  // The root file which contain the histos per layer
  TFile* theFile;
  // The root file which will contain the histos per wire (for the given layer)
  TFile* theOutputFile;

  //The event counter
  unsigned int nevents;
  //Number of events to be used for the t0 per layer histos
  unsigned int eventsForLayerT0;
  //Number of events to be used for the t0 reference per wire
  unsigned int eventsForWireT0;

  //Reject digis if far from TP peak
  unsigned int rejectDigiFromPeak;

  //Acceptance of t0 w.r.t. reference peak
  double tpPeakWidth;

  // Write relative t0's with respect to mean t0's in chamber
  bool correctByChamberMean_;

  //The wheels,sector to be calibrated (default All)
  std::string theCalibWheel;
  int selWheel;
  std::string theCalibSector;
  int selSector;

  // Map of the histos and graph by layer
  std::map<DTLayerId, TH1I*> theHistoLayerMap;
  //Histo with t0 mean per layer for all the sector
  TH1D* hT0SectorHisto;

  //Layer with histos for each wire
  std::vector<DTWireId> wireIdWithHistos;
  std::vector<std::string> cellsWithHistos;

  //Maps with t0, sigma, number of digi per wire
  std::map<DTWireId, double> theAbsoluteT0PerWire;
  std::map<DTWireId, double> theRelativeT0PerWire;
  std::map<DTWireId, double> theSigmaT0PerWire;
  std::map<DTWireId, int> nDigiPerWire;
  std::map<DTWireId, int> nDigiPerWire_ref;
  std::map<DTWireId, double> mK;
  std::map<DTWireId, double> mK_ref;
  std::map<DTWireId, double> qK;
  //Map with histo per wire for the chosen layer
  std::map<DTWireId, TH1I*> theHistoWireMap;
  std::map<DTWireId, TH1I*> theHistoWireMap_ref;
  //Map with mean and RMS of t0 per layer
  std::map<std::string, double> theT0LayerMap;
  std::map<std::string, double> theSigmaT0LayerMap;

  //DTGeometry used to loop on the SL in the endJob
  edm::ESHandle<DTGeometry> dtGeom;
  const edm::ESGetToken<DTGeometry, MuonGeometryRecord> dtGeomToken_;
};
#endif