I am reading Otto Forster's book "Lecture on Riemann surfaces" and on pages 109-110, he introduces the space $L^2(D,\mathcal{O})$ of holomorphic square-integrable functions $f:D\to \mathbb{C}$ (where $D\subset\mathbb{C}$ is open). In particular for $D=B(a,r)$ he explains that the monomials $\psi_n(z)=(z-a)^n$ form an orthogonal system with $\|\psi_n\|=\sqrt\frac{\pi}{n+1}r^{n+1}$ and that if $$f(z)=\sum_{n=0}^\infty c_n(z-a)^n\tag{1}\label{taylor}$$ is in $L^2(B(a,r),\mathcal{O})$, then by Pythagoras $$\|f\|^2_{L^2}=\sum_{n=0}^\infty \frac{\pi r^{2n+2}}{n+1}|c_n|^2.\label{pythagoras}\tag{2}$$
My only problem with all this is that it seems to me that we can only apply Parseval if the Taylor series \eqref{taylor} converges for the norm $L^2$, and it doesn't seem obvious that it does.
I know that the Taylor series converges pointwise on $B(a,r)$ and uniformly on every compact subset. I tried to apply the dominated convergence theorem to show that \eqref{taylor} also converges for the $L^2$ norm but I can't get a good integrable function to bound the differences $$\left|f-\sum_{n= 0}^Nc_n(z-a)^n\right|^2=\left|\sum_{n= N+1}^\infty c_n(z-a)^n\right|^2.$$
I also tried to use an approach similar to what is discussed in the comments of this answer, but I got stuck because I don't know that $L^2(D,\mathcal{O})$ is a Hilbert space (the proof in the book relies on \eqref{pythagoras}).
Suppose we're in the unit disc $\mathbb D$ for simplicity. Let $\sum_{n=0}^{\infty}a_nz^n$ be the Taylor series of $f$ in $\mathbb D.$ Using the orthogonality of the exponenetials, we see $$\int_{\mathbb D}|f|^2\, dA = \int_0^1 \int_0^{2\pi} |f(re^{it})|^2\, dt \, r\, dr = \int_0^1\int_0^{2\pi}|\sum_{n=0}^{\infty}a_nr^ne^{int}|^2\, dt\, r\, dr$$ $$ = \int_0^1 (2\pi \sum_{n=0}^{\infty}|a_n|^2r^{2n})\, r\, dr = 2\pi \sum_{n=0}^{\infty}|a_n|^2/(2n+2).$$ Now $|a_n|^2/(2n+2) = (\|a_nz^n\|_{L^2})^2.$ So whether $f\in L^2$ or not, the above always holds.