DTT0Calibration.h

Go to the documentation of this file.

#ifndef CalibMuon_DTT0Calibration_H
#define CalibMuon_DTT0Calibration_H
 
//#include "FWCore/Framework/interface/EDAnalyzer.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 <string>
#include <vector>
#include <map>
 
#include "TFile.h"
#include "TH1I.h"
#include "TH1D.h"
#include "TSpectrum.h"
 
//class TFile;
//class TH1I;
//class TH1D;
//class TSpectrum;
class DTT0;
 
class DTT0Calibration : public edm::one::EDAnalyzer<> {
public:
  DTT0Calibration(const edm::ParameterSet& pset);
 
  ~DTT0Calibration() 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 label used to retrieve digis from the event
  edm::EDGetTokenT<DTDigiCollection> digiToken;
 
  // The root file which contain the histos per layer
  TFile theFile;
 
  //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;
 
  //Acceptance of t0 w.r.t. reference peak
  double tpPeakWidth;
 
  //Acceptance of t0 w.r.t. reference peak
  double tpPeakWidthPerLayer;
 
  //Digi's will be rejected if too far from TP peak
  unsigned int rejectDigiFromPeak;
 
  //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;
 
  // Histogram containing position of all peaks
  TH1D hLayerPeaks;
 
  TSpectrum spectrum;
 
  //Layer with histos for each wire
  std::vector<DTWireId> wireIdWithHistos;
  std::vector<DTLayerId> layerIdWithWireHistos;
 
  //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;
  //Map with mean and RMS of t0 per layer
  std::map<std::string, double> theT0LayerMap;
  std::map<std::string, double> theSigmaT0LayerMap;
  std::map<DTLayerId, double> theTPPeakMap;
  //Ref. t0 per chamber
  std::map<DTChamberId, double> theSumT0ByChamber;
  std::map<DTChamberId, int> theCountT0ByChamber;
  std::map<DTChamberId, double> theSigmaT0ByChamber;
  std::map<DTChamberId, double> theMeanT0ByChamber;
  std::map<DTChamberId, double> theRefT0ByChamber;
 
  //DTGeometry used to loop on the SL in the endJob
  edm::ESHandle<DTGeometry> dtGeom;
};
#endif