We present our analytic results for the NLO corrections to the partial decay width t→H + b followed by b→BX for nonzero b-quark mass (m b ≠0) in the Fixed-Flavor-Number scheme (FFNs) taking the H + boson to be stable. To make our predictions for the energy distribution of the outgoing bottom-flavored hadron (B-hadron) as a function of the normalized B-energy fraction x B , we apply the General-Mass Variable-Flavor-Number scheme (GM-VFNs) in the general two Higgs doublet model (2HDM). In order to describe both the b-quark and the gluon hadronizations in top decay we use fragmentation functions extracted from data from e + e − machines. We find that the most reliable prediction for the B-hadron energy spectrum is made in the GM-VFN scheme, especially when the second model in the general 2HDM is considered.