#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 >	tokenSoAHitsDevice_

const edm::EDGetTokenT< HitsOnHost >	tokenSoAHitsHost_

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>()),
       tokenSoAHitsHost_(consumes(iConfig.getParameter<edm::InputTag>("pixelHitsSrcCPU"))),
       tokenSoAHitsDevice_(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_2022v15_cff::distance, PVValHelper::dx, PVValHelper::dy, mps_fire::i, iEvent, dqmiolumiharvest::j, or, MillePedeFileConverter_cfg::out, PixelSubdetector::PixelBarrel, and PixelSubdetector::PixelEndcap.

                                                                                              {
   const auto& rhsoaHandleHost = iEvent.getHandle(tokenSoAHitsHost_);
   const auto& rhsoaHandleDevice = iEvent.getHandle(tokenSoAHitsDevice_);
   if (not rhsoaHandleHost or not rhsoaHandleDevice) {
     edm::LogWarning out("SiPixelCompareRecHitSoA");
     if (not rhsoaHandleHost) {
       out << "reference (Host) rechits not found; ";
     }
     if (not rhsoaHandleDevice) {
       out << "target (Device) rechits not found; ";
     }
     out << "the comparison will not run.";
     return;
   }
 
   auto const& rhsoaHost = *rhsoaHandleHost;
   auto const& rhsoaDevice = *rhsoaHandleDevice;
 
   auto const& soa2dHost = rhsoaHost.const_view();
   auto const& soa2dDevice = rhsoaDevice.const_view();
 
   uint32_t nHitsHost = soa2dHost.nHits();
   uint32_t nHitsDevice = soa2dDevice.nHits();
 
   hnHits_->Fill(nHitsHost, nHitsDevice);
   auto detIds = tkGeom_->detUnitIds();
   for (uint32_t i = 0; i < nHitsHost; i++) {
     float minD = mind2cut_;
     uint32_t matchedHit = invalidHit_;
     uint16_t indHost = soa2dHost[i].detectorIndex();
     float xLocalHost = soa2dHost[i].xLocal();
     float yLocalHost = soa2dHost[i].yLocal();
     for (uint32_t j = 0; j < nHitsDevice; j++) {
       if (soa2dDevice.detectorIndex(j) == indHost) {
         float dx = xLocalHost - soa2dDevice[j].xLocal();
         float dy = yLocalHost - soa2dDevice[j].yLocal();
         float distance = dx * dx + dy * dy;
         if (distance < minD) {
           minD = distance;
           matchedHit = j;
         }
       }
     }
     DetId id = detIds[indHost];
     uint32_t chargeHost = soa2dHost[i].chargeAndStatus().charge;
     int16_t sizeXHost = std::ceil(float(std::abs(soa2dHost[i].clusterSizeX()) / 8.));
     int16_t sizeYHost = std::ceil(float(std::abs(soa2dHost[i].clusterSizeY()) / 8.));
     uint32_t chargeDevice = 0;
     int16_t sizeXDevice = -99;
     int16_t sizeYDevice = -99;
     float xLocalDevice = -999.;
     float yLocalDevice = -999.;
     if (matchedHit != invalidHit_) {
       chargeDevice = soa2dDevice[matchedHit].chargeAndStatus().charge;
       sizeXDevice = std::ceil(float(std::abs(soa2dDevice[matchedHit].clusterSizeX()) / 8.));
       sizeYDevice = std::ceil(float(std::abs(soa2dDevice[matchedHit].clusterSizeY()) / 8.));
       xLocalDevice = soa2dDevice[matchedHit].xLocal();
       yLocalDevice = soa2dDevice[matchedHit].yLocal();
     }
     switch (id.subdetId()) {
       case PixelSubdetector::PixelBarrel:
         hBchargeL_[tTopo_->pxbLayer(id) - 1]->Fill(chargeHost, chargeDevice);
         hBsizexL_[tTopo_->pxbLayer(id) - 1]->Fill(sizeXHost, sizeXDevice);
         hBsizeyL_[tTopo_->pxbLayer(id) - 1]->Fill(sizeYHost, sizeYDevice);
         hBposxL_[tTopo_->pxbLayer(id) - 1]->Fill(xLocalHost, xLocalDevice);
         hBposyL_[tTopo_->pxbLayer(id) - 1]->Fill(yLocalHost, yLocalDevice);
         hBchargeDiff_->Fill(chargeHost - chargeDevice);
         hBsizeXDiff_->Fill(sizeXHost - sizeXDevice);
         hBsizeYDiff_->Fill(sizeYHost - sizeYDevice);
         hBposXDiff_->Fill(micron_ * (xLocalHost - xLocalDevice));
         hBposYDiff_->Fill(micron_ * (yLocalHost - yLocalDevice));
         break;
       case PixelSubdetector::PixelEndcap:
         hFchargeD_[tTopo_->pxfSide(id) - 1][tTopo_->pxfDisk(id) - 1]->Fill(chargeHost, chargeDevice);
         hFsizexD_[tTopo_->pxfSide(id) - 1][tTopo_->pxfDisk(id) - 1]->Fill(sizeXHost, sizeXDevice);
         hFsizeyD_[tTopo_->pxfSide(id) - 1][tTopo_->pxfDisk(id) - 1]->Fill(sizeYHost, sizeYDevice);
         hFposxD_[tTopo_->pxfSide(id) - 1][tTopo_->pxfDisk(id) - 1]->Fill(xLocalHost, xLocalDevice);
         hFposyD_[tTopo_->pxfSide(id) - 1][tTopo_->pxfDisk(id) - 1]->Fill(yLocalHost, yLocalDevice);
         hFchargeDiff_->Fill(chargeHost - chargeDevice);
         hFsizeXDiff_->Fill(sizeXHost - sizeXDevice);
         hFsizeYDiff_->Fill(sizeYHost - sizeYDevice);
         hFposXDiff_->Fill(micron_ * (xLocalHost - xLocalDevice));
         hFposYDiff_->Fill(micron_ * (yLocalHost - yLocalDevice));
         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", "HostvsDevice RecHits per event;#Host RecHits;#Device 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("HostvsDevice RecHits Charge Barrel Layer%d;Host Charge;Device Charge",il+1), 250, 0, 100000, 250, 0, 100000);
     hBsizexL_[il] = iBook.book2I(Form("recHitsBLay%dSizex",il+1), Form("HostvsDevice RecHits SizeX Barrel Layer%d;Host SizeX;Device SizeX",il+1), 30, 0, 30, 30, 0, 30);
     hBsizeyL_[il] = iBook.book2I(Form("recHitsBLay%dSizey",il+1), Form("HostvsDevice RecHits SizeY Barrel Layer%d;Host SizeY;Device SizeY",il+1), 30, 0, 30, 30, 0, 30);
     hBposxL_[il] = iBook.book2D(Form("recHitsBLay%dPosx",il+1), Form("HostvsDevice RecHits x-pos in Barrel Layer%d;Host pos x;Device pos x",il+1), 200, -5, 5, 200,-5,5);
     hBposyL_[il] = iBook.book2D(Form("recHitsBLay%dPosy",il+1), Form("HostvsDevice RecHits y-pos in Barrel Layer%d;Host pos y;Device 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("HostvsDevice RecHits Charge Endcaps Disk%+d;Host Charge;Device Charge",id*sign+sign), 250, 0, 100000, 250, 0, 100000);
       hFsizexD_[is][id] = iBook.book2I(Form("recHitsFDisk%+dSizex",id*sign+sign), Form("HostvsDevice RecHits SizeX Endcaps Disk%+d;Host SizeX;Device SizeX",id*sign+sign), 30, 0, 30, 30, 0, 30);
       hFsizeyD_[is][id] = iBook.book2I(Form("recHitsFDisk%+dSizey",id*sign+sign), Form("HostvsDevice RecHits SizeY Endcaps Disk%+d;Host SizeY;Device SizeY",id*sign+sign), 30, 0, 30, 30, 0, 30);
       hFposxD_[is][id] = iBook.book2D(Form("recHitsFDisk%+dPosx",id*sign+sign), Form("HostvsDevice RecHits x-pos Endcaps Disk%+d;Host pos x;Device pos x",id*sign+sign), 200, -5, 5, 200, -5, 5);
       hFposyD_[is][id] = iBook.book2D(Form("recHitsFDisk%+dPosy",id*sign+sign), Form("HostvsDevice RecHits y-pos Endcaps Disk%+d;Host pos y;Device 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 (Host - Device)", 101, -50.5, 50.5);
   hFchargeDiff_ = iBook.book1D("rechitChargeDiffFpix","Charge differnce of rechits in FPix; rechit charge difference (Host - Device)", 101, -50.5, 50.5);
   hBsizeXDiff_ = iBook.book1D("rechitsizeXDiffBpix","SizeX difference of rechits in BPix; rechit sizex difference (Host - Device)", 21, -10.5, 10.5);
   hFsizeXDiff_ = iBook.book1D("rechitsizeXDiffFpix","SizeX difference of rechits in FPix; rechit sizex difference (Host - Device)", 21, -10.5, 10.5);
   hBsizeYDiff_ = iBook.book1D("rechitsizeYDiffBpix","SizeY difference of rechits in BPix; rechit sizey difference (Host - Device)", 21, -10.5, 10.5);
   hFsizeYDiff_ = iBook.book1D("rechitsizeYDiffFpix","SizeY difference of rechits in FPix; rechit sizey difference (Host - Device)", 21, -10.5, 10.5);
   hBposXDiff_ = iBook.book1D("rechitsposXDiffBpix","x-position difference of rechits in BPix; rechit x-pos difference (Host - Device)", 1000, -10, 10);
   hFposXDiff_ = iBook.book1D("rechitsposXDiffFpix","x-position difference of rechits in FPix; rechit x-pos difference (Host - Device)", 1000, -10, 10);
   hBposYDiff_ = iBook.book1D("rechitsposYDiffBpix","y-position difference of rechits in BPix; rechit y-pos difference (Host - Device)", 1000, -10, 10);
   hFposYDiff_ = iBook.book1D("rechitsposYDiffFpix","y-position difference of rechits in FPix; rechit y-pos difference (Host - Device)", 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, HLT_2022v15_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

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