In arenas of application including environmental science, economics, and medicine, it is increasingly common to consider time series of curves or functions. Many inferential procedures employed in the analysis of such data involve the long run covariance function or operator, which is analogous to the long run covariance matrix familiar to finite dimensional time series analysis and econometrics. This function may be naturally estimated using a smoothed periodogram type estimator evaluated at frequency zero that relies crucially on the choice of a bandwidth parameter. Motivated by a number of prior contributions in the finite dimensional setting, we propose a bandwidth selection method that aims to minimize the estimator's asymptotic mean squared normed error (AMSNE) in $L^2[0,1]^2$. As the AMSNE depends on unknown population quantities including the long run covariance function itself, estimates for these are plugged in in an initial step after which the estimated AMSNE can be minimized to produce an empirical optimal bandwidth. We show that the bandwidth produced in this way is asymptotically consistent with the AMSNE optimal bandwidth, with quantifiable rates, under mild stationarity and moment conditions. These results and the efficacy of the proposed methodology are evaluated by means of a comprehensive simulation study, from which we can offer practical advice on how to select the bandwidth parameter in this setting.