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.