DTCalibValidationFromMuons.h

Go to the documentation of this file.

#ifndef DTCalibValidationFromMuons_H
#define DTCalibValidationFromMuons_H

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "DataFormats/MuonDetId/interface/DTSuperLayerId.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include <FWCore/Framework/interface/EDAnalyzer.h>
#include <FWCore/Framework/interface/LuminosityBlock.h>

#include <DQMServices/Core/interface/DQMStore.h>
#include <DQMServices/Core/interface/MonitorElement.h>
#include <DQMServices/Core/interface/DQMEDAnalyzer.h>

#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "DataFormats/DTRecHit/interface/DTRecHitCollection.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment2DCollection.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"

#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "RecoMuon/TrackingTools/interface/MuonSegmentMatcher.h"


#include <FWCore/Framework/interface/ESHandle.h>

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


// To remove into CMSSW versions before 20X; removed in CMSSW 10_5_X onwards
//class DQMStore;
// To add into CMSSW versions before 20X
//class DaqMonitorBEInterface;

class MonitorElement;
class DTGeometry;
class DTChamber;

// FR class DTCalibValidationFromMuons: public edm::EDAnalyzer{
class DTCalibValidationFromMuons: public DQMEDAnalyzer{

 public:
  DTCalibValidationFromMuons(const edm::ParameterSet& pset);

  ~DTCalibValidationFromMuons() override;

  void dqmBeginRun(const edm::Run& r, const edm::EventSetup& c) override;

  // Operations
  void analyze(const edm::Event& event, const edm::EventSetup& setup) override;


 protected:
  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;


 private:

  // Switch for verbosity
  //bool debug;
  edm::ParameterSet parameters;
  int wrongSegment;
  int rightSegment;
  int nevent;
  // the geometry
  edm::ESHandle<DTGeometry> dtGeom;

  // Label of 4D segments in the event
  edm::EDGetTokenT<DTRecSegment4DCollection> segment4DToken_;

  // Label of muons in the event
  edm::EDGetTokenT<reco::MuonCollection> muonToken_;

  // Compute the distance from wire (cm) of a hits in a DTRecHit1DPair
  float recHitDistFromWire(const DTRecHit1DPair& hitPair, const DTLayer* layer);
  // Compute the distance from wire (cm) of a hits in a DTRecHit1D
  float recHitDistFromWire(const DTRecHit1D& recHit, const DTLayer* layer);
  // Compute the position with respect to the wire (cm) of a hits in a DTRecHit1DPair
  float recHitPosition(const DTRecHit1DPair& hitPair, const DTLayer* layer, const DTChamber* chamber, float segmPos, int sl);
  // Compute the position with respect to the wire (cm) of a hits in a DTRecHit1D
  float recHitPosition(const DTRecHit1D& recHit, const DTLayer* layer, const DTChamber* chamber, float segmPos, int sl);

  // Does the real job
  void compute(const DTGeometry *dtGeom, const DTRecSegment4D& segment);

  // Book a set of histograms for a give chamber
  void bookHistos(DTSuperLayerId slId, int step);
  // Fill a set of histograms for a give chamber
  void fillHistos(DTSuperLayerId slId,
          float distance,
          float residualOnDistance,
          float position,
          float residualOnPosition,
          int step);

  std::map<std::pair<DTSuperLayerId,int>, std::vector<MonitorElement*> > histosPerSL;

};
#endif


/* Local Variables: */
/* show-trailing-whitespace: t */
/* truncate-lines: t */
/* End: */