Skip to content

Commit 3bf465f

Browse files
authored
[PWGCF] Add some QA plots (#12924)
1 parent 13e76ab commit 3bf465f

File tree

1 file changed

+8
-1
lines changed

1 file changed

+8
-1
lines changed

PWGCF/Flow/Tasks/flowEseTask.cxx

Lines changed: 8 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -226,6 +226,9 @@ struct FlowEseTask {
226226
histos.add(Form("histQvecV2"), "", {HistType::kTH3F, {qvecAxis, qvecAxis, centAxis}});
227227
histos.add(Form("histMult_Cent"), "", {HistType::kTH2F, {multNumAxis, centAxis}});
228228
histos.add(Form("histQvecCent"), "", {HistType::kTH2F, {lowerQAxis, centAxis}});
229+
histos.add(Form("histPrPtCent"), "", {HistType::kTHnSparseF, {ptAxis, ptAxis, ptAxis, centAxis}});
230+
histos.add(Form("histPiPtCent"), "", {HistType::kTHnSparseF, {ptAxis, ptAxis, ptAxis, centAxis}});
231+
histos.add(Form("histPrBoostedPtCent"), "", {HistType::kTHnSparseF, {ptAxis, ptAxis, ptAxis, centAxis}});
229232

230233
for (auto i = 2; i < cfgnMods + 2; i++) {
231234
histos.add(Form("psi%d/h_lambda_cos", i), "", {HistType::kTHnSparseF, {massAxis, ptAxis, cosAxis, centAxis, epAxis}});
@@ -680,7 +683,7 @@ struct FlowEseTask {
680683
qvecRefAInd = refAId * 4 + 3 + (nmode - 2) * cfgNQvec * 4;
681684
qvecRefBInd = refBId * 4 + 3 + (nmode - 2) * cfgNQvec * 4;
682685

683-
for (auto& v0 : V0s) { // o2-linter: disable=const-ref-in-for-loop(need to modify v0)
686+
for (const auto& v0 : V0s) {
684687
auto postrack = v0.template posTrack_as<TrackCandidates>();
685688
auto negtrack = v0.template negTrack_as<TrackCandidates>();
686689

@@ -732,6 +735,10 @@ struct FlowEseTask {
732735
psi = safeATan2(collision.qvecIm()[qvecDetInd], collision.qvecRe()[qvecDetInd]) / static_cast<float>(nmode);
733736
relphi = TVector2::Phi_0_2pi(static_cast<float>(nmode) * (LambdaVec.Phi() - psi));
734737

738+
histos.fill(HIST("histPrPtCent"), protonVec.Px(), protonVec.Py(), protonVec.Pz(), collision.centFT0C());
739+
histos.fill(HIST("histPiPtCent"), pionVec.Px(), pionVec.Py(), pionVec.Pz(), collision.centFT0C());
740+
histos.fill(HIST("histPrBoostedPtCent"), protonBoostedVec.Px(), protonBoostedVec.Py(), protonBoostedVec.Pz(), collision.centFT0C());
741+
735742
if (cfgShiftCorr) {
736743
auto deltapsiFT0C = 0.0;
737744
auto deltapsiFT0A = 0.0;

0 commit comments

Comments
 (0)