  The multiplicity of heavy quarks from gluon splitting in e+e- annihilation
has now been theoretically calculated and experimentally measured at LEP.
However, the experimental measurement requires theoretical input for the shape
of the multiplicity with respect to an event shape. In this paper we calculate
the multiplicity of heavy quarks from gluon splitting in e+e- annihilation as a
function of the heavy jet mass of the event, to next-to-leading logarithmic
accuracy. We compare our result with Monte Carlo predictions.
