diff --git a/PWGHF/D2H/DataModel/ReducedDataModel.h b/PWGHF/D2H/DataModel/ReducedDataModel.h index 999be60db1b..377f1a87bb3 100644 --- a/PWGHF/D2H/DataModel/ReducedDataModel.h +++ b/PWGHF/D2H/DataModel/ReducedDataModel.h @@ -749,8 +749,8 @@ namespace hf_cand_b0_reduced { DECLARE_SOA_INDEX_COLUMN_FULL(Prong0, prong0, int, HfRed3Prongs, "_0"); //! Prong0 index DECLARE_SOA_INDEX_COLUMN_FULL(Prong1, prong1, int, HfRedTrackBases, "_1"); //! Prong1 index -DECLARE_SOA_INDEX_COLUMN_FULL(ProngD0, prongD0, int, HfRed2Prongs, "_0"); //! ProngD0 index -DECLARE_SOA_INDEX_COLUMN_FULL(ProngBachPi, prongBachPi, int, HfRedTrackBases, "_1"); //! ProngBachPi index +DECLARE_SOA_INDEX_COLUMN_FULL(ProngD0, prongD0, int, HfRed2Prongs, "_0"); //! ProngD0 index +DECLARE_SOA_INDEX_COLUMN_FULL(ProngBachPi, prongBachPi, int, HfRedTrackBases, "_1"); //! ProngBachPi index DECLARE_SOA_INDEX_COLUMN_FULL(ProngSoftPi, prongSoftPi, int, HfRedSoftPiBases, "_2"); //! ProngSoftPi index DECLARE_SOA_COLUMN(Prong0MlScoreBkg, prong0MlScoreBkg, float); //! Bkg ML score of the D daughter DECLARE_SOA_COLUMN(Prong0MlScorePrompt, prong0MlScorePrompt, float); //! Prompt ML score of the D daughter diff --git a/PWGHF/D2H/TableProducer/candidateSelectorB0ToDPiReduced.cxx b/PWGHF/D2H/TableProducer/candidateSelectorB0ToDPiReduced.cxx index 4eede1fd6c5..51ba23c8511 100644 --- a/PWGHF/D2H/TableProducer/candidateSelectorB0ToDPiReduced.cxx +++ b/PWGHF/D2H/TableProducer/candidateSelectorB0ToDPiReduced.cxx @@ -156,12 +156,13 @@ struct HfCandidateSelectorB0ToDPiReduced { } } - /// Utility function to retrieve the bach pion track - ///from the B0 candidate in the D*-pi decay channel + /// Utility function to retrieve the bach pion track + /// from the B0 candidate in the D*-pi decay channel /// \param candidate is the B0 candidate /// \return bach pion track template - auto getTrackBachPi(const T1& candidate) { + auto getTrackBachPi(const T1& candidate) + { return candidate.template prongBachPi_as(); } @@ -170,21 +171,23 @@ struct HfCandidateSelectorB0ToDPiReduced { /// \param prongBachPi is the candidate's bachelor pion prong /// \note this method is used for B0 → D*- π+ candidates with D meson ML scores template - auto getMlInputFeatures(const T1& candB0, const T2& prongBachPi) { + auto getMlInputFeatures(const T1& candB0, const T2& prongBachPi) + { auto prongSoftPi = candB0.template prongSoftPi_as(); if constexpr (withDmesMl) { - return hfMlResponse.getInputFeaturesDStarPi(candB0, prongBachPi, prongSoftPi); + return hfMlResponse.getInputFeaturesDStarPi(candB0, prongBachPi, prongSoftPi); } else { - return hfMlResponse.getInputFeaturesDStarPi(candB0, prongBachPi, prongSoftPi); + return hfMlResponse.getInputFeaturesDStarPi(candB0, prongBachPi, prongSoftPi); } } - /// Utility function to retrieve the bach pion track - ///from the B0 candidate in the D-pi decay channel + /// Utility function to retrieve the bach pion track + /// from the B0 candidate in the D-pi decay channel /// \param candidate is the B0 candidate /// \return bach pion track template - auto getTrackBachPi(const T1& candidate) { + auto getTrackBachPi(const T1& candidate) + { return candidate.template prong1_as(); } @@ -193,7 +196,8 @@ struct HfCandidateSelectorB0ToDPiReduced { /// \param prongBachPi is the candidate's bachelor pion prong /// \note this method is used for B0 → D- π+ candidates with D meson ML scores template - auto getMlInputFeatures(const T1& candB0, const T2& prongBachPi) { + auto getMlInputFeatures(const T1& candB0, const T2& prongBachPi) + { if constexpr (withDmesMl) { return hfMlResponse.getInputFeatures(candB0, prongBachPi); } else { diff --git a/PWGHF/D2H/Tasks/taskB0Reduced.cxx b/PWGHF/D2H/Tasks/taskB0Reduced.cxx index 91cb9ef6a0e..169e4fcd834 100644 --- a/PWGHF/D2H/Tasks/taskB0Reduced.cxx +++ b/PWGHF/D2H/Tasks/taskB0Reduced.cxx @@ -235,8 +235,8 @@ struct HfTaskB0Reduced { doprocessMcDplusPi, doprocessMcDplusPiWithDecayTypeCheck, doprocessMcDplusPiWithDmesMl, doprocessMcDplusPiWithDmesMlAndDecayTypeCheck, doprocessMcDplusPiWithB0Ml, doprocessMcDplusPiWithB0MlAndDecayTypeCheck}; - const AxisSpec axisMass = ((std::accumulate(processFuncDplusPi.begin(), processFuncDplusPi.end(), 0)) > 0) ? axisMassDminus : axisMassDeltaMassDStar; - std::string dMesSpecie{""}; + const AxisSpec axisMass = ((std::accumulate(processFuncDplusPi.begin(), processFuncDplusPi.end(), 0)) > 0) ? axisMassDminus : axisMassDeltaMassDStar; + std::string dMesSpecie{""}; if ((std::accumulate(processFuncDplusPi.begin(), processFuncDplusPi.end(), 0)) > 0) { dMesSpecie += "D^{#minus}"; } else {