I'm not sure I would cite this paper. You are implementing a 3D-SHORE (extension of ISMRM 2009), using spherical harmonics and Gaussian-Laguerre, as proposed in Merlet et al.