I would try polar coordinates. Your assumptions are a little strange: you integrate on z from -Infinity to +Infinity, but you assume z>=0.
z
-Infinity
+Infinity
