An equilibrium superfluid is usually modelled by a Bose gas below its transition temperature. Here, we compute rigorously the extremal equilibrium states below this transition temperature. We compute the deviation from normality of the angular momentum fluctuations both above and below this transition temperature. Finally, we compute the distribution functions of the angular momentum and its fluctuations.