In preparation of future multiobject spectrographs (MOS) whose one of the major role is to provide an extensive statistical studies of high redshifted galaxies surveyed, the demonstrator CANARY has been designed to tackle technical challenges related to open-loop adaptive optics (AO) control with jointed Natural Guide Star and Laser Guide Star tomography. We have developed a point spread function (PSF) reconstruction algorithm dedicated to multiobject adaptive optics systems using system telemetry to estimate the PSF potentially anywhere in the observed field, a prerequisite to postprocess AO-corrected observations in integral field spectroscopy. We show how to handle off-axis data to estimate the PSF using atmospheric tomography and compare it to a classical approach that uses on-axis residual phase from a truth sensor observing a natural bright source. We have reconstructed over 450 on-sky CANARY PSFs and we get standard-deviation (std) of 1.3/4.8 on the H-band Strehl ratio (SR) with 92.3% of correlation between reconstructed and sky SR. On the full-width at half-maximum, we get, respectively, 2.94 mas, 19.9 mas, and 88.3% for the bias, std, and correlation. The reference method achieves 0.4/3.5/95% on the SR and 2.71 mas/14.9 mas/92.5% on the FWHM for the bias/std/correlation.