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