Four NE213 detectors of different size have been calibrated at the accelerator facility of the PTB. The response functions were experimentally determined for 3 3 neutron energies between 1 MeV and 16 MeV and compared with Monte Carlo simulations using the NRESP7 code. The light output functions for recoil protons were found to be significantly different for all detectors even if they were of the same size. The neutron fluence determined on the basis of the response functions calculated with the corresponding light output functions agreed to better than ± 2% with reference values if energy independent adjustment factors between 0.98 and 1.03 were applied. The response matrices required for the unfolding of neutron induced pulse height spectra were therefore calculated with the NRESP7 code taking into account the adjustment factors. Similarly, the response matrices for photons were calculated with the EGS4 code, but without any adjustment. Finally, the DIFBAS code was applied for the unfolding of pure neutron- and photon-induced pulse height spectra. The resulting spectral fluences are in reasonable agreement with the results obtained by time-of-flight measurements and by spectrometry with a Ge detector.