The following is a re-post from MSE because I did not get any answer even after offering a bounty.

Towards the end of G. N. Watson's (one of the joint authors of famous book "A Course of Modern Analysis") paper "The Final Problem: An Account of the Mock Theta Functions" the following formula of Ramanujan is mentioned: \begin{align}&\int_{0}^{\infty}e^{-3\pi x^{2}}\frac{\sinh \pi x}{\sinh 3\pi x}\,dx\notag\\ &\,\,\,\,\,\,\,\,= \frac{1}{e^{2\pi/3}\sqrt{3}}\sum_{n = 0}^{\infty}\frac{e^{-2n(n + 1)\pi}}{(1 + e^{-\pi})^{2}(1 + e^{-3\pi})^{2}\dots(1 + e^{-(2n + 1)\pi})^{2}}\tag{1} \end{align} where the term corresponding to $n = 0$ in the sum on the right is $1$.

Is there way to establish this exotic integral formula? Or a reference to any existing proof of $(1)$ would be of great help.