#include <SiPixelCompareRecHitsSoA.cc>

Inheritance diagram for SiPixelCompareRecHitsSoA< T >:

Public Types
using	HitSoA = TrackingRecHitSoAView< T >

using	HitsOnHost = TrackingRecHitSoAHost< 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

void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override

	SiPixelCompareRecHitsSoA (const edm::ParameterSet &)

	~SiPixelCompareRecHitsSoA () override=default

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

Private Attributes
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomToken_

MonitorElement *	hBchargeDiff_

MonitorElement *	hBchargeL_ [4]

MonitorElement *	hBposXDiff_

MonitorElement *	hBposxL_ [4]

MonitorElement *	hBposYDiff_

MonitorElement *	hBposyL_ [4]

MonitorElement *	hBsizeXDiff_

MonitorElement *	hBsizexL_ [4]

MonitorElement *	hBsizeYDiff_

MonitorElement *	hBsizeyL_ [4]

MonitorElement *	hFchargeD_ [2][12]

MonitorElement *	hFchargeDiff_

MonitorElement *	hFposxD_ [2][12]

MonitorElement *	hFposXDiff_

MonitorElement *	hFposyD_ [2][12]

MonitorElement *	hFposYDiff_

MonitorElement *	hFsizexD_ [2][12]

MonitorElement *	hFsizeXDiff_

MonitorElement *	hFsizeyD_ [2][12]

MonitorElement *	hFsizeYDiff_

MonitorElement *	hnHits_

const float	mind2cut_

const TrackerGeometry *	tkGeom_ = nullptr

const edm::EDGetTokenT< HitsOnHost >	tokenSoAHitsCPU_

const edm::EDGetTokenT< HitsOnHost >	tokenSoAHitsGPU_

const std::string	topFolderName_

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	topoToken_

const TrackerTopology *	tTopo_ = nullptr

Static Private Attributes
static constexpr uint32_t	invalidHit_ = std::numeric_limits<uint32_t>::max()

static constexpr float	micron_ = 10000.

Detailed Description

template<typename T>
class SiPixelCompareRecHitsSoA< T >

Definition at line 32 of file SiPixelCompareRecHitsSoA.cc.

Member Typedef Documentation

◆ HitSoA

template<typename T >

using SiPixelCompareRecHitsSoA< T >::HitSoA = TrackingRecHitSoAView<T>

Definition at line 34 of file SiPixelCompareRecHitsSoA.cc.

◆ HitsOnHost

template<typename T >

using SiPixelCompareRecHitsSoA< T >::HitsOnHost = TrackingRecHitSoAHost<T>

Definition at line 35 of file SiPixelCompareRecHitsSoA.cc.

Constructor & Destructor Documentation

◆ SiPixelCompareRecHitsSoA()

template<typename T >

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

explicit

Definition at line 83 of file SiPixelCompareRecHitsSoA.cc.

     : geomToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>()),
       topoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd, edm::Transition::BeginRun>()),
       tokenSoAHitsCPU_(consumes(iConfig.getParameter<edm::InputTag>("pixelHitsSrcCPU"))),
       tokenSoAHitsGPU_(consumes(iConfig.getParameter<edm::InputTag>("pixelHitsSrcGPU"))),
       topFolderName_(iConfig.getParameter<std::string>("topFolderName")),
       mind2cut_(iConfig.getParameter<double>("minD2cut")) {}

◆ ~SiPixelCompareRecHitsSoA()

template<typename T >

SiPixelCompareRecHitsSoA< T >::~SiPixelCompareRecHitsSoA ( )

overridedefault

Member Function Documentation

◆ analyze()

template<typename T >

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

override

Definition at line 103 of file SiPixelCompareRecHitsSoA.cc.

