While playing around with random matrices and I arrived at a different formula for the mean of the limiting normal distribution for a spectral CLT for sample covariance matrices. More precisely I have the formula
\begin{align*}
& \sum\limits_{b=1}^{r-1} c^{b} \left(A(r,b) - {r \choose b}^2 \, b\right)
\end{align*}
for the expression (1.23) from Bai and Silverstein's famous paper 'CLT for Linear Spectral Statistics of Large-Dimensional Sample Covariance Matrices', where
\begin{align*}
& A(r,b) := \sum\limits_{m=1}^{b \land (r-b+1)} (2m-1) {r \choose b-m}{r \choose b+m-1} + \sum\limits_{m=1}^{b \land (r-b)} (2m+1) {r \choose b-m}{r \choose b+m} \ .
\end{align*}
In order to show that my expression is equal to that of Bai and Silverstein, I need to show the equality
$$
A(r,b) = \frac{1}{2} {2r \choose 2b} + \left( b - \frac{1}{2} \right) {r \choose b}^2
$$
for all $r \in \mathbb{N}$ and $b<r$. I am not too well versed with combinatorial identities and so far my efforts have lead nowhere. Does someone have an idea how to prove this?

Testing with the computer shows this to be true for at least all $b<r \leq 1000$.

Technically I could just use the fact that my formula and Bai and Silverstein's formula describe the same object to show this identity, which seems rather crude unless this identity is not jet known.