@@ -0,0 +1,771 @@
+# -*- coding: utf-8 -*-%0Adef highest_probability_pixel(infile):%0A%0A     %22%22%22%0A%0A     Find the highest probability pixel RA%5Bdeg%5D, DEC%5Bdeg%5D%0A%0A     %22%22%22%0A %0A     %0A     import healpy as hp%0A     import numpy as np%0A%0A     # read probability skymap%0A     hpx = hp.read_map(infile,verbose=False)%0A%0A     # number of pixels%0A     npix = len(hpx)%0A%0A     # nside: resolution for the HEALPix map%0A     nside = hp.npix2nside(npix)%0A%0A     # pixel position%0A     ipix_max = np.argmax(hpx)%0A     hpx%5Bipix_max%5D%0A%0A     # sky coordinate%0A     theta, phi = hp.pix2ang(nside, ipix_max)%0A%0A     ra_max = np.rad2deg(phi)%0A%0A     dec_max = np.rad2deg(0.5 * np.pi - theta)%0A%0A     print ' or insert the highest probability pixel located %5Cn at RA =', str(( '%25 .5f' %25 ra_max))+'%C2%B0', 'and Dec =',str(( '%25 .5f' %25 dec_max))+'%C2%B0.' %0A

Create highest_probability_pixel.py