calibPixel::CalibDigis 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, const SiPixelGainCalibrationForHLTSoAConstView gains, int numElements) const

Detailed Description

Definition at line 29 of file CalibPixel.h.

Member Function Documentation

operator()()

template<typename TAcc >

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

inline

Definition at line 31 of file CalibPixel.h.

References cuy::col, cms::alpakatools::for_each_element_in_grid_strided(), PedestalClient_cfi::gain, ecalph2::gains, SiPixelGainUtilities::getPedAndGain(), mps_fire::i, calibPixel::InvId, SiStripPI::max, phase1PixelTopology::numberOfModules, gpuClustering::numElements, hltrates_dqm_sourceclient-live_cfg::offset, cms::alpakatools::once_per_grid(), EcalCondDBWriter_cfi::pedestal, runTheMatrix::ret, SiPixelClusterThresholds::vCaltoElectronGain, HLT_2023v12_cff::VCaltoElectronGain, SiPixelClusterThresholds::vCaltoElectronGain_L1, HLT_2023v12_cff::VCaltoElectronGain_L1, SiPixelClusterThresholds::vCaltoElectronOffset, HLT_2023v12_cff::VCaltoElectronOffset, SiPixelClusterThresholds::vCaltoElectronOffset_L1, and HLT_2023v12_cff::VCaltoElectronOffset_L1.

                                                         {
      const float VCaltoElectronGain = clusterThresholds.vCaltoElectronGain;
      const float VCaltoElectronGain_L1 = clusterThresholds.vCaltoElectronGain_L1;
      const float VCaltoElectronOffset = clusterThresholds.vCaltoElectronOffset;
      const float VCaltoElectronOffset_L1 = clusterThresholds.vCaltoElectronOffset_L1;

      // 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, phase1PixelTopology::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 (dvgi.moduleId() != InvId) {
          bool isDeadColumn = false, isNoisyColumn = false;
          int row = dvgi.xx();
          int col = dvgi.yy();
          auto ret = SiPixelGainUtilities::getPedAndGain(gains, dvgi.moduleId(), col, row, isDeadColumn, isNoisyColumn);
          float pedestal = ret.first;
          float gain = ret.second;
          if (isDeadColumn | isNoisyColumn) {
            dvgi.moduleId() = InvId;
            dvgi.adc() = 0;
            printf("bad pixel at %d in %d\n", i, dvgi.moduleId());
          } else {
            float vcal = dvgi.adc() * gain - pedestal * gain;

            float conversionFactor = dvgi.moduleId() < 96 ? VCaltoElectronGain_L1 : VCaltoElectronGain;
            float offset = dvgi.moduleId() < 96 ? VCaltoElectronOffset_L1 : VCaltoElectronOffset;
#ifdef GPU_DEBUG
            auto old_adc = dvgi.adc();
#endif
            dvgi.adc() = std::max(100, int(vcal * conversionFactor + offset));
#ifdef GPU_DEBUG
            if (cms::alpakatools::once_per_grid(acc)) {
              printf(
                  "module %d pixel %d -> old_adc = %d; vcal = %.2f; conversionFactor = %.2f; offset = %.2f; new_adc = "
                  "%d \n",
                  dvgi.moduleId(),
                  i,
                  old_adc,
                  vcal,
                  conversionFactor,
                  offset,
                  dvgi.adc());
            }
#endif
          }
        }
      });
    }