def METHOD_NAME(self):
    nside = 256
    ra, dec, pix = random_sky(nside, allsky=False)
    self.assertEqual(len(ra), len(dec))
    surveypix = desimodel.footprint.tiles2pix(nside)
    theta = np.radians(90 - dec)
    phi = np.radians(ra)
    skypix = hp.ang2pix(nside, theta, phi, nest=True)
    self.assertEqual(set(surveypix), set(skypix))