We present a study into the capabilities of integrated and spatially resolved integral field spectroscopy of galaxies at z = 2-4 with the future HARMONI spectrograph for the European Extremely Large Telescope (E-ELT) using the simulation pipeline, HSIM. We focus particularly on the instrument's capabilities in stellar absorption line integral field spectroscopy, which will allow us to study the stellar kinematics and stellar population characteristics. Such measurements for star-forming and passive galaxies around the peak star formation era will provide a critical insight into the star formation, quenching and mass assembly history of high-z, and thus present-day galaxies. First, we perform a signal-to-noise study for passive galaxies at a range of stellar masses for z = 2-4, assuming different light profiles; for this population, we estimate that integrated stellar absorption line spectroscopy with HARMONI will be limited to galaxies with M-* greater than or similar to 10(10.7) M-circle dot. Secondly, we use HSIM to perform a mock observation of a typical star-forming 10(10) M-circle dot galaxy at z = 3 generated from the high-resolution cosmological simulation NUTFB. We demonstrate that the input stellar kinematics of the simulated galaxy can be accurately recovered from the integrated spectrum in a 15-h observation, using common analysis tools. Whilst spatially resolved spectroscopy is likely to remain out of reach for this particular galaxy, we estimate HARMONI's performance limits in this regime from our findings. This study demonstrates how instrument simulators such as HSIM can be used to quantify instrument performance and study observational biases on kinematics retrieval; and shows the potential of making observational predictions from cosmological simulation output data.