#include <SiPixelPhase1CompareVertexSoA.cc>

Inheritance diagram for SiPixelPhase1CompareVertexSoA:

Public Types
using	IndToEdm = std::vector< uint16_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

	SiPixelPhase1CompareVertexSoA (const edm::ParameterSet &)

	~SiPixelPhase1CompareVertexSoA () override=default

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

Private Attributes
const float	dzCut_

MonitorElement *	hchi2_

MonitorElement *	hchi2oNdof_

MonitorElement *	hntrks_

MonitorElement *	hnVertex_

MonitorElement *	hptv2_

MonitorElement *	hx_

MonitorElement *	hxdiff_

MonitorElement *	hy_

MonitorElement *	hydiff_

MonitorElement *	hz_

MonitorElement *	hzdiff_

const edm::EDGetTokenT< reco::BeamSpot >	tokenBeamSpot_

const edm::EDGetTokenT< ZVertexHeterogeneous >	tokenSoAVertexCPU_

const edm::EDGetTokenT< ZVertexHeterogeneous >	tokenSoAVertexGPU_

const std::string	topFolderName_

Detailed Description

Definition at line 24 of file SiPixelPhase1CompareVertexSoA.cc.

Member Typedef Documentation

◆ IndToEdm

using SiPixelPhase1CompareVertexSoA::IndToEdm = std::vector<uint16_t>

Definition at line 26 of file SiPixelPhase1CompareVertexSoA.cc.

Constructor & Destructor Documentation

◆ SiPixelPhase1CompareVertexSoA()

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

explicit

