SiStripRecHitsValid.h

Go to the documentation of this file.

#ifndef SiStripRecHitsValid_h
#define SiStripRecHitsValid_h
 
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
//only mine
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
//DQM services for histogram
#include "DQMServices/Core/interface/DQMStore.h"
 
#include "FWCore/ServiceRegistry/interface/Service.h"
 
//--- for SimHit association
#include "SimDataFormats/TrackingHit/interface/PSimHit.h"
#include "SimTracker/TrackerHitAssociation/interface/TrackerHitAssociator.h"
 
#include "Geometry/CommonTopologies/interface/PixelTopology.h"
#include "Geometry/CommonTopologies/interface/StripTopology.h"
#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
#include "Geometry/CommonDetUnit/interface/GeomDetType.h"
#include "Geometry/CommonDetUnit/interface/GeomDet.h"
#include "Geometry/CommonDetUnit/interface/GluedGeomDet.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/TrackerNumberingBuilder/interface/GeometricDet.h"
#include "Geometry/CommonDetUnit/interface/PixelGeomDetType.h"
#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetType.h"
 
#include <string>
#include <DQMServices/Core/interface/DQMEDAnalyzer.h>
//For RecHit
#include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2DCollection.h"
#include "DataFormats/TrackerRecHit2D/interface/SiStripMatchedRecHit2DCollection.h"
#include "SimTracker/TrackerHitAssociation/interface/TrackerHitAssociator.h"
 
class SiStripDetCabling;
class SiStripDCSStatus;
 
class SiStripRecHitsValid : public DQMEDAnalyzer {
public:
  SiStripRecHitsValid(const edm::ParameterSet& conf);
 
  ~SiStripRecHitsValid() override;
 
  struct TotalMEs {  // MEs for total detector Level
    MonitorElement* meNumTotrphi;
    MonitorElement* meNumTotStereo;
    MonitorElement* meNumTotMatched;
  };
 
  struct SubDetMEs {  // MEs for Subdetector Level
    MonitorElement* meNumrphi;
    MonitorElement* meBunchrphi;
    MonitorElement* meEventrphi;
    MonitorElement* meNumStereo;
    MonitorElement* meBunchStereo;
    MonitorElement* meEventStereo;
    MonitorElement* meNumMatched;
    MonitorElement* meBunchMatched;
    MonitorElement* meEventMatched;
  };
 
  struct LayerMEs {  // MEs for Layer Level
    MonitorElement* meWclusrphi;
    MonitorElement* meAdcrphi;
    MonitorElement* mePosxrphi;
    MonitorElement* meResolxrphi;
    MonitorElement* meResrphi;
    MonitorElement* mePullLFrphi;
    MonitorElement* mePullMFrphi;
    MonitorElement* meChi2rphi;
    MonitorElement* meNsimHitrphi;
  };
 
  struct StereoAndMatchedMEs {  // MEs for Layer Level
    MonitorElement* meWclusStereo;
    MonitorElement* meAdcStereo;
    MonitorElement* mePosxStereo;
    MonitorElement* meResolxStereo;
    MonitorElement* meResStereo;
    MonitorElement* mePullLFStereo;
    MonitorElement* mePullMFStereo;
    MonitorElement* meChi2Stereo;
    MonitorElement* meNsimHitStereo;
    MonitorElement* mePosxMatched;
    MonitorElement* mePosyMatched;
    MonitorElement* meResolxMatched;
    MonitorElement* meResolyMatched;
    MonitorElement* meResxMatched;
    MonitorElement* meResyMatched;
    MonitorElement* meChi2Matched;
    MonitorElement* meNsimHitMatched;
  };
 
  struct RecHitProperties {
    float x;
    float y;
    //    float z;
    float resolxx;
    //    float resolxy;
    float resolyy;
    float resx;
    float resy;
    float pullMF;
    int clusiz;
    float cluchg;
    float chi2;
    int NsimHit;
    int bunch;
    int event;
  };
 
protected:
  void analyze(const edm::Event& e, const edm::EventSetup& c) override;
  void bookHistograms(DQMStore::IBooker& ibooker, const edm::Run& run, const edm::EventSetup& es) override;
 
private:
  TotalMEs totalMEs;
 
  bool switchNumTotrphi;
  bool switchNumTotStereo;
  bool switchNumTotMatched;
 
