ALPAKA_ACCELERATOR_NAMESPACE::PFRecHitProducerKernelConstruct< CAL > Struct Template Reference

Public Member Functions
template<>
ALPAKA_FN_ACC bool	applyCuts (const typename HCAL::CaloRecHitSoATypeDevice::ConstView::const_element rh, const HCAL::ParameterType::ConstView params, const HCAL::TopologyTypeDevice::ConstView topology)
template<>
ALPAKA_FN_ACC bool	applyCuts (const ECAL::CaloRecHitSoATypeDevice::ConstView::const_element rh, const ECAL::ParameterType::ConstView params, const ECAL::TopologyTypeDevice::ConstView topology)
template<>
ALPAKA_FN_ACC void	constructPFRecHit (reco::PFRecHitDeviceCollection::View::element pfrh, const HCAL::CaloRecHitSoATypeDevice::ConstView::const_element rh)
template<>
ALPAKA_FN_ACC void	constructPFRecHit (reco::PFRecHitDeviceCollection::View::element pfrh, const ECAL::CaloRecHitSoATypeDevice::ConstView::const_element rh)
template<typename TAcc , typename = std::enable_if_t<alpaka::isAccelerator<TAcc>>>
ALPAKA_FN_ACC void	operator() (const TAcc &acc, const typename CAL::ParameterType::ConstView params, const typename CAL::TopologyTypeDevice::ConstView topology, const typename CAL::CaloRecHitSoATypeDevice::ConstView recHits, reco::PFRecHitDeviceCollection::View pfRecHits, uint32_t *__restrict__ denseId2pfRecHit, uint32_t *__restrict__ num_pfRecHits) const

Static Public Member Functions
static ALPAKA_FN_ACC bool	applyCuts (const typename CAL::CaloRecHitSoATypeDevice::ConstView::const_element rh, const typename CAL::ParameterType::ConstView params, const typename CAL::TopologyTypeDevice::ConstView topology)
static ALPAKA_FN_ACC void	constructPFRecHit (reco::PFRecHitDeviceCollection::View::element pfrh, const typename CAL::CaloRecHitSoATypeDevice::ConstView::const_element rh)

Detailed Description

template<typename CAL>
struct ALPAKA_ACCELERATOR_NAMESPACE::PFRecHitProducerKernelConstruct< CAL >

Definition at line 15 of file PFRecHitProducerKernel.dev.cc.

Member Function Documentation

applyCuts() [1/3]

template<typename CAL >

static ALPAKA_FN_ACC bool ALPAKA_ACCELERATOR_NAMESPACE::PFRecHitProducerKernelConstruct< CAL >::applyCuts ( const typename CAL::CaloRecHitSoATypeDevice::ConstView::const_element  rh, const typename CAL::ParameterType::ConstView  params, const typename CAL::TopologyTypeDevice::ConstView  topology  )

static

applyCuts() [2/3]

template<>

ALPAKA_FN_ACC bool ALPAKA_ACCELERATOR_NAMESPACE::PFRecHitProducerKernelConstruct< HCAL >::applyCuts	(	const typename HCAL::CaloRecHitSoATypeDevice::ConstView::const_element	rh,
		const HCAL::ParameterType::ConstView	params,
		const HCAL::TopologyTypeDevice::ConstView	topology
	)

Definition at line 53 of file PFRecHitProducerKernel.dev.cc.

