calibPixel::CalibDigis< debug > Struct Template 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

template<bool debug = false>
struct calibPixel::CalibDigis< debug >

Definition at line 26 of file CalibPixel.h.

Member Function Documentation

operator()()

template<bool debug = false>

template<typename TAcc >

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

inline

Definition at line 28 of file CalibPixel.h.

References cuy::col, ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), debug, Hcal_Conditions_forGlobalTag_cff::gain, ecalph2::gains, SiPixelGainUtilities::getPedAndGain(), mps_fire::i, gpuClustering::invalidModuleId, WZElectronSkims53X_cff::max, phase1PixelTopology::numberOfModules, gpuClustering::numElements, HLT_IsoTrack_cff::offset, cms::alpakatools::once_per_grid(), Hcal_Conditions_forGlobalTag_cff::pedestal, runTheMatrix::ret, cms::alpakatools::uniform_elements(), SiPixelClusterThresholds::vCaltoElectronGain, HLT_2024v14_cff::VCaltoElectronGain, SiPixelClusterThresholds::vCaltoElectronGain_L1, HLT_2024v14_cff::VCaltoElectronGain_L1, SiPixelClusterThresholds::vCaltoElectronOffset, HLT_2024v14_cff::VCaltoElectronOffset, SiPixelClusterThresholds::vCaltoElectronOffset_L1, and HLT_2024v14_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() = 0;
        clus_view[0].moduleStart() = 0;
      }
      for (auto i : cms::alpakatools::uniform_elements(acc, phase1PixelTopology::numberOfModules)) {
        clus_view[i].clusInModule() = 0;
      }

      for (auto i : cms::alpakatools::uniform_elements(acc, numElements)) {
        auto dvgi = view[i];
        if (dvgi.moduleId() == ::pixelClustering::invalidModuleId)
          continue;

        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) {
          if constexpr (debug)
            printf("bad pixel at %d in %d\n", i, dvgi.moduleId());
          dvgi.moduleId() = ::pixelClustering::invalidModuleId;
          dvgi.adc() = 0;
        } 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
        }
      }
    }