OK, so if I try to map just the square roots to the range (-inf, 0) I get:
1) z^2-1=-R where R ranges from (0, inf).
Solving for Z yields z = \pm \sqrt{1-R}, which for R<1 gives (-1,1), and for R>1 gives (-i inf, +i inf) so I get branch cuts from (-1,1) and the imaginary axis.
2) solving...