We study the heat semigroup maximal operator associated with a well-known orthonormal system in the
$d$-dimensional ball. The corresponding heat kernel is shown to satisfy Gaussian bounds. As a consequence, we can prove weighted
$L^{p}$ estimates, as well as some weighted inequalities in mixed norm spaces, for this maximal operator.