GPU: fix bug when filling subthreshold clusters

davidrohr · davidrohr · commit 4f7d71b64ef2 · 2024-11-25T07:59:57.000+01:00
diff --git a/GPU/GPUTracking/Merger/GPUTPCGMTrackParam.cxx b/GPU/GPUTracking/Merger/GPUTPCGMTrackParam.cxx
@@ -223,7 +223,7 @@ GPUd() bool GPUTPCGMTrackParam::Fit(GPUTPCGMMerger* GPUrestrict() merger, int32_
         bool dodEdx = param.par.dodEdx && param.dodEdxDownscaled && param.rec.tpc.adddEdxSubThresholdClusters && iWay == nWays - 1 && CAMath::Abs(cluster.row - lastRow) == 2 && cluster.leg == clusters[maxN - 1].leg;
         dodEdx = AttachClustersPropagate(merger, cluster.slice, lastRow, cluster.row, iTrk, cluster.leg == clusters[maxN - 1].leg, prop, inFlyDirection, GPUCA_MAX_SIN_PHI, dodEdx);
         if (dodEdx) {
-          dEdx.fillSubThreshold(lastRow - 1, param);
+          dEdx.fillSubThreshold(lastRow - wayDirection, param);
         }
       }
 

Original file line number	Diff line number	Diff line change
`@@ -223,7 +223,7 @@ GPUd() bool GPUTPCGMTrackParam::Fit(GPUTPCGMMerger* GPUrestrict() merger, int32_`
`223`	`223`	`bool dodEdx = param.par.dodEdx && param.dodEdxDownscaled && param.rec.tpc.adddEdxSubThresholdClusters && iWay == nWays - 1 && CAMath::Abs(cluster.row - lastRow) == 2 && cluster.leg == clusters[maxN - 1].leg;`
`224`	`224`	`dodEdx = AttachClustersPropagate(merger, cluster.slice, lastRow, cluster.row, iTrk, cluster.leg == clusters[maxN - 1].leg, prop, inFlyDirection, GPUCA_MAX_SIN_PHI, dodEdx);`
`225`	`225`	`if (dodEdx) {`
`226`		`- dEdx.fillSubThreshold(lastRow - 1, param);`
	`226`	`+ dEdx.fillSubThreshold(lastRow - wayDirection, param);`
`227`	`227`	`}`
`228`	`228`	`}`
`229`	`229`