I've already see a proof by Marko Riedel which I list it follows:
The standard way to treat these sums is to integrate $$ f(z) = \frac{1}{(z+\alpha)^2} \pi \cot(\pi z)$$ along a contour consisting of a circle of radius $R$ and with $R$ going to infinity and hence being larger than $\alpha$, where the circle does not pass through the poles on the real axis. Now along the semicircle in the upper half plane we have $$|f(z)| \le \frac{1}{(R-|\alpha|)^2}\pi \left|\frac{e^{i\pi R\exp(i\theta)} + e^{-i\pi R\exp(i\theta)}} {e^{i\pi R\exp(i\theta)} - e^{-i\pi R\exp(i\theta)}}\right|= \frac{1}{(R-|\alpha|)^2} \pi \left|\frac{e^{2i\pi R\exp(i\theta)}+1}{e^{2i\pi R\exp(i\theta)}-1}\right| < \frac{1}{(R-|\alpha|)^2} \pi \frac{1+e^{-2\pi R\sin(\theta)}e^{2i\pi R\cos(\theta)}} {1-e^{-2\pi R\sin(\theta)}e^{2i\pi R\cos(\theta)}}$$ This last term is clearly $O(1/R^2)$ as $R$ goes to infinity as the quotient of the two exponentials goes to one since $\exp(-R)$ vanishes and there is no singularity when $\theta = 0$ or $\theta = \pi$ as $R\cos\theta = \pm R$, which is not an integer by the assumption that the circle avoids the poles and hence cannot be one.
My question:
It is clear that for each $z=Re^{i \theta}$ , $f(z)=O(1/R^2)$ ,but can this guarantee that the last term on RHS is uniform bounded ? Since $\theta$ varies on a compact interval $[0, \pi]$ , I want to show that for each $ \theta \in [0, \pi]$ , and a fixed positive number $M \gt 1 $ there exist an open ball contains $\theta$ and a fixed $R_{\theta}\gt 0$ , for every $R \ge R_{\theta}$ and $x \in $ the open ball $$\frac{1+e^{-2\pi R\sin(x)}e^{2i\pi R\cos(x)}}
{1-e^{-2\pi R\sin(x)}e^{2i\pi R\cos(x)}} \le M$$ For $\theta \neq 0,\pi$ , it is easy to find the desired $R_{\theta}$ and the open ball , but if $\theta = 0$ or $ \theta = \pi$ , for every open ball containing $\theta$ , it behave erratically near $\theta$ , I have no idea how to deal with this .