Skip to content

[PWGLF] Added QA histograms for phi meson #18180

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
mpuccio merged 1 commit into from
Oct 7, 2026
+62 −0
Merged

[PWGLF] Added QA histograms for phi meson #18180

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
62 changes: 62 additions & 0 deletions PWGLF/Tasks/Resonances/kstarqa.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -173,6 +173,7 @@ struct Kstarqa {
HistogramRegistry hInvMass{"hInvMass", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
HistogramRegistry hPID{"hPID", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
HistogramRegistry hOthers{"hOthers", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
HistogramRegistry hShift{"hShift", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};

// Confugrable for QA histograms
Configurable<bool> calcLikeSign{"calcLikeSign", true, "Calculate Like Sign"};
Expand Down Expand Up @@ -359,6 +360,21 @@ struct Kstarqa {

if (doprocessRecPhi) {
hInvMass.add("PhiMCChecks", "PhiMCChecks", kTHnSparseF, {ptAxis, etaAxis, phiAxis, subAxis, subAxis, subAxis}); // pT, eta, phi, deltaEta, deltaPhi, deltaR

hShift.add("h2D_dpx_vs_px_KaonFromPhi", "Kaon from Phi: dpx vs px", kTH2F, {{subAxis}, {ptAxis}});
hShift.add("h2D_dpy_vs_py_KaonFromPhi", "Kaon from Phi: dpy vs py", kTH2F, {{subAxis}, {ptAxis}});
hShift.add("h2D_dpz_vs_pz_KaonFromPhi", "Kaon from Phi: dpz vs pz", kTH2F, {{subAxis}, {ptAxis}});
hShift.add("h2D_dpt_vs_pt_KaonFromPhi", "Kaon from Phi: dpt vs pt", kTH2F, {{subAxis}, {ptAxis}});

hShift.add("h2D_dpx_vs_px_AllKaons", "Kaon from Phi: deta vs eta", kTH2F, {{subAxis}, {etaAxis}});
hShift.add("h2D_dpy_vs_py_AllKaons", "Kaon from Phi: dphi vs phi", kTH2F, {{subAxis}, {phiAxis}});
hShift.add("h2D_dpz_vs_pz_AllKaons", "Kaon from Phi: dR vs eta", kTH2F, {{subAxis}, {etaAxis}});
hShift.add("h2D_dpt_vs_pt_AllKaons", "Kaon from Phi: dR vs phi", kTH2F, {{subAxis}, {phiAxis}});

hShift.add("h2D_dpx_vs_px_Phi", "Pion from Phi: deta vs eta", kTH2F, {{subAxis}, {etaAxis}});
hShift.add("h2D_dpy_vs_py_Phi", "Pion from Phi: dphi vs phi", kTH2F, {{subAxis}, {phiAxis}});
hShift.add("h2D_dpz_vs_pz_Phi", "Pion from Phi: dR vs eta", kTH2F, {{subAxis}, {etaAxis}});
hShift.add("h2D_dpt_vs_pt_Phi", "Pion from Phi: dR vs phi", kTH2F, {{subAxis}, {phiAxis}});
}

// Signal Loss & Event Loss in Light Ion Collisions
Expand Down Expand Up @@ -3055,6 +3071,20 @@ struct Kstarqa {
continue;
}

// All Kaons (Filling single-track kaon resolution histograms)
const auto mctrack1_all = track1.mcParticle();
if (std::abs(mctrack1_all.pdgCode()) == PDG_t::kKPlus) {
double dpx_all = mctrack1_all.px() - track1.px();
double dpy_all = mctrack1_all.py() - track1.py();
double dpz_all = mctrack1_all.pz() - track1.pz();
double dpt_all = mctrack1_all.pt() - track1.pt();

hShift.fill(HIST("h2D_dpx_vs_px_AllKaons"), dpx_all, track1.px());
hShift.fill(HIST("h2D_dpy_vs_py_AllKaons"), dpy_all, track1.py());
hShift.fill(HIST("h2D_dpz_vs_pz_AllKaons"), dpz_all, track1.pz());
hShift.fill(HIST("h2D_dpt_vs_pt_AllKaons"), dpt_all, track1.pt());
}

auto track1ID = track1.index();
for (const auto& track2 : tracks) {
rEventSelection.fill(HIST("recMCparticles"), 0.5);
Expand Down Expand Up @@ -3169,6 +3199,38 @@ struct Kstarqa {
auto genPhiPhi = mothertrack1.phi();

hInvMass.fill(HIST("PhiMCChecks"), mother.pt(), mother.eta(), mother.phi(), genpTPhi - mother.pt(), genEtaPhi - mother.eta(), genPhiPhi - mother.phi());

// Kaons from Phi daughter (Filling for daughter track1 & daughter track2)
double dpx_k1 = mctrack1.px() - track1.px();
double dpy_k1 = mctrack1.py() - track1.py();
double dpz_k1 = mctrack1.pz() - track1.pz();
double dpt_k1 = mctrack1.pt() - track1.pt();

hShift.fill(HIST("h2D_dpx_vs_px_KaonFromPhi"), dpx_k1, track1.px());
hShift.fill(HIST("h2D_dpy_vs_py_KaonFromPhi"), dpy_k1, track1.py());
hShift.fill(HIST("h2D_dpz_vs_pz_KaonFromPhi"), dpz_k1, track1.pz());
hShift.fill(HIST("h2D_dpt_vs_pt_KaonFromPhi"), dpt_k1, track1.pt());

double dpx_k2 = mctrack2.px() - track2.px();
double dpy_k2 = mctrack2.py() - track2.py();
double dpz_k2 = mctrack2.pz() - track2.pz();
double dpt_k2 = mctrack2.pt() - track2.pt();

hShift.fill(HIST("h2D_dpx_vs_px_KaonFromPhi"), dpx_k2, track2.px());
hShift.fill(HIST("h2D_dpy_vs_py_KaonFromPhi"), dpy_k2, track2.py());
hShift.fill(HIST("h2D_dpz_vs_pz_KaonFromPhi"), dpz_k2, track2.pz());
hShift.fill(HIST("h2D_dpt_vs_pt_KaonFromPhi"), dpt_k2, track2.pt());

// Phi Meson
double dpx_phi = mothertrack1.px() - mother.Px();
double dpy_phi = mothertrack1.py() - mother.Py();
double dpz_phi = mothertrack1.pz() - mother.Pz();
double dpt_phi = mothertrack1.pt() - mother.Pt();

hShift.fill(HIST("h2D_dpx_vs_px_Phi"), dpx_phi, mother.Px());
hShift.fill(HIST("h2D_dpy_vs_py_Phi"), dpy_phi, mother.Py());
hShift.fill(HIST("h2D_dpz_vs_pz_Phi"), dpz_phi, mother.Pz());
hShift.fill(HIST("h2D_dpt_vs_pt_Phi"), dpt_phi, mother.Pt());
}
}
}
Expand Down
Loading