HEALpy - Getting spherical harmonics from as function of pixel, l and m

Question

I am in need to calculate an integral of spherical harmonics and their products of different orders (l,m) over a fraction of a sphere (let\'s say for simplicity over a half a sp