CSCSegmentValidation.h

Go to the documentation of this file.

#ifndef CSCSegmentValidation_h
#define CSCSegmentValidation_h
 
#include "DQMServices/Core/interface/DQMStore.h"
#include "DataFormats/CSCRecHit/interface/CSCSegment.h"
#include "DataFormats/CSCRecHit/interface/CSCSegmentCollection.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "Geometry/CSCGeometry/interface/CSCLayer.h"
#include "Validation/MuonCSCDigis/interface/CSCBaseValidation.h"
 
class CSCSegmentValidation : public CSCBaseValidation {
public:
  CSCSegmentValidation(const edm::InputTag &inputTag, 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;
  static int whatChamberType(int detId);
 
  edm::EDGetTokenT<CSCSegmentCollection> segments_Token_;
 
  // 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 *theThetaPullPlots[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