Spectral functions that are used in neutrino event, generators to model quasielastic (QE) scattering from nuclear targets include Fermi gas, Local Thomas Fermi gas (LTF), Bodek-Ritchie Fermi gas with high momentum tail, and the Benhar-Fantoni two dimensional spectral function. We find that the ν dependence of predictions of these spectral functions for the QE differential cross sections (d 2 σ /d Q 2 dν) are in disagreement with the prediction of the ψ superscaling function which is extracted from fits to quasielastic electron scattering data on nuclear targets. It is known that spectral functions do not fully describe quasielastic scattering because they only model the initial state. Final state interactions distort the shape of the differential cross section at the peak and increase the cross section at the tails of the distribution. We show that the kinematic distributions predicted by the ψ superscaling formalism can be well described with a modified effective spectral function (ESF). By construction, models using ESF in combination with the transverse enhancement contribution correctly predict electron QE scattering data.