Proof of the transformation formula of the Dedekind eta function using the Euler-Maclaurin formula

Question

Using the Euler-Maclaurin summation formula, we can prove that$$\log η(τ)=\frac{iπτ}{12}-\frac{iπ}{12τ}-\log\sqrt{-iτ}-\int_0^∞p(y)\left[\frac{2πiτ}{e^{-2πiτy}-1}+\frac1y\right]\,\mathrm dy,$$where $\eta(\tau)$ is the Dedekind eta function, and $p(y)=\left\{y\right\}-\frac{1}{2}$ ($\{y\}$ being the fractional part of $y$). My question: How do we use this formula to deduce the transformation formula: $$\eta\left(\frac{a\tau+b}{c\tau+d}\right)=\epsilon(a,b,c,d)[-i(c\tau+d)]^{\frac{1}{2}}\eta(\tau),$$ where the integers $a,b,c$ and $d$ satisfy $c>0$ and $ad-bc=1$, and $\epsilon(a,b,c,d)$ is given by: $$\epsilon(a,b,c,d)=\exp \left[i\pi\left(\frac{a+b}{12c}-s(d,c)-\frac{1}{4}\right)\right],$$ where $s(d,c)$ is the Dedekind sum, and is given by: $$s(d,c)=\sum_{n=1}^{c-1}\frac{n}{c}\left(\frac{dn}{c}-\left \lfloor \frac{dn}{c} \right \rfloor -\frac{1}{2}\right).$$ My gut feeling tells me that the Mittag-Leffler expansion: $$\frac{1}{e^{z}-1}=\frac{1}{z}-\frac{1}{2}+2z\sum_{n=1}^{\infty}\frac{1}{z^{2}+4\pi^{2}n^{2}}$$ is key here, but I have no idea how to proceed!