Additional symmetries in a theta-like function

Question

Let $\Theta : \mathfrak{h}\times \mathfrak{h} \to \mathbb{R}$ be defined as follows $$ \Theta(z,\tau) = \sum_{\omega_1, \omega_2 \ \in \ \mathbb{Z}\tau + \mathbb{Z}} \exp\Big(-2\pi\frac{| \omega_1 z + \omega_2|^2 + | \omega_1\bar{z} + \omega_2|^2}{\Im(z)\Im(\tau)}\Big) $$
This function is clearly invariant under the action of $\operatorname{SL}_2(\mathbb{Z})\times\operatorname{SL}_2(\mathbb{Z})$ acting on $(z,\tau)$ by linear fractional transformations. However it also satisfies the symmetry $$ \Theta(z,\tau) = \Theta(\tau,z)$$ This last symmetry can be proven by noticing that $$ |(m_1\tau+n_1)z + (m_2\tau+n_2)|^2 = |(m_1z+m_2)\tau + (n_1z+n_2)|^2.$$ I am curious to know whether there is are any more symmetries satified by this function. A way to determine this might be to find a more conceptual proof of this final identity - perhaps a proof which treats all symmetries on an equal footing, and places the function above within a larger framework. Any ideas for a more conceptual proof which would determine the full group of symmetries would be most welcome.