diff --git a/PWGHF/DataModel/AliasTables.h b/PWGHF/DataModel/AliasTables.h
index 98594b546db..bb7cce309a3 100644
--- a/PWGHF/DataModel/AliasTables.h
+++ b/PWGHF/DataModel/AliasTables.h
@@ -47,6 +47,7 @@ using TracksPidPr = soa::Join<aod::pidTPCFullPr, aod::pidTOFFullPr>;
 using TracksPidDe = soa::Join<aod::pidTPCFullDe, aod::pidTOFFullDe>;
 using TracksPidTr = soa::Join<aod::pidTPCFullTr, aod::pidTOFFullTr>;
 using TracksPidHe = soa::Join<aod::pidTPCFullHe, aod::pidTOFFullHe>;
+using TracksPidAl = soa::Join<aod::pidTPCFullAl, aod::pidTOFFullAl>;
 
 using TracksPidTinyEl = soa::Join<aod::pidTPCEl, aod::pidTOFEl>;
 using TracksPidTinyMu = soa::Join<aod::pidTPCMu, aod::pidTOFMu>;
@@ -56,6 +57,7 @@ using TracksPidTinyPr = soa::Join<aod::pidTPCPr, aod::pidTOFPr>;
 using TracksPidTinyDe = soa::Join<aod::pidTPCDe, aod::pidTOFDe>;
 using TracksPidTinyTr = soa::Join<aod::pidTPCTr, aod::pidTOFTr>;
 using TracksPidTinyHe = soa::Join<aod::pidTPCHe, aod::pidTOFHe>;
+using TracksPidTinyAl = soa::Join<aod::pidTPCAl, aod::pidTOFAl>;
 } // namespace o2::aod
 
 #endif // PWGHF_DATAMODEL_ALIASTABLES_H_
diff --git a/PWGHF/DataModel/CandidateReconstructionTables.h b/PWGHF/DataModel/CandidateReconstructionTables.h
index 3af629eb0e4..50c51aec9ef 100644
--- a/PWGHF/DataModel/CandidateReconstructionTables.h
+++ b/PWGHF/DataModel/CandidateReconstructionTables.h
@@ -46,6 +46,7 @@ DECLARE_SOA_COLUMN(TpcTofNSigmaPr, tpcTofNSigmaPr, float); //! Combined NSigma s
 DECLARE_SOA_COLUMN(TpcTofNSigmaDe, tpcTofNSigmaDe, float); //! Combined NSigma separation with the TPC & TOF detectors for deuteron
 DECLARE_SOA_COLUMN(TpcTofNSigmaTr, tpcTofNSigmaTr, float); //! Combined NSigma separation with the TPC & TOF detectors for triton
 DECLARE_SOA_COLUMN(TpcTofNSigmaHe, tpcTofNSigmaHe, float); //! Combined NSigma separation with the TPC & TOF detectors for helium
+DECLARE_SOA_COLUMN(TpcTofNSigmaAl, tpcTofNSigmaAl, float); //! Combined NSigma separation with the TPC & TOF detectors for alpha
 } // namespace pid_tpc_tof_static_full
 
 namespace pid_tpc_tof_static_tiny
@@ -59,6 +60,7 @@ DECLARE_SOA_COLUMN(TpcTofNSigmaPr, tpcTofNSigmaPr, float); //! Combined NSigma s
 DECLARE_SOA_COLUMN(TpcTofNSigmaDe, tpcTofNSigmaDe, float); //! Combined NSigma separation with the TPC & TOF detectors for deuteron
 DECLARE_SOA_COLUMN(TpcTofNSigmaTr, tpcTofNSigmaTr, float); //! Combined NSigma separation with the TPC & TOF detectors for triton
 DECLARE_SOA_COLUMN(TpcTofNSigmaHe, tpcTofNSigmaHe, float); //! Combined NSigma separation with the TPC & TOF detectors for helium
+DECLARE_SOA_COLUMN(TpcTofNSigmaAl, tpcTofNSigmaAl, float); //! Combined NSigma separation with the TPC & TOF detectors for alpha
 } // namespace pid_tpc_tof_static_tiny
 
 // Extension of per particle tables
@@ -72,12 +74,14 @@ DECLARE_SOA_TABLE(PidTpcTofFullKa, "AOD", "PIDTPCTOFFULLKA", //! Table of the TP
                   pid_tpc_tof_static_full::TpcTofNSigmaKa);
 DECLARE_SOA_TABLE(PidTpcTofFullPr, "AOD", "PIDTPCTOFFULLPR", //! Table of the TPC & TOF Combined NSigma for proton
                   pid_tpc_tof_static_full::TpcTofNSigmaPr);
