Skip to content

[PWGHF] [xic0omegac0] Update codes to fix PID assignment in KFParticle #16381

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
fcatalan92 merged 13 commits into from
May 27, 2026
+53 −26
Merged

[PWGHF] [xic0omegac0] Update codes to fix PID assignment in KFParticle #16381

Changes from all commits
Commits
Show all changes
13 commits
Select commit Hold shift + click to select a range
a9467ce
Assigning o2::track::PID value for getTrackParCovFromKFP function
May 7, 2026
8cd6195
Merge branch 'AliceO2Group:master' into kfupdate
jinhyunni May 26, 2026
9226eb3
Fix cpplint issue
May 26, 2026
706420b
Please consider the following formatting changes
alibuild May 26, 2026
2c4b562
Merge pull request #22 from alibuild/alibot-cleanup-16381
jinhyunni May 26, 2026
65566ca
Remove unnecessary lines
May 26, 2026
09a6a5c
Remove unnecessary lines
May 26, 2026
c621847
Remove unnecessary printout for debugging
May 26, 2026
600df9b
Merge branch 'master' of https://github.com/AliceO2Group/O2Physics in…
May 27, 2026
c8b1aa9
Added helper function to give right input for getTrackParCovFromKFP()
May 27, 2026
d20de6e
Please consider the following formatting changes
alibuild May 27, 2026
4fc0b25
Merge pull request #23 from alibuild/alibot-cleanup-16381
jinhyunni May 27, 2026
ab89135
Update return value for default case in switch statement
May 27, 2026
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
79 changes: 53 additions & 26 deletions PWGHF/TableProducer/candidateCreatorXic0Omegac0Qa.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -261,7 +261,7 @@ struct HfCandidateCreatorXic0Omegac0Qa {
int pdgIdOfV0DauPos{}, pdgIdOfV0DauNeg{}, pdgIdOfBach{}, pdgIdOfCharmBach{};
int pdgIdOfV0{}, pdgIdOfCascade{}, pdgIdOfCharmBaryon{};

// Track PID - Used in DCAFitter
// Track PID: PID value of tracks defined under o2::track::PID namespace
int trackPidOfCascade{};

// Mass of daughter tracks & V0s & cascades & charm baryons;
Expand Down Expand Up @@ -362,20 +362,21 @@ struct HfCandidateCreatorXic0Omegac0Qa {
massOfV0 = o2::constants::physics::MassLambda;
massOfCascade = o2::constants::physics::MassOmegaMinus;
}
LOGF(info, "~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~");
LOGF(info, "PDG ID of V0 positive daughter: %d", pdgIdOfV0DauPos);
LOGF(info, "PDG ID of V0 negative daughter: %d", pdgIdOfV0DauNeg);
LOGF(info, "PDG ID of Bachelor: %d", pdgIdOfBach);
LOGF(info, "PDG ID of Charm Bachelor: %d", pdgIdOfCharmBach);
LOGF(info, "-------------------------------------------");
LOGF(info, "PDG ID of V0: %d", pdgIdOfV0);
LOGF(info, "PDG ID of Cascade: %d", pdgIdOfCascade);
LOGF(info, "PDG ID of Charm Baryon: %d", pdgIdOfCharmBaryon);
LOGF(info, "-------------------------------------------");
LOGF(info, "Mass of V0 set as: %f", massOfV0);
LOGF(info, "Mass of CharmBach set as: %f", massOfCharmBach);
LOGF(info, "Mass of Casc as: %f", massOfCascade);
LOGF(info, "~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~");
LOGF(info, "~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~");
LOGF(info, " Chosen reconstruction info: ToXiPi(%d)/ ToOmegaPi(%d)/ ToOmegaKa(%d)", xipiEnabledDca + xipiEnabledKf, omegapiEnabledDca + omegapiEnabledKf, omegakaEnabledDca + omegakaEnabledKf);
LOGF(info, "~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~");
LOGF(info, " PDG ID of V0 positive daughter: %d", pdgIdOfV0DauPos);
LOGF(info, " PDG ID of V0 negative daughter: %d", pdgIdOfV0DauNeg);
LOGF(info, " PDG ID of V0: %d", pdgIdOfV0);
LOGF(info, " PDG ID of Bachelor: %d", pdgIdOfBach);
LOGF(info, " PDG ID of Cascade: %d", pdgIdOfCascade);
LOGF(info, " PDG ID of Charm Bachelor: %d", pdgIdOfCharmBach);
LOGF(info, " PDG ID of Charm Baryon: %d", pdgIdOfCharmBaryon);
LOGF(info, "-----------------------------------------------------------------------");
LOGF(info, " Mass of V0 set as: %f", massOfV0);
LOGF(info, " Mass of CharmBach set as: %f", massOfCharmBach);
LOGF(info, " Mass of Casc as: %f", massOfCascade);
LOGF(info, "~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~");

// Add histogram to indicate which sv method was used
registry.add("hVertexerType", "Use KF or DCAFitterN;Vertexer type;entries", {kTH1F, {{2, 0.0, 2.0}}});
Expand Down Expand Up @@ -491,6 +492,30 @@ struct HfCandidateCreatorXic0Omegac0Qa {

} // end of initialization

// helper function to convert PDG of KFParticle object into appropriate track PID.
// exclusively used to make input of getTrackParCovFromKFP()
// \brief kfPdg is Pdg value of KFParticle object
o2::track::PID::ID convertPDGIntoPID(const int kfPdg)
{
switch (std::abs(kfPdg)) {
case kPiPlus:
return o2::track::PID::Pion;
case kKPlus:
return o2::track::PID::Kaon;
case kProton:
return o2::track::PID::Proton;
case kLambda0:
return o2::track::PID::Lambda;
case kXiMinus:
return o2::track::PID::XiMinus;
case kOmegaMinus:
return o2::track::PID::OmegaMinus;
default:
LOGF(fatal, "Undefined PDG value from KFParticle given for conversion(%d). Replace it with PID of Pion", kfPdg);
return o2::track::PID::Pion;
}
}

// template function for running charm baryon reconstruction with DCAFitter
/// \brief centEstimator is for different centrality estimators
/// \brief decayChannel is for different decay channels. 0 for XiczeroOmegaczeroToXiPi, 1 for OmegaczeroToOmegaPi, 2 for OmegaczeroToOmeagaK
Expand Down Expand Up @@ -934,6 +959,7 @@ struct HfCandidateCreatorXic0Omegac0Qa {
const KFPTrack kfTrackBach = createKFPTrackFromTrack(bachTrack);

bool isAnti = (bachTrack.signed1Pt() > 0 ? true : false);
int bachCharge = (bachTrack.signed1Pt() > 0 ? +1 : -1);

KFParticle kfPos(kfTrack0, (isAnti ? -pdgIdOfV0DauNeg : pdgIdOfV0DauPos));
KFParticle kfNeg(kfTrack1, (isAnti ? -pdgIdOfV0DauPos : pdgIdOfV0DauNeg));
Expand Down Expand Up @@ -1082,6 +1108,18 @@ struct HfCandidateCreatorXic0Omegac0Qa {
kfCharmBaryonToPv.SetProductionVertex(kfPv);

//----------Reconstruct information after vertex fit----------
auto trackParCovV0DauPos = getTrackParCovFromKFP(kfPos, convertPDGIntoPID(kfPos.GetPDG()), 1);
auto trackParCovV0DauNeg = getTrackParCovFromKFP(kfNeg, convertPDGIntoPID(kfNeg.GetPDG()), -1);
auto trackParCovBach = getTrackParCovFromKFP(kfBachToCasc, convertPDGIntoPID(kfBachToCasc.GetPDG()), bachCharge);
auto trackParCovCharmBach = getTrackParCovFromKFP(kfCharmBachToCharmBaryon, convertPDGIntoPID(kfCharmBachToCharmBaryon.GetPDG()), -bachCharge);
auto trackParCovCasc = getTrackParCovFromKFP(kfCascToCharmBaryon, convertPDGIntoPID(kfCascToCharmBaryon.GetPDG()), bachCharge);

trackParCovV0DauPos.setAbsCharge(1);
trackParCovV0DauNeg.setAbsCharge(1);
trackParCovBach.setAbsCharge(1);
trackParCovCharmBach.setAbsCharge(1);
trackParCovCasc.setAbsCharge(1);

std::array<float, 3> vertexV0 = {kfV0.GetX(), kfV0.GetY(), kfV0.GetZ()};
std::array<float, 3> vertexCasc = {kfCasc.GetX(), kfCasc.GetY(), kfCasc.GetZ()};

Expand All @@ -1098,17 +1136,6 @@ struct HfCandidateCreatorXic0Omegac0Qa {

std::array<float, 3> pvCoord = {collision.posX(), collision.posY(), collision.posZ()};

auto trackParCovV0DauPos = getTrackParCovFromKFP(kfPos, kfPos.GetPDG(), 1);
auto trackParCovV0DauNeg = getTrackParCovFromKFP(kfNeg, kfNeg.GetPDG(), -1);
auto trackParCovBach = getTrackParCovFromKFP(kfBachToCasc, kfBachToCasc.GetPDG(), (isAnti ? 1 : -1));
auto trackParCovCharmBach = getTrackParCovFromKFP(kfCharmBachToCharmBaryon, kfCharmBachToCharmBaryon.GetPDG(), (isAnti ? -1 : 1));
auto trackParCovCasc = getTrackParCovFromKFP(kfCascToCharmBaryon, kfCascToCharmBaryon.GetPDG(), (isAnti ? 1 : -1));
trackParCovV0DauPos.setAbsCharge(1);
trackParCovV0DauNeg.setAbsCharge(1);
trackParCovBach.setAbsCharge(1);
trackParCovCharmBach.setAbsCharge(1);
trackParCovCasc.setAbsCharge(1);

//----------Calculate physical quantities and fill candidate table----------

// impact parameters
Expand Down
Loading