Skip to content

Conversation

@FDUEnrich
Copy link
Contributor

No description provided.

FDUEnrich and others added 24 commits October 12, 2025 15:07
[PWGHF] Please consider the following formatting changes to AliceO2Group#13353
Added new mass columns for Omega and Xi particles and updated event selection criteria.
Please consider the following formatting changes to AliceO2Group#14252
@FDUEnrich FDUEnrich marked this pull request as ready for review December 24, 2025 10:23
@vkucera vkucera marked this pull request as draft December 24, 2025 16:19
@FDUEnrich FDUEnrich marked this pull request as ready for review December 26, 2025 06:32
@vkucera vkucera marked this pull request as draft December 26, 2025 13:11
@FDUEnrich FDUEnrich marked this pull request as ready for review December 29, 2025 13:29
@alibuild
Copy link
Collaborator

Error while checking build/O2Physics/o2 for 9d4f31c at 2025-12-30 15:40:

## sw/BUILD/O2Physics-latest/log
c++: fatal error: Killed signal terminated program cc1plus
ninja: build stopped: subcommand failed.

Full log here.

Comment on lines 355 to 356
{"hMassXiPi", "inv. mass #Xi + #pi;inv. mass (GeV/#it{c}^{2})", {HistType::kTH1D, {{400, 1.5, 3.}}}},
{"hMassXiPiVsPt", "inv. mass #Xi + #pi;inv. mass (GeV/#it{c}^{2});p_{T} (GeV/#it{c})", {HistType::kTH2D, {{400, 1.5, 3.}, {10, 0., 10.}}}},
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

hMassXiPi is redundant

Comment on lines 808 to 810
default:
passSelectedChannel = true; // unexpected code -> do not reject
break;
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Why do you silently accept unexpected cases?

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

What is the point of doing this?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I want to either distinguish different decay channel reconstructions or simultaneously reconstruct two or more decay channels. (The logical error has been fixed.)

Copy link
Collaborator

@vkucera vkucera Jan 12, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

According to your documentation, the case with all channels is matched to the value AnyChannel, so how can the code do anything useful if it receives an unexpected value?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

If I want AnyChannels or other channel, and neither is satisfied, it will enter the default option and continue skipping the loop.

@vkucera vkucera marked this pull request as draft January 8, 2026 00:36
Updated variable names for clarity and consistency in the treeCreatorOmegacSt.cxx file. Removed unused histogram for mass Xi + pi.
@FDUEnrich FDUEnrich marked this pull request as ready for review January 8, 2026 13:43
@vkucera vkucera marked this pull request as draft January 12, 2026 06:02
Comment on lines +769 to +770
const auto ptCascade = std::hypot(momenta[0][0], momenta[0][1]);
const auto ptPionOrKaon = std::hypot(momenta[1][0], momenta[1][1]);
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

You could just use RecoDecay::pt(momenta[0]).

std::array<double, NDaughters> const massesOmegacToOmegaPi{o2::constants::physics::MassOmegaMinus, o2::constants::physics::MassPiPlus};
std::array<double, NDaughters> const massesOmegacToOmegaK{o2::constants::physics::MassOmegaMinus, o2::constants::physics::MassKPlus};
std::array<double, NDaughters> const massesXicDaughters{o2::constants::physics::MassXiMinus, o2::constants::physics::MassPiPlus};
std::array<double, NDaughters> const massOmegaK{o2::constants::physics::MassOmegaMinus, o2::constants::physics::MassKPlus};
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Wrong name. This is not a single mass.

Comment on lines +780 to +784
const bool massOmegacToOmegaPiPass = std::abs(massOmegaPion - o2::constants::physics::MassOmegaC0) < massWindowOmegaC;
const bool massOmegacToOmegaKPass = std::abs(massOmegaKaon - o2::constants::physics::MassOmegaC0) < massWindowOmegaC;
const bool massOmegacToXiPiPass = std::abs(massXiPion - o2::constants::physics::MassOmegaC0) < massWindowOmegaC;
const bool massXicToOmegaKaPass = std::abs(massOmegaKaon - o2::constants::physics::MassXiC0) < massWindowXiC;
const bool massXicToXiPiPass = std::abs(massXiPion - o2::constants::physics::MassXiC0) < massWindowXiC;
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

These are not masses.

Comment on lines +791 to +800
case XicToXiPi:
passSelectedChannel = massXicToXiPiPass;
break;
case XicToOmegaKa:
passSelectedChannel = massXicToOmegaKaPass;
break;
case OmegacToOmegaPi:
passSelectedChannel = massOmegacToOmegaPiPass;
break;
case OmegacToOmegaKa:
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Why don't you calculate the conditions here (only when needed)?

Comment on lines +869 to +874
if (cpaCharmedBaryon < cpaCharmedBaryonMin)
continue;
if (cpaXYCharmedBaryon < cpaXYCharmedBaryonMin)
continue;
if (cpaCasc < cpaCascMin)
continue;
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Missing braces everywhere.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

pwghf PWG-HF

Development

Successfully merging this pull request may close these issues.

3 participants