Skip to content

[PWGLF] Fix xi1530Analysisqa.cxx #12673

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
dmallick2 merged 7 commits into from
Aug 21, 2025
Merged

[PWGLF] Fix xi1530Analysisqa.cxx #12673

Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
6e64f20
added initializerevetqa.cxx and updated xi1530Analysisqa.cxx
mjkim525 Aug 20, 2025
a4821a9
Please consider the following formatting changes
alibuild Aug 20, 2025
3bb699b
Merge pull request #6 from alibuild/alibot-cleanup-12657
mjkim525 Aug 20, 2025
78a29da
Merge branch 'AliceO2Group:master' into master
mjkim525 Aug 21, 2025
aeba5a8
Fix xi1530Analysisqa.cxx
mjkim525 Aug 21, 2025
f90f3e2
Please consider the following formatting changes
alibuild Aug 21, 2025
6ae7ebd
Merge pull request #7 from alibuild/alibot-cleanup-12673
mjkim525 Aug 21, 2025
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
119 changes: 82 additions & 37 deletions PWGLF/Tasks/Resonances/xi1530Analysisqa.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1058,13 +1058,15 @@ struct Xi1530Analysisqa {
aod::ResoTracks const& resoTracks,
aod::ResoCascades const& cascTracks)
{
auto linkRow = collisionIndex.iteratorAt(resoCollision.globalIndex());
auto collId = linkRow.collisionId(); // Take original collision global index matched with resoCollision
if (cRecoINELgt0) {
auto linkRow = collisionIndex.iteratorAt(resoCollision.globalIndex());
auto collId = linkRow.collisionId(); // Take original collision global index matched with resoCollision

auto coll = collisions.iteratorAt(collId); // Take original collision matched with resoCollision
auto coll = collisions.iteratorAt(collId); // Take original collision matched with resoCollision

if (cRecoINELgt0 && !coll.isInelGt0()) // Check reco INELgt0 (at least one PV track in |eta| < 1) about the collision
return;
if (!coll.isInelGt0()) // Check reco INELgt0 (at least one PV track in |eta| < 1) about the collision
return;
}

histos.fill(HIST("QAevent/hEvtCounterSameE"), 1.0);
auto multiplicity = resoCollision.cent();
Expand All @@ -1073,52 +1075,92 @@ struct Xi1530Analysisqa {

// Reconstructed level MC for the track
void processMC(ResoMCCols::iterator const& resoCollision,
aod::ResoCollisionColls const& resoCollisionIndex,
aod::ResoCollisionColls const& collisionIndex,
soa::Join<aod::ResoCollisionCandidatesMC, aod::PVMults> const& collisionsMC,
soa::Join<aod::ResoCascades, aod::ResoMCCascades> const& cascTracks,
soa::Join<aod::ResoTracks, aod::ResoMCTracks> const& resoTracks,
soa::Join<aod::McCollisions, aod::McCentFT0Ms> const&)
{
if (!resoCollision.isInAfterAllCuts() || (std::abs(resoCollision.posZ()) > cZvertCutMC)) // MC event selection, all cuts missing vtx cut
return;
float multiplicity;
if (cMCCent && cRecoINELgt0) {
auto linkRow = collisionIndex.iteratorAt(resoCollision.globalIndex());
auto collId = linkRow.collisionId(); // Take original collision global index matched with resoCollision

auto linkRow = resoCollisionIndex.iteratorAt(resoCollision.globalIndex());
const int collId = linkRow.collisionId();
auto coll = collisionsMC.iteratorAt(collId); // Take original collision matched with resoCollision

auto coll = collisionsMC.iteratorAt(collId);
if (!coll.isInelGt0()) // Check reco INELgt0 (at least one PV track in |eta| < 1) about the collision
return;

if (cRecoINELgt0 && !coll.isInelGt0())
return;
auto mcColl = coll.mcCollision_as<soa::Join<aod::McCollisions, aod::McCentFT0Ms>>();
multiplicity = mcColl.centFT0M();
} else if (!cMCCent && cRecoINELgt0) {
auto linkRow = collisionIndex.iteratorAt(resoCollision.globalIndex());
auto collId = linkRow.collisionId(); // Take original collision global index matched with resoCollision

auto coll = collisionsMC.iteratorAt(collId); // Take original collision matched with resoCollision

auto mcColl = coll.mcCollision_as<soa::Join<aod::McCollisions, aod::McCentFT0Ms>>();
if (!coll.isInelGt0()) // Check reco INELgt0 (at least one PV track in |eta| < 1) about the collision
return;

auto multiplicityReco = resoCollision.cent(); // Reco level multiplicity per.
auto multiplicityGen = mcColl.centFT0M(); // Gen level multiplicity per.
multiplicity = resoCollision.cent();
} else if (cMCCent && !cRecoINELgt0) {
auto linkRow = collisionIndex.iteratorAt(resoCollision.globalIndex());
auto collId = linkRow.collisionId(); // Take original collision global index matched with resoCollision

float multiplicity = cMCCent ? multiplicityGen : multiplicityReco;
auto coll = collisionsMC.iteratorAt(collId); // Take original collision matched with resoCollision

auto mcColl = coll.mcCollision_as<soa::Join<aod::McCollisions, aod::McCentFT0Ms>>();
multiplicity = mcColl.centFT0M();
} else {
multiplicity = resoCollision.cent();
}

if (!resoCollision.isInAfterAllCuts() || (std::abs(resoCollision.posZ()) > cZvertCutMC)) // MC event selection, all cuts missing vtx cut
return;

fillHistograms<false, true, false>(resoCollision, multiplicity, resoTracks, cascTracks);
}

// Truth level MC for the track with reco event
void processMCTrue(ResoMCCols::iterator const& resoCollision,
aod::ResoCollisionColls const& resoCollisionIndex,
aod::ResoCollisionColls const& collisionIndex,
aod::ResoMCParents const& resoParents,
aod::ResoCollisionCandidatesMC const& collisionsMC,
soa::Join<aod::McCollisions, aod::McCentFT0Ms> const&)
{
float multiplicity;
if (cMCCent && cRecoINELgt0) {
auto linkRow = collisionIndex.iteratorAt(resoCollision.globalIndex());
auto collId = linkRow.collisionId(); // Take original collision global index matched with resoCollision

auto linkRow = resoCollisionIndex.iteratorAt(resoCollision.globalIndex());
const int collId = linkRow.collisionId();
auto coll = collisionsMC.iteratorAt(collId); // Take original collision matched with resoCollision

auto coll = collisionsMC.iteratorAt(collId);
if (!coll.isInelGt0()) // Check reco INELgt0 (at least one PV track in |eta| < 1) about the collision
return;

auto mcColl = coll.mcCollision_as<soa::Join<aod::McCollisions, aod::McCentFT0Ms>>();
auto mcColl = coll.mcCollision_as<soa::Join<aod::McCollisions, aod::McCentFT0Ms>>();
multiplicity = mcColl.centFT0M();
} else if (!cMCCent && cRecoINELgt0) {
auto linkRow = collisionIndex.iteratorAt(resoCollision.globalIndex());
auto collId = linkRow.collisionId(); // Take original collision global index matched with resoCollision

auto multiplicityReco = resoCollision.cent(); // Reco level multiplicity per.
auto multiplicityGen = mcColl.centFT0M(); // Gen level multiplicity per.
auto coll = collisionsMC.iteratorAt(collId); // Take original collision matched with resoCollision

float multiplicity = cMCCent ? multiplicityGen : multiplicityReco;
if (!coll.isInelGt0()) // Check reco INELgt0 (at least one PV track in |eta| < 1) about the collision
return;

multiplicity = resoCollision.cent();
} else if (cMCCent && !cRecoINELgt0) {
auto linkRow = collisionIndex.iteratorAt(resoCollision.globalIndex());
auto collId = linkRow.collisionId(); // Take original collision global index matched with resoCollision

auto coll = collisionsMC.iteratorAt(collId); // Take original collision matched with resoCollision

auto mcColl = coll.mcCollision_as<soa::Join<aod::McCollisions, aod::McCentFT0Ms>>();
multiplicity = mcColl.centFT0M();
} else {
multiplicity = resoCollision.cent();
}

for (const auto& part : resoParents) { // loop over all pre-filtered MC particles
if (std::abs(part.pdgCode()) != kXiStar || std::abs(part.y()) >= cfgRapidityCut)
Expand Down Expand Up @@ -1171,13 +1213,15 @@ struct Xi1530Analysisqa {
aod::ResoMicroTracks const& resomicrotracks,
aod::ResoCascades const& cascTracks)
{
auto linkRow = collisionIndex.iteratorAt(resoCollision.globalIndex());
auto collId = linkRow.collisionId(); // Take original collision global index matched with resoCollision
if (cRecoINELgt0) {
auto linkRow = collisionIndex.iteratorAt(resoCollision.globalIndex());
auto collId = linkRow.collisionId(); // Take original collision global index matched with resoCollision

auto coll = collisions.iteratorAt(collId); // Take original collision matched with resoCollision
auto coll = collisions.iteratorAt(collId); // Take original collision matched with resoCollision

if (cRecoINELgt0 && !coll.isInelGt0()) // Check reco INELgt0 (at least one PV track in |eta| < 1) about the collision
return;
if (!coll.isInelGt0()) // Check reco INELgt0 (at least one PV track in |eta| < 1) about the collision
return;
}

histos.fill(HIST("QAevent/hEvtCounterSameE"), 1.0);
auto multiplicity = resoCollision.cent();
Expand Down Expand Up @@ -1222,15 +1266,16 @@ struct Xi1530Analysisqa {

for (const auto& [collision1, tracks1, collision2, tracks2] : pairs) {

const auto rcIdx = collision1.globalIndex();

const auto linkRow = collisionIndex.iteratorAt(rcIdx);
const auto collId = linkRow.collisionId();
if (cRecoINELgt0) {
const auto rcIdx = collision1.globalIndex();
auto linkRow = collisionIndex.iteratorAt(rcIdx);
auto collId = linkRow.collisionId(); // Take original collision global index matched with resoCollision

auto coll = collisions.iteratorAt(collId);
if (cRecoINELgt0 && !coll.isInelGt0()) // Check reco INELgt0 (at least one PV track in |eta| < 1) about the collision
continue;
auto coll = collisions.iteratorAt(collId); // Take original collision matched with resoCollision

if (!coll.isInelGt0()) // Check reco INELgt0 (at least one PV track in |eta| < 1) about the collision
continue;
}
histos.fill(HIST("QAevent/hEvtCounterMixedE"), 1.0);
auto multiplicity = collision1.cent();
fillHistograms<true, false, true>(collision1, multiplicity, tracks1, tracks2);
Expand Down
Loading