-DECLARE_SOA_TABLE(PidTpcTofFullDe, "AOD", "PIDTPCTOFFULLDe", //! Table of the TPC & TOF Combined NSigma for deuteron
+DECLARE_SOA_TABLE(PidTpcTofFullDe, "AOD", "PIDTPCTOFFULLDE", //! Table of the TPC & TOF Combined NSigma for deuteron
                   pid_tpc_tof_static_full::TpcTofNSigmaDe);
-DECLARE_SOA_TABLE(PidTpcTofFullTr, "AOD", "PIDTPCTOFFULLTr", //! Table of the TPC & TOF Combined NSigma for triton
+DECLARE_SOA_TABLE(PidTpcTofFullTr, "AOD", "PIDTPCTOFFULLTR", //! Table of the TPC & TOF Combined NSigma for triton
                   pid_tpc_tof_static_full::TpcTofNSigmaTr);
-DECLARE_SOA_TABLE(PidTpcTofFullHe, "AOD", "PIDTPCTOFFULLHe", //! Table of the TPC & TOF Combined NSigma for helium
+DECLARE_SOA_TABLE(PidTpcTofFullHe, "AOD", "PIDTPCTOFFULLHE", //! Table of the TPC & TOF Combined NSigma for helium
                   pid_tpc_tof_static_full::TpcTofNSigmaHe);
+DECLARE_SOA_TABLE(PidTpcTofFullAl, "AOD", "PIDTPCTOFFULLAL", //! Table of the TPC & TOF Combined NSigma for alpha
+                  pid_tpc_tof_static_full::TpcTofNSigmaAl);
 // Extension of per particle tables
 DECLARE_SOA_TABLE(PidTpcTofTinyEl, "AOD", "PIDTPCTOFTINYEL", //! Table of the TPC & TOF Combined NSigma for electron
                   pid_tpc_tof_static_tiny::TpcTofNSigmaEl);
@@ -91,10 +95,12 @@ DECLARE_SOA_TABLE(PidTpcTofTinyPr, "AOD", "PIDTPCTOFTINYPR", //! Table of the TP
                   pid_tpc_tof_static_tiny::TpcTofNSigmaPr);
 DECLARE_SOA_TABLE(PidTpcTofTinyDe, "AOD", "PIDTPCTOFTINYDE", //! Table of the TPC & TOF Combined NSigma for deuteron
                   pid_tpc_tof_static_tiny::TpcTofNSigmaDe);
-DECLARE_SOA_TABLE(PidTpcTofTinyTr, "AOD", "PIDTPCTOFTINYTr", //! Table of the TPC & TOF Combined NSigma for triton
+DECLARE_SOA_TABLE(PidTpcTofTinyTr, "AOD", "PIDTPCTOFTINYTR", //! Table of the TPC & TOF Combined NSigma for triton
                   pid_tpc_tof_static_tiny::TpcTofNSigmaTr);
-DECLARE_SOA_TABLE(PidTpcTofTinyHe, "AOD", "PIDTPCTOFTINYHe", //! Table of the TPC & TOF Combined NSigma for helium
+DECLARE_SOA_TABLE(PidTpcTofTinyHe, "AOD", "PIDTPCTOFTINYHE", //! Table of the TPC & TOF Combined NSigma for helium
                   pid_tpc_tof_static_tiny::TpcTofNSigmaHe);
+DECLARE_SOA_TABLE(PidTpcTofTinyAl, "AOD", "PIDTPCTOFTINYAL", //! Table of the TPC & TOF Combined NSigma for alpha
+                  pid_tpc_tof_static_tiny::TpcTofNSigmaAl);
 // general decay properties
 namespace hf_cand
 {
@@ -197,6 +203,9 @@ DECLARE_SOA_COLUMN(NSigTpcTr2, nSigTpcTr2, float);           //! TPC nSigma for
 DECLARE_SOA_COLUMN(NSigTpcHe0, nSigTpcHe0, float);           //! TPC nSigma for helium hypothesis - prong 0
 DECLARE_SOA_COLUMN(NSigTpcHe1, nSigTpcHe1, float);           //! TPC nSigma for helium hypothesis - prong 1
 DECLARE_SOA_COLUMN(NSigTpcHe2, nSigTpcHe2, float);           //! TPC nSigma for helium hypothesis - prong 2
+DECLARE_SOA_COLUMN(NSigTpcAl0, nSigTpcAl0, float);           //! TPC nSigma for alpha hypothesis - prong 0
+DECLARE_SOA_COLUMN(NSigTpcAl1, nSigTpcAl1, float);           //! TPC nSigma for alpha hypothesis - prong 1
+DECLARE_SOA_COLUMN(NSigTpcAl2, nSigTpcAl2, float);           //! TPC nSigma for alpha hypothesis - prong 2
 DECLARE_SOA_COLUMN(NSigTofPi0, nSigTofPi0, float);           //! TOF nSigma for pion hypothesis - prong 0
 DECLARE_SOA_COLUMN(NSigTofPi1, nSigTofPi1, float);           //! TOF nSigma for pion hypothesis - prong 1
 DECLARE_SOA_COLUMN(NSigTofPi2, nSigTofPi2, float);           //! TOF nSigma for pion hypothesis - prong 2
@@ -215,6 +224,9 @@ DECLARE_SOA_COLUMN(NSigTofTr2, nSigTofTr2, float);           //! TOF nSigma for
 DECLARE_SOA_COLUMN(NSigTofHe0, nSigTofHe0, float);           //! TOF nSigma for helium hypothesis - prong 0
 DECLARE_SOA_COLUMN(NSigTofHe1, nSigTofHe1, float);           //! TOF nSigma for helium hypothesis - prong 1
 DECLARE_SOA_COLUMN(NSigTofHe2, nSigTofHe2, float);           //! TOF nSigma for helium hypothesis - prong 2
+DECLARE_SOA_COLUMN(NSigTofAl0, nSigTofAl0, float);           //! TOF nSigma for alpha hypothesis - prong 0
+DECLARE_SOA_COLUMN(NSigTofAl1, nSigTofAl1, float);           //! TOF nSigma for alpha hypothesis - prong 1
+DECLARE_SOA_COLUMN(NSigTofAl2, nSigTofAl2, float);           //! TOF nSigma for alpha hypothesis - prong 2
 DECLARE_SOA_DYNAMIC_COLUMN(TpcTofNSigmaPi0, tpcTofNSigmaPi0, //! Combined NSigma separation with the TPC & TOF detectors for pion - prong 0
                            [](float tpcNSigmaPi0, float tofNSigmaPi0) -> float { return pid_tpc_tof_utils::combineNSigma<false /*tiny*/>(tpcNSigmaPi0, tofNSigmaPi0); });
 DECLARE_SOA_DYNAMIC_COLUMN(TpcTofNSigmaPi1, tpcTofNSigmaPi1, //! Combined NSigma separation with the TPC & TOF detectors for pion - prong 1
@@ -251,6 +263,12 @@ DECLARE_SOA_DYNAMIC_COLUMN(TpcTofNSigmaHe1, tpcTofNSigmaHe1, //! Combined NSigma
                            [](float tpcNSigmaHe1, float tofNSigmaHe1) -> float { return pid_tpc_tof_utils::combineNSigma<false /*tiny*/>(tpcNSigmaHe1, tofNSigmaHe1); });
 DECLARE_SOA_DYNAMIC_COLUMN(TpcTofNSigmaHe2, tpcTofNSigmaHe2, //! Combined NSigma separation with the TPC & TOF detectors for helium - prong 2
                            [](float tpcNSigmaHe2, float tofNSigmaHe2) -> float { return pid_tpc_tof_utils::combineNSigma<false /*tiny*/>(tpcNSigmaHe2, tofNSigmaHe2); });
+DECLARE_SOA_DYNAMIC_COLUMN(TpcTofNSigmaAl0, tpcTofNSigmaAl0, //! Combined NSigma separation with the TPC & TOF detectors for alpha - prong 0
+                           [](float tpcNSigmaAl0, float tofNSigmaAl0) -> float { return pid_tpc_tof_utils::combineNSigma<false /*tiny*/>(tpcNSigmaAl0, tofNSigmaAl0); });
+DECLARE_SOA_DYNAMIC_COLUMN(TpcTofNSigmaAl1, tpcTofNSigmaAl1, //! Combined NSigma separation with the TPC & TOF detectors for alpha - prong 1
+                           [](float tpcNSigmaAl1, float tofNSigmaAl1) -> float { return pid_tpc_tof_utils::combineNSigma<false /*tiny*/>(tpcNSigmaAl1, tofNSigmaAl1); });
+DECLARE_SOA_DYNAMIC_COLUMN(TpcTofNSigmaAl2, tpcTofNSigmaAl2, //! Combined NSigma separation with the TPC & TOF detectors for alpha - prong 2
+                           [](float tpcNSigmaAl2, float tofNSigmaAl2) -> float { return pid_tpc_tof_utils::combineNSigma<false /*tiny*/>(tpcNSigmaAl2, tofNSigmaAl2); });
 // tiny (binned) option
 DECLARE_SOA_DYNAMIC_COLUMN(TpcTofNSigmaTinyPi0, tpcTofNSigmaTinyPi0, //! Combined NSigma separation with the TPC & TOF detectors for pion - prong 0
                            [](float tpcNSigmaPi0, float tofNSigmaPi0) -> float { return pid_tpc_tof_utils::combineNSigma<true /*tiny*/>(tpcNSigmaPi0, tofNSigmaPi0); });
@@ -809,15 +827,24 @@ DECLARE_SOA_TABLE(HfCand3Prong1PidTr, "AOD", "HFCAND3P1PIDTR", //!
 DECLARE_SOA_TABLE(HfCand3Prong2PidTr, "AOD", "HFCAND3P2PIDTR", //!
                   hf_cand::NSigTpcTr2, hf_cand::NSigTofTr2,
                   hf_cand::TpcTofNSigmaTr2<hf_cand::NSigTpcTr2, hf_cand::NSigTofTr2>);
-DECLARE_SOA_TABLE(HfCand3Prong0PidHe, "AOD", "HFCAND3P0PIDHe", //!
+DECLARE_SOA_TABLE(HfCand3Prong0PidHe, "AOD", "HFCAND3P0PIDHE", //!
                   hf_cand::NSigTpcHe0, hf_cand::NSigTofHe0,
-                  hf_cand::TpcTofNSigmaTr0<hf_cand::NSigTpcHe0, hf_cand::NSigTofHe0>);
-DECLARE_SOA_TABLE(HfCand3Prong1PidHe, "AOD", "HFCAND3P1PIDHe", //!
+                  hf_cand::TpcTofNSigmaHe0<hf_cand::NSigTpcHe0, hf_cand::NSigTofHe0>);
+DECLARE_SOA_TABLE(HfCand3Prong1PidHe, "AOD", "HFCAND3P1PIDHE", //!
                   hf_cand::NSigTpcHe1, hf_cand::NSigTofHe1,
                   hf_cand::TpcTofNSigmaHe1<hf_cand::NSigTpcHe1, hf_cand::NSigTofHe1>);
-DECLARE_SOA_TABLE(HfCand3Prong2PidHe, "AOD", "HFCAND3P2PIDHe", //!
+DECLARE_SOA_TABLE(HfCand3Prong2PidHe, "AOD", "HFCAND3P2PIDHE", //!
                   hf_cand::NSigTpcHe2, hf_cand::NSigTofHe2,
                   hf_cand::TpcTofNSigmaHe2<hf_cand::NSigTpcHe2, hf_cand::NSigTofHe2>);
+DECLARE_SOA_TABLE(HfCand3Prong0PidAl, "AOD", "HFCAND3P0PIDAL", //!
+                  hf_cand::NSigTpcAl0, hf_cand::NSigTofAl0,
+                  hf_cand::TpcTofNSigmaAl0<hf_cand::NSigTpcAl0, hf_cand::NSigTofAl0>);
+DECLARE_SOA_TABLE(HfCand3Prong1PidAl, "AOD", "HFCAND3P1PIDAL", //!
+                  hf_cand::NSigTpcAl1, hf_cand::NSigTofAl1,
+                  hf_cand::TpcTofNSigmaAl1<hf_cand::NSigTpcAl1, hf_cand::NSigTofAl1>);
+DECLARE_SOA_TABLE(HfCand3Prong2PidAl, "AOD", "HFCAND3P2PIDAL", //!
+                  hf_cand::NSigTpcAl2, hf_cand::NSigTofAl2,
+                  hf_cand::TpcTofNSigmaAl2<hf_cand::NSigTpcAl2, hf_cand::NSigTofAl2>);
 
 using HfCand3Prong = HfCand3ProngExt;
 using HfCand3ProngWPidPiKaPr = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3Prong0PidPr, HfCand3Prong0PidKa, HfCand3Prong1PidPi, HfCand3Prong1PidPr, HfCand3Prong1PidKa, HfCand3Prong2PidPi, HfCand3Prong2PidPr, HfCand3Prong2PidKa>;
@@ -825,6 +852,7 @@ using HfCand3ProngWPidPiKa = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3
 using HfCand3ProngWPidPiKaDe = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3Prong0PidDe, HfCand3Prong0PidKa, HfCand3Prong1PidPi, HfCand3Prong1PidDe, HfCand3Prong1PidKa, HfCand3Prong2PidPi, HfCand3Prong2PidDe, HfCand3Prong2PidKa>;
 using HfCand3ProngWPidPiKaTr = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3Prong0PidTr, HfCand3Prong0PidKa, HfCand3Prong1PidPi, HfCand3Prong1PidTr, HfCand3Prong1PidKa, HfCand3Prong2PidPi, HfCand3Prong2PidTr, HfCand3Prong2PidKa>;
 using HfCand3ProngWPidPiKaHe = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3Prong0PidHe, HfCand3Prong0PidKa, HfCand3Prong1PidPi, HfCand3Prong1PidHe, HfCand3Prong1PidKa, HfCand3Prong2PidPi, HfCand3Prong2PidHe, HfCand3Prong2PidKa>;
+using HfCand3ProngWPidPiKaAl = soa::Join<HfCand3Prong, HfCand3Prong0PidPi, HfCand3Prong0PidAl, HfCand3Prong0PidKa, HfCand3Prong1PidPi, HfCand3Prong1PidAl, HfCand3Prong1PidKa, HfCand3Prong2PidPi, HfCand3Prong2PidAl, HfCand3Prong2PidKa>;
 
 DECLARE_SOA_TABLE(HfCand3ProngKF, "AOD", "HFCAND3PKF",
                   hf_cand_3prong::KfXError, hf_cand_3prong::KfYError, hf_cand_3prong::KfZError,
diff --git a/PWGHF/TableProducer/candidateCreator3Prong.cxx b/PWGHF/TableProducer/candidateCreator3Prong.cxx
index b88ea949b3e..334f185a169 100644
--- a/PWGHF/TableProducer/candidateCreator3Prong.cxx
+++ b/PWGHF/TableProducer/candidateCreator3Prong.cxx
@@ -21,6 +21,7 @@
 
 #include "PWGHF/Core/CentralityEstimation.h"
 #include "PWGHF/Core/DecayChannels.h"
+#include "PWGHF/Core/SelectorCuts.h"
 #include "PWGHF/DataModel/AliasTables.h"
 #include "PWGHF/DataModel/CandidateReconstructionTables.h"
 #include "PWGHF/DataModel/TrackIndexSkimmingTables.h"
@@ -78,6 +79,7 @@
 #include <vector>
 
 using namespace o2;
+using namespace o2::analysis;
 using namespace o2::hf_evsel;
 using namespace o2::hf_trkcandsel;
 using namespace o2::aod::hf_cand_3prong;
@@ -98,14 +100,23 @@ struct HfCandidateCreator3Prong {
   Produces<aod::HfCand3Prong0PidKa> rowProng0PidKa;
   Produces<aod::HfCand3Prong0PidPr> rowProng0PidPr;
   Produces<aod::HfCand3Prong0PidDe> rowProng0PidDe;
+  Produces<aod::HfCand3Prong0PidTr> rowProng0PidTr;
+  Produces<aod::HfCand3Prong0PidHe> rowProng0PidHe;
+  Produces<aod::HfCand3Prong0PidAl> rowProng0PidAl;
   Produces<aod::HfCand3Prong1PidPi> rowProng1PidPi;
   Produces<aod::HfCand3Prong1PidKa> rowProng1PidKa;
   Produces<aod::HfCand3Prong1PidPr> rowProng1PidPr;
   Produces<aod::HfCand3Prong1PidDe> rowProng1PidDe;
+  Produces<aod::HfCand3Prong1PidTr> rowProng1PidTr;
+  Produces<aod::HfCand3Prong1PidHe> rowProng1PidHe;
+  Produces<aod::HfCand3Prong1PidAl> rowProng1PidAl;
   Produces<aod::HfCand3Prong2PidPi> rowProng2PidPi;
   Produces<aod::HfCand3Prong2PidKa> rowProng2PidKa;
   Produces<aod::HfCand3Prong2PidPr> rowProng2PidPr;
   Produces<aod::HfCand3Prong2PidDe> rowProng2PidDe;
+  Produces<aod::HfCand3Prong2PidTr> rowProng2PidTr;
+  Produces<aod::HfCand3Prong2PidHe> rowProng2PidHe;
+  Produces<aod::HfCand3Prong2PidAl> rowProng2PidAl;
 
   // vertexing
   Configurable<bool> propagateToPCA{"propagateToPCA", true, "create tracks version propagated to PCA"};
@@ -126,11 +137,13 @@ struct HfCandidateCreator3Prong {
   Configurable<bool> createDs{"createDs", false, "enable Ds+/- candidate creation"};
   Configurable<bool> createLc{"createLc", false, "enable Lc+/- candidate creation"};
   Configurable<bool> createXic{"createXic", false, "enable Xic+/- candidate creation"};
-  Configurable<bool> createCd{"createCd", false, "enable Cd candidate creation"};
+  Configurable<bool> createCharmNuclei{"createCharmNuclei", false, "enable createCharmNuclei candidate creation"};
   // KF
   Configurable<bool> applyTopoConstraint{"applyTopoConstraint", false, "apply origin from PV hypothesis for created candidate, works only in KF mode"};
   Configurable<bool> applyInvMassConstraint{"applyInvMassConstraint", false, "apply particle type hypothesis to recalculate created candidate's momentum, works only in KF mode"};
 
+  Configurable<LabeledArray<float>> tpcPidBBParamsLightNuclei{"tpcPidBBParamsLightNuclei", {hf_presel_lightnuclei::BetheBlochParams[0], hf_presel_lightnuclei::NParticleRows, hf_presel_lightnuclei::NBetheBlochParams, hf_presel_lightnuclei::labelsRowsNucleiType, hf_presel_lightnuclei::labelsBetheBlochParams}, "TPC PID Bethe–Bloch parameter configurations for light nuclei (deuteron, triton, helium-3)"};
+
   HfEventSelection hfEvSel;        // event selection and monitoring
   o2::vertexing::DCAFitterN<3> df; // 3-prong vertex fitter
   Service<o2::ccdb::BasicCCDBManager> ccdb{};
@@ -143,10 +156,10 @@ struct HfCandidateCreator3Prong {
 
   using FilteredHf3Prongs = soa::Filtered<aod::Hf3Prongs>;
   using FilteredPvRefitHf3Prongs = soa::Filtered<soa::Join<aod::Hf3Prongs, aod::HfPvRefit3Prong>>;
-  using TracksWCovExtraPidPiKaPrDe = soa::Join<aod::TracksWCovExtra, aod::TracksPidPi, aod::PidTpcTofFullPi, aod::TracksPidKa, aod::PidTpcTofFullKa, aod::TracksPidPr, aod::PidTpcTofFullPr, aod::TracksPidDe, aod::PidTpcTofFullDe>;
+  using TracksWCovExtraPidPiKaPrLightNuclei = soa::Join<aod::TracksWCovExtra, aod::TracksPidPi, aod::PidTpcTofFullPi, aod::TracksPidKa, aod::PidTpcTofFullKa, aod::TracksPidPr, aod::PidTpcTofFullPr, aod::TracksPidDe, aod::PidTpcTofFullDe, aod::TracksPidHe, aod::PidTpcTofFullHe, aod::TracksPidTr, aod::PidTpcTofFullTr, aod::TracksPidAl, aod::PidTpcTofFullAl>;
 
   // filter candidates
-  Filter filterSelected3Prongs = (createDplus && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::DplusToPiKPi))) != static_cast<uint8_t>(0)) || (createDs && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::DsToKKPi))) != static_cast<uint8_t>(0)) || (createLc && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::LcToPKPi))) != static_cast<uint8_t>(0)) || (createXic && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::XicToPKPi))) != static_cast<uint8_t>(0)) || (createCd && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::CdToDeKPi))) != static_cast<uint8_t>(0));
+  Filter filterSelected3Prongs = (createDplus && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::DplusToPiKPi))) != static_cast<uint8_t>(0)) || (createDs && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::DsToKKPi))) != static_cast<uint8_t>(0)) || (createLc && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::LcToPKPi))) != static_cast<uint8_t>(0)) || (createXic && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::XicToPKPi))) != static_cast<uint8_t>(0)) || (createCharmNuclei && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::CdToDeKPi))) != static_cast<uint8_t>(0)) || (createCharmNuclei && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::CtToTrKPi))) != static_cast<uint8_t>(0)) || (createCharmNuclei && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::ChToHeKPi))) != static_cast<uint8_t>(0)) || (createCharmNuclei && (o2::aod::hf_track_index::hfflag & static_cast<uint8_t>(BIT(DecayType::CaToAlKPi))) != static_cast<uint8_t>(0));
 
   std::shared_ptr<TH1> hCandidates;
   HistogramRegistry registry{"registry"};