References funct::abs(), reco::ceil(), HLT_2023v12_cff::distance, PVValHelper::dx, PVValHelper::dy, mps_fire::i, iEvent, dqmiolumiharvest::j, or, MillePedeFileConverter_cfg::out, PixelSubdetector::PixelBarrel, and PixelSubdetector::PixelEndcap.

                                                                                              {
   const auto& rhsoaHandleCPU = iEvent.getHandle(tokenSoAHitsCPU_);
   const auto& rhsoaHandleGPU = iEvent.getHandle(tokenSoAHitsGPU_);
   if (not rhsoaHandleCPU or not rhsoaHandleGPU) {
     edm::LogWarning out("SiPixelCompareRecHitSoA");
     if (not rhsoaHandleCPU) {
       out << "reference (CPU) rechits not found; ";
     }
     if (not rhsoaHandleGPU) {
       out << "target (GPU) rechits not found; ";
     }
     out << "the comparison will not run.";
     return;
   }
 
   auto const& rhsoaCPU = *rhsoaHandleCPU;
   auto const& rhsoaGPU = *rhsoaHandleGPU;
 
   auto const& soa2dCPU = rhsoaCPU.const_view();
   auto const& soa2dGPU = rhsoaGPU.const_view();
 
   uint32_t nHitsCPU = soa2dCPU.nHits();
   uint32_t nHitsGPU = soa2dGPU.nHits();
 
   hnHits_->Fill(nHitsCPU, nHitsGPU);
   auto detIds = tkGeom_->detUnitIds();
   for (uint32_t i = 0; i < nHitsCPU; i++) {
     float minD = mind2cut_;
     uint32_t matchedHit = invalidHit_;
     uint16_t indCPU = soa2dCPU[i].detectorIndex();
     float xLocalCPU = soa2dCPU[i].xLocal();
     float yLocalCPU = soa2dCPU[i].yLocal();
     for (uint32_t j = 0; j < nHitsGPU; j++) {
       if (soa2dGPU.detectorIndex(j) == indCPU) {
         float dx = xLocalCPU - soa2dGPU[j].xLocal();
         float dy = yLocalCPU - soa2dGPU[j].yLocal();
         float distance = dx * dx + dy * dy;
         if (distance < minD) {
           minD = distance;
           matchedHit = j;
         }
       }
     }
     DetId id = detIds[indCPU];
     uint32_t chargeCPU = soa2dCPU[i].chargeAndStatus().charge;
     int16_t sizeXCPU = std::ceil(float(std::abs(soa2dCPU[i].clusterSizeX()) / 8.));
     int16_t sizeYCPU = std::ceil(float(std::abs(soa2dCPU[i].clusterSizeY()) / 8.));
     uint32_t chargeGPU = 0;
     int16_t sizeXGPU = -99;
     int16_t sizeYGPU = -99;
     float xLocalGPU = -999.;
     float yLocalGPU = -999.;
     if (matchedHit != invalidHit_) {
       chargeGPU = soa2dGPU[matchedHit].chargeAndStatus().charge;
       sizeXGPU = std::ceil(float(std::abs(soa2dGPU[matchedHit].clusterSizeX()) / 8.));
       sizeYGPU = std::ceil(float(std::abs(soa2dGPU[matchedHit].clusterSizeY()) / 8.));
       xLocalGPU = soa2dGPU[matchedHit].xLocal();
       yLocalGPU = soa2dGPU[matchedHit].yLocal();
     }
     switch (id.subdetId()) {
       case PixelSubdetector::PixelBarrel:
         hBchargeL_[tTopo_->pxbLayer(id) - 1]->Fill(chargeCPU, chargeGPU);
         hBsizexL_[tTopo_->pxbLayer(id) - 1]->Fill(sizeXCPU, sizeXGPU);
         hBsizeyL_[tTopo_->pxbLayer(id) - 1]->Fill(sizeYCPU, sizeYGPU);
         hBposxL_[tTopo_->pxbLayer(id) - 1]->Fill(xLocalCPU, xLocalGPU);
         hBposyL_[tTopo_->pxbLayer(id) - 1]->Fill(yLocalCPU, yLocalGPU);
         hBchargeDiff_->Fill(chargeCPU - chargeGPU);
         hBsizeXDiff_->Fill(sizeXCPU - sizeXGPU);
         hBsizeYDiff_->Fill(sizeYCPU - sizeYGPU);
         hBposXDiff_->Fill(micron_ * (xLocalCPU - xLocalGPU));
         hBposYDiff_->Fill(micron_ * (yLocalCPU - yLocalGPU));
         break;
       case PixelSubdetector::PixelEndcap:
         hFchargeD_[tTopo_->pxfSide(id) - 1][tTopo_->pxfDisk(id) - 1]->Fill(chargeCPU, chargeGPU);
         hFsizexD_[tTopo_->pxfSide(id) - 1][tTopo_->pxfDisk(id) - 1]->Fill(sizeXCPU, sizeXGPU);
         hFsizeyD_[tTopo_->pxfSide(id) - 1][tTopo_->pxfDisk(id) - 1]->Fill(sizeYCPU, sizeYGPU);
         hFposxD_[tTopo_->pxfSide(id) - 1][tTopo_->pxfDisk(id) - 1]->Fill(xLocalCPU, xLocalGPU);
         hFposyD_[tTopo_->pxfSide(id) - 1][tTopo_->pxfDisk(id) - 1]->Fill(yLocalCPU, yLocalGPU);
         hFchargeDiff_->Fill(chargeCPU - chargeGPU);
         hFsizeXDiff_->Fill(sizeXCPU - sizeXGPU);
         hFsizeYDiff_->Fill(sizeYCPU - sizeYGPU);
         hFposXDiff_->Fill(micron_ * (xLocalCPU - xLocalGPU));
         hFposYDiff_->Fill(micron_ * (yLocalCPU - yLocalGPU));
         break;
     }
   }
 }

