Phase2TrackerValidateDigi.h

Go to the documentation of this file.

#ifndef Phase2TrackerValidateDigi_h
#define Phase2TrackerValidateDigi_h

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"

#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "SimDataFormats/Track/interface/SimTrackContainer.h"
#include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
#include "SimDataFormats/Vertex/interface/SimVertexContainer.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DQMServices/Core/interface/DQMStore.h"

// DQM Histograming
class PixelDigiSimLink;
class SimTrack;
class SimHit;
class TrackerTopology;
class PixelDigi;
class Phase2TrackerDigi;
class TrackerGeometry;
class TrackerDigiGeometryRecord;
class TrackerTopologyRcd;
class Phase2TrackerValidateDigi : public DQMEDAnalyzer {
public:
  explicit Phase2TrackerValidateDigi(const edm::ParameterSet&);
  ~Phase2TrackerValidateDigi() override;
  void bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iRun, edm::EventSetup const& iSetup) override;
  void analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) override;
  void dqmBeginRun(const edm::Run& iRun, const edm::EventSetup& iSetup) override;
  std::string getHistoId(uint32_t det_id, bool flag);

  struct DigiMEs {
    MonitorElement* SimTrackPt;
    MonitorElement* SimTrackEta;
    MonitorElement* SimTrackPhi;
    MonitorElement* MatchedTrackPt;
    MonitorElement* MatchedTrackPhi;
    MonitorElement* MatchedTrackEta;
    MonitorElement* MissedHitTrackPt;
    MonitorElement* MissedHitTrackPhi;
    MonitorElement* MissedHitTrackEta;
    MonitorElement* MissedDigiTrackPt;
    MonitorElement* MissedDigiTrackPhi;
    MonitorElement* MissedDigiTrackEta;
    MonitorElement* MissedDigiSimHitElossP;
    MonitorElement* MissedDigiSimHitElossS;
    MonitorElement* MatchedSimHitElossP;
    MonitorElement* MatchedSimHitElossS;
    MonitorElement* SimHitDx;
    MonitorElement* SimHitDy;
    MonitorElement* SimHitDz;
    MonitorElement* BunchXTimeBin;
    MonitorElement* FractionOfOOTDigis;
    MonitorElement* MissedDigiLocalXposVsYPos;
    MonitorElement* MissedDigiTimeWindow;
    int nHits;
    int nDigis;
  };

private:
  MonitorElement* nSimulatedTracks;
  MonitorElement* nSimulatedTracksP;
  MonitorElement* nSimulatedTracksS;

  MonitorElement* nSimVertices;

  MonitorElement* SimulatedTrackPt;
  MonitorElement* SimulatedTrackEta;
  MonitorElement* SimulatedTrackPhi;

  MonitorElement* SimulatedTrackPtP;
  MonitorElement* SimulatedTrackEtaP;
  MonitorElement* SimulatedTrackPhiP;

  MonitorElement* SimulatedTrackPtS;
  MonitorElement* SimulatedTrackEtaS;
  MonitorElement* SimulatedTrackPhiS;

  MonitorElement* SimulatedXYPositionMap;
  MonitorElement* SimulatedRZPositionMap;

  MonitorElement* MatchedXYPositionMap;
  MonitorElement* MatchedRZPositionMap;

  MonitorElement* SimulatedTOFEtaMap;
  MonitorElement* SimulatedTOFPhiMap;
  MonitorElement* SimulatedTOFRMap;
  MonitorElement* SimulatedTOFZMap;

  MonitorElement* nSimHitsPerTrack;
  float etaCut_;
  float ptCut_;
  float tofUpperCut_;
  float tofLowerCut_;

  void bookLayerHistos(DQMStore::IBooker& ibooker, unsigned int det_id, bool flag);
  unsigned int getSimTrackId(const edm::DetSetVector<PixelDigiSimLink>* simLinks,
                             const DetId& detId,
                             unsigned int& channel);
  int matchedSimTrack(edm::Handle<edm::SimTrackContainer>& SimTk, unsigned int simTrkId);
  bool isPrimary(const SimTrack& simTrk, const PSimHit& simHit);

  void fillHistogram(MonitorElement* th1, MonitorElement* th2, MonitorElement* th3, float val, int primary);
  int fillSimHitInfo(const edm::Event& iEvent, const SimTrack simTrk);

  bool findOTDigi(unsigned int detid, unsigned int id);

  bool findITPixelDigi(unsigned int detid, unsigned int id);
  void fillOTBXInfo();
  void fillITPixelBXInfo();
  void fillHitsPerTrack();

  edm::ParameterSet config_;
  std::map<std::string, DigiMEs> layerMEs;

  bool pixelFlag_;
  std::string geomType_;

  edm::InputTag otDigiSrc_;
  edm::InputTag otDigiSimLinkSrc_;
  edm::InputTag itPixelDigiSrc_;
  edm::InputTag itPixelDigiSimLinkSrc_;
  std::vector<edm::InputTag> pSimHitSrc_;
  edm::InputTag simTrackSrc_;
  edm::InputTag simVertexSrc_;

  const edm::EDGetTokenT<edm::DetSetVector<Phase2TrackerDigi>> otDigiToken_;
  const edm::EDGetTokenT<edm::DetSetVector<PixelDigiSimLink>> otDigiSimLinkToken_;
  const edm::EDGetTokenT<edm::DetSetVector<PixelDigi>> itPixelDigiToken_;
  const edm::EDGetTokenT<edm::DetSetVector<PixelDigiSimLink>> itPixelDigiSimLinkToken_;
  const edm::EDGetTokenT<edm::SimTrackContainer> simTrackToken_;
  const edm::EDGetTokenT<edm::SimVertexContainer> simVertexToken_;
  const edm::DetSetVector<PixelDigiSimLink>* otSimLink_;
  const edm::DetSetVector<PixelDigiSimLink>* itSimLink_;
  const edm::DetSetVector<Phase2TrackerDigi>* otdigis_;
  const edm::DetSetVector<PixelDigi>* itdigis_;
  std::vector<edm::EDGetTokenT<edm::PSimHitContainer>> simHitTokens_;
  edm::Handle<edm::PSimHitContainer> simHits;
  edm::Handle<edm::SimTrackContainer> simTracks;
  edm::Handle<edm::SimVertexContainer> simVertices;
  const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> geomToken_;
  const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> topoToken_;
  const TrackerGeometry* tkGeom_ = nullptr;
  const TrackerTopology* tTopo_ = nullptr;

  const float GeVperElectron;  // 3.7E-09
  const float cval;            // cm/ns
};
#endif