|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
|
#ifndef CondFormats_SiPixelObjects_interface_alpaka_SiPixelGainCalibrationForHLTUtilities_h
#define CondFormats_SiPixelObjects_interface_alpaka_SiPixelGainCalibrationForHLTUtilities_h
#include <cstdint>
#include <alpaka/alpaka.hpp>
#include "CondFormats/SiPixelObjects/interface/SiPixelGainCalibrationForHLTLayout.h"
struct SiPixelGainUtilities {
ALPAKA_FN_HOST_ACC ALPAKA_FN_ACC ALPAKA_FN_INLINE static std::pair<float, float> getPedAndGain(
const SiPixelGainCalibrationForHLTSoAConstView& view,
uint32_t moduleInd,
int col,
int row,
bool& isDeadColumn,
bool& isNoisyColumn) {
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);
};
};
#endif // CondFormats_SiPixelObjects_interface_alpaka_SiPixelGainCalibrationForHLTUtilities_h
|