Neutron peripheral contamination in patients undergoing high-energy photon radiotherapy is considered as a risk factor for secondary cancer induction. Organ-specific neutron-equivalent dose estimation is therefore essential for a reasonable assessment of these associated risks. This work aimed to develop a method to estimate neutron-equivalent doses in multiple organs of radiotherapy patients. The method involved the convolution, at 16 reference points in an anthropomorphic phantom, of the normalized Monte Carlo neutron fluence energy spectra with the kerma and energy-dependent radiation weighting factor. This was then scaled with the total neutron fluence measured with passive detectors, at the same reference points, in order to obtain the equivalent doses in organs. The latter were correlated with the readings of a neutron digital detector located inside the treatment room during phantom irradiation. This digital detector, designed and developed by our group, integrates the thermal neutron fluence. The correlation model, applied to the digital detector readings during patient irradiation, enables the online estimation of neutron-equivalent doses in organs. The model takes into account the specific irradiation site, the field parameters (energy, field size, angle incidence, etc) and the installation (linac and bunker geometry). This method, which is suitable for routine clinical use, will help to systematically generate the dosimetric data essential for the improvement of current risk-estimation models.