◆ bookHistograms()

template<typename T >

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

override

Definition at line 195 of file SiPixelCompareRecHitsSoA.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::book2I(), dqm::implementation::NavigatorBase::cd(), l1ctLayer2EG_cff::id, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, dqm::implementation::NavigatorBase::setCurrentFolder(), and Validation_hcalonly_cfi::sign.

                                                                               {
   iBook.cd();
   iBook.setCurrentFolder(topFolderName_);
 
   // clang-format off
   //Global
   hnHits_ = iBook.book2I("nHits", "CPUvsGPU RecHits per event;#CPU RecHits;#GPU RecHits", 200, 0, 5000,200, 0, 5000);
   //Barrel Layer
   for(unsigned int il=0;il<tkGeom_->numberOfLayers(PixelSubdetector::PixelBarrel);il++){
     hBchargeL_[il] = iBook.book2I(Form("recHitsBLay%dCharge",il+1), Form("CPUvsGPU RecHits Charge Barrel Layer%d;CPU Charge;GPU Charge",il+1), 250, 0, 100000, 250, 0, 100000);
     hBsizexL_[il] = iBook.book2I(Form("recHitsBLay%dSizex",il+1), Form("CPUvsGPU RecHits SizeX Barrel Layer%d;CPU SizeX;GPU SizeX",il+1), 30, 0, 30, 30, 0, 30);
     hBsizeyL_[il] = iBook.book2I(Form("recHitsBLay%dSizey",il+1), Form("CPUvsGPU RecHits SizeY Barrel Layer%d;CPU SizeY;GPU SizeY",il+1), 30, 0, 30, 30, 0, 30);
     hBposxL_[il] = iBook.book2D(Form("recHitsBLay%dPosx",il+1), Form("CPUvsGPU RecHits x-pos in Barrel Layer%d;CPU pos x;GPU pos x",il+1), 200, -5, 5, 200,-5,5);
     hBposyL_[il] = iBook.book2D(Form("recHitsBLay%dPosy",il+1), Form("CPUvsGPU RecHits y-pos in Barrel Layer%d;CPU pos y;GPU pos y",il+1), 200, -5, 5, 200,-5,5);
   }
   //Endcaps
   //Endcaps Disk
   for(int is=0;is<2;is++){
     int sign=is==0? -1:1;
     for(unsigned int id=0;id<tkGeom_->numberOfLayers(PixelSubdetector::PixelEndcap);id++){
       hFchargeD_[is][id] = iBook.book2I(Form("recHitsFDisk%+dCharge",id*sign+sign), Form("CPUvsGPU RecHits Charge Endcaps Disk%+d;CPU Charge;GPU Charge",id*sign+sign), 250, 0, 100000, 250, 0, 100000);
       hFsizexD_[is][id] = iBook.book2I(Form("recHitsFDisk%+dSizex",id*sign+sign), Form("CPUvsGPU RecHits SizeX Endcaps Disk%+d;CPU SizeX;GPU SizeX",id*sign+sign), 30, 0, 30, 30, 0, 30);
       hFsizeyD_[is][id] = iBook.book2I(Form("recHitsFDisk%+dSizey",id*sign+sign), Form("CPUvsGPU RecHits SizeY Endcaps Disk%+d;CPU SizeY;GPU SizeY",id*sign+sign), 30, 0, 30, 30, 0, 30);
       hFposxD_[is][id] = iBook.book2D(Form("recHitsFDisk%+dPosx",id*sign+sign), Form("CPUvsGPU RecHits x-pos Endcaps Disk%+d;CPU pos x;GPU pos x",id*sign+sign), 200, -5, 5, 200, -5, 5);
       hFposyD_[is][id] = iBook.book2D(Form("recHitsFDisk%+dPosy",id*sign+sign), Form("CPUvsGPU RecHits y-pos Endcaps Disk%+d;CPU pos y;GPU pos y",id*sign+sign), 200, -5, 5, 200, -5, 5);
     }
   }
   //1D differences
   hBchargeDiff_ = iBook.book1D("rechitChargeDiffBpix","Charge differnce of rechits in BPix; rechit charge difference (CPU - GPU)", 101, -50.5, 50.5);
   hFchargeDiff_ = iBook.book1D("rechitChargeDiffFpix","Charge differnce of rechits in FPix; rechit charge difference (CPU - GPU)", 101, -50.5, 50.5);
   hBsizeXDiff_ = iBook.book1D("rechitsizeXDiffBpix","SizeX difference of rechits in BPix; rechit sizex difference (CPU - GPU)", 21, -10.5, 10.5);
   hFsizeXDiff_ = iBook.book1D("rechitsizeXDiffFpix","SizeX difference of rechits in FPix; rechit sizex difference (CPU - GPU)", 21, -10.5, 10.5);
   hBsizeYDiff_ = iBook.book1D("rechitsizeYDiffBpix","SizeY difference of rechits in BPix; rechit sizey difference (CPU - GPU)", 21, -10.5, 10.5);
   hFsizeYDiff_ = iBook.book1D("rechitsizeYDiffFpix","SizeY difference of rechits in FPix; rechit sizey difference (CPU - GPU)", 21, -10.5, 10.5);
   hBposXDiff_ = iBook.book1D("rechitsposXDiffBpix","x-position difference of rechits in BPix; rechit x-pos difference (CPU - GPU)", 1000, -10, 10);
   hFposXDiff_ = iBook.book1D("rechitsposXDiffFpix","x-position difference of rechits in FPix; rechit x-pos difference (CPU - GPU)", 1000, -10, 10);
   hBposYDiff_ = iBook.book1D("rechitsposYDiffBpix","y-position difference of rechits in BPix; rechit y-pos difference (CPU - GPU)", 1000, -10, 10);
   hFposYDiff_ = iBook.book1D("rechitsposYDiffFpix","y-position difference of rechits in FPix; rechit y-pos difference (CPU - GPU)", 1000, -10, 10);
 }

◆ dqmBeginRun()

template<typename T >

void SiPixelCompareRecHitsSoA< T >::dqmBeginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 94 of file SiPixelCompareRecHitsSoA.cc.

References edm::EventSetup::getData().

                                                                                              {
   tkGeom_ = &iSetup.getData(geomToken_);
   tTopo_ = &iSetup.getData(topoToken_);
 }

◆ fillDescriptions()

template<typename T >

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

static

Definition at line 238 of file SiPixelCompareRecHitsSoA.cc.

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

                                                                                              {
   // monitorpixelRecHitsSoA
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("pixelHitsSrcCPU", edm::InputTag("siPixelRecHitsPreSplittingSoA@cpu"));
   desc.add<edm::InputTag>("pixelHitsSrcGPU", edm::InputTag("siPixelRecHitsPreSplittingSoA@cuda"));
   desc.add<std::string>("topFolderName", "SiPixelHeterogeneous/PixelRecHitsCompareGPUvsCPU");
   desc.add<double>("minD2cut", 0.0001);
   descriptions.addWithDefaultLabel(desc);
 }