I am probably making some very annoying mistake but for whatever it is worth let me ask it here.

Consider a hypersurface $S$ in $\mathbb{R}^{n+1}$. Let $\nu$ be the unit outward normal. Let $f:\mathbb{R}^{n+1}\rightarrow\mathbb{R}$ be defined as $f(p)=r^2$, where $r$ is the radial distance of $p$ to the origin. Let $z=f|_S$, $\sigma$ be the second fundamental tensor, and for an orthonormal local framefield $\{E_i\}$ let $h=\frac{1}{n}\sum\sigma(E_i,E_i)$. Then we have the following -

$\Delta z = 2n(1+\langle h,\vec{r}\rangle)$, This comes from the (possibly preprint) paper of Oscar Garay titled "application of reilly's formula"

Where $\Delta$ is defined as $\Delta g=-\textrm{Trace}(X\mapsto\nabla_X\textrm{grad}g)$

My problem is I am getting the wrong sign always ... I get $\Delta z = 2n(-1+\langle h,\vec{r}\rangle)$