calibPixel::CalibDigisPhase2 Struct Reference

#include <CalibPixel.h>

Public Member Functions
template<typename TAcc >
ALPAKA_FN_ACC void	operator() (const TAcc &acc, SiPixelClusterThresholds clusterThresholds, SiPixelDigisSoAView view, SiPixelClustersSoAView clus_view, int numElements) const

Detailed Description

Definition at line 90 of file CalibPixel.h.

Member Function Documentation

operator()()

template<typename TAcc >

ALPAKA_FN_ACC void calibPixel::CalibDigisPhase2::operator() ( const TAcc &  acc, SiPixelClusterThresholds  clusterThresholds, SiPixelDigisSoAView  view, SiPixelClustersSoAView  clus_view, int  numElements  ) const

inline

Definition at line 92 of file CalibPixel.h.

References SiPixelClusterThresholds::electronPerADCGain, HLT_2023v12_cff::ElectronPerADCGain, cms::alpakatools::for_each_element_in_grid_strided(), mps_fire::i, createfilelist::int, pixelClustering::invalidModuleId, SiStripPI::max, SiStripPI::min, ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, phase2PixelTopology::numberOfModules, gpuClustering::numElements, cms::alpakatools::once_per_grid(), SiPixelClusterThresholds::phase2DigiBaseline, HLT_2023v12_cff::Phase2DigiBaseline, SiPixelClusterThresholds::phase2KinkADC, HLT_2023v12_cff::Phase2KinkADC, SiPixelClusterThresholds::phase2ReadoutMode, and HLT_2023v12_cff::Phase2ReadoutMode.

                                                         {
      const float ElectronPerADCGain = clusterThresholds.electronPerADCGain;
      const int8_t Phase2ReadoutMode = clusterThresholds.phase2ReadoutMode;
      const uint16_t Phase2DigiBaseline = clusterThresholds.phase2DigiBaseline;
      const uint8_t Phase2KinkADC = clusterThresholds.phase2KinkADC;

      // zero for next kernels...
      if (cms::alpakatools::once_per_grid(acc)) {
        clus_view[0].clusModuleStart() = clus_view[0].moduleStart() = 0;
      }

      cms::alpakatools::for_each_element_in_grid_strided(
          acc, phase2PixelTopology::numberOfModules, [&](uint32_t i) { clus_view[i].clusInModule() = 0; });
      cms::alpakatools::for_each_element_in_grid_strided(acc, numElements, [&](uint32_t i) {
        auto dvgi = view[i];
        if (pixelClustering::invalidModuleId != dvgi.moduleId()) {
          const int mode = (Phase2ReadoutMode < -1 ? -1 : Phase2ReadoutMode);
          int adc_int = dvgi.adc();
          if (mode < 0)
            adc_int = int(adc_int * ElectronPerADCGain);
          else {
            if (adc_int < Phase2KinkADC)
              adc_int = int((adc_int + 0.5) * ElectronPerADCGain);
            else {
              const int8_t dspp = (Phase2ReadoutMode < 10 ? Phase2ReadoutMode : 10);
              const int8_t ds = int8_t(dspp <= 1 ? 1 : (dspp - 1) * (dspp - 1));
              adc_int -= Phase2KinkADC;
              adc_int *= ds;
              adc_int += Phase2KinkADC;
              adc_int = ((adc_int + 0.5 * ds) * ElectronPerADCGain);
            }
            adc_int += int(Phase2DigiBaseline);
          }
          dvgi.adc() = std::min(adc_int, int(std::numeric_limits<uint16_t>::max()));
        }
      });
    }