-
Notifications
You must be signed in to change notification settings - Fork 622
[PWGHF] Add preselection and centrality to st tree creator #14252
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.
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
base: master
Are you sure you want to change the base?
[PWGHF] Add preselection and centrality to st tree creator #14252
Conversation
Fix a name bug.
[PWGHF] Please consider the following formatting changes to AliceO2Group#13353
Added new mass columns for Omega and Xi particles and updated event selection criteria.
…ty-to-STTree-creatore
…ty-to-STTree-creatore
Please consider the following formatting changes to AliceO2Group#14252
…ty-to-STTree-creatore
| {"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.}}}}, |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
hMassXiPi is redundant
| default: | ||
| passSelectedChannel = true; // unexpected code -> do not reject | ||
| break; |
There was a problem hiding this comment.
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?
There was a problem hiding this comment.
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?
There was a problem hiding this comment.
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.)
There was a problem hiding this comment.
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?
There was a problem hiding this comment.
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.
…ty-to-STTree-creatore
Updated variable names for clarity and consistency in the treeCreatorOmegacSt.cxx file. Removed unused histogram for mass Xi + pi.
| const auto ptCascade = std::hypot(momenta[0][0], momenta[0][1]); | ||
| const auto ptPionOrKaon = std::hypot(momenta[1][0], momenta[1][1]); |
There was a problem hiding this comment.
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}; |
There was a problem hiding this comment.
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.
| 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; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
These are not masses.
| case XicToXiPi: | ||
| passSelectedChannel = massXicToXiPiPass; | ||
| break; | ||
| case XicToOmegaKa: | ||
| passSelectedChannel = massXicToOmegaKaPass; | ||
| break; | ||
| case OmegacToOmegaPi: | ||
| passSelectedChannel = massOmegacToOmegaPiPass; | ||
| break; | ||
| case OmegacToOmegaKa: |
There was a problem hiding this comment.
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)?
| if (cpaCharmedBaryon < cpaCharmedBaryonMin) | ||
| continue; | ||
| if (cpaXYCharmedBaryon < cpaXYCharmedBaryonMin) | ||
| continue; | ||
| if (cpaCasc < cpaCascMin) | ||
| continue; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Missing braces everywhere.
No description provided.