@@ -133,7 +133,7 @@ GPUd() bool GPUTPCGMTrackParam::Fit(GPUTPCGMMerger& GPUrestrict() merger, int32_
133133 const float clAlpha = param.Alpha (clusters[ihit].sector );
134134 float xx, yy, zz;
135135 {
136- const ClusterNative& GPUrestrict () cl = merger.GetConstantMem ()-> ioPtrs .clustersNative ->clustersLinear [clusters[ihit].num ];
136+ const ClusterNative& GPUrestrict () cl = merger.GetConstantMem () -> ioPtrs.clustersNative ->clustersLinear [clusters[ihit].num ];
137137 merger.GetConstantMem ()->calibObjects .fastTransformHelper ->Transform (clusters[ihit].sector , clusters[ihit].row , cl.getPad (), cl.getTime (), xx, yy, zz, mTOffset );
138138 }
139139 // clang-format off
@@ -367,7 +367,7 @@ GPUdii() int32_t GPUTPCGMTrackParam::FitHit(GPUTPCGMMerger& GPUrestrict() merger
367367GPUdii () float GPUTPCGMTrackParam::FindBestInterpolatedHit(GPUTPCGMMerger& GPUrestrict () merger, gputpcgmmergertypes::InterpolationErrorHit& GPUrestrict() inter, uint8_t sector, uint8_t row, const float deltaZ, const float sumInvSqrtCharge, const int nAvgCharge, const GPUTPCGMPropagator& GPUrestrict() prop, int32_t iTrk)
368368{
369369 const GPUParam& GPUrestrict () param = merger.Param ();
370- const GPUTPCTracker& GPUrestrict () tracker = *(merger.GetConstantMem ()-> tpcTrackers + sector);
370+ const GPUTPCTracker& GPUrestrict () tracker = *(merger.GetConstantMem () -> tpcTrackers + sector);
371371 const GPUTPCRow& GPUrestrict () rowData = tracker.Row (row);
372372 GPUglobalref () const cahit2* hits = tracker.HitData (rowData);
373373 GPUglobalref () const calink* firsthit = tracker.FirstHitInBin (rowData);
@@ -421,7 +421,7 @@ GPUdii() float GPUTPCGMTrackParam::FindBestInterpolatedHit(GPUTPCGMMerger& GPUre
421421
422422 if (dy * dy < sy2 && dz * dz < sz2) {
423423 float err2YA, err2ZA;
424- const ClusterNative& GPUrestrict () cl = merger.GetConstantMem ()-> ioPtrs .clustersNative ->clustersLinear [idOffset + ids[ih]];
424+ const ClusterNative& GPUrestrict () cl = merger.GetConstantMem () -> ioPtrs.clustersNative ->clustersLinear [idOffset + ids[ih]];
425425 const auto clflags = cl.getFlags () & GPUTPCGMMergedTrackHit::clustererAndSharedFlags;
426426 const float time = cl.getTime ();
427427 const float invSqrtCharge = CAMath::InvSqrt (cl.qMax );
@@ -455,13 +455,18 @@ GPUdii() float GPUTPCGMTrackParam::FindBestInterpolatedHit(GPUTPCGMMerger& GPUre
455455 }
456456 }
457457 }
458- CADEBUG (const auto * dbgCand = &merger.ClusterCandidates ()[(iTrk * GPUCA_ROW_COUNT + row) * param.rec .tpc .rebuildTrackInFitClusterCandidates ]; for (int dbgi = 0 ; dbgi < nCandidates; dbgi++) { if (dbgCand[dbgi].id > 1 ) printf (" iTrk %d Row %d Hit %d Candidate %d hit %d err %f\n " , iTrk, (int )row, cluster.num , dbgi, dbgCand[dbgi].id - 2 , dbgCand[dbgi].error ); else break ; }});
459- }
460- if (nCandidates == 0 ) {
461- merger.ClusterCandidates ()[(iTrk * GPUCA_ROW_COUNT + row) * param.rec .tpc .rebuildTrackInFitClusterCandidates + 0 ].id = 1 ;
462- }
458+ CADEBUG (const auto * dbgCand = &merger.ClusterCandidates ()[(iTrk * GPUCA_ROW_COUNT + row) * param.rec .tpc .rebuildTrackInFitClusterCandidates ]; for (int dbgi = 0 ; dbgi < nCandidates; dbgi++) {
459+ if (dbgCand[dbgi].id > 1 )
460+ printf (" iTrk %d Row %d Hit %d Candidate %d hit %d err %f\n " , iTrk, (int )row, cluster.num , dbgi, dbgCand[dbgi].id - 2 , dbgCand[dbgi].error );
461+ else
462+ break ; }
463+ });
463464 }
464- return uncorrectedY;
465+ if (nCandidates == 0 ) {
466+ merger.ClusterCandidates ()[(iTrk * GPUCA_ROW_COUNT + row) * param.rec .tpc .rebuildTrackInFitClusterCandidates + 0 ].id = 1 ;
467+ }
468+ }
469+ return uncorrectedY;
465470}
466471
467472GPUdii () void GPUTPCGMTrackParam::FitdEdx(GPUdEdx& GPUrestrict () dEdx, GPUdEdx& GPUrestrict() dEdxAlt, const GPUParam& GPUrestrict() param, bool finalFit, const GPUCalibObjectsConst& GPUrestrict() calib, int ihit, int ihitMergeFirst, int wayDirection, const ClusterNative* GPUrestrict() clustersArray, const GPUTPCGMMergedTrackHit* GPUrestrict() clusters, uint8_t clusterState, float zz)
@@ -573,7 +578,7 @@ GPUd() int32_t GPUTPCGMTrackParam::MergeDoubleRowClusters(int32_t& ihit, int32_t
573578 xx = yy = zz = 0 .f ;
574579 clusterState = 0 ;
575580 while (true ) {
576- const ClusterNative& GPUrestrict () cl = merger.GetConstantMem ()-> ioPtrs .clustersNative ->clustersLinear [clusters[ihit].num ];
581+ const ClusterNative& GPUrestrict () cl = merger.GetConstantMem () -> ioPtrs.clustersNative ->clustersLinear [clusters[ihit].num ];
577582 float clamp = cl.qTot ;
578583 float clx, cly, clz;
579584 merger.GetConstantMem ()->calibObjects .fastTransformHelper ->Transform (clusters[ihit].sector , clusters[ihit].row , cl.getPad (), cl.getTime (), clx, cly, clz, mTOffset );
@@ -625,7 +630,7 @@ GPUd() float GPUTPCGMTrackParam::AttachClusters(const GPUTPCGMMerger& GPUrestric
625630 if (merger.Param ().rec .tpc .disableRefitAttachment & 1 ) {
626631 return -1e6f;
627632 }
628- const GPUTPCTracker& GPUrestrict () tracker = *(merger.GetConstantMem ()-> tpcTrackers + sector);
633+ const GPUTPCTracker& GPUrestrict () tracker = *(merger.GetConstantMem () -> tpcTrackers + sector);
629634 const GPUTPCRow& GPUrestrict () row = tracker.Row (iRow);
630635 GPUglobalref () const cahit2* hits = tracker.HitData (row);
631636 GPUglobalref () const calink* firsthit = tracker.FirstHitInBin (row);
@@ -920,7 +925,7 @@ GPUd() float GPUTPCGMTrackParam::ShiftZ(const GPUTPCGMMergedTrackHit* GPUrestric
920925 if (N == 0 ) {
921926 N = 1 ;
922927 }
923- const auto & GPUrestrict () cls = merger.GetConstantMem ()-> ioPtrs .clustersNative ->clustersLinear ;
928+ const auto & GPUrestrict () cls = merger.GetConstantMem () -> ioPtrs.clustersNative ->clustersLinear ;
924929 float z0 = cls[clusters[0 ].num ].getTime (), zn = cls[clusters[N - 1 ].num ].getTime ();
925930 const auto tmp = zn > z0 ? std::array<float , 3 >{zn, z0, GPUTPCGeometry::Row2X (clusters[N - 1 ].row )} : std::array<float , 3 >{z0, zn, GPUTPCGeometry::Row2X (clusters[0 ].row )};
926931 return ShiftZ (merger, clusters[0 ].sector , tmp[0 ], tmp[1 ], tmp[2 ]);
0 commit comments