DTt0DBValidation.h

Go to the documentation of this file.

#ifndef DTt0DBValidation_H
#define DTt0DBValidation_H

#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "RecoMuon/TrackingTools/interface/MuonServiceProxy.h"

#include "DataFormats/MuonDetId/interface/DTLayerId.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "CondFormats/DataRecord/interface/DTT0Rcd.h"

#include <fstream>
#include <string>
#include <vector>

class DTT0;
class TFile;

class DTt0DBValidation : public edm::EDAnalyzer {
public:
  typedef dqm::legacy::MonitorElement MonitorElement;
  typedef dqm::legacy::DQMStore DQMStore;
  DTt0DBValidation(const edm::ParameterSet &pset);

  ~DTt0DBValidation() override;

  // Read the DTGeometry and the t0 DB
  void beginRun(const edm::Run &run, const edm::EventSetup &setup) override;
  void endRun(edm::Run const &, edm::EventSetup const &) override;
  void endJob() override;
  void analyze(const edm::Event &event, const edm::EventSetup &setup) override {}

private:
  void bookHistos(DTLayerId lId, int firstWire, int lastWire);
  void bookHistos(int wheel);

  DQMStore *dbe_;
  // Switch for verbosity
  std::string metname_;
  // The DB label
  edm::ESGetToken<DTT0, DTT0Rcd> labelDBRef_;
  edm::ESGetToken<DTT0, DTT0Rcd> labelDB_;

  // The file which will contain the difference plot
  bool outputMEsInRootFile_;
  std::string outputFileName_;

  std::string t0TestName_;

  // The DTGeometry
  edm::ESGetToken<DTGeometry, MuonGeometryRecord> muonGeomToken_;
  const DTGeometry *dtGeom;

  // The t0 map
  const DTT0 *tZeroRefMap_;
  const DTT0 *tZeroMap_;

  // Map of the t0 and sigma per wire
  std::map<DTWireId, std::vector<float>> t0RefMap_;
  std::map<DTWireId, std::vector<float>> t0Map_;

  // Map of the t0 difference histos per layer
  std::map<DTLayerId, MonitorElement *> t0DiffHistos_;

  // Summary histos
  std::map<int, MonitorElement *> wheelSummary_;
};
#endif