This paper introduces the exact, analytic integration of all final state variables for the process of nonlinear Compton scattering in an intense plane wave laser pulse, improving upon a previously slow and challenging numerical approach. Computationally simple and insightful formulae are derived for the total scattering probability and mean energy-momentum of the emitted radiation. The general form of the effective mass appears explicitely. We consider several limiting cases, and present a quantum correction to Larmor's formula. Numerical results are plotted and analysed in detail.a. Introduction: The recent and predicted progress in laser technology leading to very high peak intensities justify the need for a better understanding of so-called nonlinear QED, describing phenomena occurring in fields so strong that their effects cannot be treated perturbatively.Unfortunately, the complexity of the processes inside these ultra-intense laser beams has meant that several simplifications have had to be used to make practical computations feasible. The laser beam is usually supposed to be in a coherent state, which can be well approximated by a classical field. Due to the relatively small frequencies, unless massive particles of very high Lorentz factor are involved, quantum effects are small, so a fully classical description may often be justified. For instance, in discussing the scattering of an electron in a laser beam, one may consider Thompson scattering instead of its quantum counterpart, nonlinear Compton scattering (NLCS) [1]. The classical approximation allows for a realistic description of the laser field and the inclusion of radiation reaction (RR) [2], but is unable to describe important quantum effects, such as nonperturbative pair creation from vacuum [3], the trident process [4], or vacuum birefringence [5].A treatment of these processes in the framework of nonlinear QED, even in a semiclassical approach, has not yet been performed without further approximations, such as replacing the laser field by an idealized plane wave, thus allowing for analytical (Volkov) solutions to the Dirac equation. This disregards the strong spatial focusing of the beam needed to attain high intensities. In addition, for a long time, results were restricted to infinite, monochromatic plane waves, or even, giving up periodicity, to a crossed field model [6,7]. Only recently, the more realistic short pulse plane waves came into use [8][9][10].In [8,11,12], the photonic and electronic distributions resulting from NLCS were described in detail for some model pulses. In principle, by integrating these distributions, the total probability and expectation values, such as for the emitted radiation's energy-momentum, can be obtained. However, previous papers did never plot these quantities, because of the great numerical challenge posed by this task. By a change in integration order and a different regularization, allowing for all final state integrals to be performed analytically, we obtain formulae that are not only very easy to com...