def search_mocs(mocname, mocfile): conn = get_conn() cur = conn.cursor() moc = MOCFinder(mocfile) conn.execute("""insert into mocs (moc_name, moc_file) values ('%s', '%s')""" % (mocname, mocfile)) for uid, ra, dec in cur.execute(""" select uid, ra, dec from clusters order by uid""").fetchall(): if moc.is_in(ra, dec): conn.execute("""insert into cluster_in_moc (uid, moc_name) values (%s, '%s')""" % (uid, mocname)) print mocname, uid conn.commit() cur.close()
cosphi = np.cos(phirange) sinphi = np.sin(phirange) zrange = np.linspace(-20., 20., num=Zsize) for ir, r in enumerate(np.linspace(0.1, 20., num=Rsize)): x0 = r*cosphi - R0 y0 = r*sinphi r0 = np.sqrt(r**2 + R0**2 - 2*r*R0*cosphi) l = np.arctan2(y0/r0, x0/r0) print ir for iphi, phi in enumerate(phirange): bval = np.zeros(Zsize) r1 = np.sqrt(r0[iphi]**2 + zrange**2) b = np.degrees(np.arctan2(-zrange/r1, r0[iphi]/r1)) ll = np.degrees(np.ones(Zsize)*l[iphi]) for iz, coo in enumerate(np.array(gal2eq(ll, b)).T): ra, de = coo #print ra, de if moc.is_in(ra, de): incmap[ir, iz] += 1 incmap = incmap / float(PhiSize) np.savetxt('%s.map' % sys.argv[1], incmap) fig, ax = plt.subplots() ax.imshow(incmap.T, cmap=cm.jet, interpolation='none') plt.yticks(xrange(0, 40, 5), ['%.2f' % -z for z in zrange[::5]]) plt.xlabel('R (kpc)') plt.ylabel('Z (kpc)') plt.savefig('%s.map.png' % sys.argv[1], format='png')