SiPixelGainUtilities Struct Reference

#include <SiPixelGainCalibrationForHLTUtilities.h>

Static Public Member Functions
ALPAKA_FN_HOST_ACC ALPAKA_FN_ACC static ALPAKA_FN_INLINE std::pair< float, float >	getPedAndGain (const SiPixelGainCalibrationForHLTSoAConstView &view, uint32_t moduleInd, int col, int row, bool &isDeadColumn, bool &isNoisyColumn)

Detailed Description

Definition at line 10 of file SiPixelGainCalibrationForHLTUtilities.h.

Member Function Documentation

getPedAndGain()

ALPAKA_FN_HOST_ACC ALPAKA_FN_ACC static ALPAKA_FN_INLINE std::pair<float, float> SiPixelGainUtilities::getPedAndGain ( const SiPixelGainCalibrationForHLTSoAConstView &  view, uint32_t  moduleInd, int  col, int  row, bool &  isDeadColumn, bool &  isNoisyColumn  )

inlinestatic

Definition at line 11 of file SiPixelGainCalibrationForHLTUtilities.h.

References ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::ALPAKA_ASSERT_ACC(), cuy::col, mps_fire::end, ALCARECOEcalPhiSym_cff::float, HLT_IsoTrack_cff::offset, alignCSCRings::s, and isotrackNtupler::start.

Referenced by calibPixel::CalibDigis< debug >::operator()().

                           {
    auto start = view.modStarts()[moduleInd];
    auto end = view.modEnds()[moduleInd];
    auto nCols = view.modCols()[moduleInd];
    // determine what averaged data block we are in (there should be 1 or 2 of these depending on if plaquette is 1 by X or 2 by X
    unsigned int lengthOfColumnData = (end - start) / nCols;
    unsigned int lengthOfAveragedDataInEachColumn = 2;  // we always only have two values per column averaged block
    unsigned int numberOfDataBlocksToSkip = row / view.numberOfRowsAveragedOver();

    auto offset = start + col * lengthOfColumnData + lengthOfAveragedDataInEachColumn * numberOfDataBlocksToSkip;
    ALPAKA_ASSERT_ACC(offset < end);
    ALPAKA_ASSERT_ACC(offset < 3088384);
    ALPAKA_ASSERT_ACC(0 == offset % 2);

    auto lp = view.v_pedestals();
    auto s = lp[offset / 2];

    isDeadColumn = (s.ped & 0xFF) == view.deadFlag();
    isNoisyColumn = (s.ped & 0xFF) == view.noisyFlag();
    float decodeGain = float(s.gain & 0xFF) * view.gainPrecision() + view.minGain();
    float decodePed = float(s.ped & 0xFF) * view.pedPrecision() + view.minPed();

    return std::make_pair(decodePed, decodeGain);
  };