The purpose of the question is to check which eigenfunctions of ##L^2## and ##L_z## may appear when the original wavefunction is expanded as a sum of those eigenfunctions, which are the spherical harmonics ##Y_{lm}##. Now if you remember the characteristic dependency of ##Y_{lm}## w.r.t to the three coordinates x,y, and z is governed by the indices ##l## and ##m##. Using this fact you can get the first rough idea of which ##l## and ##m## are allowed to appear because they conform with the dependency of the original wavefunction to three coordinates. First of all, do you know how the spherical harmonics looks like and its dependency on x,y,z?
For example, ##Y_{2,-2} \propto (x-iy)^2##, by comparing this with the manner x,y, and z appear in the original wavefunction, do you think ##Y_{2,-2}## will appear in the expansion?