DTT0CalibrationNew.h

Go to the documentation of this file.

#ifndef CalibMuon_DTT0CalibrationNew_H
#define CalibMuon_DTT0CalibrationNew_H

#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "FWCore/Framework/interface/ESHandle.h"

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

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

class DTT0CalibrationNew : public edm::EDAnalyzer {
public:
  DTT0CalibrationNew(const edm::ParameterSet& pset);

  virtual ~DTT0CalibrationNew();

  // Operations

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

  void endJob();


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
  std::string digiLabel;

  std::string dbLabel;

  // 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;

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

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

  //Size of timeBox 
  unsigned int timeBoxWidth; 

  //Digi's will be rejected if too far from TP peak 
  unsigned int rejectDigiFromPeak;

  unsigned int retryForLayerT0;

  //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;
  
  TSpectrum* spectrum;

  //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;
  //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> theRefT0ByChamber;

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