@@ -189,7 +202,7 @@ struct HfCandidateCreator3Prong {
     if (nProcessesUpc > 0 && isRun2) {
       LOGP(fatal, "Process function for UPC is only available in Run 3!");
     }
-    std::array<bool, 5> creationFlags = {createDplus, createDs, createLc, createXic, createCd};
+    std::array<bool, 5> creationFlags = {createDplus, createDs, createLc, createXic, createCharmNuclei};
     if (std::accumulate(creationFlags.begin(), creationFlags.end(), 0) == 0) {
       LOGP(fatal, "At least one particle specie should be enabled for the creation.");
     }
@@ -258,17 +271,26 @@ struct HfCandidateCreator3Prong {
       fillProngPid<HfProngSpecies::Proton>(track1, rowProng1PidPr);
       fillProngPid<HfProngSpecies::Proton>(track2, rowProng2PidPr);
     }
-    if (createCd) {
-      fillProngPid<HfProngSpecies::Deuteron>(track0, rowProng0PidDe);
-      fillProngPid<HfProngSpecies::Deuteron>(track1, rowProng1PidDe);
-      fillProngPid<HfProngSpecies::Deuteron>(track2, rowProng2PidDe);
+    if (createCharmNuclei) {
+      fillProngPidLightNuclei<HfProngSpecies::Deuteron>(track0, rowProng0PidDe, tpcPidBBParamsLightNuclei);
+      fillProngPidLightNuclei<HfProngSpecies::Deuteron>(track1, rowProng1PidDe, tpcPidBBParamsLightNuclei);
+      fillProngPidLightNuclei<HfProngSpecies::Deuteron>(track2, rowProng2PidDe, tpcPidBBParamsLightNuclei);
+      fillProngPidLightNuclei<HfProngSpecies::Helium3>(track0, rowProng0PidHe, tpcPidBBParamsLightNuclei);
+      fillProngPidLightNuclei<HfProngSpecies::Helium3>(track1, rowProng1PidHe, tpcPidBBParamsLightNuclei);
+      fillProngPidLightNuclei<HfProngSpecies::Helium3>(track2, rowProng2PidHe, tpcPidBBParamsLightNuclei);
+      fillProngPidLightNuclei<HfProngSpecies::Triton>(track0, rowProng0PidTr, tpcPidBBParamsLightNuclei);
+      fillProngPidLightNuclei<HfProngSpecies::Triton>(track1, rowProng1PidTr, tpcPidBBParamsLightNuclei);
+      fillProngPidLightNuclei<HfProngSpecies::Triton>(track2, rowProng2PidTr, tpcPidBBParamsLightNuclei);
+      fillProngPidLightNuclei<HfProngSpecies::Alpha>(track0, rowProng0PidAl, tpcPidBBParamsLightNuclei);
+      fillProngPidLightNuclei<HfProngSpecies::Alpha>(track1, rowProng1PidAl, tpcPidBBParamsLightNuclei);
+      fillProngPidLightNuclei<HfProngSpecies::Alpha>(track2, rowProng2PidAl, tpcPidBBParamsLightNuclei);
     }
   }
 
   template <bool DoPvRefit, bool ApplyUpcSel, o2::hf_centrality::CentralityEstimator CentEstimator, typename Coll, typename Cand, typename BCsType>
   void runCreator3ProngWithDCAFitterN(Coll const&,
                                       Cand const& rowsTrackIndexProng3,
-                                      TracksWCovExtraPidPiKaPrDe const&,
+                                      TracksWCovExtraPidPiKaPrLightNuclei const&,
                                       BCsType const& bcs)
   {
     // loop over triplets of track indices
@@ -288,9 +310,9 @@ struct HfCandidateCreator3Prong {
         continue;
       }
 
-      auto track0 = rowTrackIndexProng3.template prong0_as<TracksWCovExtraPidPiKaPrDe>();
-      auto track1 = rowTrackIndexProng3.template prong1_as<TracksWCovExtraPidPiKaPrDe>();
-      auto track2 = rowTrackIndexProng3.template prong2_as<TracksWCovExtraPidPiKaPrDe>();
+      auto track0 = rowTrackIndexProng3.template prong0_as<TracksWCovExtraPidPiKaPrLightNuclei>();
+      auto track1 = rowTrackIndexProng3.template prong1_as<TracksWCovExtraPidPiKaPrLightNuclei>();
+      auto track2 = rowTrackIndexProng3.template prong2_as<TracksWCovExtraPidPiKaPrLightNuclei>();
       auto trackParVar0 = getTrackParCov(track0);
       auto trackParVar1 = getTrackParCov(track1);
       auto trackParVar2 = getTrackParCov(track2);
@@ -445,7 +467,7 @@ struct HfCandidateCreator3Prong {
   template <bool DoPvRefit, bool ApplyUpcSel, o2::hf_centrality::CentralityEstimator CentEstimator, typename Coll, typename Cand, typename BCsType>
   void runCreator3ProngWithKFParticle(Coll const&,
                                       Cand const& rowsTrackIndexProng3,
-                                      TracksWCovExtraPidPiKaPrDe const&,
+                                      TracksWCovExtraPidPiKaPrLightNuclei const&,
                                       BCsType const& bcs)
   {
     for (const auto& rowTrackIndexProng3 : rowsTrackIndexProng3) {
@@ -463,9 +485,9 @@ struct HfCandidateCreator3Prong {
         continue;
       }
 
-      auto track0 = rowTrackIndexProng3.template prong0_as<TracksWCovExtraPidPiKaPrDe>();
-      auto track1 = rowTrackIndexProng3.template prong1_as<TracksWCovExtraPidPiKaPrDe>();
-      auto track2 = rowTrackIndexProng3.template prong2_as<TracksWCovExtraPidPiKaPrDe>();
+      auto track0 = rowTrackIndexProng3.template prong0_as<TracksWCovExtraPidPiKaPrLightNuclei>();
+      auto track1 = rowTrackIndexProng3.template prong1_as<TracksWCovExtraPidPiKaPrLightNuclei>();
+      auto track2 = rowTrackIndexProng3.template prong2_as<TracksWCovExtraPidPiKaPrLightNuclei>();
 
       /// Set the magnetic field from ccdb.
       /// The static instance of the propagator was already modified in the HFTrackIndexSkimCreator,
@@ -694,7 +716,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using DCA fitter  w/ PV refit and w/o centrality selections
   void processPvRefitWithDCAFitterN(soa::Join<aod::Collisions, aod::EvSels> const& collisions,
                                     FilteredPvRefitHf3Prongs const& rowsTrackIndexProng3,
-                                    TracksWCovExtraPidPiKaPrDe const& tracks,
+                                    TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                     aod::BCsWithTimestamps const& bcWithTimeStamps)
   {
     runCreator3ProngWithDCAFitterN</*doPvRefit*/ true, false, CentralityEstimator::None>(collisions, rowsTrackIndexProng3, tracks, bcWithTimeStamps);
@@ -704,7 +726,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using DCA fitter  w/o PV refit and w/o centrality selections
   void processNoPvRefitWithDCAFitterN(soa::Join<aod::Collisions, aod::EvSels> const& collisions,
                                       FilteredHf3Prongs const& rowsTrackIndexProng3,
-                                      TracksWCovExtraPidPiKaPrDe const& tracks,
+                                      TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                       aod::BCsWithTimestamps const& bcWithTimeStamps)
   {
     runCreator3ProngWithDCAFitterN</*doPvRefit*/ false, false, CentralityEstimator::None>(collisions, rowsTrackIndexProng3, tracks, bcWithTimeStamps);
@@ -714,7 +736,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using KFParticle package  w/ PV refit and w/o centrality selections
   void processPvRefitWithKFParticle(soa::Join<aod::Collisions, aod::EvSels> const& collisions,
                                     FilteredPvRefitHf3Prongs const& rowsTrackIndexProng3,
-                                    TracksWCovExtraPidPiKaPrDe const& tracks,
+                                    TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                     aod::BCsWithTimestamps const& bcWithTimeStamps)
   {
     runCreator3ProngWithKFParticle</*doPvRefit*/ true, false, CentralityEstimator::None>(collisions, rowsTrackIndexProng3, tracks, bcWithTimeStamps);
@@ -724,7 +746,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using KFParticle package  w/o PV refit and w/o centrality selections
   void processNoPvRefitWithKFParticle(soa::Join<aod::Collisions, aod::EvSels> const& collisions,
                                       FilteredHf3Prongs const& rowsTrackIndexProng3,
-                                      TracksWCovExtraPidPiKaPrDe const& tracks,
+                                      TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                       aod::BCsWithTimestamps const& bcWithTimeStamps)
   {
     runCreator3ProngWithKFParticle</*doPvRefit*/ false, false, CentralityEstimator::None>(collisions, rowsTrackIndexProng3, tracks, bcWithTimeStamps);
@@ -740,7 +762,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using DCA fitter  w/ PV refit and w/ centrality selection on FT0C
   void processPvRefitWithDCAFitterNCentFT0C(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs> const& collisions,
                                             FilteredPvRefitHf3Prongs const& rowsTrackIndexProng3,
-                                            TracksWCovExtraPidPiKaPrDe const& tracks,
+                                            TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                             aod::BCsWithTimestamps const& bcWithTimeStamps)
   {
     runCreator3ProngWithDCAFitterN</*doPvRefit*/ true, false, CentralityEstimator::FT0C>(collisions, rowsTrackIndexProng3, tracks, bcWithTimeStamps);
@@ -750,7 +772,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using DCA fitter  w/o PV refit and  w/ centrality selection on FT0C
   void processNoPvRefitWithDCAFitterNCentFT0C(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs> const& collisions,
                                               FilteredHf3Prongs const& rowsTrackIndexProng3,
-                                              TracksWCovExtraPidPiKaPrDe const& tracks,
+                                              TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                               aod::BCsWithTimestamps const& bcWithTimeStamps)
   {
     runCreator3ProngWithDCAFitterN</*doPvRefit*/ false, false, CentralityEstimator::FT0C>(collisions, rowsTrackIndexProng3, tracks, bcWithTimeStamps);
@@ -760,7 +782,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using KFParticle package  w/ PV refit and w/ centrality selection on FT0C
   void processPvRefitWithKFParticleCentFT0C(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs> const& collisions,
                                             FilteredPvRefitHf3Prongs const& rowsTrackIndexProng3,
-                                            TracksWCovExtraPidPiKaPrDe const& tracks,
+                                            TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                             aod::BCsWithTimestamps const& bcWithTimeStamps)
   {
     runCreator3ProngWithKFParticle</*doPvRefit*/ true, false, CentralityEstimator::FT0C>(collisions, rowsTrackIndexProng3, tracks, bcWithTimeStamps);
@@ -770,7 +792,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using KFParticle package  w/o PV refit and  w/ centrality selection on FT0C
   void processNoPvRefitWithKFParticleCentFT0C(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs> const& collisions,
                                               FilteredHf3Prongs const& rowsTrackIndexProng3,
-                                              TracksWCovExtraPidPiKaPrDe const& tracks,
+                                              TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                               aod::BCsWithTimestamps const& bcWithTimeStamps)
   {
     runCreator3ProngWithKFParticle</*doPvRefit*/ false, false, CentralityEstimator::FT0C>(collisions, rowsTrackIndexProng3, tracks, bcWithTimeStamps);
@@ -786,7 +808,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using DCA fitter  w/ PV refit and w/ centrality selection on FT0M
   void processPvRefitWithDCAFitterNCentFT0M(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
                                             FilteredPvRefitHf3Prongs const& rowsTrackIndexProng3,
-                                            TracksWCovExtraPidPiKaPrDe const& tracks,
+                                            TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                             aod::BCsWithTimestamps const& bcWithTimeStamps)
   {
     runCreator3ProngWithDCAFitterN</*doPvRefit*/ true, false, CentralityEstimator::FT0M>(collisions, rowsTrackIndexProng3, tracks, bcWithTimeStamps);
@@ -796,7 +818,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using DCA fitter  w/o PV refit and  w/ centrality selection on FT0M
   void processNoPvRefitWithDCAFitterNCentFT0M(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
                                               FilteredHf3Prongs const& rowsTrackIndexProng3,
-                                              TracksWCovExtraPidPiKaPrDe const& tracks,
+                                              TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                               aod::BCsWithTimestamps const& bcWithTimeStamps)
   {
     runCreator3ProngWithDCAFitterN</*doPvRefit*/ false, false, CentralityEstimator::FT0M>(collisions, rowsTrackIndexProng3, tracks, bcWithTimeStamps);
@@ -806,7 +828,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using KFParticle package  w/ PV refit and w/ centrality selection on FT0M
   void processPvRefitWithKFParticleCentFT0M(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
                                             FilteredPvRefitHf3Prongs const& rowsTrackIndexProng3,
-                                            TracksWCovExtraPidPiKaPrDe const& tracks,
+                                            TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                             aod::BCsWithTimestamps const& bcWithTimeStamps)
   {
     runCreator3ProngWithKFParticle</*doPvRefit*/ true, false, CentralityEstimator::FT0M>(collisions, rowsTrackIndexProng3, tracks, bcWithTimeStamps);
@@ -816,7 +838,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using KFParticle package  w/o PV refit and  w/ centrality selection on FT0M
   void processNoPvRefitWithKFParticleCentFT0M(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
                                               FilteredHf3Prongs const& rowsTrackIndexProng3,
-                                              TracksWCovExtraPidPiKaPrDe const& tracks,
+                                              TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                               aod::BCsWithTimestamps const& bcWithTimeStamps)
   {
     runCreator3ProngWithKFParticle</*doPvRefit*/ false, false, CentralityEstimator::FT0M>(collisions, rowsTrackIndexProng3, tracks, bcWithTimeStamps);
@@ -832,7 +854,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using DCA fitter  w/ PV refit and w/ centrality selection on UPC
   void processPvRefitWithDCAFitterNUpc(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
                                        FilteredPvRefitHf3Prongs const& rowsTrackIndexProng3,
-                                       TracksWCovExtraPidPiKaPrDe const& tracks,
+                                       TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                        aod::BcFullInfos const& bcWithTimeStamps,
                                        aod::FT0s const& /*ft0s*/,
                                        aod::FV0As const& /*fv0as*/,
@@ -846,7 +868,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using DCA fitter  w/o PV refit and  w/ centrality selection on UPC
   void processNoPvRefitWithDCAFitterNUpc(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
                                          FilteredHf3Prongs const& rowsTrackIndexProng3,
-                                         TracksWCovExtraPidPiKaPrDe const& tracks,
+                                         TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                          aod::BcFullInfos const& bcWithTimeStamps,
                                          aod::FT0s const& /*ft0s*/,
                                          aod::FV0As const& /*fv0as*/,
@@ -860,7 +882,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using KFParticle package  w/ PV refit and w/ centrality selection on UPC
   void processPvRefitWithKFParticleUpc(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
                                        FilteredPvRefitHf3Prongs const& rowsTrackIndexProng3,
-                                       TracksWCovExtraPidPiKaPrDe const& tracks,
+                                       TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                        aod::BcFullInfos const& bcWithTimeStamps,
                                        aod::FT0s const& /*ft0s*/,
                                        aod::FV0As const& /*fv0as*/,
@@ -874,7 +896,7 @@ struct HfCandidateCreator3Prong {
   /// @brief process function using KFParticle package  w/o PV refit and  w/ centrality selection on UPC
   void processNoPvRefitWithKFParticleUpc(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
                                          FilteredHf3Prongs const& rowsTrackIndexProng3,
-                                         TracksWCovExtraPidPiKaPrDe const& tracks,
+                                         TracksWCovExtraPidPiKaPrLightNuclei const& tracks,
                                          aod::BcFullInfos const& bcWithTimeStamps,
                                          aod::FT0s const& /*ft0s*/,
                                          aod::FV0As const& /*fv0as*/,
diff --git a/PWGHF/TableProducer/mcPidTof.cxx b/PWGHF/TableProducer/mcPidTof.cxx
index b08e912e4c2..dba9eb2753f 100644
--- a/PWGHF/TableProducer/mcPidTof.cxx
+++ b/PWGHF/TableProducer/mcPidTof.cxx
@@ -12,7 +12,7 @@
 ///
 /// \file   mcPidTof.cxx
 /// \author Fabrizio Grosa fabrizio.grosa@cern.ch
-/// \brief  Task to produce PID tables for TOF split for pi, K, p, de, tr, he3 copied from https://github.com/AliceO2Group/O2Physics/blob/master/Common/TableProducer/PID/pidTofMerge.cxx
+/// \brief  Task to produce PID tables for TOF split for pi, K, p, de, tr, he3, alpha copied from https://github.com/AliceO2Group/O2Physics/blob/master/Common/TableProducer/PID/pidTofMerge.cxx
 ///         It works only for MC and adds the possibility to apply postcalibrations for MC.
 ///
 
@@ -722,6 +722,7 @@ static constexpr int IdxPr = 4;
 static constexpr int IdxDe = 5;
 static constexpr int IdxTr = 6;
 static constexpr int IdxHe = 7;
+static constexpr int IdxAl = 8;
 
 /// Task to produce the response table
 struct McPidTof {
@@ -732,6 +733,7 @@ struct McPidTof {
   Produces<o2::aod::pidTOFDe> tablePIDDe;
   Produces<o2::aod::pidTOFTr> tablePIDTr;
   Produces<o2::aod::pidTOFHe> tablePIDHe;
+  Produces<o2::aod::pidTOFAl> tablePIDAl;
 
   // Tables to produce (full)
   Produces<o2::aod::pidTOFFullPi> tablePIDFullPi;
@@ -740,6 +742,7 @@ struct McPidTof {
   Produces<o2::aod::pidTOFFullDe> tablePIDFullDe;
   Produces<o2::aod::pidTOFFullTr> tablePIDFullTr;
   Produces<o2::aod::pidTOFFullHe> tablePIDFullHe;
+  Produces<o2::aod::pidTOFFullAl> tablePIDFullAl;
 
   // Detector response parameters
   o2::pid::tof::TOFResoParamsV3 mRespParamsV3;
@@ -877,6 +880,14 @@ struct McPidTof {
         }
         break;
       }
+      case IdxAl: {
+        if (fullTable) {
+          tablePIDFullAl.reserve(size);
+        } else {
+          tablePIDAl.reserve(size);
+        }
+        break;
+      }
       default:
         LOG(fatal) << "Wrong particle ID in reserveTable() for " << (fullTable ? "full" : "tiny") << " tables";
         break;
@@ -929,6 +940,13 @@ struct McPidTof {
           aod::pidtof_tiny::binning::packInTable(-999.f, tablePIDHe);
         }
         break;
+      case IdxAl:
+        if (fullTable) {
+          tablePIDFullAl(-999.f, -999.f);
+        } else {
+          aod::pidtof_tiny::binning::packInTable(-999.f, tablePIDAl);
+        }
+        break;
       default:
         LOG(fatal) << "Wrong particle ID in makeTableEmpty() for " << (fullTable ? "full" : "tiny") << " tables";
         break;
@@ -1012,6 +1030,7 @@ struct McPidTof {
     constexpr auto ResponseDe = ResponseImplementation<PID::Deuteron>();
     constexpr auto ResponseTr = ResponseImplementation<PID::Triton>();
     constexpr auto ResponseHe = ResponseImplementation<PID::Helium3>();
+    constexpr auto ResponseAl = ResponseImplementation<PID::Alpha>();
 
     mTOFCalibConfig.processSetup(mRespParamsV3, ccdb, bcs.iteratorAt(0)); // Update the calibration parameters
 
@@ -1092,7 +1111,7 @@ struct McPidTof {
             nSigma = ResponseTr.GetSeparation(mRespParamsV3, trk);
             if (enableMcRecalib && trk.has_mcParticle()) {
               if (std::abs(trk.mcParticle().pdgCode()) == o2::constants::physics::kTriton) { // we rescale only true signal
-                nSigma = applyMcRecalib(IdxPr, trk.pt(), nSigma);                            // FIXME: currently postcalibrations for protons applied to deuterons, to be checked
+                nSigma = applyMcRecalib(IdxPr, trk.pt(), nSigma);                            // FIXME: currently postcalibrations for protons applied to tritons, to be checked
               }
             }
             aod::pidtof_tiny::binning::packInTable(nSigma, tablePIDTr);
@@ -1102,12 +1121,22 @@ struct McPidTof {
             nSigma = ResponseHe.GetSeparation(mRespParamsV3, trk);
             if (enableMcRecalib && trk.has_mcParticle()) {
               if (std::abs(trk.mcParticle().pdgCode()) == o2::constants::physics::kHelium3) { // we rescale only true signal
-                nSigma = applyMcRecalib(IdxPr, trk.pt(), nSigma);                             // FIXME: currently postcalibrations for protons applied to deuterons, to be checked
+                nSigma = applyMcRecalib(IdxPr, trk.pt(), nSigma);                             // FIXME: currently postcalibrations for protons applied to helium3, to be checked
               }
             }
             aod::pidtof_tiny::binning::packInTable(nSigma, tablePIDHe);
             break;
           }
+          case IdxAl: {
+            nSigma = ResponseAl.GetSeparation(mRespParamsV3, trk);
+            if (enableMcRecalib && trk.has_mcParticle()) {
+              if (std::abs(trk.mcParticle().pdgCode()) == o2::constants::physics::kAlpha) { // we rescale only true signal
+                nSigma = applyMcRecalib(IdxPr, trk.pt(), nSigma);                           // FIXME: currently postcalibrations for protons applied to alpha, to be checked
+              }
+            }
+            aod::pidtof_tiny::binning::packInTable(nSigma, tablePIDAl);
+            break;
+          }
           default:
             LOG(fatal) << "Wrong particle ID for standard tables";
             break;
@@ -1168,7 +1197,7 @@ struct McPidTof {
             nSigma = ResponseTr.GetSeparation(mRespParamsV3, trk, resolution);
             if (enableMcRecalib && trk.has_mcParticle()) {
               if (std::abs(trk.mcParticle().pdgCode()) == o2::constants::physics::kTriton) { // we rescale only true signal
-                nSigma = applyMcRecalib(IdxPr, trk.pt(), nSigma);                            // FIXME: currently postcalibrations for protons applied to deuterons, to be checked
+                nSigma = applyMcRecalib(IdxPr, trk.pt(), nSigma);                            // FIXME: currently postcalibrations for protons applied to tritons, to be checked
               }
             }
             tablePIDFullTr(resolution, nSigma);
@@ -1179,12 +1208,23 @@ struct McPidTof {
             nSigma = ResponseHe.GetSeparation(mRespParamsV3, trk, resolution);
             if (enableMcRecalib && trk.has_mcParticle()) {
               if (std::abs(trk.mcParticle().pdgCode()) == o2::constants::physics::kHelium3) { // we rescale only true signal
-                nSigma = applyMcRecalib(IdxPr, trk.pt(), nSigma);                             // FIXME: currently postcalibrations for protons applied to deuterons, to be checked
+                nSigma = applyMcRecalib(IdxPr, trk.pt(), nSigma);                             // FIXME: currently postcalibrations for protons applied to helium3, to be checked
               }
             }
             tablePIDFullHe(resolution, nSigma);
             break;
           }
+          case IdxAl: {
+            resolution = ResponseAl.GetExpectedSigma(mRespParamsV3, trk);
+            nSigma = ResponseAl.GetSeparation(mRespParamsV3, trk, resolution);
+            if (enableMcRecalib && trk.has_mcParticle()) {
+              if (std::abs(trk.mcParticle().pdgCode()) == o2::constants::physics::kAlpha) { // we rescale only true signal
+                nSigma = applyMcRecalib(IdxPr, trk.pt(), nSigma);                           // FIXME: currently postcalibrations for protons applied to alpha, to be checked
+              }
+            }
+            tablePIDFullAl(resolution, nSigma);
+            break;
+          }
           default:
             LOG(fatal) << "Wrong particle ID for full tables";
             break;
diff --git a/PWGHF/TableProducer/pidCreator.cxx b/PWGHF/TableProducer/pidCreator.cxx
index 1b5d7ca73f9..d108a443539 100644
--- a/PWGHF/TableProducer/pidCreator.cxx
+++ b/PWGHF/TableProducer/pidCreator.cxx
@@ -51,6 +51,8 @@ struct HfPidCreator {
   Produces<aod::PidTpcTofTinyTr> trackPidTinyTr;
   Produces<aod::PidTpcTofFullHe> trackPidFullHe;
   Produces<aod::PidTpcTofTinyHe> trackPidTinyHe;
+  Produces<aod::PidTpcTofFullAl> trackPidFullAl;
+  Produces<aod::PidTpcTofTinyAl> trackPidTinyAl;
 
   /// Function to check whether the process function flag matches the need for filling the table
   /// \param initContext  workflow context (argument of the init function)
@@ -89,6 +91,8 @@ struct HfPidCreator {
     checkTableSwitch(initContext, "PidTpcTofTinyTr", doprocessTinyTr);
     checkTableSwitch(initContext, "PidTpcTofFullHe", doprocessFullHe);
     checkTableSwitch(initContext, "PidTpcTofTinyHe", doprocessTinyHe);
+    checkTableSwitch(initContext, "PidTpcTofFullAl", doprocessFullAl);
+    checkTableSwitch(initContext, "PidTpcTofTinyAl", doprocessTinyAl);
   }
 
   void processDummy(aod::Collisions const&) {}
@@ -121,6 +125,7 @@ struct HfPidCreator {
   PROCESS_PID(De)
   PROCESS_PID(Tr)
   PROCESS_PID(He)
+  PROCESS_PID(Al)
 
 #undef PROCESS_PID
 };
diff --git a/PWGHF/Utils/utilsPid.h b/PWGHF/Utils/utilsPid.h
index 8be0035453d..3f7aa28f5cc 100644
--- a/PWGHF/Utils/utilsPid.h
+++ b/PWGHF/Utils/utilsPid.h
@@ -20,7 +20,9 @@
 #include "Common/DataModel/PIDResponseTOF.h"
 #include "Common/DataModel/PIDResponseTPC.h"
 
+#include <CommonConstants/PhysicsConstants.h>
 #include <Framework/Logger.h>
+#include <MathUtils/BetheBlochAleph.h>
 
 #include <cstdint>
 
@@ -33,7 +35,8 @@ enum HfProngSpecies : uint8_t {
   Proton,
   Deuteron,
   Triton,
-  Helium,
+  Helium3,
+  Alpha,
   NHfProngSpecies
 };
 
@@ -45,6 +48,131 @@ enum PidMethod {
   NPidMethods
 };
 
+/// Compute TPC nσ for light nuclei (De/Tr/He3/Al) using a Bethe–Bloch parameter configuration (BB-based PID).
+///
+/// \tparam TrackType   Track/ASoA row type providing TPC accessors.
+/// \param track        Track to be tested.
+/// \param species  Species selector: 3=Deuteron, 4=Triton, 5=Helium3, 6=Alpha.
+/// \param bbParams is Bethe–Bloch  parameters
+/// \return             TPC nσ for the chosen nucleus hypothesis (or -999 if not applicable).
+template <typename TrackType>
+float getNSigmaTpcLightNucleiBetheBloch(const TrackType& track,
+                                        HfProngSpecies species,
+                                        const o2::framework::Configurable<o2::framework::LabeledArray<float>>& bbParams)
+{
+  constexpr float RetValNotApplicable = -999.f;
+
+  if (!track.hasTPC()) {
+    return RetValNotApplicable;
+  }
+  const int row = static_cast<int>(species) - static_cast<int>(HfProngSpecies::Deuteron);
+
+  if (row < 0 || row >= HfProngSpecies::NHfProngSpecies) {
+    return RetValNotApplicable;
+  }
+
+  // Columns: [0..4] BB params, [5] relative resolution (sigma/mean)
+
+  const double relRes = bbParams->get(row, 5u);
+  if (relRes <= 0.) {
+    return RetValNotApplicable;
+  }
+
+  const double bb0 = bbParams->get(row, 0u);
+  const double bb1 = bbParams->get(row, 1u);
+  const double bb2 = bbParams->get(row, 2u);
+  const double bb3 = bbParams->get(row, 3u);
+  const double bb4 = bbParams->get(row, 4u);
+
+  double mass = 0.;
+  switch (species) {
+    case HfProngSpecies::Deuteron:
+      mass = o2::constants::physics::MassDeuteron;
+      break;
+    case HfProngSpecies::Triton:
+      mass = o2::constants::physics::MassTriton;
+      break;
+    case HfProngSpecies::Helium3:
+      mass = o2::constants::physics::MassHelium3;
+      break;
+    case HfProngSpecies::Alpha:
+      mass = o2::constants::physics::MassAlpha;
+      break;
+    default:
+      LOG(fatal) << "Unhandled HfProngSpecies " << static_cast<int>(species);
+  }
+
+  const int charge = (species == HfProngSpecies::Helium3 || species == HfProngSpecies::Alpha) ? 2 : 1;
+  const float rigidity = track.tpcInnerParam(); // p/|q|
+
+  const double x = static_cast<double>(charge) * static_cast<double>(rigidity) / mass;
+  const double expBethe = o2::common::BetheBlochAleph(x, bb0, bb1, bb2, bb3, bb4);
+  const double expSigma = expBethe * relRes;
+
+  if (expSigma <= 0.) {
+    return RetValNotApplicable;
+  }
+
+  return static_cast<float>((track.tpcSignal() - expBethe) / expSigma);
+}
+
+/// \brief Function to fill tables with HF prong PID information
+/// \tparam specPid particle species
+/// \tparam TTrack datatype of the prong track
+/// \tparam TCursor datatype of the cursor of the prong PID table to fill
+/// \param track prong track
+/// \param bbParams is Bethe–Bloch  parameters (only for light nuclei)
+/// \param rowPid cursor of the prong PID table to fill
+template <HfProngSpecies SpecPid, typename TTrack, typename TCursor>
+void fillProngPidLightNuclei(TTrack const& track,
+                             TCursor& rowPid,
+                             const o2::framework::Configurable<o2::framework::LabeledArray<float>>& bbParams)
+{
+
+  // get PID information for the daughter tracks
+  // TODO: add here the code for a possible PID post-calibrations in MC
+  float nSigTpc = -999.f;
+  float nSigTof = -999.f;
+  if constexpr (SpecPid == HfProngSpecies::Deuteron) {
+    // deuteron PID
+    if (track.hasTPC()) {
+      nSigTpc = getNSigmaTpcLightNucleiBetheBloch(track, SpecPid, bbParams);
+    }
+    if (track.hasTOF()) {
+      nSigTof = track.tofNSigmaDe();
+    }
+  } else if constexpr (SpecPid == HfProngSpecies::Triton) {
+    // triton PID
+    if (track.hasTPC()) {
+      nSigTpc = getNSigmaTpcLightNucleiBetheBloch(track, SpecPid, bbParams);
+    }
+    if (track.hasTOF()) {
+      nSigTof = track.tofNSigmaTr();
+    }
+  } else if constexpr (SpecPid == HfProngSpecies::Helium3) {
+    // helium3 PID
+    if (track.hasTPC()) {
+      nSigTpc = getNSigmaTpcLightNucleiBetheBloch(track, SpecPid, bbParams);
+    }
+    if (track.hasTOF()) {
+      nSigTof = track.tofNSigmaHe();
+    }
+  } else if constexpr (SpecPid == HfProngSpecies::Alpha) {
+    // helium4 PID
+    if (track.hasTPC()) {
+      nSigTpc = getNSigmaTpcLightNucleiBetheBloch(track, SpecPid, bbParams);
+    }
+    if (track.hasTOF()) {
+      nSigTof = track.tofNSigmaAl();
+    }
+  } else {
+    LOG(fatal) << "Unsupported PID. Supported species in HF framework:  HfProngSpecies::Deuteron, HfProngSpecies::Triton, HfProngSpecies::Helium3, HfProngSpecies::Alpha";
+  }
+
+  // fill candidate prong PID rows
+  rowPid(nSigTpc, nSigTof);
+}
+
 /// Function to combine TPC and TOF nSigma
 /// \tparam Tiny switch between full and tiny (binned) PID tables
 /// \param nSigmaTpc is the (binned) nSigma separation in TPC (if Tiny = true)
@@ -111,32 +239,8 @@ void fillProngPid(TTrack const& track, TCursor& rowPid)
     if (track.hasTOF()) {
       nSigTof = track.tofNSigmaPr();
     }
-  } else if constexpr (SpecPid == HfProngSpecies::Deuteron) {
-    // deuteron PID
-    if (track.hasTPC()) {
-      nSigTpc = track.tpcNSigmaDe();
-    }
-    if (track.hasTOF()) {
-      nSigTof = track.tofNSigmaDe();
-    }
-  } else if constexpr (SpecPid == HfProngSpecies::Triton) {
-    // triton PID
-    if (track.hasTPC()) {
-      nSigTpc = track.tpcNSigmaTr();
-    }
-    if (track.hasTOF()) {
-      nSigTof = track.tofNSigmaTr();
-    }
-  } else if constexpr (SpecPid == HfProngSpecies::Helium) {
-    // triton PID
-    if (track.hasTPC()) {
-      nSigTpc = track.tpcNSigmaHe();
-    }
-    if (track.hasTOF()) {
-      nSigTof = track.tofNSigmaHe();
-    }
   } else {
-    LOG(fatal) << "Unsupported PID. Supported species in HF framework: HfProngSpecies::Pion, HfProngSpecies::Kaon, HfProngSpecies::Proton, HfProngSpecies::Deuteron, HfProngSpecies::Triton, HfProngSpecies::Helium";
+    LOG(fatal) << "Unsupported PID. Supported species in HF framework: HfProngSpecies::Pion, HfProngSpecies::Kaon, HfProngSpecies::Proton";
   }
 
   // fill candidate prong PID rows