SiPixelCompareTrackSoA< T > Class Template Reference

#include <SiPixelCompareTrackSoA.cc>

Inheritance diagram for SiPixelCompareTrackSoA< T >:

Public Types
using	PixelTrackSoA = TrackSoAHeterogeneousHost< T >

Public Member Functions
void	analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup) override
void	bookHistograms (DQMStore::IBooker &ibooker, edm::Run const &iRun, edm::EventSetup const &iSetup) override
	SiPixelCompareTrackSoA (const edm::ParameterSet &)
	~SiPixelCompareTrackSoA () override=default

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Attributes
const float	dr2cut_
MonitorElement *	hCharge_
MonitorElement *	hchi2_
MonitorElement *	hChi2VsEta_
MonitorElement *	hChi2VsPhi_
MonitorElement *	hCurvdiffMatched_
MonitorElement *	heta_
MonitorElement *	hetadiffMatched_
MonitorElement *	hnHits_
MonitorElement *	hnHitsVsEta_
MonitorElement *	hnHitsVsPhi_
MonitorElement *	hnLayers_
MonitorElement *	hnLayersVsEta_
MonitorElement *	hnLayersVsPhi_
MonitorElement *	hnLooseAndAboveTracks_
MonitorElement *	hnLooseAndAboveTracks_matched_
MonitorElement *	hnTracks_
MonitorElement *	hphi_
MonitorElement *	hphi_z_tkAllCPU_
MonitorElement *	hphi_z_tkAllCPUMatched_
MonitorElement *	hphidiffMatched_
MonitorElement *	hpt_
MonitorElement *	hpt_eta_tkAllCPU_
MonitorElement *	hpt_eta_tkAllCPUMatched_
MonitorElement *	hptdiffMatched_
MonitorElement *	hptLogLog_
MonitorElement *	hquality_
MonitorElement *	htip_
MonitorElement *	htipdiffMatched_
MonitorElement *	hz_
MonitorElement *	hzdiffMatched_
const pixelTrack::Quality	minQuality_
const edm::EDGetTokenT< PixelTrackSoA >	tokenSoATrackCPU_
const edm::EDGetTokenT< PixelTrackSoA >	tokenSoATrackGPU_
const std::string	topFolderName_
const bool	useQualityCut_

Detailed Description

template<typename T>
class SiPixelCompareTrackSoA< T >

Definition at line 67 of file SiPixelCompareTrackSoA.cc.

Member Typedef Documentation

PixelTrackSoA

template<typename T >

using SiPixelCompareTrackSoA< T >::PixelTrackSoA = TrackSoAHeterogeneousHost<T>

Definition at line 69 of file SiPixelCompareTrackSoA.cc.

Constructor & Destructor Documentation

SiPixelCompareTrackSoA()

template<typename T >

SiPixelCompareTrackSoA< T >::SiPixelCompareTrackSoA ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 124 of file SiPixelCompareTrackSoA.cc.

    : tokenSoATrackCPU_(consumes<PixelTrackSoA>(iConfig.getParameter<edm::InputTag>("pixelTrackSrcCPU"))),
      tokenSoATrackGPU_(consumes<PixelTrackSoA>(iConfig.getParameter<edm::InputTag>("pixelTrackSrcGPU"))),
      topFolderName_(iConfig.getParameter<std::string>("topFolderName")),
      useQualityCut_(iConfig.getParameter<bool>("useQualityCut")),
      minQuality_(pixelTrack::qualityByName(iConfig.getParameter<std::string>("minQuality"))),
      dr2cut_(iConfig.getParameter<double>("deltaR2cut")) {}

~SiPixelCompareTrackSoA()

template<typename T >

SiPixelCompareTrackSoA< T >::~SiPixelCompareTrackSoA ( )

overridedefault

Member Function Documentation

analyze()

template<typename T >

