Invariant-mass distribution of $\mu^+ \mu^-$ pairs in Pb-Pb collisions at $\sNN$ = 2.76 TeV with the $\Upsilon(1S)$, $\Upsilon(2S)$, and $\Upsilon(3S)$ peaks fitted by the sum of three extended Crystal Ball functions with identical relative widths and identical relative displacements from the PDG mass values. The tail shape is fixed by the embedding-MC simulation and the combinatorial background is parametrized by an exponential.