References hcalRecHitTable_cff::depth, hcalRecHitTable_cff::detId, ALPAKA_ACCELERATOR_NAMESPACE::particleFlowRecHitProducer::HCAL::detId2denseId(), ALPAKA_ACCELERATOR_NAMESPACE::particleFlowRecHitProducer::HCAL::getDepth(), ALPAKA_ACCELERATOR_NAMESPACE::particleFlowRecHitProducer::getSubdet(), HcalBarrel, HcalEndcap, ALPAKA_ACCELERATOR_NAMESPACE::particleFlowRecHitProducer::HCAL::kInvalidDenseId, ALPAKA_ACCELERATOR_NAMESPACE::particleFlowRecHitProducer::HCAL::kMaxDepthHB, submitPVValidationJobs::params, DiMuonV_cfg::threshold, and HLT_2024v14_cff::topology.

                                                      {
    // Reject HCAL recHits below enery threshold
    float threshold = 9999.;
    const uint32_t detId = rh.detId();
    const uint32_t depth = HCAL::getDepth(detId);
    const uint32_t subdet = getSubdet(detId);

    // skip bad channels
    if (rh.chi2() < 0)
      return false;

    if (topology.cutsFromDB()) {
      const auto& denseId = HCAL::detId2denseId(detId);
      if (denseId != HCAL::kInvalidDenseId) {
        threshold = topology.noiseThreshold()[denseId];
      } else {
        printf("Encountered invalid denseId for detId %u (subdetector %u)!", detId, subdet);
        return false;
      }
    } else {
      if (subdet == HcalBarrel) {
        threshold = params.energyThresholds()[depth - 1];
      } else if (subdet == HcalEndcap) {
        threshold = params.energyThresholds()[depth - 1 + HCAL::kMaxDepthHB];
      } else {
        printf("Rechit with detId %u has invalid subdetector %u!\n", detId, subdet);
        return false;
      }
    }
    return rh.energy() >= threshold;
  }

applyCuts() [3/3]

template<>

ALPAKA_FN_ACC bool ALPAKA_ACCELERATOR_NAMESPACE::PFRecHitProducerKernelConstruct< ECAL >::applyCuts	(	const ECAL::CaloRecHitSoATypeDevice::ConstView::const_element	rh,
		const ECAL::ParameterType::ConstView	params,
		const ECAL::TopologyTypeDevice::ConstView	topology
	)

Definition at line 89 of file PFRecHitProducerKernel.dev.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::particleFlowRecHitProducer::ECAL::checkFlag(), ALPAKA_ACCELERATOR_NAMESPACE::particleFlowRecHitProducer::ECAL::detId2denseId(), HLT_2024v14_cff::kTowerRecovered, EcalSeverityLevel::kWeird, and submitPVValidationJobs::params.

                                                      {
    // Reject ECAL recHits below energy threshold
    if (rh.energy() < params.energyThresholds()[ECAL::detId2denseId(rh.detId())])
      return false;

    // Reject ECAL recHits of bad quality
    if ((ECAL::checkFlag(rh.flags(), ECAL::Flags::kOutOfTime) && rh.energy() > params.cleaningThreshold()) ||
        ECAL::checkFlag(rh.flags(), ECAL::Flags::kTowerRecovered) || ECAL::checkFlag(rh.flags(), ECAL::Flags::kWeird) ||
        ECAL::checkFlag(rh.flags(), ECAL::Flags::kDiWeird))
      return false;

    return true;
  }

constructPFRecHit() [1/3]

template<typename CAL >

static ALPAKA_FN_ACC void ALPAKA_ACCELERATOR_NAMESPACE::PFRecHitProducerKernelConstruct< CAL >::constructPFRecHit ( reco::PFRecHitDeviceCollection::View::element  pfrh, const typename CAL::CaloRecHitSoATypeDevice::ConstView::const_element  rh  )

static

constructPFRecHit() [2/3]

template<>

ALPAKA_FN_ACC void ALPAKA_ACCELERATOR_NAMESPACE::PFRecHitProducerKernelConstruct< HCAL >::constructPFRecHit	(	reco::PFRecHitDeviceCollection::View::element	pfrh,
		const HCAL::CaloRecHitSoATypeDevice::ConstView::const_element	rh
	)