  bool switchNumrphi;
  bool switchBunchrphi;
  bool switchEventrphi;
  bool switchNumStereo;
  bool switchBunchStereo;
  bool switchEventStereo;
  bool switchNumMatched;
  bool switchBunchMatched;
  bool switchEventMatched;
 
  bool switchWclusrphi;
  bool switchAdcrphi;
  bool switchPosxrphi;
  bool switchResolxrphi;
  bool switchResrphi;
  bool switchPullLFrphi;
  bool switchPullMFrphi;
  bool switchChi2rphi;
  bool switchNsimHitrphi;
  bool switchWclusStereo;
  bool switchAdcStereo;
  bool switchPosxStereo;
  bool switchResolxStereo;
  bool switchResStereo;
  bool switchPullLFStereo;
  bool switchPullMFStereo;
  bool switchChi2Stereo;
  bool switchNsimHitStereo;
  bool switchPosxMatched;
  bool switchPosyMatched;
  bool switchResolxMatched;
  bool switchResolyMatched;
  bool switchResxMatched;
  bool switchResyMatched;
  bool switchChi2Matched;
  bool switchNsimHitMatched;
 
  std::string topFolderName_;
  std::vector<std::string> SubDetList_;
 
  std::map<std::string, LayerMEs> LayerMEsMap;
  std::map<std::string, StereoAndMatchedMEs> StereoAndMatchedMEsMap;
  std::map<std::string, SubDetMEs> SubDetMEsMap;
  std::map<std::string, std::vector<uint32_t> > LayerDetMap;
  std::map<std::string, std::vector<uint32_t> > StereoAndMatchedDetMap;
 
  edm::ESHandle<SiStripDetCabling> SiStripDetCabling_;
 
  std::pair<LocalPoint, LocalVector> projectHit(const PSimHit& hit,
                                                const StripGeomDetUnit* stripDet,
                                                const BoundPlane& plane);
  void createMEs(DQMStore::IBooker& ibooker, const edm::EventSetup& es);
  void createTotalMEs(DQMStore::IBooker& ibooker);
  void createLayerMEs(DQMStore::IBooker& ibooker, std::string label);
  void createSubDetMEs(DQMStore::IBooker& ibooker, std::string label);
  void createStereoAndMatchedMEs(DQMStore::IBooker& ibooker, std::string label);
 
  MonitorElement* bookME1D(DQMStore::IBooker& ibooker,
                           const char* ParameterSetLabel,
                           const char* HistoName,
                           const char* HistoTitle);
 
  inline void fillME(MonitorElement* ME, float value1) {
    if (ME != nullptr)
      ME->Fill(value1);
  }
  inline void fillME(MonitorElement* ME, float value1, float value2) {
    if (ME != nullptr)
      ME->Fill(value1, value2);
  }
  inline void fillME(MonitorElement* ME, float value1, float value2, float value3) {
    if (ME != nullptr)
      ME->Fill(value1, value2, value3);
  }
  inline void fillME(MonitorElement* ME, float value1, float value2, float value3, float value4) {
    if (ME != nullptr)
      ME->Fill(value1, value2, value3, value4);
  }
 
  edm::ParameterSet conf_;
  TrackerHitAssociator::Config trackerHitAssociatorConfig_;
  unsigned long long m_cacheID_;
  //const StripTopology* topol;
 
  /* static const int MAXHIT = 1000; */
 
  //  std::vector<RecHitProperties> rechitrphi;
  //  std::vector<RecHitProperties> rechitstereo;
  //  std::vector<RecHitProperties> rechitmatched;
  RecHitProperties rechitpro;
 
  void rechitanalysis(SiStripRecHit2D const rechit, const StripTopology& topol, TrackerHitAssociator& associate);
  void rechitanalysis_matched(SiStripMatchedRecHit2D const rechit,
                              const GluedGeomDet* gluedDet,
                              TrackerHitAssociator& associate);
 
  edm::EDGetTokenT<SiStripMatchedRecHit2DCollection> matchedRecHitsToken_;
  edm::EDGetTokenT<SiStripRecHit2DCollection> rphiRecHitsToken_;
  edm::EDGetTokenT<SiStripRecHit2DCollection> stereoRecHitsToken_;
};
 
#endif