From Rademacher's book (Topics in Analytic Number Theory) I'm using the functional equation of $\vartheta_2(0|\tau) = 2\sum_{m=0}^\infty q^{\left(m+\frac12\right)^2} = \vartheta_2(\tau)$ and the fact that it vanishes at the cusps $\tau = \frac{a}{b}$, $a$ odd, $b$ even.
The order of vanishing at infinity is just $\frac14$. I'd like to find the order of vanishing at the cusp $\tau = \frac12$ (no particular reason for this choice). Using say $A = \left(\begin{array}{lr} 1 & 0 \\ 2 & 1 \end{array}\right) \in \Gamma_0(2)$ and plugging into the functional equation I get $$\vartheta_2\left(\frac{\tau}{2\tau+1}\right) = i e^{-\pi i/4} \sqrt{-i(2\tau+1)}\vartheta_2(\tau)$$ As I understand it I'd want to plug $\tau = i\infty$ into this, which would give me $\vartheta_2\left(\frac12\right)$ in terms of $\vartheta_2(i\infty) = q^\frac{1}{4}+...$ However, I don't know how to handle the $\sqrt{2\tau+1}$ term if I was to put $\tau = i\infty$, so I feel like there's a big hole in my understanding.
The equation from Rademacher is given as $$\vartheta_2\left(\frac{\upsilon}{c\tau+d}\bigg|\frac{a\tau+b}{c\tau+d}\right) = i^{3-d-(a+1)/2}e^{\pi i (dc+ab-ac-1)/4}\sqrt{\frac{c\tau+d}{i}}e^{\pi i c\upsilon^2/(c\tau+d)}\vartheta_{1-c, 1-d}(\upsilon|\tau)$$ on page 182, for $c>0$, $a$ odd.