From a61dcdf7894201ffcf917dd9a9d396929ede4aa9 Mon Sep 17 00:00:00 2001 From: jikim1290 Date: Wed, 18 Feb 2026 19:34:17 +0900 Subject: [PATCH 1/4] adding track selection and correction --- PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx | 7 +++++++ 1 file changed, 7 insertions(+) diff --git a/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx b/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx index 3e495a3bc42..454cf399142 100644 --- a/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx +++ b/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx @@ -334,6 +334,13 @@ struct jEPFlowAnalysis { epFlowHistograms.fill(HIST("EpResQvecRefARefBxy"), i + 2, cent, qx_shifted[2] * qy_shifted[1] - qx_shifted[1] * qy_shifted[2]); for (const auto& track : tracks) { + if (trackSel(track)) + continue; + + if (cfgEffCor) { + weight /= getEfficiencyCorrection(effMap, track.eta(), track.pt(), cent, coll.posZ()); + } + float vn = std::cos((i + 2) * (track.phi() - eps[0])); float vnSin = std::sin((i + 2) * (track.phi() - eps[0])); From 10e37705bb8f023b67b6dd2752d258fa85f05d88 Mon Sep 17 00:00:00 2001 From: jikim1290 Date: Wed, 18 Feb 2026 20:01:05 +0900 Subject: [PATCH 2/4] adding track selection and correction --- PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx b/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx index 454cf399142..7dbdda44d6d 100644 --- a/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx +++ b/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx @@ -344,11 +344,11 @@ struct jEPFlowAnalysis { float vn = std::cos((i + 2) * (track.phi() - eps[0])); float vnSin = std::sin((i + 2) * (track.phi() - eps[0])); - epFlowHistograms.fill(HIST("vncos"), i + 2, cent, track.pt(), vn * weight); - epFlowHistograms.fill(HIST("vnsin"), i + 2, cent, track.pt(), vnSin * weight); + epFlowHistograms.fill(HIST("vncos"), i + 2, cent, track.pt(), vn, weight); + epFlowHistograms.fill(HIST("vnsin"), i + 2, cent, track.pt(), vnSin, weight); - epFlowHistograms.fill(HIST("SPvnxx"), i + 2, cent, track.pt(), (std::cos(track.phi() * static_cast(i + 2)) * qx_shifted[0] + std::sin(track.phi() * static_cast(i + 2)) * qy_shifted[0]) * weight); - epFlowHistograms.fill(HIST("SPvnxy"), i + 2, cent, track.pt(), (std::sin(track.phi() * static_cast(i + 2)) * qx_shifted[0] - std::cos(track.phi() * static_cast(i + 2)) * qy_shifted[0]) * weight); + epFlowHistograms.fill(HIST("SPvnxx"), i + 2, cent, track.pt(), (std::cos(track.phi() * static_cast(i + 2)) * qx_shifted[0] + std::sin(track.phi() * static_cast(i + 2)) * qy_shifted[0]), weight); + epFlowHistograms.fill(HIST("SPvnxy"), i + 2, cent, track.pt(), (std::sin(track.phi() * static_cast(i + 2)) * qx_shifted[0] - std::cos(track.phi() * static_cast(i + 2)) * qy_shifted[0]), weight); } } } From f32a0d0051621230b48d00afe71d38561ca56bf2 Mon Sep 17 00:00:00 2001 From: jikim1290 Date: Wed, 18 Feb 2026 22:37:51 +0900 Subject: [PATCH 3/4] adding track selection and correction --- PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx b/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx index 7dbdda44d6d..3b1bf3ab8ec 100644 --- a/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx +++ b/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx @@ -249,6 +249,11 @@ struct jEPFlowAnalysis { return; } + if (cfgEffCor) { + auto bc = coll.bc_as(); + effMap = ccdb->getForTimeStamp>(cfgEffCorDir, bc.timestamp()); + } + float cent = coll.cent(); epFlowHistograms.fill(HIST("hCentrality"), cent); epFlowHistograms.fill(HIST("hVertex"), coll.posZ()); From fb994ea4048d84ece390f84b8c650e673af3270a Mon Sep 17 00:00:00 2001 From: jikim1290 Date: Thu, 19 Feb 2026 00:51:40 +0900 Subject: [PATCH 4/4] adding track selection and correction --- PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx | 7 +++++++ 1 file changed, 7 insertions(+) diff --git a/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx b/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx index 3b1bf3ab8ec..79c5c4e2df6 100644 --- a/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx +++ b/PWGCF/JCorran/Tasks/jEPFlowAnalysis.cxx @@ -292,6 +292,13 @@ struct jEPFlowAnalysis { float weight = 1.0; + qx_shifted[0] = coll.qvecRe()[4 * detId + harmInd]; + qy_shifted[0] = coll.qvecIm()[4 * detId + harmInd]; + qx_shifted[1] = coll.qvecRe()[4 * refAId + harmInd]; + qy_shifted[1] = coll.qvecIm()[4 * refAId + harmInd]; + qx_shifted[2] = coll.qvecRe()[4 * refBId + harmInd]; + qy_shifted[2] = coll.qvecIm()[4 * refBId + harmInd]; + if (cfgManShiftCorr) { constexpr int kShiftBins = 10; for (int ishift = 1; ishift <= kShiftBins; ishift++) {