From 535544a347c4f096e8624f7f213627fa9e0018f8 Mon Sep 17 00:00:00 2001 From: huinaibing Date: Sun, 14 Jun 2026 22:51:51 +0800 Subject: [PATCH] [PWGCF] Add more graph with different weight --- PWGCF/Flow/Tasks/pidFlowPtCorr.cxx | 15 +++++++++++++++ 1 file changed, 15 insertions(+) diff --git a/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx b/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx index d779826c5fd..c09a098c529 100644 --- a/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx +++ b/PWGCF/Flow/Tasks/pidFlowPtCorr.cxx @@ -557,6 +557,8 @@ struct PidFlowPtCorr { oba4Ch->Add(new TNamed("hMeanPtWeightOne", "hMeanPtWeightOne")); oba4Ch->Add(new TNamed("ptAveWeightOne", "ptAveWeightOne")); oba4Ch->Add(new TNamed("ptSquareAveWeightOne", "ptSquareAveWeightOne")); + + oba4Ch->Add(new TNamed("hMeanPtWeightFull", "hMeanPtWeightFull")); // end fill TObjArray for charged // init fFCCh @@ -572,6 +574,8 @@ struct PidFlowPtCorr { oba4PID->Add(new TNamed("covV2PtPID", "covV2PtPID")); oba4PID->Add(new TNamed("c22TrackWeightPID", "c22TrackWeightPID")); + oba4PID->Add(new TNamed("hMeanPtWeightCharged", "hMeanPtWeightCharged")); + fFCPi->SetName("FlowContainerPi"); fFCPi->Initialize(oba4PID, axisMultiplicity, cfgFlowNbootstrap); @@ -1082,6 +1086,8 @@ struct PidFlowPtCorr { registry.fill(HIST("meanptCentNbs/hCharged"), ptSum / nch, cent, rndm * cfgFlowNbootstrap, val, nch * dnx); registry.fill(HIST("meanptCentNbs/hChargedMeanpt"), ptSum / nch, cent, rndm * cfgFlowNbootstrap, ptSum / nch, 1.); + + fFCCh->FillProfile("hMeanPtWeightFull", cent, (ptSum / nch), nch * dnx, rndm); } /** @@ -1124,6 +1130,9 @@ struct PidFlowPtCorr { registry.fill(HIST("meanptCentNbs/hChargedPionWithNpair"), pidPtSum / nPid, cent, rndm * cfgFlowNbootstrap, val, dnx); registry.fill(HIST("meanptCentNbs/hPionMeanpt"), pidPtSum / nPid, cent, rndm * cfgFlowNbootstrap, pidPtSum / nPid, 1.); + fFCPi->FillProfile("hMeanPtWeightFull", cent, (pidPtSum / nPid), nPid * npairPid, rndm); + fFCPi->FillProfile("hMeanPtWeightCharged", cent, (pidPtSum / nPid), dnx * nPid, rndm); + if (switchsOpts.cfgClosureTest.value != 0) { double npair4c22pure = fGFW->Calculate(corrconfigs.at(29), 0, true).real(); if (npair4c22pure > minVal4Float) @@ -1159,6 +1168,9 @@ struct PidFlowPtCorr { registry.fill(HIST("meanptCentNbs/hChargedKaonWithNpair"), pidPtSum / nPid, cent, rndm * cfgFlowNbootstrap, val, dnx); registry.fill(HIST("meanptCentNbs/hKaonMeanpt"), pidPtSum / nPid, cent, rndm * cfgFlowNbootstrap, pidPtSum / nPid, 1.); + fFCKa->FillProfile("hMeanPtWeightFull", cent, (pidPtSum / nPid), nPid * npairPid, rndm); + fFCKa->FillProfile("hMeanPtWeightCharged", cent, (pidPtSum / nPid), dnx * nPid, rndm); + break; // end kaon @@ -1176,6 +1188,9 @@ struct PidFlowPtCorr { registry.fill(HIST("meanptCentNbs/hChargedProtonWithNpair"), pidPtSum / nPid, cent, rndm * cfgFlowNbootstrap, val, dnx); registry.fill(HIST("meanptCentNbs/hProtonMeanpt"), pidPtSum / nPid, cent, rndm * cfgFlowNbootstrap, pidPtSum / nPid, 1.); + fFCPr->FillProfile("hMeanPtWeightFull", cent, (pidPtSum / nPid), nPid * npairPid, rndm); + fFCPr->FillProfile("hMeanPtWeightCharged", cent, (pidPtSum / nPid), dnx * nPid, rndm); + break; // end proton