@@ -365,11 +365,11 @@ void* GPUTPCGMMerger::SetPointersRefitScratch(void* mem)
365365
366366void * GPUTPCGMMerger::SetPointersOutput (void * mem)
367367{
368- computePointerWithAlignment (mem, mOutputTracks , mNMaxTracks );
368+ computePointerWithAlignment (mem, mMergedTracks , mNMaxTracks );
369369 if (mRec ->GetParam ().dodEdxEnabled ) {
370- computePointerWithAlignment (mem, mOutputTracksdEdx , mNMaxTracks );
370+ computePointerWithAlignment (mem, mMergedTracksdEdx , mNMaxTracks );
371371 if (mRec ->GetParam ().rec .tpc .dEdxClusterRejectionFlagMask != mRec ->GetParam ().rec .tpc .dEdxClusterRejectionFlagMaskAlt ) {
372- computePointerWithAlignment (mem, mOutputTracksdEdxAlt , mNMaxTracks );
372+ computePointerWithAlignment (mem, mMergedTracksdEdxAlt , mNMaxTracks );
373373 }
374374 }
375375 computePointerWithAlignment (mem, mClusters , mNMaxOutputTrackClusters );
@@ -1318,7 +1318,7 @@ GPUd() void GPUTPCGMMerger::MergeCEFill(const GPUTPCGMSectorTrack* track, const
13181318 const float x0 = GPUTPCGeometry::Row2X (attempt == 0 ? 63 : cls.row );
13191319 if (track->TransportToX (this , x0, Param ().bzCLight , b, GPUCA_MAX_SIN_PHI_LOW)) {
13201320 b.SetTrackID (itr);
1321- b.SetNClusters (mOutputTracks [itr].NClusters ());
1321+ b.SetNClusters (mMergedTracks [itr].NClusters ());
13221322 if (CAMath::Abs (b.Cov ()[4 ]) >= 0 .5f ) {
13231323 b.SetCov (4 , 0 .5f ); // TODO: Is this needed and better than the cut in BorderTrack?
13241324 }
@@ -1339,11 +1339,11 @@ GPUd() void GPUTPCGMMerger::MergeCE(int32_t nBlocks, int32_t nThreads, int32_t i
13391339{
13401340 const ClusterNative* cls = Param ().par .earlyTpcTransform ? nullptr : mConstantMem ->ioPtrs .clustersNative ->clustersLinear ;
13411341 for (uint32_t i = iBlock * nThreads + iThread; i < mMemory ->nMergedTracks ; i += nThreads * nBlocks) {
1342- if (mOutputTracks [i].CSide () == 0 && mTrackLinks [i] >= 0 ) {
1342+ if (mMergedTracks [i].CSide () == 0 && mTrackLinks [i] >= 0 ) {
13431343 if (mTrackLinks [mTrackLinks [i]] != (int32_t )i) {
13441344 continue ;
13451345 }
1346- GPUTPCGMMergedTrack* trk[2 ] = {&mOutputTracks [i], &mOutputTracks [mTrackLinks [i]]};
1346+ GPUTPCGMMergedTrack* trk[2 ] = {&mMergedTracks [i], &mMergedTracks [mTrackLinks [i]]};
13471347
13481348 if (!trk[1 ]->OK () || trk[1 ]->CCE ()) {
13491349 continue ;
@@ -1459,7 +1459,7 @@ GPUd() void GPUTPCGMMerger::MergeCE(int32_t nBlocks, int32_t nThreads, int32_t i
14591459 }
14601460 }
14611461
1462- // for (int32_t i = 0;i < mMemory->nMergedTracks;i++) {if (mOutputTracks [i].CCE() == false) {mOutputTracks [i].SetNClusters(0);mOutputTracks [i].SetOK(false);}} //Remove all non-CE tracks
1462+ // for (int32_t i = 0;i < mMemory->nMergedTracks;i++) {if (mMergedTracks [i].CCE() == false) {mMergedTracks [i].SetNClusters(0);mMergedTracks [i].SetOK(false);}} //Remove all non-CE tracks
14631463}
14641464
14651465namespace o2 ::gpu::internal
@@ -1752,7 +1752,7 @@ GPUd() void GPUTPCGMMerger::CollectMergedTracks(int32_t nBlocks, int32_t nThread
17521752 continue ;
17531753 }
17541754
1755- GPUTPCGMMergedTrack& mergedTrack = mOutputTracks [iOutputTrack];
1755+ GPUTPCGMMergedTrack& mergedTrack = mMergedTracks [iOutputTrack];
17561756
17571757 mergedTrack.SetFlags (0 );
17581758 mergedTrack.SetOK (1 );
@@ -1825,7 +1825,7 @@ GPUd() void GPUTPCGMMerger::SortTracks(int32_t nBlocks, int32_t nThreads, int32_
18251825{
18261826#ifndef GPUCA_SPECIALIZE_THRUST_SORTS
18271827 if (iThread == 0 && iBlock == 0 ) {
1828- GPUCommonAlgorithm::sortDeviceDynamic (mTrackOrderProcess , mTrackOrderProcess + mMemory ->nMergedTracks , GPUTPCGMMergerSortTracks_comp (mOutputTracks ));
1828+ GPUCommonAlgorithm::sortDeviceDynamic (mTrackOrderProcess , mTrackOrderProcess + mMemory ->nMergedTracks , GPUTPCGMMergerSortTracks_comp (mMergedTracks ));
18291829 }
18301830#endif
18311831}
@@ -1834,7 +1834,7 @@ GPUd() void GPUTPCGMMerger::SortTracksQPt(int32_t nBlocks, int32_t nThreads, int
18341834{
18351835#ifndef GPUCA_SPECIALIZE_THRUST_SORTS
18361836 if (iThread == 0 && iBlock == 0 ) {
1837- GPUCommonAlgorithm::sortDeviceDynamic (mTrackSort , mTrackSort + mMemory ->nMergedTracks , GPUTPCGMMergerSortTracksQPt_comp (mOutputTracks ));
1837+ GPUCommonAlgorithm::sortDeviceDynamic (mTrackSort , mTrackSort + mMemory ->nMergedTracks , GPUTPCGMMergerSortTracksQPt_comp (mMergedTracks ));
18381838 }
18391839#endif
18401840}
@@ -1843,7 +1843,7 @@ GPUd() void GPUTPCGMMerger::PrepareClustersForFit1(int32_t nBlocks, int32_t nThr
18431843{
18441844 for (uint32_t i = iBlock * nThreads + iThread; i < mMemory ->nMergedTracks ; i += nBlocks * nThreads) {
18451845 mTrackOrderAttach [mTrackSort [i]] = i;
1846- const GPUTPCGMMergedTrack& trk = mOutputTracks [i];
1846+ const GPUTPCGMMergedTrack& trk = mMergedTracks [i];
18471847 if (trk.OK ()) {
18481848 for (uint32_t j = 0 ; j < trk.NClusters (); j++) {
18491849 mClusterAttachment [mClusters [trk.FirstClusterRef () + j].num ] = attachAttached | attachGood;
@@ -1884,7 +1884,7 @@ GPUd() void GPUTPCGMMerger::Finalize0(int32_t nBlocks, int32_t nThreads, int32_t
18841884GPUd () void GPUTPCGMMerger::Finalize1(int32_t nBlocks, int32_t nThreads, int32_t iBlock, int32_t iThread)
18851885{
18861886 for (uint32_t i = iBlock * nThreads + iThread; i < mMemory ->nMergedTracks ; i += nThreads * nBlocks) {
1887- const GPUTPCGMMergedTrack& trk = mOutputTracks [i];
1887+ const GPUTPCGMMergedTrack& trk = mMergedTracks [i];
18881888 if (!trk.OK () || trk.NClusters () == 0 ) {
18891889 continue ;
18901890 }
@@ -1919,7 +1919,7 @@ GPUd() void GPUTPCGMMerger::MergeLoopersInit(int32_t nBlocks, int32_t nThreads,
19191919{
19201920 const float lowPtThresh = Param ().rec .tpc .rejectQPtB5 * 1 .1f ; // Might need to merge tracks above the threshold with parts below the threshold
19211921 for (uint32_t i = get_global_id (0 ); i < mMemory ->nMergedTracks ; i += get_global_size (0 )) {
1922- const auto & trk = mOutputTracks [i];
1922+ const auto & trk = mMergedTracks [i];
19231923 const auto & p = trk.GetParam ();
19241924 const float qptabs = CAMath::Abs (p.GetQPt ());
19251925 if (trk.NClusters () && qptabs * Param ().qptB5Scaler > 5 .f && qptabs * Param ().qptB5Scaler <= lowPtThresh) {
@@ -1983,7 +1983,7 @@ GPUd() void GPUTPCGMMerger::MergeLoopersMain(int32_t nBlocks, int32_t nThreads,
19831983#if GPUCA_MERGE_LOOPER_MC && !defined(GPUCA_GPUCODE)
19841984 std::vector<int64_t > paramLabels (mMemory ->nLooperMatchCandidates );
19851985 for (uint32_t i = 0 ; i < mMemory ->nLooperMatchCandidates ; i++) {
1986- paramLabels[i] = GetTrackLabel (mOutputTracks [params[i].id ]);
1986+ paramLabels[i] = GetTrackLabel (mMergedTracks [params[i].id ]);
19871987 }
19881988 /* std::vector<bool> dropped(mMemory->nLooperMatchCandidates);
19891989 std::vector<bool> droppedMC(mMemory->nLooperMatchCandidates);
@@ -2005,8 +2005,8 @@ GPUd() void GPUTPCGMMerger::MergeLoopersMain(int32_t nBlocks, int32_t nThreads,
20052005 // bs |= 1;
20062006 continue ;
20072007 }
2008- const auto & trk1 = mOutputTracks [params[i].id ];
2009- const auto & trk2 = mOutputTracks [params[j].id ];
2008+ const auto & trk1 = mMergedTracks [params[i].id ];
2009+ const auto & trk2 = mMergedTracks [params[j].id ];
20102010 const auto & param1 = trk1.GetParam ();
20112011 const auto & param2 = trk2.GetParam ();
20122012 if (CAMath::Abs (param1.GetDzDs ()) > 0 .03f && CAMath::Abs (param2.GetDzDs ()) > 0 .03f && param1.GetDzDs () * param2.GetDzDs () * param1.GetQPt () * param2.GetQPt () < 0 ) {
@@ -2045,7 +2045,7 @@ GPUd() void GPUTPCGMMerger::MergeLoopersMain(int32_t nBlocks, int32_t nThreads,
20452045 const int64_t label2 = paramLabels[j];
20462046 bool labelEQ = label1 != -1 && label1 == label2;
20472047 if (1 || EQ || labelEQ) {
2048- // printf("Matching track %d/%d %u-%u (%ld/%ld): dist %f side %d %d, tgl %f %f, qpt %f %f, x %f %f, y %f %f\n", (int32_t)EQ, (int32_t)labelEQ, i, j, label1, label2, d, (int32_t)mOutputTracks [params[i].id].CSide(), (int32_t)mOutputTracks [params[j].id].CSide(), params[i].tgl, params[j].tgl, params[i].qpt, params[j].qpt, params[i].x, params[j].x, params[i].y, params[j].y);
2048+ // printf("Matching track %d/%d %u-%u (%ld/%ld): dist %f side %d %d, tgl %f %f, qpt %f %f, x %f %f, y %f %f\n", (int32_t)EQ, (int32_t)labelEQ, i, j, label1, label2, d, (int32_t)mMergedTracks [params[i].id].CSide(), (int32_t)mMergedTracks [params[j].id].CSide(), params[i].tgl, params[j].tgl, params[i].qpt, params[j].qpt, params[i].x, params[j].x, params[i].y, params[j].y);
20492049 static auto & tup = GPUROOTDump<TNtuple>::get (" mergeloopers" , " labeleq:sides:d2xy:tgl1:tgl2:qpt1:qpt2:dz:dzcorr:dtgl:dqpt:dznorm:bs" );
20502050 tup.Fill ((float )labelEQ, (trk1.CSide () ? 1 : 0 ) | (trk2.CSide () ? 2 : 0 ), d2xy, param1.GetDzDs (), param2.GetDzDs (), param1.GetQPt (), param2.GetQPt (), CAMath::Abs (params[j].refz ) - CAMath::Abs (params[i].refz ), dzcorr, dtgl, dqpt, dznorm, bs);
20512051 static auto tup2 = GPUROOTDump<TNtuple>::getNew (" mergeloopers2" , " labeleq:refz1:refz2:tgl1:tgl2:qpt1:qpt2:snp1:snp2:a1:a2:dzn:phasecor:phasedir:dzcorr" );
@@ -2063,9 +2063,9 @@ GPUd() void GPUTPCGMMerger::MergeLoopersMain(int32_t nBlocks, int32_t nThreads,
20632063 }*/
20642064#endif
20652065 if (EQ) {
2066- mOutputTracks [params[j].id ].SetMergedLooper (true );
2066+ mMergedTracks [params[j].id ].SetMergedLooper (true );
20672067 if (CAMath::Abs (param2.GetQPt () * Param ().qptB5Scaler ) >= Param ().rec .tpc .rejectQPtB5 ) {
2068- mOutputTracks [params[i].id ].SetMergedLooper (true );
2068+ mMergedTracks [params[i].id ].SetMergedLooper (true );
20692069 }
20702070 }
20712071 }
0 commit comments