From 96e41889866f27476804d9ce6b1d89b0e48ab75b Mon Sep 17 00:00:00 2001
From: David Rohr <drohr@jwdt.org>
Date: Wed, 27 Nov 2024 14:08:25 +0100
Subject: [PATCH] GPU TPC: Fix TPC clusterizer qMax cut

---
 GPU/GPUTracking/TPCClusterFinder/ClusterAccumulator.cxx | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/GPU/GPUTracking/TPCClusterFinder/ClusterAccumulator.cxx b/GPU/GPUTracking/TPCClusterFinder/ClusterAccumulator.cxx
index 8988126f7a15e..e8176ecb60d78 100644
--- a/GPU/GPUTracking/TPCClusterFinder/ClusterAccumulator.cxx
+++ b/GPU/GPUTracking/TPCClusterFinder/ClusterAccumulator.cxx
@@ -27,6 +27,10 @@ GPUd() bool ClusterAccumulator::toNative(const ChargePos& pos, Charge q, tpc::Cl
   if (cn.qTot <= param.rec.tpc.cfQTotCutoff) {
     return false;
   }
+  cn.qMax = q;
+  if (cn.qMax <= param.rec.tpc.cfQMaxCutoff) {
+    return false;
+  }
   if (mTimeMean < param.rec.tpc.clustersShiftTimebinsClusterizer) {
     return false;
   }
@@ -48,7 +52,6 @@ GPUd() bool ClusterAccumulator::toNative(const ChargePos& pos, Charge q, tpc::Cl
   flags |= (wasSplitInPad) ? tpc::ClusterNative::flagSplitPad : 0;
   flags |= (isSingleCluster) ? tpc::ClusterNative::flagSingle : 0;
 
-  cn.qMax = q;
   cn.setTimeFlags(mTimeMean - param.rec.tpc.clustersShiftTimebinsClusterizer, flags);
   cn.setPad(mPadMean);
   cn.setSigmaTime(mTimeSigma);