@@ -348,30 +348,35 @@ struct FactorialMomentsTask {
348348 fqEvent = {{{{{0 , 0 , 0 , 0 , 0 , 0 }}}}};
349349 binConEvent = {{{0 , 0 , 0 , 0 , 0 }}};
350350 for (auto const & track : tracks) {
351- if (track.hasTPC ()) {
352- histos.fill (HIST (" mCollID" ), track.collisionId ());
353- histos.fill (HIST (" mEta" ), track.eta ());
354- histos.fill (HIST (" mPt" ), track.pt ());
355- histos.fill (HIST (" mPhi" ), track.phi ());
356- histos.fill (HIST (" mNFindableClsTPC" ), track.tpcNClsFindable ());
357- histos.fill (HIST (" mNClsTPC" ), track.tpcNClsFound ());
358- histos.fill (HIST (" mNClsITS" ), track.itsNCls ());
359- histos.fill (HIST (" mChi2TPC" ), track.tpcChi2NCl ());
360- histos.fill (HIST (" mChi2ITS" ), track.itsChi2NCl ());
361- histos.fill (HIST (" mChi2TRD" ), track.trdChi2 ());
362- histos.fill (HIST (" mDCAxy" ), track.dcaXY ());
363- histos.fill (HIST (" mDCAx" ), track.dcaZ ());
364- histos.fill (HIST (" mDCAxyPt" ), track.pt (), track.dcaXY ());
365- histos.fill (HIST (" mDCAzPt" ), track.pt (), track.dcaZ ());
366- histos.fill (HIST (" mNSharedClsTPC" ), track.tpcNClsShared ());
367- histos.fill (HIST (" mCrossedRowsTPC" ), track.tpcNClsCrossedRows ());
368- histos.fill (HIST (" mNFinClsminusCRows" ), track.tpcNClsFindableMinusCrossedRows ());
369- histos.fill (HIST (" mNFractionShClsTPC" ), track.tpcFractionSharedCls ());
370- histos.fill (HIST (" mSharedClsvsPt" ), track.pt (), track.tpcNClsShared ());
371- histos.fill (HIST (" mSharedClsProbvsPt" ), track.pt (), track.tpcFractionSharedCls () / track.tpcNClsCrossedRows ());
372- checkpT (track);
373- }
351+ if (useITS && !track.hasITS ())
352+ continue ;
353+ if (useTPC && !track.hasTPC ())
354+ continue ;
355+ if (useGlobal && !track.isGlobalTrack ())
356+ continue ;
357+ histos.fill (HIST (" mCollID" ), track.collisionId ());
358+ histos.fill (HIST (" mEta" ), track.eta ());
359+ histos.fill (HIST (" mPt" ), track.pt ());
360+ histos.fill (HIST (" mPhi" ), track.phi ());
361+ histos.fill (HIST (" mNFindableClsTPC" ), track.tpcNClsFindable ());
362+ histos.fill (HIST (" mNClsTPC" ), track.tpcNClsFound ());
363+ histos.fill (HIST (" mNClsITS" ), track.itsNCls ());
364+ histos.fill (HIST (" mChi2TPC" ), track.tpcChi2NCl ());
365+ histos.fill (HIST (" mChi2ITS" ), track.itsChi2NCl ());
366+ histos.fill (HIST (" mChi2TRD" ), track.trdChi2 ());
367+ histos.fill (HIST (" mDCAxy" ), track.dcaXY ());
368+ histos.fill (HIST (" mDCAx" ), track.dcaZ ());
369+ histos.fill (HIST (" mDCAxyPt" ), track.pt (), track.dcaXY ());
370+ histos.fill (HIST (" mDCAzPt" ), track.pt (), track.dcaZ ());
371+ histos.fill (HIST (" mNSharedClsTPC" ), track.tpcNClsShared ());
372+ histos.fill (HIST (" mCrossedRowsTPC" ), track.tpcNClsCrossedRows ());
373+ histos.fill (HIST (" mNFinClsminusCRows" ), track.tpcNClsFindableMinusCrossedRows ());
374+ histos.fill (HIST (" mNFractionShClsTPC" ), track.tpcFractionSharedCls ());
375+ histos.fill (HIST (" mSharedClsvsPt" ), track.pt (), track.tpcNClsShared ());
376+ histos.fill (HIST (" mSharedClsProbvsPt" ), track.pt (), track.tpcFractionSharedCls () / track.tpcNClsCrossedRows ());
377+ checkpT (track);
374378 }
379+
375380 for (int iPt = 0 ; iPt < numPt; ++iPt) {
376381 if (countTracks[iPt] > 0 ) {
377382 mHistArrQA [iPt * 4 + 3 ]->Fill (countTracks[iPt]);
@@ -470,13 +475,11 @@ struct FactorialMomentsTask {
470475 }
471476 }
472477 for (auto iPt = 0 ; iPt < numPt; ++iPt) {
473- // if (countTracks[iPt] > 0)countTracks = {0, 0, 0, 0, 0};
474478 if (countTracks[iPt] > 0 ) {
475479 mHistArrQA [iPt * 4 + 3 ]->Fill (countTracks[iPt]);
476480 }
477481 }
478482 histos.fill (HIST (" mEventSelected" ), 6 );
479- // Calculate the normalized factorial moments
480483 calculateMoments (mHistArrReset );
481484 }
482485 PROCESS_SWITCH (FactorialMomentsTask, processMCRec, " main process function" , false );
0 commit comments