Write particle positions using functionObjects

I've searched extensively within the source code and online and looks like there isn't any method available to write the position of Lagrangian particles, perhaps appending some Lagrangian data to each position, in the same way as cutplanes and isosurfaces can be written using function objects, i.e. having a file written with a given frequency in the corresponding time folder within the postProcessing folder in vtk or vtp format.

It would be great to have a confirmation for this and, if so, to consider the development of such function object and spare the user from writing the whole solution in order to reconstruct and view the dynamics of particles with a decent frequency, especially for large meshes.