EXPERIMENTAL! GPU: Process dEdx with full qTot range

Author: davidrohr

DataFormats/Detectors/TPC/include/DataFormatsTPC/ClusterNative.h

@@ -87,6 +87,13 @@ struct ClusterNative {
   }
 
   GPUd() uint16_t getQmax() const { return qMax; }
+  GPUd() float getFullQtot() const
+  {
+    if (isSaturated()) [[unlikely]] {
+      return getSaturatedQtot();
+    }
+    return qTot;
+  }
   GPUd() uint16_t getQtot() const
   {
     if (isSaturated()) [[unlikely]] {

GPU/GPUTracking/Definitions/GPUDefParametersWrapper.h

@@ -37,7 +37,7 @@ namespace o2::gpu
 #endif
 
 #define GPUCA_PAR_MERGER_INTERPOLATION_ERROR_TYPE_A GPUCA_DETERMINISTIC_CODE(float, GPUCA_PAR_MERGER_INTERPOLATION_ERROR_TYPE)
-#define GPUCA_PAR_DEDX_STORAGE_TYPE_A GPUCA_DETERMINISTIC_CODE(float, GPUCA_PAR_DEDX_STORAGE_TYPE)
+#define GPUCA_PAR_DEDX_STORAGE_TYPE_A float
 
 // #define GPUCA_TRACKLET_CONSTRUCTOR_DO_PROFILE                       // Output Profiling Data for Tracklet Constructor Tracklet Scheduling
 

GPU/GPUTracking/Merger/GPUTPCGMTrackParam.cxx

@@ -294,7 +294,7 @@ GPUd() bool GPUTPCGMTrackParam::Fit(GPUTPCGMMerger* GPUrestrict() merger, int32_
               const int32_t clusterCount = (ihit - ihitMergeFirst) * wayDirection + 1;
               for (int32_t iTmp = ihitMergeFirst; iTmp != ihit + wayDirection; iTmp += wayDirection) {
                 const ClusterNative& cl = merger->GetConstantMem()->ioPtrs.clustersNative->clustersLinear[cluster.num];
-                qtot += cl.qTot;
+                qtot += cl.getFullQtot();
                 qmax = CAMath::Max<float>(qmax, cl.qMax);
                 pad += cl.getPad();
                 relTime += cl.getTime();