CSCSegmentValidation.h

Go to the documentation of this file.

#ifndef Validation_CSCRecHits_CSCSegmentValidation_h
#define Validation_CSCRecHits_CSCSegmentValidation_h

#include "Validation/MuonCSCDigis/interface/CSCBaseValidation.h"
#include "DataFormats/CSCRecHit/interface/CSCSegmentCollection.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"

class CSCSegmentValidation : public CSCBaseValidation {
public:
  CSCSegmentValidation(const edm::ParameterSet &, edm::ConsumesCollector &&iC);
  ~CSCSegmentValidation() override;
  void bookHistograms(DQMStore::IBooker &);
  void analyze(const edm::Event &, const edm::EventSetup &) override;

private:
  void plotResolution(const PSimHit &simHit, const CSCSegment &recHit, const CSCLayer *layer, int chamberType);

  bool hasSegment(int chamberId) const;

  edm::EDGetTokenT<CSCSegmentCollection> segments_Token_;
  edm::InputTag inputTag_;

  // map to count how many layers are hit.
  // First index is chamber detId
  // second is layers that have hits
  typedef std::map<int, std::vector<int>> ChamberHitMap;
  ChamberHitMap theLayerHitsPerChamber;
  void fillLayerHitsPerChamber();
  void fillEfficiencyPlots();
  // decide which SimHit to use for the direction
  const PSimHit *keyHit(int chamberId) const;

  typedef std::map<int, std::vector<CSCSegment>> ChamberSegmentMap;
  ChamberSegmentMap theChamberSegmentMap;
  // the number of hits in a chamber to make it a shower
  int theShowerThreshold;

  MonitorElement *theNPerEventPlot;
  MonitorElement *theNRecHitsPlot;
  MonitorElement *theNPerChamberTypePlot;
  MonitorElement *theRdPhiResolutionPlots[10];
  MonitorElement *theRdPhiPullPlots[10];
  MonitorElement *theThetaResolutionPlots[10];
  MonitorElement *thedXdZResolutionPlots[10];
  MonitorElement *thedXdZPullPlots[10];
  MonitorElement *thedYdZResolutionPlots[10];
  MonitorElement *thedYdZPullPlots[10];

  MonitorElement *theTypePlot4HitsNoShower;
  MonitorElement *theTypePlot4HitsNoShowerSeg;
  MonitorElement *theTypePlot4HitsShower;
  MonitorElement *theTypePlot4HitsShowerSeg;
  MonitorElement *theTypePlot5HitsNoShower;
  MonitorElement *theTypePlot5HitsNoShowerSeg;
  MonitorElement *theTypePlot5HitsShower;
  MonitorElement *theTypePlot5HitsShowerSeg;
  MonitorElement *theTypePlot6HitsNoShower;
  MonitorElement *theTypePlot6HitsNoShowerSeg;
  MonitorElement *theTypePlot6HitsShower;
  MonitorElement *theTypePlot6HitsShowerSeg;
};

#endif