MuonDTDigis.h

Go to the documentation of this file.

#ifndef Validation_MuonDTDigis_h
#define Validation_MuonDTDigis_h

#include <FWCore/Framework/interface/Event.h>
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
#include "SimDataFormats/TrackingHit/interface/PSimHit.h"

#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "Geometry/DTGeometry/interface/DTLayer.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"

#include "DataFormats/DTDigi/interface/DTDigi.h"
#include "DataFormats/DTDigi/interface/DTDigiCollection.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "DataFormats/MuonDetId/interface/DTLayerId.h"

#include "SimMuon/DTDigitizer/test/Histograms.h"

#include <vector>

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

class TH1F;
class TFile;
class PSimHit;
class hDigis;

namespace edm
{
  class ParameterSet;
  class Event;
  class EventSetup;
}

class MuonDTDigis : public DQMEDAnalyzer{

 public:
  // Constructor
  explicit MuonDTDigis(const edm::ParameterSet& pset);

  // Destructor
  ~MuonDTDigis() override;

 protected:
  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
  // Analysis
  void analyze(const edm::Event & event, const edm::EventSetup& eventSetup) override;

  hDigis* WheelHistos(int wheel);

 private:
  typedef std::map<DTWireId, std::vector<const PSimHit*> > DTWireIdMap;

  edm::EDGetTokenT<edm::PSimHitContainer> SimHitToken_;
  edm::EDGetTokenT<DTDigiCollection> DigiToken_;

  // Switch for debug output
  bool verbose_;

  // Monitor elements
  MonitorElement* meDigiTimeBox_;
  MonitorElement* meDigiTimeBox_wheel2m_;
  MonitorElement* meDigiTimeBox_wheel1m_;
  MonitorElement* meDigiTimeBox_wheel0_;
  MonitorElement* meDigiTimeBox_wheel1p_;
  MonitorElement* meDigiTimeBox_wheel2p_;
  MonitorElement* meDigiEfficiency_;
  MonitorElement* meDigiEfficiencyMu_;
  MonitorElement* meDoubleDigi_;
  MonitorElement* meSimvsDigi_;
  MonitorElement* meWire_DoubleDigi_;

  MonitorElement* meMB1_sim_occup_;
  MonitorElement* meMB1_digi_occup_;
  MonitorElement* meMB2_sim_occup_;
  MonitorElement* meMB2_digi_occup_;
  MonitorElement* meMB3_sim_occup_;
  MonitorElement* meMB3_digi_occup_;
  MonitorElement* meMB4_sim_occup_;
  MonitorElement* meMB4_digi_occup_;

  std::vector<MonitorElement*> meDigiTimeBox_SL_;
  MonitorElement* meDigiHisto_;

  //  TH1F *DigiTimeBox;
  //  TFile *file_more_plots;

  std::unique_ptr<hDigis> hDigis_global;
  std::unique_ptr<hDigis> hDigis_W0;
  std::unique_ptr<hDigis> hDigis_W1;
  std::unique_ptr<hDigis> hDigis_W2;
  std::unique_ptr<hHits> hAllHits;

};

#endif