Skip to content

[PWGHF] Implement TPC pid for light nuclei based on Bethe-Bloch parametrization in HF track skimming task #14729

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

Draft
zhangbiao-phy wants to merge 9 commits into
base: master
Choose a base branch
from
+187 −146
Draft

[PWGHF] Implement TPC pid for light nuclei based on Bethe-Bloch parametrization in HF track skimming task #14729

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
027e597
enabled TPC pid for light nuclei based on Bethe-Bloch parametrization
zhangbiao-phy Jan 31, 2026
590e563
Add PID QA plots
zhangbiao-phy Jan 31, 2026
5200e72
fix typo
zhangbiao-phy Jan 31, 2026
129256a
fix warning
zhangbiao-phy Jan 31, 2026
01c6f00
fix warning
zhangbiao-phy Jan 31, 2026
309a35b
fix error
zhangbiao-phy Feb 1, 2026
457ae60
fix comments from vit
zhangbiao-phy Feb 1, 2026
79aeaaa
fix comments from vit and move the track quality and pid selection to…
zhangbiao-phy Feb 8, 2026
e117d65
fix the changes related to default setting
zhangbiao-phy Feb 8, 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
25 changes: 16 additions & 9 deletions PWGHF/Core/SelectorCuts.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -67,30 +67,37 @@ static const std::vector<std::string> labelsCutVarTrack = {"min_dcaxytoprimary",

namespace hf_presel_pid
{
static constexpr int NPidRows = 7; // number of PID channels / rows
static constexpr int NPidRows = 4; // number of PID channels / rows
static constexpr int NPidCuts = 6; // number of cuts per PID (TPC+TOF)
// default values for the PID cuts for protons in the track-index-skim-creator
constexpr float CutsPid[NPidRows][NPidCuts] = {{0.f, 1000.f, 5.f, 0.f, 1000.f, 5.f},
{0.f, 1000.f, 5.f, 0.f, 1000.f, 5.f},
{0.f, 1000.f, 5.f, 0.f, 1000.f, 5.f},
{0.f, 1000.f, 5.f, 0.f, 1000.f, 5.f},
{0.f, 1000.f, 5.f, 0.f, 1000.f, 5.f},
{0.f, 1000.f, 5.f, 0.f, 1000.f, 5.f},
{0.f, 1000.f, 5.f, 0.f, 1000.f, 5.f}};
static const std::vector<std::string> labelsCutsPid = {"minPtTpc", "maxPtTpc", "nSigmaMaxTpc", "minPtTof", "maxPtTof", "nSigmaMaxTof"};
static const std::vector<std::string> labelsRowsPid = {"ProtonInLcToPKPi", "ProtonInXicToPKPi", "ProtonInLcToPK0S", "KaonIn3Prongs", "DeuteronInCdToDeKPi", "TritonInCtToTrKPi", "HeliumInChToHeKPi"};
static const std::vector<std::string> labelsRowsPid = {"ProtonInLcToPKPi", "ProtonInXicToPKPi", "ProtonInLcToPK0S", "KaonIn3Prongs"};
} // namespace hf_presel_pid

namespace hf_presel_lightnuclei
{

static constexpr int NParticleRows = 3; // number of particles / rows
static constexpr int NVarCuts = 10; // number of cuts for each particles
static constexpr int NBetheBlochParams = 6; // number of parameters for Bethe-Bloch

// default values for the track cuts for lightnuclei in the track-index-skim-creator
constexpr float CutsTrackQuality[3][9] = {{-4, 3, 5., 0., 100, 100, 0.83, 160., 1.},
{-4, 3, 5., 0., 100, 100, 0.83, 160., 1.},
{-4, 3, 5., 0., 100, 100, 0.83, 160., 1.}};
static const std::vector<std::string> labelsCutsTrack = {"nSigmaMinIts", "minItsClusterSizes", "minItsCluster", "minItsIbCluster", "minTpcCluster", "minTpcRow", "minTpcCrossedOverFound", "maxTpcShared", "maxTpcFracShared"};
constexpr float CutsTrackQuality[NParticleRows][NVarCuts] = {{-4.f, 3.f, 5.f, 0.f, 100.f, 100.f, 0.83, 160.f, 1.f, 5.f},
{-4.f, 3.f, 5.f, 0.f, 100.f, 100.f, 0.83, 160.f, 1.f, 5.f},
{-4.f, 3.f, 5.f, 0.f, 100.f, 100.f, 0.83, 160.f, 1.f, 5.f}};
static const std::vector<std::string> labelsCutsTrack = {"nSigmaMinIts", "minItsClusterSizes", "minItsCluster", "minItsIbCluster", "minTpcCluster", "minTpcRow", "minTpcCrossedOverFound", "maxTpcShared", "maxTpcFracShared", "maxTPCnSigmaBB"};
static const std::vector<std::string> labelsRowsNucleiType = {"Deutron", "Triton", "Helium3"};

constexpr float BetheBlochParams[NParticleRows][NBetheBlochParams] = {{5.39302, 7.859534, 0.004048, 2.323197, 1.609307, 0.09},
{5.39302, 7.859534, 0.004048, 2.323197, 1.609307, 0.09},
{-126.55736, -0.858569, 1.11164, 1.21032, 2.656374, 0.09}};

static const std::vector<std::string> labelsBetheBlochParams = {"p0", "p1", "p2", "p3", "p4", "resolution"};

} // namespace hf_presel_lightnuclei

namespace hf_cuts_bdt_multiclass
Expand Down
Loading
Loading