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"))