I think r is a function of and the solution of:
That's the domain in the x-y plane where the integration is being done (just complete the square when you equate both equations.
Doing that I get (in cylindrical coordinates):
and is one of the roots to the solution of that equation above. I don't know. Looks pretty messy to me. Maybe another easier way though.