Using Rosenblatt transformation, How can I generate multivariate Cauchy distribution with the density:

f(x)=\frac{\Gamma (b)}{\pi(1+x^Tx)^b}

x=(x_1,...,x_d)^T\in \mathbb{R} ^d

where b=(d+1)/2