#include <SiPixelPhase1CompareTrackSoA.cc>

Inheritance diagram for SiPixelPhase1CompareTrackSoA:

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

	SiPixelPhase1CompareTrackSoA (const edm::ParameterSet &)

	~SiPixelPhase1CompareTrackSoA () override=default

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

Private Attributes
const float	dr2cut_

MonitorElement *	hchi2_

MonitorElement *	hChi2VsEta_

MonitorElement *	hChi2VsPhi_

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 *	hz_

MonitorElement *	hzdiffMatched_

const pixelTrack::Quality	minQuality_

const edm::EDGetTokenT< PixelTrackHeterogeneous >	tokenSoATrackCPU_

const edm::EDGetTokenT< PixelTrackHeterogeneous >	tokenSoATrackGPU_

const std::string	topFolderName_

const bool	useQualityCut_

Detailed Description

Definition at line 65 of file SiPixelPhase1CompareTrackSoA.cc.

Constructor & Destructor Documentation

◆ SiPixelPhase1CompareTrackSoA()

SiPixelPhase1CompareTrackSoA::SiPixelPhase1CompareTrackSoA ( const edm::ParameterSet & iConfig )

explicit

Definition at line 115 of file SiPixelPhase1CompareTrackSoA.cc.

     : tokenSoATrackCPU_(consumes<PixelTrackHeterogeneous>(iConfig.getParameter<edm::InputTag>("pixelTrackSrcCPU"))),
       tokenSoATrackGPU_(consumes<PixelTrackHeterogeneous>(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")) {}

◆ ~SiPixelPhase1CompareTrackSoA()

SiPixelPhase1CompareTrackSoA::~SiPixelPhase1CompareTrackSoA ( )

overridedefault

Member Function Documentation

◆ analyze()

void SiPixelPhase1CompareTrackSoA::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 126 of file SiPixelPhase1CompareTrackSoA.cc.

References reco::deltaPhi(), reco::deltaR2(), dr2cut_, dqm::impl::MonitorElement::Fill(), hchi2_, heta_, hetadiffMatched_, hnHits_, hnLayers_, hnLooseAndAboveTracks_, hnLooseAndAboveTracks_matched_, hnTracks_, hphi_, hphi_z_tkAllCPU_, hphi_z_tkAllCPUMatched_, hphidiffMatched_, hpt_, hpt_eta_tkAllCPU_, hpt_eta_tkAllCPUMatched_, hptdiffMatched_, hptLogLog_, htip_, hz_, hzdiffMatched_, iEvent, minQuality_, notFound, or, MillePedeFileConverter_cfg::out, tokenSoATrackCPU_, tokenSoATrackGPU_, and useQualityCut_.

                                                                                               {
   const auto& tsoaHandleCPU = iEvent.getHandle(tokenSoATrackCPU_);
   const auto& tsoaHandleGPU = iEvent.getHandle(tokenSoATrackGPU_);
   if (not tsoaHandleCPU or not tsoaHandleGPU) {
     edm::LogWarning out("SiPixelPhase1CompareTrackSoA");
     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->get();
   auto const& tsoaGPU = *tsoaHandleGPU->get();
   auto maxTracksCPU = tsoaCPU.stride();  //this should be same for both?
   auto maxTracksGPU = tsoaGPU.stride();  //this should be same for both?
   auto const* qualityCPU = tsoaCPU.qualityData();
   auto const* qualityGPU = tsoaGPU.qualityData();
   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 (tsoaGPU.nHits(jt) == 0)
       break;  // this is a guard
     if (!(tsoaGPU.pt(jt) > 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) {
     if (tsoaCPU.nHits(it) == 0)
       break;  // this is a guard
     if (!(tsoaCPU.pt(it) > 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_;
     float etacpu = tsoaCPU.eta(it);
     float phicpu = tsoaCPU.phi(it);
     for (auto gid : looseTrkidxGPU) {
       float etagpu = tsoaGPU.eta(gid);
       float phigpu = tsoaGPU.phi(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, tsoaCPU.pt(it));  //all CPU tk
     hphi_z_tkAllCPU_->Fill(phicpu, tsoaCPU.zip(it));
     if (closestTkidx == notFound)
       continue;
     nLooseAndAboveTracksCPU_matchedGPU++;
 
     hchi2_->Fill(tsoaCPU.chi2(it), tsoaGPU.chi2(closestTkidx));
     hnHits_->Fill(tsoaCPU.nHits(it), tsoaGPU.nHits(closestTkidx));
     hnLayers_->Fill(tsoaCPU.nLayers(it), tsoaGPU.nLayers(closestTkidx));
     hpt_->Fill(tsoaCPU.pt(it), tsoaGPU.pt(closestTkidx));
     hptLogLog_->Fill(tsoaCPU.pt(it), tsoaGPU.pt(closestTkidx));
     heta_->Fill(etacpu, tsoaGPU.eta(closestTkidx));
     hphi_->Fill(phicpu, tsoaGPU.phi(closestTkidx));
     hz_->Fill(tsoaCPU.zip(it), tsoaGPU.zip(closestTkidx));
     htip_->Fill(tsoaCPU.tip(it), tsoaGPU.tip(closestTkidx));
     hptdiffMatched_->Fill(tsoaCPU.pt(it) - tsoaGPU.pt(closestTkidx));
     hetadiffMatched_->Fill(etacpu - tsoaGPU.eta(closestTkidx));
     hphidiffMatched_->Fill(reco::deltaPhi(phicpu, tsoaGPU.phi(closestTkidx)));
     hzdiffMatched_->Fill(tsoaCPU.zip(it) - tsoaGPU.zip(closestTkidx));
     hpt_eta_tkAllCPUMatched_->Fill(etacpu, tsoaCPU.pt(it));  //matched to gpu
     hphi_z_tkAllCPUMatched_->Fill(phicpu, tsoaCPU.zip(it));
   }
   hnTracks_->Fill(nTracksCPU, nTracksGPU);
   hnLooseAndAboveTracks_->Fill(nLooseAndAboveTracksCPU, nLooseAndAboveTracksGPU);
   hnLooseAndAboveTracks_matched_->Fill(nLooseAndAboveTracksCPU, nLooseAndAboveTracksCPU_matchedGPU);
 }

◆ bookHistograms()

void SiPixelPhase1CompareTrackSoA::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

override

Definition at line 225 of file SiPixelPhase1CompareTrackSoA.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2I(), dqm::implementation::NavigatorBase::cd(), hchi2_, heta_, hetadiffMatched_, hnHits_, hnLayers_, hnLooseAndAboveTracks_, hnLooseAndAboveTracks_matched_, hnTracks_, hphi_, hphi_z_tkAllCPU_, hphi_z_tkAllCPUMatched_, hphidiffMatched_, hpt_, hpt_eta_tkAllCPU_, hpt_eta_tkAllCPUMatched_, hptdiffMatched_, hptLogLog_, htip_, hz_, hzdiffMatched_, M_PI, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and topFolderName_.

                                                                                {
   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 either TH2I or 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.);
 
   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.);
   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);
   //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()

void SiPixelPhase1CompareTrackSoA::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 268 of file SiPixelPhase1CompareTrackSoA.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);
 }