From 8b3152e6557b76147b03af4bb81cd556338fb022 Mon Sep 17 00:00:00 2001 From: huinaibing Date: Tue, 28 Apr 2026 19:12:54 +0800 Subject: [PATCH 1/3] [PWGCF] Add additional hists --- PWGCF/Flow/Tasks/pidFlowPtCorr.cxx | 21 ++++++++++++++++----- 1 file changed, 16 insertions(+), 5 deletions(-) diff --git a/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx b/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx index a19e1bc506e..e954abc3105 100644 --- a/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx +++ b/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx @@ -58,10 +58,11 @@ #include +#include + #include #include #include -#include #include #include #include @@ -556,6 +557,8 @@ struct PidFlowPtCorr { TObjArray* oba4PID = reinterpret_cast(oba4Ch->Clone()); oba4PID->Add(new TNamed("c22pure", "c22pure")); oba4PID->Add(new TNamed("c32pure", "c32pure")); + oba4PID->Add(new TNamed("covV2PtPID", "covV2PtPID")); + oba4PID->Add(new TNamed("c22TrackWeightPID", "c22TrackWeightPID")); fFCPi->SetName("FlowContainerPi"); fFCPi->Initialize(oba4PID, axisMultiplicity, cfgFlowNbootstrap); @@ -889,7 +892,6 @@ struct PidFlowPtCorr { const float ptMax = pidPtRangeOpts.cfgPtMax4TOFPiKa.value; if (pt > ptMin && pt < ptMax) { - isPionBool = std::hypot(tpcNsigma, tofNsigma) < circleCutOpts.cfgCircleCutSigmaPi.value; } else { // Fallback: TPC only cut @@ -2241,14 +2243,23 @@ struct PidFlowPtCorr { fFCCh->FillProfile("hMeanPt", cent, (ptSum / nch), nch, rndm); - if (nPionWeighted > 0) + if (nPionWeighted > 0) { fFCPi->FillProfile("hMeanPt", cent, (pionPtSum / nPionWeighted), nPionWeighted, rndm); + fillFCvnpt(MyParticleType::kPion, corrconfigs.at(29), cent, rndm, pionPtSum, nPionWeighted, "covV2PtPID"); + fillFCvnpt(MyParticleType::kPion, corrconfigs.at(29), cent, rndm, nPionWeighted, nPionWeighted, "c22TrackWeightPID"); + } - if (nKaonWeighted > 0) + if (nKaonWeighted > 0) { fFCKa->FillProfile("hMeanPt", cent, (kaonPtSum / nKaonWeighted), nKaonWeighted, rndm); + fillFCvnpt(MyParticleType::kKaon, corrconfigs.at(30), cent, rndm, kaonPtSum, nKaonWeighted, "covV2PtPID"); + fillFCvnpt(MyParticleType::kKaon, corrconfigs.at(30), cent, rndm, nKaonWeighted, nKaonWeighted, "c22TrackWeightPID"); + } - if (nProtonWeighted > 0) + if (nProtonWeighted > 0) { fFCPr->FillProfile("hMeanPt", cent, (protonPtSum / nProtonWeighted), nProtonWeighted, rndm); + fillFCvnpt(MyParticleType::kProton, corrconfigs.at(31), cent, rndm, protonPtSum, nProtonWeighted, "covV2PtPID"); + fillFCvnpt(MyParticleType::kProton, corrconfigs.at(31), cent, rndm, nProtonWeighted, nProtonWeighted, "c22TrackWeightPID"); + } double nchDiff = nch * nch - nchSquare; if (nchDiff > minVal4Float) { From e040241f5a5090c5ad06c5d5f8e2b9adb4b1d489 Mon Sep 17 00:00:00 2001 From: huinaibing Date: Tue, 28 Apr 2026 19:17:46 +0800 Subject: [PATCH 2/3] Fix extra modifications caused by overlooking the updated header files --- PWGCF/Flow/Tasks/pidFlowPtCorr.cxx | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx b/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx index e954abc3105..7676b4ecb9c 100644 --- a/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx +++ b/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx @@ -58,11 +58,10 @@ #include -#include - #include #include #include +#include #include #include #include From 4c1904dcf3c8309c2795ddbbd1f61e8654317e55 Mon Sep 17 00:00:00 2001 From: huinaibing Date: Tue, 28 Apr 2026 21:56:19 +0800 Subject: [PATCH 3/3] replace reinterpret_cast --- PWGCF/Flow/Tasks/pidFlowPtCorr.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx b/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx index 7676b4ecb9c..69cf73a8a16 100644 --- a/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx +++ b/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx @@ -553,7 +553,7 @@ struct PidFlowPtCorr { // init fFCPID // note that need to add c22pure and c32pure - TObjArray* oba4PID = reinterpret_cast(oba4Ch->Clone()); + TObjArray* oba4PID = dynamic_cast(oba4Ch->Clone()); oba4PID->Add(new TNamed("c22pure", "c22pure")); oba4PID->Add(new TNamed("c32pure", "c32pure")); oba4PID->Add(new TNamed("covV2PtPID", "covV2PtPID"));