void SiPixelCompareTrackSoA< T >::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Definition at line 136 of file SiPixelCompareTrackSoA.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, reco::deltaPhi(), reco::deltaR2(), PV3DBase< T, PVType, FrameType >::eta(), iEvent, nHits, notFound, or, MillePedeFileConverter_cfg::out, qcdUeDQM_cfi::tip, and ComparisonHelper::zip().

                                                                                           {
  using helper = TracksUtilities<T>;
  const auto& tsoaHandleCPU = iEvent.getHandle(tokenSoATrackCPU_);
  const auto& tsoaHandleGPU = iEvent.getHandle(tokenSoATrackGPU_);
  if (not tsoaHandleCPU or not tsoaHandleGPU) {
    edm::LogWarning out("SiPixelCompareTrackSoA");
    if (not tsoaHandleCPU) {
      out << "reference (cpu) tracks not found; ";
    }
    if (not tsoaHandleGPU) {
      out << "target (gpu) tracks not found; ";
    }
    out << "the comparison will not run.";
    return;
  }

  auto const& tsoaCPU = *tsoaHandleCPU;
  auto const& tsoaGPU = *tsoaHandleGPU;
  auto maxTracksCPU = tsoaCPU.view().metadata().size();  //this should be same for both?
  auto maxTracksGPU = tsoaGPU.view().metadata().size();  //this should be same for both?
  auto const* qualityCPU = tsoaCPU.view().quality();
  auto const* qualityGPU = tsoaGPU.view().quality();
  int32_t nTracksCPU = 0;
  int32_t nTracksGPU = 0;
  int32_t nLooseAndAboveTracksCPU = 0;
  int32_t nLooseAndAboveTracksCPU_matchedGPU = 0;
  int32_t nLooseAndAboveTracksGPU = 0;

  //Loop over GPU tracks and store the indices of the loose tracks. Whats happens if useQualityCut_ is false?
  std::vector<int32_t> looseTrkidxGPU;
  for (int32_t jt = 0; jt < maxTracksGPU; ++jt) {
    if (helper::nHits(tsoaGPU.view(), jt) == 0)
      break;  // this is a guard
    if (!(tsoaGPU.view()[jt].pt() > 0.))
      continue;
    nTracksGPU++;
    if (useQualityCut_ && qualityGPU[jt] < minQuality_)
      continue;
    nLooseAndAboveTracksGPU++;
    looseTrkidxGPU.emplace_back(jt);
  }

  //Now loop over CPU tracks//nested loop for loose gPU tracks
  for (int32_t it = 0; it < maxTracksCPU; ++it) {
    int nHitsCPU = helper::nHits(tsoaCPU.view(), it);

    if (nHitsCPU == 0)
      break;  // this is a guard

    float ptCPU = tsoaCPU.view()[it].pt();
    float etaCPU = tsoaCPU.view()[it].eta();
    float phiCPU = helper::phi(tsoaCPU.view(), it);
    float zipCPU = helper::zip(tsoaCPU.view(), it);
    float tipCPU = helper::tip(tsoaCPU.view(), it);

    if (!(ptCPU > 0.))
      continue;
    nTracksCPU++;
    if (useQualityCut_ && qualityCPU[it] < minQuality_)
      continue;
    nLooseAndAboveTracksCPU++;
    //Now loop over loose GPU trk and find the closest in DeltaR//do we need pt cut?
    const int32_t notFound = -1;
    int32_t closestTkidx = notFound;
    float mindr2 = dr2cut_;

    for (auto gid : looseTrkidxGPU) {
      float etaGPU = tsoaGPU.view()[gid].eta();
      float phiGPU = helper::phi(tsoaGPU.view(), gid);
      float dr2 = reco::deltaR2(etaCPU, phiCPU, etaGPU, phiGPU);
      if (dr2 > dr2cut_)
        continue;  // this is arbitrary
      if (mindr2 > dr2) {
        mindr2 = dr2;
        closestTkidx = gid;
      }
    }

    hpt_eta_tkAllCPU_->Fill(etaCPU, ptCPU);  //all CPU tk
    hphi_z_tkAllCPU_->Fill(phiCPU, zipCPU);
    if (closestTkidx == notFound)
      continue;
    nLooseAndAboveTracksCPU_matchedGPU++;

    hchi2_->Fill(tsoaCPU.view()[it].chi2(), tsoaGPU.view()[closestTkidx].chi2());
    hCharge_->Fill(helper::charge(tsoaCPU.view(), it), helper::charge(tsoaGPU.view(), closestTkidx));
    hnHits_->Fill(helper::nHits(tsoaCPU.view(), it), helper::nHits(tsoaGPU.view(), closestTkidx));
    hnLayers_->Fill(tsoaCPU.view()[it].nLayers(), tsoaGPU.view()[closestTkidx].nLayers());
    hpt_->Fill(tsoaCPU.view()[it].pt(), tsoaGPU.view()[closestTkidx].pt());
    hptLogLog_->Fill(tsoaCPU.view()[it].pt(), tsoaGPU.view()[closestTkidx].pt());
    heta_->Fill(etaCPU, tsoaGPU.view()[closestTkidx].eta());
    hphi_->Fill(etaCPU, helper::phi(tsoaGPU.view(), closestTkidx));
    hz_->Fill(zipCPU, helper::zip(tsoaGPU.view(), closestTkidx));
    htip_->Fill(tipCPU, helper::tip(tsoaGPU.view(), closestTkidx));
    hptdiffMatched_->Fill(tsoaCPU.view()[it].pt() - tsoaGPU.view()[closestTkidx].pt());
    hCurvdiffMatched_->Fill((helper::charge(tsoaCPU.view(), it) / tsoaCPU.view()[it].pt()) -
                            (helper::charge(tsoaGPU.view(), closestTkidx) / tsoaGPU.view()[closestTkidx].pt()));
    hetadiffMatched_->Fill(etaCPU - tsoaGPU.view()[closestTkidx].eta());
    hphidiffMatched_->Fill(reco::deltaPhi(etaCPU, helper::phi(tsoaGPU.view(), closestTkidx)));
    hzdiffMatched_->Fill(zipCPU - helper::zip(tsoaGPU.view(), closestTkidx));
    htipdiffMatched_->Fill(tipCPU - helper::tip(tsoaGPU.view(), closestTkidx));
    hpt_eta_tkAllCPUMatched_->Fill(etaCPU, tsoaCPU.view()[it].pt());  //matched to gpu
    hphi_z_tkAllCPUMatched_->Fill(etaCPU, zipCPU);
  }
  hnTracks_->Fill(nTracksCPU, nTracksGPU);
  hnLooseAndAboveTracks_->Fill(nLooseAndAboveTracksCPU, nLooseAndAboveTracksGPU);
  hnLooseAndAboveTracks_matched_->Fill(nLooseAndAboveTracksCPU, nLooseAndAboveTracksCPU_matchedGPU);
}

