The following are targets for this run:
This time, using a priority list with cuts from Risa. Another difference is a change from slits with 6" on each side to 3"
In [345]:
# This changes the current directory to the base saga directory - make sure to run this first!
# This is necessary to be able to import the py files and use the right directories,
# while keeping all the notebooks in their own directory.
import os
import sys
if 'saga_base_dir' not in locals():
saga_base_dir = os.path.abspath('..')
if saga_base_dir not in sys.path:
os.chdir(saga_base_dir)
In [346]:
for module in ['hosts', 'targeting', 'magellan']:
if module in globals():
reload(globals()[module])
else:
globals()[module] = __import__(module)
#g = targeting.get_gama() #re-caches the gama catalog
from astropy import units as u
from astropy.io import ascii
In [347]:
%matplotlib inline
from matplotlib.pyplot import *
from matplotlib import rcParams
rcParams['image.interpolation'] = 'none'
rcParams['figure.figsize'] = (10,10)
In [6]:
goog_hosts = hosts.get_saga_hosts_from_google('etollerud@gmail.com')
Out[6]:
In [99]:
targ_hosts = [h for h in goog_hosts if h.name.lower() in 'gilgamesh,aeneid,anak,iliad,ulysses'.split(',')]
targ_hosts.append(hosts.NSAHost(132339))
targ_hosts.append(hosts.NSAHost(131531))
For Gilgamesh, we design 2 masks, because of airmass limitations and the need to get the two -2 pri targets
In [483]:
currhost = [h for h in targ_hosts if h.name.lower()=='gilgamesh'][0]
currhost.shortname = 'Gilg'
In [484]:
currhost.physical_to_projected(300*u.kpc)
Out[484]:
In [487]:
target_catalog = ascii.read('catalogs/2014Jul_imacs_pris/imacs_{0}.dat'.format(currhost.nsaid),
names=['objID', 'ra', 'dec', 'r', 'imacs_pri'])
magellan.build_imacs_targeting_files(currhost, 'Marchi/Munoz', '2014-08-24', targs=target_catalog,
inclhost=False, weakdlimit=True, pdecide=1, onlygals=False, overwrite=False)
#onlygals=False because this target catalog already has been cleaned of stars
In [513]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True,skipnums=[])
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')
For AnaK we design 4 masks, to try to get near-completeness in pri1
In [471]:
currhost = [h for h in targ_hosts if h.name=='AnaK'][0]
currhost.shortname = 'AnaK'
In [362]:
currhost.physical_to_projected(300*u.kpc)
Out[362]:
In [369]:
target_catalog = ascii.read('catalogs/2014Jul_imacs_pris/imacs_{0}.dat'.format(currhost.nsaid),
names=['objID', 'ra', 'dec', 'r', 'imacs_pri'])
magellan.build_imacs_targeting_files(currhost, 'Marchi/Munoz', '2014-08-24', targs=target_catalog,
inclhost=False, weakdlimit=True, pdecide=1, onlygals=False, overwrite=False)
#onlygals=False because this target catalog already has been cleaned of stars
Now actually do the selection w/ intgui
In [472]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True)
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')
For all others, 2 masks/host
In [473]:
currhost = [h for h in targ_hosts if h.name.lower()=='aeneid'][0]
currhost.shortname = 'Aen2'
In [376]:
currhost.physical_to_projected(300*u.kpc)
Out[376]:
In [377]:
target_catalog = ascii.read('catalogs/2014Jul_imacs_pris/imacs_{0}.dat'.format(currhost.nsaid),
names=['objID', 'ra', 'dec', 'r', 'imacs_pri'])
magellan.build_imacs_targeting_files(currhost, 'Marchi/Munoz', '2014-08-24', targs=target_catalog,
inclhost=False, weakdlimit=True, pdecide=1, onlygals=False, overwrite=False)
#onlygals=False because this target catalog already has been cleaned of stars
In [474]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True)
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')
For all others, 2 masks/host
In [475]:
currhost = [h for h in targ_hosts if h.name.lower()=='iliad'][0]
currhost.shortname = 'Ilid'
In [390]:
currhost.physical_to_projected(300*u.kpc)
Out[390]:
In [393]:
target_catalog = ascii.read('catalogs/2014Jul_imacs_pris/imacs_{0}.dat'.format(currhost.nsaid),
names=['objID', 'ra', 'dec', 'r', 'imacs_pri'])
magellan.build_imacs_targeting_files(currhost, 'Marchi/Munoz', '2014-08-24', targs=target_catalog,
inclhost=False, weakdlimit=True, pdecide=1, onlygals=False, overwrite=False)
#onlygals=False because this target catalog already has been cleaned of stars
In [476]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True)
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')
For all others, 2 masks/host
In [477]:
currhost = [h for h in targ_hosts if h.name.lower()=='ulysses'][0]
currhost.shortname = 'Ulys'
In [399]:
currhost.physical_to_projected(300*u.kpc)
Out[399]:
In [400]:
target_catalog = ascii.read('catalogs/2014Jul_imacs_pris/imacs_{0}.dat'.format(currhost.nsaid),
names=['objID', 'ra', 'dec', 'r', 'imacs_pri'])
magellan.build_imacs_targeting_files(currhost, 'Marchi/Munoz', '2014-08-24', targs=target_catalog,
inclhost=False, weakdlimit=True, pdecide=1, onlygals=False, overwrite=False)
#onlygals=False because this target catalog already has been cleaned of stars
In [478]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True)
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')
For all others, 2 masks/host
In [479]:
currhost = [h for h in targ_hosts if h.nsaid==132339][0]
currhost.shortname = 'NS39'
In [406]:
currhost.physical_to_projected(300*u.kpc)
Out[406]:
In [407]:
target_catalog = ascii.read('catalogs/2014Jul_imacs_pris/imacs_{0}.dat'.format(currhost.nsaid),
names=['objID', 'ra', 'dec', 'r', 'imacs_pri'])
magellan.build_imacs_targeting_files(currhost, 'Marchi/Munoz', '2014-08-24', targs=target_catalog,
inclhost=False, weakdlimit=True, pdecide=1, onlygals=False, overwrite=False)
#onlygals=False because this target catalog already has been cleaned of stars
In [480]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True)
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')
For all others, 2 masks/host
In [481]:
currhost = [h for h in targ_hosts if h.nsaid==132339][0]
currhost.shortname = 'NS31'
In [438]:
currhost.physical_to_projected(300*u.kpc)
Out[438]:
In [439]:
target_catalog = ascii.read('catalogs/2014Jul_imacs_pris/imacs_{0}.dat'.format(currhost.nsaid),
names=['objID', 'ra', 'dec', 'r', 'imacs_pri'])
magellan.build_imacs_targeting_files(currhost, 'Marchi/Munoz', '2014-08-24', targs=target_catalog,
inclhost=False, weakdlimit=True, pdecide=1, onlygals=False, overwrite=False)
#onlygals=False because this target catalog already has been cleaned of stars
In [482]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True)
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')