diff --git a/PWGCF/EbyEFluctuations/Tasks/FactorialMomentsTask.cxx b/PWGCF/EbyEFluctuations/Tasks/FactorialMomentsTask.cxx index 6e39fc150ff..f91bfb86088 100644 --- a/PWGCF/EbyEFluctuations/Tasks/FactorialMomentsTask.cxx +++ b/PWGCF/EbyEFluctuations/Tasks/FactorialMomentsTask.cxx @@ -41,6 +41,9 @@ using namespace o2::framework::expressions; TH1D* tmpFqErr[6][5][52]; struct FactorialMomentsTask { + Configurable useITS{"useITS", false, "Select tracks with ITS"}; + Configurable useTPC{"useTPC", false, "Select tracks with TPC"}; + Configurable useGlobal{"useGlobal", true, "Select global tracks"}; Configurable applyCheckPtForRec{"applyCheckPtForRec", false, "Apply checkpT for reconstructed tracks"}; Configurable applyCheckPtForMC{"applyCheckPtForMC", true, "Apply checkpT for MC-generated tracks"}; Configurable centralEta{"centralEta", 0.9, "eta limit for tracks"}; @@ -402,9 +405,12 @@ struct FactorialMomentsTask { fqEvent = {{{{{0, 0, 0, 0, 0, 0}}}}}; binConEvent = {{{0, 0, 0, 0, 0}}}; for (auto const& track : colltracks) { - // if (track.hasITS()) - // if (track.hasTPC()) - // if (track.isGlobalTrack()) { + if (useITS && !track.hasITS()) + continue; + if (useTPC && !track.hasTPC()) + continue; + if (useGlobal && !track.isGlobalTrack()) + continue; histos.fill(HIST("mCollID"), track.collisionId()); histos.fill(HIST("mEta"), track.eta()); histos.fill(HIST("mPt"), track.pt()); @@ -532,8 +538,6 @@ struct FactorialMomentsTask { for (int iPt = 0; iPt < numPt; ++iPt) { if (countTracks[iPt] > 0) { mHistArrQA[iPt * 4 + 3]->Fill(countTracks[iPt]); - } else { - return; } } @@ -589,8 +593,6 @@ struct FactorialMomentsTask { for (int iPt = 0; iPt < numPt; ++iPt) { if (countTracks[iPt] > 0) { mHistArrQA[iPt * 4 + 3]->Fill(countTracks[iPt]); - } else { - return; } } // Calculate the normalized factorial moments