Square Function Estimates and Local Smoothing for Fourier Integral Operators