Inheritance diagram for SiPixelCompareRecHitsSoAAlpaka< T >:

Public Types
using	HitsOnHost = TrackingRecHitHost< T >

Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

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

	SiPixelCompareRecHitsSoAAlpaka (const edm::ParameterSet &)

	~SiPixelCompareRecHitsSoAAlpaka () override=default

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

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.

Additional Inherited Members
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

template<typename T>
class SiPixelCompareRecHitsSoAAlpaka< T >

Definition at line 21 of file SiPixelCompareRecHitsSoAAlpaka.cc.

Member Typedef Documentation

◆ HitsOnHost

template<typename T >

using SiPixelCompareRecHitsSoAAlpaka< T >::HitsOnHost = TrackingRecHitHost<T>

Definition at line 23 of file SiPixelCompareRecHitsSoAAlpaka.cc.

Constructor & Destructor Documentation

◆ SiPixelCompareRecHitsSoAAlpaka()

template<typename T >

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

explicit

Definition at line 71 of file SiPixelCompareRecHitsSoAAlpaka.cc.

     : geomToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>()),
       topoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd, edm::Transition::BeginRun>()),
       tokenSoAHitsHost_(consumes(iConfig.getParameter<edm::InputTag>("pixelHitsSrcHost"))),
       tokenSoAHitsDevice_(consumes(iConfig.getParameter<edm::InputTag>("pixelHitsSrcDevice"))),
       topFolderName_(iConfig.getParameter<std::string>("topFolderName")),
       mind2cut_(iConfig.getParameter<double>("minD2cut")) {}

◆ ~SiPixelCompareRecHitsSoAAlpaka()

template<typename T >

SiPixelCompareRecHitsSoAAlpaka< T >::~SiPixelCompareRecHitsSoAAlpaka ( )

overridedefault

Member Function Documentation

◆ analyze()

template<typename T >

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 92 of file SiPixelCompareRecHitsSoAAlpaka.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& rhsoaHandleHost = iEvent.getHandle(tokenSoAHitsHost_);
   const auto& rhsoaHandleDevice = iEvent.getHandle(tokenSoAHitsDevice_);
   if (not rhsoaHandleHost or not rhsoaHandleDevice) {
     edm::LogWarning out("SiPixelCompareRecHitsSoAAlpaka");
     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.metadata().size();
   uint32_t nHitsDevice = soa2dDevice.metadata().size();
 
   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 SiPixelCompareRecHitsSoAAlpaka< T >::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 184 of file SiPixelCompareRecHitsSoAAlpaka.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 SiPixelCompareRecHitsSoAAlpaka< T >::dqmBeginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 83 of file SiPixelCompareRecHitsSoAAlpaka.cc.

References edm::EventSetup::getData().

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

◆ fillDescriptions()

template<typename T >

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

static

Definition at line 227 of file SiPixelCompareRecHitsSoAAlpaka.cc.

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

                                                                                                    {
   // monitorpixelRecHitsSoAAlpaka
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("pixelHitsSrcHost", edm::InputTag("siPixelRecHitsPreSplittingAlpakaSerial"));
   desc.add<edm::InputTag>("pixelHitsSrcDevice", edm::InputTag("siPixelRecHitsPreSplittingAlpaka"));
   desc.add<std::string>("topFolderName", "SiPixelHeterogeneous/PixelRecHitsCompareDeviceVSHost");
   desc.add<double>("minD2cut", 0.0001);
   descriptions.addWithDefaultLabel(desc);
 }