bookHistograms()

template<typename T >

void SiPixelCompareTrackSoA< T >::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  iRun, edm::EventSetup const &  iSetup  )

override

Definition at line 249 of file SiPixelCompareTrackSoA.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2I(), dqm::implementation::NavigatorBase::cd(), M_PI, dqm::implementation::NavigatorBase::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                            {
  iBook.cd();
  iBook.setCurrentFolder(topFolderName_);

  // clang-format off
  std::string toRep = "Number of tracks";
  // FIXME: all the 2D correlation plots are quite heavy in terms of memory consumption, so a as soon as DQM supports THnSparse
  // these should be moved to a less resource consuming format
  hnTracks_ = iBook.book2I("nTracks", fmt::sprintf("%s per event; CPU; GPU",toRep), 501, -0.5, 500.5, 501, -0.5, 500.5);
  hnLooseAndAboveTracks_ = iBook.book2I("nLooseAndAboveTracks", fmt::sprintf("%s (quality #geq loose) per event; CPU; GPU",toRep), 501, -0.5, 500.5, 501, -0.5, 500.5);
  hnLooseAndAboveTracks_matched_ = iBook.book2I("nLooseAndAboveTracks_matched", fmt::sprintf("%s (quality #geq loose) per event; CPU; GPU",toRep), 501, -0.5, 500.5, 501, -0.5, 500.5);

  toRep = "Number of all RecHits per track (quality #geq loose)";
  hnHits_ = iBook.book2I("nRecHits", fmt::sprintf("%s;CPU;GPU",toRep), 15, -0.5, 14.5, 15, -0.5, 14.5);

  toRep = "Number of all layers per track (quality #geq loose)";
  hnLayers_ = iBook.book2I("nLayers", fmt::sprintf("%s;CPU;GPU",toRep), 15, -0.5, 14.5, 15, -0.5, 14.5);

  toRep = "Track (quality #geq loose) #chi^{2}/ndof";
  hchi2_ = iBook.book2I("nChi2ndof", fmt::sprintf("%s;CPU;GPU",toRep), 40, 0., 20., 40, 0., 20.);

  toRep = "Track (quality #geq loose) charge";
  hCharge_ = iBook.book2I("charge",fmt::sprintf("%s;CPU;GPU",toRep),3, -1.5, 1.5, 3, -1.5, 1.5);

  hpt_ = iBook.book2I("pt", "Track (quality #geq loose) p_{T} [GeV];CPU;GPU", 200, 0., 200., 200, 0., 200.);
  hptLogLog_ = make2DIfLog(iBook, true, true, "ptLogLog", "Track (quality #geq loose) p_{T} [GeV];CPU;GPU", 200, log10(0.5), log10(200.), 200, log10(0.5), log10(200.));
  heta_ = iBook.book2I("eta", "Track (quality #geq loose) #eta;CPU;GPU", 30, -3., 3., 30, -3., 3.);
  hphi_ = iBook.book2I("phi", "Track (quality #geq loose) #phi;CPU;GPU", 30, -M_PI, M_PI, 30, -M_PI, M_PI);
  hz_ = iBook.book2I("z", "Track (quality #geq loose) z [cm];CPU;GPU", 30, -30., 30., 30, -30., 30.);
  htip_ = iBook.book2I("tip", "Track (quality #geq loose) TIP [cm];CPU;GPU", 100, -0.5, 0.5, 100, -0.5, 0.5);
  //1D difference plots
  hptdiffMatched_ = iBook.book1D("ptdiffmatched", " p_{T} diff [GeV] between matched tracks; #Delta p_{T} [GeV]", 60, -30., 30.);
  hCurvdiffMatched_ = iBook.book1D("curvdiffmatched", "q/p_{T} diff [GeV] between matched tracks; #Delta q/p_{T} [GeV]", 60, -30., 30.);
  hetadiffMatched_ = iBook.book1D("etadiffmatched", " #eta diff between matched tracks; #Delta #eta", 160, -0.04 ,0.04);
  hphidiffMatched_ = iBook.book1D("phidiffmatched", " #phi diff between matched tracks; #Delta #phi",  160, -0.04 ,0.04);
  hzdiffMatched_ = iBook.book1D("zdiffmatched", " z diff between matched tracks; #Delta z [cm]", 300, -1.5, 1.5);
  htipdiffMatched_ = iBook.book1D("tipdiffmatched", " TIP diff between matched tracks; #Delta TIP [cm]", 300, -1.5, 1.5);
  //2D plots for eff
  hpt_eta_tkAllCPU_ = iBook.book2I("ptetatrkAllCPU", "Track (quality #geq loose) on CPU; #eta; p_{T} [GeV];", 30, -M_PI, M_PI, 200, 0., 200.);
  hpt_eta_tkAllCPUMatched_ = iBook.book2I("ptetatrkAllCPUmatched", "Track (quality #geq loose) on CPU matched to GPU track; #eta; p_{T} [GeV];", 30, -M_PI, M_PI, 200, 0., 200.);

  hphi_z_tkAllCPU_ = iBook.book2I("phiztrkAllCPU", "Track (quality #geq loose) on CPU; #phi; z [cm];",  30, -M_PI, M_PI, 30, -30., 30.);
  hphi_z_tkAllCPUMatched_ = iBook.book2I("phiztrkAllCPUmatched", "Track (quality #geq loose) on CPU; #phi; z [cm];", 30, -M_PI, M_PI, 30, -30., 30.);

}

