def allsky_survey(pointing_file): dir = os.path.dirname(__file__) savedir = os.path.join(dir, '../../data/allsky_test/') # FOV Dimensions FOV_Dim = (2048 * 6.2 / 3600) * (pi / 180) # Base Dimension FOV_phi = FOV_Dim * 2 FOV_theta = FOV_Dim skymap = SkyMap2(nside=2**8, LVF_theta=FOV_theta, LVF_phi=FOV_phi, cap_theta=pi) #print skymap.skydata.sky_df #skymap.skydata.sky_df.save("testdf") pointings1 = read_in_date(pointing_file) pointings3 = pointings1[(pointings1[:, 0] <= 200)] pointings1 = pointings3[(pointings3[:, 2] <= pi + 1.5 * FOV_phi) & (pointings3[:, 2] >= pi - 1.5 * FOV_phi)] print "Length of pointings", pointings1.shape[0] days = np.unique(pointings1[:, 0]) print "Number of Days", len(days) # Initializations fig_sub, ax_sub = plt.subplots(1, 4, subplot_kw=dict(projection="mollweide")) #sum_plot_days=[] least_plot_days = [182, 183, 184, 185, 186, 187, 188, 189] for day in days: points_in = map(tuple, pointings1[pointings1[:, 0] == day, 1:4]) for i, tupler in enumerate(points_in): print "day:%s, step:%s" % (day, i) skymap.make_dicts(day, tupler) if (day in least_plot_days) & (i == 0): least_hits_arr = np.copy(skymap.least_hits()) theta_in = least_hits_arr[:, 0] phi_in = least_hits_arr[:, 1] hits_in = least_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in, "least_allsky: Day %s" % (day - 1), "least_allsky_7_day_%s" % (day - 1)) least_hits_arr = np.copy(skymap.least_hits()) theta_in = least_hits_arr[:, 0] phi_in = least_hits_arr[:, 1] hits_in = least_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in, "least_allsky: Day %s" % ("FINAL"), "least_allsky_7_day_%s" % ("FINAL")) np.save("check180_3", least_hits_arr) skymap.skydata.sky_df.save("testdf_3")
def allsky_survey(pointing_file): dir = os.path.dirname(__file__) savedir = os.path.join(dir, '../../data/allsky_test/') # FOV Dimensions FOV_Dim = (2048 * 6.2 / 3600) * (pi / 180) # Base Dimension FOV_phi = FOV_Dim * 2 FOV_theta = FOV_Dim skymap = SkyMap2(nside=2**8, LVF_theta=FOV_theta, LVF_phi=FOV_phi, cap_theta=pi) pointings1 = read_in_date(pointing_file) days = np.unique(pointings1[:, 0]) # Initializations fig_sub, ax_sub = plt.subplots(1, 4, subplot_kw=dict(projection="mollweide")) sum_plot_days = [2, 8, 9, 10, 11, 12, 13, 31] least_plot_days = [2, 8, 9, 10, 11, 12, 13, 31] for day in days: points_in = map(tuple, pointings1[pointings1[:, 0] == day, 1:4]) for i, tupler in enumerate(points_in): print "day:%s, step:%s" % (day, i) skymap.make_dicts(day, tupler) if (day in least_plot_days) & (i == 0): least_hits_arr = np.copy(skymap.least_hits()) theta_in = least_hits_arr[:, 0] phi_in = least_hits_arr[:, 1] hits_in = least_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in, "least_allsky: Day %s" % (day - 1), "least_allsky_8_day_%s" % (day - 1)) if (day in sum_plot_days) & (i == 0): sum_hits_arr = np.copy(skymap.sum_hits()) theta_in = sum_hits_arr[:, 0] phi_in = sum_hits_arr[:, 1] hits_in = sum_hits_arr[:, 2] moll_plot2(theta_in, phi_in, hits_in, "sum_allsky: Day %s" % (day - 1), "sum_allsky_8_day_%s" % (day - 1))
def allsky_survey(pointing_file): dir= os.path.dirname(__file__) savedir = os.path.join(dir, '../../data/allsky_test/') # FOV Dimensions FOV_Dim=(2048*6.2/3600)*(pi/180) # Base Dimension FOV_phi=FOV_Dim*2 FOV_theta=FOV_Dim skymap = SkyMap2(nside=2**8,LVF_theta=FOV_theta, LVF_phi=FOV_phi, cap_theta=pi) pointings1 = read_in_date(pointing_file) days = np.unique(pointings1[:,0]) # Initializations fig_sub, ax_sub= plt.subplots(1, 4, subplot_kw=dict(projection="mollweide")) sum_plot_days=[2,8,10,31, 45, 92,183,184,274] least_plot_days=[92,93,183,184,274,275] for day in days: points_in=map(tuple, pointings1[ pointings1[:,0]==day, 1:4] ) for i, tupler in enumerate(points_in): print "day:%s, step:%s" %(day, i) skymap.make_dicts(day, tupler ) if (day in least_plot_days) & (i==0): least_hits_arr=skymap.least_hits() theta_in=least_hits_arr[:, 0] phi_in=least_hits_arr[:, 1] hits_in=least_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in , "least_allsky: Day %s"%(day-1), "least_allsky_6_day_%s" %(day-1)) #if (day in sum_plot_days) & (i==0): # sum_hits_arr=skymap.sum_hits() # theta_in=sum_hits_arr[:, 0] # phi_in=sum_hits_arr[:, 1] # hits_in=sum_hits_arr[:, 2] # moll_plot(theta_in, phi_in, hits_in , "sum_allsky: Day %s"%(day-1), "sum_allsky_4_day_%s" %(day-1)) least_hits_arr=skymap.least_hits() theta_in=least_hits_arr[:, 0] phi_in=least_hits_arr[:, 1] hits_in=least_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in , "least_allsky: Day %s"%("FINAL"), "least_allsky_6_day_%s" %("FINAL"))
def allsky_survey(pointing_file): dir= os.path.dirname(__file__) savedir = os.path.join(dir, '../../data/allsky_test/') # FOV Dimensions FOV_Dim=(2048*6.2/3600)*(pi/180) # Base Dimension FOV_phi=FOV_Dim*2 FOV_theta=FOV_Dim skymap = SkyMap2(nside=2**8,LVF_theta=FOV_theta, LVF_phi=FOV_phi, cap_theta=pi) #print skymap.skydata.sky_df #skymap.skydata.sky_df.save("testdf") pointings1 = read_in_date(pointing_file) pointings3=pointings1[(pointings1[:,0]<= 200)] pointings1=pointings3[(pointings3[:,2]<=pi+1.5*FOV_phi) & (pointings3[:,2]>=pi-1.5*FOV_phi)] print "Length of pointings", pointings1.shape[0] days = np.unique(pointings1[:,0]) print "Number of Days", len(days) # Initializations fig_sub, ax_sub= plt.subplots(1, 4, subplot_kw=dict(projection="mollweide")) #sum_plot_days=[] least_plot_days=[182,183,184,185,186,187,188,189] for day in days: points_in=map(tuple, pointings1[ pointings1[:,0]==day, 1:4] ) for i, tupler in enumerate(points_in): print "day:%s, step:%s" %(day, i) skymap.make_dicts(day, tupler ) if (day in least_plot_days) & (i==0): least_hits_arr=np.copy(skymap.least_hits()) theta_in=least_hits_arr[:, 0] phi_in=least_hits_arr[:, 1] hits_in=least_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in , "least_allsky: Day %s"%(day-1), "least_allsky_7_day_%s" %(day-1)) least_hits_arr=np.copy(skymap.least_hits()) theta_in=least_hits_arr[:, 0] phi_in=least_hits_arr[:, 1] hits_in=least_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in , "least_allsky: Day %s"%("FINAL"), "least_allsky_7_day_%s" %("FINAL")) np.save("check180_3", least_hits_arr) skymap.skydata.sky_df.save("testdf_3")
for j in range(5): for i, lambdai in enumerate(lambda_in1): tupleri=(theta_in1[j], phi_in1[j], lambda_in1[i]) interim_list.append(tupleri) ex_list=[interim_list[0]] print sky_dat1.get_hits(ex_list) time3=time.time() sky_dat1.increment_hit(ex_list) time4=time.time() print sky_dat1.get_hits(ex_list) print "Time of Execution Ex: %s" %(time4-time3) print sky_dat1.get_hits(interim_list) time3=time.time() sky_dat1.increment_hit(interim_list) time4=time.time() print "Time of Execution Interim: %s" %(time4-time3) print sky_dat1.get_hits(interim_list) least_array1= sky_dat1.least_hits_array() print least_array1 thetas=least_array1[:,0] phis=least_array1[:,1] hits= least_array1[:,2] moll_plot(thetas, phis, hits, "least_allsky", "least_allsky_ex")
def allsky_survey(pointing_file, save_suffix): dir= os.path.dirname(__file__) savedir = os.path.join(dir, '../../data/allsky_test/') # FOV Dimensions FOV_Dim=(2048*6.2/3600)*(pi/180) # Base Dimension FOV_phi=FOV_Dim*2 FOV_theta=FOV_Dim Nstrip=21 # Number of strips on each band of FOV skymap = SkyMap2(nside=2**8,LVF_theta=FOV_theta, LVF_phi=FOV_phi, cap_theta=pi) pointings1 = read_in_date(pointing_file) days = np.unique(pointings1[:,0]) # Initializations skyplot1=SkyPlots2() fig_sub, ax_sub= plt.subplots(1, 4, subplot_kw=dict(projection="mollweide")) sub_ind=0 sum_plot_days=[2,91,182,273,365] least_plot_days=[91,181,273,365] plot_lambda_ranges=[(.75,2.34),(.75,.76),(1.30,1.33),(2.26,2.34)] for day in days: points_in=map(tuple, pointings1[ pointings1[:,0]==day, 1:4] ) for i, tupler in enumerate(points_in): print "day:%s, step:%s" %(day, i) skymap.make_dicts(day, tupler ) if (day in least_plot_days) & (i==0): least_hits_arr=skymap.least_hits() theta_in=least_hits_arr[:, 0] phi_in=least_hits_arr[:, 1] hits_in=least_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in , "least_allsky: Day %s"%(day), "least_allsky_4_day_%s" %(day)) if (day in sum_plot_days) & (i==0): sum_hits_arr=skymap.sum_hits() theta_in=sum_hits_arr[:, 0] phi_in=sum_hits_arr[:, 1] hits_in=sum_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in , "sum_allsky: Day %s"%(day), "sum_allsky_4_day_%s" %(day)) least_hits_arr=skymap.least_hits() theta_in=least_hits_arr[:, 0] phi_in=least_hits_arr[:, 1] hits_in=least_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in , "least_allsky: Day %s"%("FINAL"), "least_allsky_4_day_%s" %("FINAL")) sum_hits_arr=skymap.sum_hits() theta_in=sum_hits_arr[:, 0] phi_in=sum_hits_arr[:, 1] hits_in=sum_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in , "sum_allsky: Day %s"%("FINAL"), "sum_allsky_4_day_%s" %("FINAL"))
def allsky_survey(pointing_file, save_suffix): dir = os.path.dirname(__file__) savedir = os.path.join(dir, '../../data/allsky_test/') # FOV Dimensions FOV_Dim = (2048 * 6.2 / 3600) * (pi / 180) # Base Dimension FOV_phi = FOV_Dim * 2 FOV_theta = FOV_Dim Nstrip = 21 # Number of strips on each band of FOV skymap = SkyMap2(nside=2**8, LVF_theta=FOV_theta, LVF_phi=FOV_phi, cap_theta=pi) pointings1 = read_in_date(pointing_file) days = np.unique(pointings1[:, 0]) # Initializations skyplot1 = SkyPlots2() fig_sub, ax_sub = plt.subplots(1, 4, subplot_kw=dict(projection="mollweide")) sub_ind = 0 sum_plot_days = [2, 91, 182, 273, 365] least_plot_days = [91, 181, 273, 365] plot_lambda_ranges = [(.75, 2.34), (.75, .76), (1.30, 1.33), (2.26, 2.34)] for day in days: points_in = map(tuple, pointings1[pointings1[:, 0] == day, 1:4]) for i, tupler in enumerate(points_in): print "day:%s, step:%s" % (day, i) skymap.make_dicts(day, tupler) if (day in least_plot_days) & (i == 0): least_hits_arr = skymap.least_hits() theta_in = least_hits_arr[:, 0] phi_in = least_hits_arr[:, 1] hits_in = least_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in, "least_allsky: Day %s" % (day), "least_allsky_4_day_%s" % (day)) if (day in sum_plot_days) & (i == 0): sum_hits_arr = skymap.sum_hits() theta_in = sum_hits_arr[:, 0] phi_in = sum_hits_arr[:, 1] hits_in = sum_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in, "sum_allsky: Day %s" % (day), "sum_allsky_4_day_%s" % (day)) least_hits_arr = skymap.least_hits() theta_in = least_hits_arr[:, 0] phi_in = least_hits_arr[:, 1] hits_in = least_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in, "least_allsky: Day %s" % ("FINAL"), "least_allsky_4_day_%s" % ("FINAL")) sum_hits_arr = skymap.sum_hits() theta_in = sum_hits_arr[:, 0] phi_in = sum_hits_arr[:, 1] hits_in = sum_hits_arr[:, 2] moll_plot(theta_in, phi_in, hits_in, "sum_allsky: Day %s" % ("FINAL"), "sum_allsky_4_day_%s" % ("FINAL"))