m_2d = hp.norm() ** hp.uniform()

data = m_2d.rvs(1_000)
plt.scatter(data[:,0], data[:,1], c=m_2d.pdf(data), vmin=0)
plt.colorbar(label='PDF')