SiPixelCompareVertexSoA Class Reference

#include <SiPixelCompareVertexSoA.cc>

Inheritance diagram for SiPixelCompareVertexSoA:

Public Types
using	IndToEdm = std::vector< uint16_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
	SiPixelCompareVertexSoA (const edm::ParameterSet &)
	~SiPixelCompareVertexSoA () 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
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	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 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< ZVertexSoAHost >	tokenSoAVertexCPU_
const edm::EDGetTokenT< ZVertexSoAHost >	tokenSoAVertexGPU_
const std::string	topFolderName_

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

Definition at line 24 of file SiPixelCompareVertexSoA.cc.

Member Typedef Documentation

IndToEdm

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

Definition at line 26 of file SiPixelCompareVertexSoA.cc.

Constructor & Destructor Documentation

SiPixelCompareVertexSoA()

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

explicit

Definition at line 57 of file SiPixelCompareVertexSoA.cc.

    : tokenSoAVertexCPU_(consumes<ZVertexSoAHost>(iConfig.getParameter<edm::InputTag>("pixelVertexSrcCPU"))),
      tokenSoAVertexGPU_(consumes<ZVertexSoAHost>(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")) {}

~SiPixelCompareVertexSoA()

SiPixelCompareVertexSoA::~SiPixelCompareVertexSoA ( )

overridedefault

Member Function Documentation

analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 67 of file SiPixelCompareVertexSoA.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("SiPixelCompareVertexSoA");
    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;
  int nVerticesCPU = vsoaCPU.view().nvFinal();
  auto const& vsoaGPU = *vsoaHandleGPU;
  int nVerticesGPU = vsoaGPU.view().nvFinal();

  auto bsHandle = iEvent.getHandle(tokenBeamSpot_);
  float x0 = 0., y0 = 0., z0 = 0., dxdz = 0., dydz = 0.;
  if (!bsHandle.isValid()) {
    edm::LogWarning("SiPixelCompareVertexSoA") << "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.view()[ivc].sortInd();
    auto zc = vsoaCPU.view()[sic].zv();
    auto xc = x0 + dxdz * zc;
    auto yc = y0 + dydz * zc;
    zc += z0;

    auto ndofCPU = vsoaCPU.view()[sic].ndof();
    auto chi2CPU = vsoaCPU.view()[sic].chi2();

    const int32_t notFound = -1;
    int32_t closestVtxidx = notFound;
    float mindz = dzCut_;

    for (int ivg = 0; ivg < nVerticesGPU; ivg++) {
      auto sig = vsoaGPU.view()[ivg].sortInd();
      auto zgc = vsoaGPU.view()[sig].zv() + 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.view()[closestVtxidx].zv();
    auto xg = x0 + dxdz * zg;
    auto yg = y0 + dydz * zg;
    zg += z0;
    auto ndofGPU = vsoaGPU.view()[closestVtxidx].ndof();
    auto chi2GPU = vsoaGPU.view()[closestVtxidx].chi2();

    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.view()[sic].ptv2(), vsoaGPU.view()[closestVtxidx].ptv2());
    hntrks_->Fill(ndofCPU + 1, ndofGPU + 1);
  }
  hnVertex_->Fill(nVerticesCPU, nVerticesGPU);
}

bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 153 of file SiPixelCompareVertexSoA.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 Vertices;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 #sum (p_{T})^{2};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 SiPixelCompareVertexSoA::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 175 of file SiPixelCompareVertexSoA.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, HLT_2022v15_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);
}

Member Data Documentation

dzCut_

const float SiPixelCompareVertexSoA::dzCut_

private

Definition at line 38 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze().

hchi2_

MonitorElement* SiPixelCompareVertexSoA::hchi2_

private

Definition at line 43 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze(), and bookHistograms().

hchi2oNdof_

MonitorElement* SiPixelCompareVertexSoA::hchi2oNdof_

private

Definition at line 44 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze(), and bookHistograms().

hntrks_

MonitorElement* SiPixelCompareVertexSoA::hntrks_

private

Definition at line 46 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze(), and bookHistograms().

hnVertex_

MonitorElement* SiPixelCompareVertexSoA::hnVertex_

private

Definition at line 39 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze(), and bookHistograms().

hptv2_

MonitorElement* SiPixelCompareVertexSoA::hptv2_

private

Definition at line 45 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze(), and bookHistograms().

hx_

MonitorElement* SiPixelCompareVertexSoA::hx_

private

Definition at line 40 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze(), and bookHistograms().

hxdiff_

MonitorElement* SiPixelCompareVertexSoA::hxdiff_

private

Definition at line 47 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze(), and bookHistograms().

hy_

MonitorElement* SiPixelCompareVertexSoA::hy_

private

Definition at line 41 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze(), and bookHistograms().

hydiff_

MonitorElement* SiPixelCompareVertexSoA::hydiff_

private

Definition at line 48 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze(), and bookHistograms().

hz_

MonitorElement* SiPixelCompareVertexSoA::hz_

private

Definition at line 42 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze(), and bookHistograms().

hzdiff_

MonitorElement* SiPixelCompareVertexSoA::hzdiff_

private

Definition at line 49 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze(), and bookHistograms().

tokenBeamSpot_

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

private

Definition at line 36 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze().

tokenSoAVertexCPU_

const edm::EDGetTokenT<ZVertexSoAHost> SiPixelCompareVertexSoA::tokenSoAVertexCPU_

private

Definition at line 34 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze().

tokenSoAVertexGPU_

const edm::EDGetTokenT<ZVertexSoAHost> SiPixelCompareVertexSoA::tokenSoAVertexGPU_

private

Definition at line 35 of file SiPixelCompareVertexSoA.cc.

Referenced by analyze().

topFolderName_

const std::string SiPixelCompareVertexSoA::topFolderName_

private

Definition at line 37 of file SiPixelCompareVertexSoA.cc.

Referenced by bookHistograms().