The Institute for Radiological protection and Nuclear Safety was engaged in the EC funded EVIDOS project to provide reference spectrometry data using its Bonner sphere system. The data were processed by means of two unfolding codes, NUBAY and GRAVEL, both provided by the Physikalisch-Technische Bundesanstalt. The NUBAY program, based on Bayesian parameter estimation methods, assumes a parameterised spectrum and provides posterior probability distributions for the parameters. The code GRAVEL, an iterative algorithm based on SAND-II, was used with various default spectra, among them the NUBAY solution. The BS measurements were used to establish the neutron fluence energy distributions and reference values for the neutron ambient dose equivalent. As this quantity depends strongly on the high energy neutrons, a sensitivity analysis was done by unfolding the BS data with GRAVEL using the NUBAY solution spectrum as default with various changes in the parameters of the high energy peak. This new method of analysing Bonner sphere data allowed the determination of reliable neutron spectra, as well as a very good estimate of the corresponding integral quantities with small associated uncertainties.