You look to make same confusion as your previous post about Fourier transformation. Let $S_N(f)$ the partial Fourier series. The convergence of $S_N(f)$ (if it converge to something, it converge to $f$) needs more condition that only $L^2$. In my memories, there is a something with Holder continuity. But, by a theorem, $S_n(f)$ will converge to $f$ in the $L^2-$sense. This mean that $$\lim_{N\to \infty }\int_0^1|S_N(f)-f|^2=0,$$
which is not the same. In particular, using Minkovski inequality, you will indeed have that the limit of $S_N(f)\in L^2$, since Minkovski tels you that $$\int_0^1|S_N(f)|^2\leq \int_0^1|S_N(f)-f|^2+\int_0^1|f|^2,$$
(notice that it's not an obvious result). So, all we can say it that the limit of $S_N(f)$ (in the $L^2$ sense) will be $L^2$ (but this is nothing surprising since it's $f$ which is in $L^2$). But you can't say more a priori.
Here there is an interesting trap
If you can find the subtlety, I think you will understand more those $L^2$ space. We know that $L^2$ are Hilbert space. An orthonormal basis of $L^2(0,1)$ is given by $\{e^{inx}\}_{n\in\mathbb Z}$. Then, $$f(x)=\sum_{k=-\infty }^\infty \left<f,e^{inx}\right>e^{inx}=\sum_{k=-\infty }^\infty \int_0^1 f(x)e^{-iny}\mathrm d ye^{inx}=\sum_{n=-\infty }^\infty c_ne^{-inx}$$
$$=\lim_{N\to \infty }\sum_{n=-N}^N c_ne^{inx}=\lim_{N\to \infty }S_N(x),$$
where $c_n=\int_0^1 f(x)e^{-inx}\mathrm d x$. As we can see, if $S_N$ converge in the $L^2-$sense to $f$, then it converge in the usual sense to $f$.
Question : What's wrong here ?