Definition at line 107 of file PFRecHitProducerKernel.dev.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::particleFlowRecHitProducer::HCAL::detId2denseId(), ALPAKA_ACCELERATOR_NAMESPACE::particleFlowRecHitProducer::HCAL::getDepth(), ALPAKA_ACCELERATOR_NAMESPACE::particleFlowRecHitProducer::getSubdet(), PFLayer::HCAL_BARREL1, PFLayer::HCAL_ENDCAP, HcalBarrel, HcalEndcap, and PFLayer::NONE.

                                                                  {
    pfrh.detId() = rh.detId();
    pfrh.denseId() = HCAL::detId2denseId(rh.detId());
    pfrh.energy() = rh.energy();
    pfrh.time() = rh.timeM0();
    pfrh.depth() = HCAL::getDepth(pfrh.detId());
    const uint32_t subdet = getSubdet(pfrh.detId());
    if (subdet == HcalBarrel)
      pfrh.layer() = PFLayer::HCAL_BARREL1;
    else if (subdet == HcalEndcap)
      pfrh.layer() = PFLayer::HCAL_ENDCAP;
    else
      pfrh.layer() = PFLayer::NONE;
  }

constructPFRecHit() [3/3]

template<>

ALPAKA_FN_ACC void ALPAKA_ACCELERATOR_NAMESPACE::PFRecHitProducerKernelConstruct< ECAL >::constructPFRecHit	(	reco::PFRecHitDeviceCollection::View::element	pfrh,
		const ECAL::CaloRecHitSoATypeDevice::ConstView::const_element	rh
	)

Definition at line 125 of file PFRecHitProducerKernel.dev.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::particleFlowRecHitProducer::ECAL::detId2denseId(), PFLayer::ECAL_BARREL, PFLayer::ECAL_ENDCAP, EcalBarrel, EcalEndcap, ALPAKA_ACCELERATOR_NAMESPACE::particleFlowRecHitProducer::getSubdet(), and PFLayer::NONE.

                                                                  {
    pfrh.detId() = rh.detId();
    pfrh.denseId() = ECAL::detId2denseId(rh.detId());
    pfrh.energy() = rh.energy();
    pfrh.time() = rh.time();
    pfrh.depth() = 1;
    const uint32_t subdet = getSubdet(pfrh.detId());
    if (subdet == EcalBarrel)
      pfrh.layer() = PFLayer::ECAL_BARREL;
    else if (subdet == EcalEndcap)
      pfrh.layer() = PFLayer::ECAL_ENDCAP;
    else
      pfrh.layer() = PFLayer::NONE;
  }

operator()()

template<typename CAL >

template<typename TAcc , typename = std::enable_if_t<alpaka::isAccelerator<TAcc>>>

ALPAKA_FN_ACC void ALPAKA_ACCELERATOR_NAMESPACE::PFRecHitProducerKernelConstruct< CAL >::operator() ( const TAcc &  acc, const typename CAL::ParameterType::ConstView  params, const typename CAL::TopologyTypeDevice::ConstView  topology, const typename CAL::CaloRecHitSoATypeDevice::ConstView  recHits, reco::PFRecHitDeviceCollection::View  pfRecHits, uint32_t *__restrict__  denseId2pfRecHit, uint32_t *__restrict__  num_pfRecHits  ) const

inline

Definition at line 17 of file PFRecHitProducerKernel.dev.cc.

References cms::cudacompat::atomicInc(), mps_fire::i, dqmiolumiharvest::j, submitPVValidationJobs::params, HLT_2024v14_cff::pfRecHits, FastTrackerRecHitMaskProducer_cfi::recHits, HLT_2024v14_cff::topology, and cms::alpakatools::uniform_elements().

                                                                              {
      // Strided loop over CaloRecHits
      for (int32_t i : cms::alpakatools::uniform_elements(acc, recHits.metadata().size())) {
        // Check energy thresholds/quality cuts (specialised for HCAL/ECAL)
        if (!applyCuts(recHits[i], params, topology))
          continue;

        // Use atomic operation to determine index of the PFRecHit to be constructed
        // The index needs to be unique and consequtive across all threads in all blocks.
        // This is achieved using the alpaka::hierarchy::Blocks argument.
        const uint32_t j = alpaka::atomicInc(acc, num_pfRecHits, 0xffffffff, alpaka::hierarchy::Blocks{});

        // Construct PFRecHit from CAL recHit (specialised for HCAL/ECAL)
        constructPFRecHit(pfRecHits[j], recHits[i]);

        // Fill denseId -> pfRecHit index map
        denseId2pfRecHit[CAL::detId2denseId(pfRecHits.detId(j))] = j;
      }
    }