GPU TPC: For constraining track time0, use all clusters, not only those that were not reject in current refit round

davidrohr · davidrohr · commit 7980c6282bad · 2023-01-20T11:00:42.000+01:00
diff --git a/GPU/GPUTracking/Merger/GPUTPCGMTrackParam.cxx b/GPU/GPUTracking/Merger/GPUTPCGMTrackParam.cxx
@@ -79,8 +79,8 @@ GPUd() bool GPUTPCGMTrackParam::Fit(GPUTPCGMMerger* GPUrestrict() merger, int iT
     }
   }
 
-  int nWays = param.rec.tpc.nWays;
-  int maxN = N;
+  const int nWays = param.rec.tpc.nWays;
+  const int maxN = N;
   int ihitStart = 0;
   float covYYUpd = 0.f;
   float lastUpdateX = -1.f;
@@ -362,7 +362,7 @@ GPUd() bool GPUTPCGMTrackParam::Fit(GPUTPCGMMerger* GPUrestrict() merger, int iT
       }
     }
     if (((nWays - iWay) & 1) && (clusters[0].slice < 18) == (clusters[maxN - 1].slice < 18)) {
-      ShiftZ2(clusters, clustersXYZ, merger, N);
+      ShiftZ2(clusters, clustersXYZ, merger, maxN);
     }
   }
   ConstrainSinPhi();

Original file line number	Diff line number	Diff line change
`@@ -79,8 +79,8 @@ GPUd() bool GPUTPCGMTrackParam::Fit(GPUTPCGMMerger* GPUrestrict() merger, int iT`
`79`	`79`	`}`
`80`	`80`	`}`
`81`	`81`
`82`		`- int nWays = param.rec.tpc.nWays;`
`83`		`- int maxN = N;`
	`82`	`+ const int nWays = param.rec.tpc.nWays;`
	`83`	`+ const int maxN = N;`
`84`	`84`	`int ihitStart = 0;`
`85`	`85`	`float covYYUpd = 0.f;`
`86`	`86`	`float lastUpdateX = -1.f;`
`@@ -362,7 +362,7 @@ GPUd() bool GPUTPCGMTrackParam::Fit(GPUTPCGMMerger* GPUrestrict() merger, int iT`
`362`	`362`	`}`
`363`	`363`	`}`
`364`	`364`	`if (((nWays - iWay) & 1) && (clusters[0].slice < 18) == (clusters[maxN - 1].slice < 18)) {`
`365`		`- ShiftZ2(clusters, clustersXYZ, merger, N);`
	`365`	`+ ShiftZ2(clusters, clustersXYZ, merger, maxN);`
`366`	`366`	`}`
`367`	`367`	`}`
`368`	`368`	`ConstrainSinPhi();`