The following are targets for this run:

  • gilgamesh
  • aeneid
  • anak
  • iliad
  • ulysses
  • nsa 132339
  • nsa 131531

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


Password for "etollerud@gmail.com":········
Out[6]:
[<hosts.NSAHost object w/ name 'Odyssey' AKA: ['NGC6181', 'NSA147100']>,
 <hosts.NSAHost object w/ name 'Iliad' AKA: ['NGC7393', 'NSA150238']>,
 <hosts.NSAHost object w/ name 'Alice' AKA: ['NGC4030', 'NSA140594']>,
 <hosts.NSAHost object w/ name 'Beowulf' AKA: ['NGC2750', 'NSA135667']>,
 <hosts.NSAHost object w/ name 'Gilgamesh' AKA: ['NGC5962', 'NSA166313']>,
 <hosts.NSAHost object w/ name 'Hamlet' AKA: ['NGC5899', 'NSA166035']>,
 <hosts.NSAHost object w/ name 'Aeneid' AKA: ['NGC-99', 'NSA148734']>,
 <hosts.NSAHost object w/ name 'AnaK' AKA: ['NGC7716', 'NSA61945']>,
 <hosts.NSAHost object w/ name 'Ulysses' AKA: ['NGC7428', 'NSA150340']>,
 <hosts.NSAHost object w/ name 'Sopranos' AKA: ['NGC5750', 'NSA165536']>,
 <hosts.NSAHost object w/ name 'Dune' AKA: ['NGC4045', 'NSA13927']>]

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

Gilgamesh

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]:
$37.1883 \; \mathrm{{}^{\prime}}$

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


Wrote catalog to imacs_targets/Gilg.cat
Wrote obs file to imacs_targets/Gilg_ini.obs

In [513]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True,skipnums=[])
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')


['imacs_targets/Gilg_1.SMF', 'imacs_targets/Gilg_2.SMF']
Total targets already observed= 240
Priority distribution of slitted targets:
Priority -2.0 : 2 targeted, 0 remaining
Priority 1.0 : 4 targeted, 13 remaining
Priority 2.0 : 139 targeted, 800 remaining
Priority 3.0 : 2 targeted, 54 remaining
Priority 4.0 : 2 targeted, 2907 remaining
Priority 5.0 : 91 targeted, 4206 remaining

AnaK

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]:
$29.9346 \; \mathrm{{}^{\prime}}$

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


Wrote catalog to imacs_targets/AnaK.cat
Wrote obs file to imacs_targets/AnaK_ini.obs

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


['imacs_targets/AnaK_1.SMF', 'imacs_targets/AnaK_2.SMF', 'imacs_targets/AnaK_3.SMF', 'imacs_targets/AnaK_4.SMF']
Total targets already observed= 406
Priority distribution of slitted targets:
Priority 1.0 : 36 targeted, 33 remaining
Priority 2.0 : 140 targeted, 166 remaining
Priority 3.0 : 7 targeted, 438 remaining
Priority 4.0 : 10 targeted, 2037 remaining
Priority 5.0 : 213 targeted, 2864 remaining

Aeneid

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]:
$19.8766 \; \mathrm{{}^{\prime}}$

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


Wrote catalog to imacs_targets/Aen2.cat
Wrote obs file to imacs_targets/Aen2_ini.obs

In [474]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True)
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')


['imacs_targets/Aen2_1.SMF', 'imacs_targets/Aen2_2.SMF']
Total targets already observed= 185
Priority distribution of slitted targets:
Priority 1.0 : 23 targeted, 15 remaining
Priority 2.0 : 54 targeted, 57 remaining
Priority 3.0 : 0 targeted, 619 remaining
Priority 5.0 : 108 targeted, 3743 remaining

Iliad

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]:
$19.2734 \; \mathrm{{}^{\prime}}$

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


Wrote catalog to imacs_targets/Ilid.cat
Wrote obs file to imacs_targets/Ilid_ini.obs

In [476]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True)
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')


['imacs_targets/Ilid_1.SMF', 'imacs_targets/Ilid_2.SMF']
Total targets already observed= 230
Priority distribution of slitted targets:
Priority 1.0 : 45 targeted, 20 remaining
Priority 2.0 : 126 targeted, 218 remaining
Priority 3.0 : 1 targeted, 511 remaining
Priority 5.0 : 58 targeted, 5317 remaining

Ulysses

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]:
$24.6749 \; \mathrm{{}^{\prime}}$

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


Wrote catalog to imacs_targets/Ulys.cat
Wrote obs file to imacs_targets/Ulys_ini.obs

In [478]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True)
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')


['imacs_targets/Ulys_1.SMF', 'imacs_targets/Ulys_2.SMF']
Total targets already observed= 230
Priority distribution of slitted targets:
Priority 1.0 : 37 targeted, 50 remaining
Priority 2.0 : 121 targeted, 334 remaining
Priority 3.0 : 0 targeted, 444 remaining
Priority 5.0 : 72 targeted, 5221 remaining

NSA 132339

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]:
$27.9989 \; \mathrm{{}^{\prime}}$

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


Wrote catalog to imacs_targets/NS39.cat
Wrote obs file to imacs_targets/NS39_ini.obs

In [480]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True)
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')


['imacs_targets/NS39_1.SMF', 'imacs_targets/NS39_2.SMF']
Total targets already observed= 242
Priority distribution of slitted targets:
Priority 1.0 : 51 targeted, 73 remaining
Priority 2.0 : 122 targeted, 560 remaining
Priority 3.0 : 0 targeted, 498 remaining
Priority 5.0 : 69 targeted, 6028 remaining

NSA 131531

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]:
$27.9989 \; \mathrm{{}^{\prime}}$

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


Wrote catalog to imacs_targets/NS31.cat
Wrote obs file to imacs_targets/NS31_ini.obs

In [482]:
magellan.plot_imacs_masks(currhost, eastleft=True, plotpris=True, showrvir=True, showfootprint=True)
plt.savefig('imacs_targets/'+currhost.shortname+'_imacs.png')


['imacs_targets/NS31_1.SMF', 'imacs_targets/NS31_2.SMF']
Total targets already observed= 241
Priority distribution of slitted targets:
Priority 1.0 : 41 targeted, 83 remaining
Priority 2.0 : 131 targeted, 551 remaining
Priority 3.0 : 0 targeted, 498 remaining
Priority 5.0 : 69 targeted, 6028 remaining