fillDescriptions()

template<typename T >

void SiPixelCompareTrackSoA< T >::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 298 of file SiPixelCompareTrackSoA.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, HLT_2022v15_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                           {
  // monitorpixelTrackSoA
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("pixelTrackSrcCPU", edm::InputTag("pixelTracksSoA@cpu"));
  desc.add<edm::InputTag>("pixelTrackSrcGPU", edm::InputTag("pixelTracksSoA@cuda"));
  desc.add<std::string>("topFolderName", "SiPixelHeterogeneous/PixelTrackCompareGPUvsCPU");
  desc.add<bool>("useQualityCut", true);
  desc.add<std::string>("minQuality", "loose");
  desc.add<double>("deltaR2cut", 0.04);
  descriptions.addWithDefaultLabel(desc);
}

Member Data Documentation

dr2cut_

template<typename T >

const float SiPixelCompareTrackSoA< T >::dr2cut_

private

Definition at line 83 of file SiPixelCompareTrackSoA.cc.

hCharge_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hCharge_

private

Definition at line 93 of file SiPixelCompareTrackSoA.cc.

hchi2_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hchi2_

private

Definition at line 94 of file SiPixelCompareTrackSoA.cc.

hChi2VsEta_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hChi2VsEta_

private

Definition at line 96 of file SiPixelCompareTrackSoA.cc.

hChi2VsPhi_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hChi2VsPhi_

private

Definition at line 95 of file SiPixelCompareTrackSoA.cc.

hCurvdiffMatched_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hCurvdiffMatched_

private