Definition at line 56 of file SiPixelPhase1CompareVertexSoA.cc.

     : tokenSoAVertexCPU_(consumes<ZVertexHeterogeneous>(iConfig.getParameter<edm::InputTag>("pixelVertexSrcCPU"))),
       tokenSoAVertexGPU_(consumes<ZVertexHeterogeneous>(iConfig.getParameter<edm::InputTag>("pixelVertexSrcGPU"))),
       tokenBeamSpot_(consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamSpotSrc"))),
       topFolderName_(iConfig.getParameter<std::string>("topFolderName")),
       dzCut_(iConfig.getParameter<double>("dzCut")) {}

◆ ~SiPixelPhase1CompareVertexSoA()

SiPixelPhase1CompareVertexSoA::~SiPixelPhase1CompareVertexSoA ( )

overridedefault

Member Function Documentation

◆ analyze()

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

override

Definition at line 66 of file SiPixelPhase1CompareVertexSoA.cc.

References funct::abs(), cms::cuda::bs, dxdz, dydz, dzCut_, dqm::impl::MonitorElement::Fill(), hchi2_, hchi2oNdof_, hntrks_, hnVertex_, hptv2_, hx_, hxdiff_, hy_, hydiff_, hz_, hzdiff_, iEvent, notFound, or, MillePedeFileConverter_cfg::out, tokenBeamSpot_, tokenSoAVertexCPU_, tokenSoAVertexGPU_, and HLTMuonOfflineAnalyzer_cfi::z0.

                                                                                                {
   const auto& vsoaHandleCPU = iEvent.getHandle(tokenSoAVertexCPU_);
   const auto& vsoaHandleGPU = iEvent.getHandle(tokenSoAVertexGPU_);
   if (not vsoaHandleCPU or not vsoaHandleGPU) {
     edm::LogWarning out("SiPixelPhase1CompareTrackSoA");
     if (not vsoaHandleCPU) {
       out << "reference (cpu) tracks not found; ";
     }
     if (not vsoaHandleGPU) {
       out << "target (gpu) tracks not found; ";
     }
     out << "the comparison will not run.";
     return;
   }
 
   auto const& vsoaCPU = *vsoaHandleCPU->get();
   int nVerticesCPU = vsoaCPU.nvFinal;
   auto const& vsoaGPU = *vsoaHandleGPU->get();
   int nVerticesGPU = vsoaGPU.nvFinal;
 
   auto bsHandle = iEvent.getHandle(tokenBeamSpot_);
   float x0 = 0., y0 = 0., z0 = 0., dxdz = 0., dydz = 0.;
   if (!bsHandle.isValid()) {
     edm::LogWarning("PixelVertexProducer") << "No beamspot found. returning vertexes with (0,0,Z) ";
   } else {
     const reco::BeamSpot& bs = *bsHandle;
     x0 = bs.x0();
     y0 = bs.y0();
     z0 = bs.z0();
     dxdz = bs.dxdz();
     dydz = bs.dydz();
   }
 
   for (int ivc = 0; ivc < nVerticesCPU; ivc++) {
     auto sic = vsoaCPU.sortInd[ivc];
     auto zc = vsoaCPU.zv[sic];
     auto xc = x0 + dxdz * zc;
     auto yc = y0 + dydz * zc;
     zc += z0;
 
     auto ndofCPU = vsoaCPU.ndof[sic];
     auto chi2CPU = vsoaCPU.chi2[sic];
 
     const int32_t notFound = -1;
     int32_t closestVtxidx = notFound;
     float mindz = dzCut_;
 
     for (int ivg = 0; ivg < nVerticesGPU; ivg++) {
       auto sig = vsoaGPU.sortInd[ivg];
       auto zgc = vsoaGPU.zv[sig] + z0;
       auto zDist = std::abs(zc - zgc);
       //insert some matching condition
       if (zDist > dzCut_)
         continue;
       if (mindz > zDist) {
         mindz = zDist;
         closestVtxidx = sig;
       }
     }
     if (closestVtxidx == notFound)
       continue;
 
     auto zg = vsoaGPU.zv[closestVtxidx];
     auto xg = x0 + dxdz * zg;
     auto yg = y0 + dydz * zg;
     zg += z0;
     auto ndofGPU = vsoaGPU.ndof[closestVtxidx];
     auto chi2GPU = vsoaGPU.chi2[closestVtxidx];
 
     hx_->Fill(xc - x0, xg - x0);
     hy_->Fill(yc - y0, yg - y0);
     hz_->Fill(zc, zg);
     hxdiff_->Fill(xc - xg);
     hydiff_->Fill(yc - yg);
     hzdiff_->Fill(zc - zg);
     hchi2_->Fill(chi2CPU, chi2GPU);
     hchi2oNdof_->Fill(chi2CPU / ndofCPU, chi2GPU / ndofGPU);
     hptv2_->Fill(vsoaCPU.ptv2[sic], vsoaGPU.ptv2[closestVtxidx]);
     hntrks_->Fill(ndofCPU + 1, ndofGPU + 1);
   }
   hnVertex_->Fill(nVerticesCPU, nVerticesGPU);
 }

◆ bookHistograms()

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

override

Definition at line 152 of file SiPixelPhase1CompareVertexSoA.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2I(), dqm::implementation::NavigatorBase::cd(), hchi2_, hchi2oNdof_, hntrks_, hnVertex_, hptv2_, hx_, hxdiff_, hy_, hydiff_, hz_, hzdiff_, dqm::implementation::NavigatorBase::setCurrentFolder(), and topFolderName_.

                                                                                 {
   ibooker.cd();
   ibooker.setCurrentFolder(topFolderName_);
 
   // 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
   hnVertex_ = ibooker.book2I("nVertex", "# of Vertex;CPU;GPU", 101, -0.5, 100.5, 101, -0.5, 100.5);
   hx_ = ibooker.book2I("vx", "Vertez x - Beamspot x;CPU;GPU", 50, -0.1, 0.1, 50, -0.1, 0.1);
   hy_ = ibooker.book2I("vy", "Vertez y - Beamspot y;CPU;GPU", 50, -0.1, 0.1, 50, -0.1, 0.1);
   hz_ = ibooker.book2I("vz", "Vertez z;CPU;GPU", 30, -30., 30., 30, -30., 30.);
   hchi2_ = ibooker.book2I("chi2", "Vertex chi-squared;CPU;GPU", 40, 0., 20., 40, 0., 20.);
   hchi2oNdof_ = ibooker.book2I("chi2oNdof", "Vertex chi-squared/Ndof;CPU;GPU", 40, 0., 20., 40, 0., 20.);
   hptv2_ = ibooker.book2I("ptsq", "Vertex p_T squared;CPU;GPU", 200, 0., 200., 200, 0., 200.);
   hntrks_ = ibooker.book2I("ntrk", "#tracks associated;CPU;GPU", 100, -0.5, 99.5, 100, -0.5, 99.5);
   hntrks_ = ibooker.book2I("ntrk", "#tracks associated;CPU;GPU", 100, -0.5, 99.5, 100, -0.5, 99.5);
   hxdiff_ = ibooker.book1D("vxdiff", ";Vertex x difference (CPU - GPU);#entries", 100, -0.001, 0.001);
   hydiff_ = ibooker.book1D("vydiff", ";Vertex y difference (CPU - GPU);#entries", 100, -0.001, 0.001);
   hzdiff_ = ibooker.book1D("vzdiff", ";Vertex z difference (CPU - GPU);#entries", 100, -2.5, 2.5);
 }

◆ fillDescriptions()

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

static

Definition at line 174 of file SiPixelPhase1CompareVertexSoA.cc.

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

                                                                                                {
   // monitorpixelVertexSoA
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("pixelVertexSrcCPU", edm::InputTag("pixelVerticesSoA@cpu"));
   desc.add<edm::InputTag>("pixelVertexSrcGPU", edm::InputTag("pixelVerticesSoA@cuda"));
   desc.add<edm::InputTag>("beamSpotSrc", edm::InputTag("offlineBeamSpot"));
   desc.add<std::string>("topFolderName", "SiPixelHeterogeneous/PixelVertexCompareSoAGPUvsCPU");
   desc.add<double>("dzCut", 1.);
   descriptions.addWithDefaultLabel(desc);
 }