Several Collisional-Radiative (CR) models have been developed to calculate the attenuation and the population of excited states of hydrogen or deuterium beams injected into tokamak plasmas. The datasets generated by these CR models are needed for the modelling of beam ion deposition and (excited) beam densities in current experiments, and the reliability of this data will be crucial to obtain helium ash densities on ITER combining charge exchange and beam emission spectroscopy. Good agreement between the dierent CR models for the Neutral Beam (NB) is found, if corrections to the fundamental cross sections are taken into account. First the Ha and Hb beam emission spectra from JET are compared with the expected intensities. Second, the line ratios within the Stark multiplet are compared with the predictions of a sublevel resolved model. The measured intensity of the full multiplet is 30% lower than expected on the basis of beam attenuation codes and the updated beam emission rates, but apart from the atomic data this could also be due to the characterization of the NB path and line of sight integration and the absolute calibration of the optics. The modelled n = 3 to n = 4 population agrees very well with the ratio of the measured H to H beam emission intensities. Good agreement is found as well between the neutral beam power fractions measured with beam emission in plasma and on the JET Neutral Beam Test Bed. The Stark line ratios and s/p intensity ratio deviate from a statistical distribution, in agreement with the CR model in parabolic states from Marchuk et al.