Definition at line 106 of file SiPixelCompareTrackSoA.cc.

heta_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::heta_

private

Definition at line 99 of file SiPixelCompareTrackSoA.cc.

hetadiffMatched_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hetadiffMatched_

private

Definition at line 107 of file SiPixelCompareTrackSoA.cc.

hnHits_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hnHits_

private

Definition at line 87 of file SiPixelCompareTrackSoA.cc.

hnHitsVsEta_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hnHitsVsEta_

private

Definition at line 89 of file SiPixelCompareTrackSoA.cc.

hnHitsVsPhi_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hnHitsVsPhi_

private

Definition at line 88 of file SiPixelCompareTrackSoA.cc.

hnLayers_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hnLayers_

private

Definition at line 90 of file SiPixelCompareTrackSoA.cc.

hnLayersVsEta_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hnLayersVsEta_

private

Definition at line 92 of file SiPixelCompareTrackSoA.cc.

hnLayersVsPhi_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hnLayersVsPhi_

private

Definition at line 91 of file SiPixelCompareTrackSoA.cc.

hnLooseAndAboveTracks_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hnLooseAndAboveTracks_

private

Definition at line 85 of file SiPixelCompareTrackSoA.cc.

hnLooseAndAboveTracks_matched_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hnLooseAndAboveTracks_matched_

private

Definition at line 86 of file SiPixelCompareTrackSoA.cc.

hnTracks_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hnTracks_

private

Definition at line 84 of file SiPixelCompareTrackSoA.cc.

hphi_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hphi_

private

Definition at line 100 of file SiPixelCompareTrackSoA.cc.

hphi_z_tkAllCPU_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hphi_z_tkAllCPU_

private

Definition at line 115 of file SiPixelCompareTrackSoA.cc.

hphi_z_tkAllCPUMatched_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hphi_z_tkAllCPUMatched_

private

Definition at line 116 of file SiPixelCompareTrackSoA.cc.

hphidiffMatched_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hphidiffMatched_

private

Definition at line 108 of file SiPixelCompareTrackSoA.cc.

hpt_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hpt_

private

Definition at line 97 of file SiPixelCompareTrackSoA.cc.

hpt_eta_tkAllCPU_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hpt_eta_tkAllCPU_

private

Definition at line 113 of file SiPixelCompareTrackSoA.cc.

hpt_eta_tkAllCPUMatched_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hpt_eta_tkAllCPUMatched_

private

Definition at line 114 of file SiPixelCompareTrackSoA.cc.

hptdiffMatched_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hptdiffMatched_

private

Definition at line 105 of file SiPixelCompareTrackSoA.cc.

hptLogLog_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hptLogLog_

private

Definition at line 98 of file SiPixelCompareTrackSoA.cc.

hquality_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hquality_

private

Definition at line 103 of file SiPixelCompareTrackSoA.cc.

htip_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::htip_

private

Definition at line 102 of file SiPixelCompareTrackSoA.cc.

htipdiffMatched_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::htipdiffMatched_

private

Definition at line 110 of file SiPixelCompareTrackSoA.cc.

hz_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hz_

private

Definition at line 101 of file SiPixelCompareTrackSoA.cc.

hzdiffMatched_

template<typename T >

MonitorElement* SiPixelCompareTrackSoA< T >::hzdiffMatched_

private

Definition at line 109 of file SiPixelCompareTrackSoA.cc.

minQuality_

template<typename T >

const pixelTrack::Quality SiPixelCompareTrackSoA< T >::minQuality_

private

Definition at line 82 of file SiPixelCompareTrackSoA.cc.

tokenSoATrackCPU_

template<typename T >

const edm::EDGetTokenT<PixelTrackSoA> SiPixelCompareTrackSoA< T >::tokenSoATrackCPU_

private

Definition at line 78 of file SiPixelCompareTrackSoA.cc.

tokenSoATrackGPU_

template<typename T >

const edm::EDGetTokenT<PixelTrackSoA> SiPixelCompareTrackSoA< T >::tokenSoATrackGPU_

private

Definition at line 79 of file SiPixelCompareTrackSoA.cc.

topFolderName_

template<typename T >

const std::string SiPixelCompareTrackSoA< T >::topFolderName_

private

Definition at line 80 of file SiPixelCompareTrackSoA.cc.

useQualityCut_

template<typename T >

const bool SiPixelCompareTrackSoA< T >::useQualityCut_

private

Definition at line 81 of file SiPixelCompareTrackSoA.cc.