-
Notifications
You must be signed in to change notification settings - Fork 621
[Common,PWGHF] Complete the charmed hypernuclei search and add light-nuclei track-quality cut to the skimming task #14390
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?
Conversation
… of light nuclei in the skimming task
|
O2 linter results: ❌ 25 errors, |
fgrosa
left a comment
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.
Hi @zhangbiao-phy thanks a lot for the development! I still didn't check the full PR in detail, but I have a couple of important comments to be kept in mind
Common/Core/TrackSelectorPID.h
Outdated
| // ITS | ||
| float mPtItsMin = 0.; ///< minimum pT for TPC PID [GeV/c] | ||
| float mPtItsMax = 100.; ///< maximum pT for TPC PID [GeV/c] | ||
| float mNSigmaItsMin = -3.; ///< minimum number of TPC σ | ||
| float mNSigmaItsMax = 3.; ///< maximum number of TPC σ | ||
| float mNSigmaItsMinCondTpc = 0.; ///< minimum number of TPC σ if combined with TOF | ||
| float mNSigmaItsMaxCondTpc = 0.; ///< maximum number of TPC σ if combined with TOF |
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 is this here?
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.
Thanks for spotting this! I had it in my first version, where I added all the functions related to ITS also to the PID selector class. However, I realized that they are not commonly used in HF, so I decided not to include them. I then forgot to remove it.
|
|
||
| int8_t flag = BIT(ChannelsProtonPid::NChannelsProtonPid + 2) - 1; // all bits on (including the kaon one) | ||
| for (auto iChannel{0u}; iChannel < ChannelsProtonPid::NChannelsProtonPid + 2; ++iChannel) { | ||
| int8_t flag = BIT(ChannelsProtonPid::NChannelsProtonPid + 4) - 1; // all bits on (including the kaon one) |
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 is this hard-coded?
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 added two decay channels. but I fixed it now
|
|
||
| using SelectedCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::HfSelCollision>>; | ||
| using TracksWithPVRefitAndDCA = soa::Join<aod::TracksWCovDcaExtra, aod::HfPvRefitTrack>; | ||
| using TracksWithPVRefitAndDCA = soa::Join<aod::FullTracks, aod::TracksWCovDcaExtra, aod::HfPvRefitTrack>; |
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?
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.
remove it
| #include "Common/DataModel/Centrality.h" | ||
| #include "Common/DataModel/CollisionAssociationTables.h" | ||
| #include "Common/DataModel/EventSelection.h" | ||
| #include "Common/DataModel/PIDResponseITS.h" |
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?
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.
which need it for ITS nsigma PID in:
auto tracksWithItsPid = soa::Attach<TTracks, aod::pidits::ITSNSigmaDe, aod::pidits::ITSNSigmaTr, aod::pidits::ITSNSigmaHe>(tracks);
| const std::vector<double> ptBinsMl{0., 1.e10}; | ||
| const std::vector<int> cutDirMl{o2::cuts_ml::CutDirection::CutGreater, o2::cuts_ml::CutDirection::CutSmaller, o2::cuts_ml::CutDirection::CutSmaller}; | ||
| const std::array<LabeledArray<double>, kN3ProngDecays - 1> thresholdMlScore3Prongs{config.thresholdMlScoreDplusToPiKPi, config.thresholdMlScoreLcToPiKP, config.thresholdMlScoreDsToPiKK, config.thresholdMlScoreXicToPiKP}; | ||
| const std::array<LabeledArray<double>, kN3ProngDecays - 3> thresholdMlScore3Prongs{config.thresholdMlScoreDplusToPiKPi, config.thresholdMlScoreLcToPiKP, config.thresholdMlScoreDsToPiKK, config.thresholdMlScoreXicToPiKP}; |
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 is this hard-coded?
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.
fixed
|
|
||
| // initialise 3-prong ML responses | ||
| for (int iDecay3P{0}; iDecay3P < kN3ProngDecays - 1; ++iDecay3P) { | ||
| for (int iDecay3P{0}; iDecay3P < kN3ProngDecays - 3; ++iDecay3P) { |
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 is this hard-coded?
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.
fixed
[PWGHF] Please consider the following formatting changes to AliceO2Group#14390
Tag @fgrosa for the informations