notebook.community

Edit and run

Targets are SW/NGC5485 and Beowulf, but targeting only those that don't have successful redshifts. 




In [ ]:



    


# 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 [1]:



    


import wiyn
import mmthecto
import magellan
import targeting
import hosts

from astropy import coordinates as coords
from astropy import units as u
from astropy.io import ascii



    


















    






WARNING:astropy:AstropyBackwardsIncompatibleChangeWarning: Using the `ICRSCoordinates` name for this class is deprecated in v0.3, and will be removed in the next version. Use `ICRS` instead.










    






WARNING: AstropyBackwardsIncompatibleChangeWarning: Using the `ICRSCoordinates` name for this class is deprecated in v0.3, and will be removed in the next version. Use `ICRS` instead. [astropy.coordinates.old_builtin_systems_names]

















In [2]:



    


h1 = hosts.beowulf
h2 = hosts.starwars

Beowulf





In [5]:



    


print h1.distmod,h1.projected_to_physical(30*u.arcmin)



    


















    






33.0456834686 354.8003713 kpc

















In [6]:



    


h1.usnob_environs_query(dl=True)



    


















    






File catalogs/Beowulf_usnob.dat exists - not downloading anything.

















In [7]:



    


h1.sdss_environs_query(dl=True, inclphotzs=True)



    


















    






File catalogs/Beowulf_sdss.dat exists - not downloading anything.

















In [8]:



    


fullcat = h1.get_sdss_catalog()



    











In [95]:



    


# a master catalog already exists from Feb, with bad fops removed.
# need to remove these already observed, or those overly-clustered
tomatch = ascii.read('hydra_targets/Beowulf_successful_premar2014.dat')
gotzcoos = coords.ICRS(tomatch['RA']*u.deg, tomatch['DEC']*u.deg)
removed, reason = wiyn.update_master_catalog('hydra_targets/{h1.name}.cat'.format(**locals()),
                           'hydra_targets/{h1.name}_take2.cat'.format(**locals()),
                           consolidategrps=None, removecoords=gotzcoos)



    


















    






191 objects in master catalog matched to be removed of 191 possible

















In [96]:



    


#only useful if consolidategrps above is not None
removedduetogrp = [a for a,b in zip(removed, reason) if b!='matched']
targeting.sampled_imagelist([coords.Angle(line.split()[2], u.hourangle).deg for line in removedduetogrp], 
                            [coords.Angle(line.split()[3], u.deg).deg for line in removedduetogrp]);



    











In [ ]:



    


wiyn.imagelist_from_master('hydra_targets/{h1.name}_take2.cat'.format(**locals()),'target')



    











In [114]:



    


h1.sdss_image_cutout(targets=wiyn.parse_master('hydra_targets/{h1.name}.cat'.format(**locals()),'target'),
                     scale=30*u.arcmin, imagesize=(1024, 1024),raoffset=0*u.arcmin)
#use this with the next two to figure out what to manually remove because it's the same object twice



    


















    
Out[114]:
























In [113]:



    


h1.sdss_image_cutout(targets=wiyn.parse_master('hydra_targets/{h1.name}_take2.cat'.format(**locals()),'target'),
                     scale=30*u.arcmin, imagesize=(1024, 1024),raoffset=0*u.arcmin)
#use this with the previous  and next to figure out what to manually remove because it's the same object twice



    


















    
Out[113]:
























In [122]:



    


coo1 = coords.ICRS(136.35975*u.deg, 25.55069*u.deg)
ang1 = 1*u.arcmin
print wiyn.imagelist_from_master('hydra_targets/{h1.name}_take2.cat'.format(**locals()),'target',nearloc=(coo1, ang1))

coo2 = coords.ICRS(136.57173*u.deg, 25.58107*u.deg)
ang2 = 30*u.arcsec
print wiyn.imagelist_from_master('hydra_targets/{h1.name}_take2.cat'.format(**locals()),'target',nearloc=(coo2, ang2))
#use these with the previous two to figure out what to manually remove because it's the same object twice



    


















    






name ra dec
J090525.93+253304.3   136.358045833 25.5512083333
J090526.11+253249.8   136.3588125 25.5471916667
name ra dec
J090617.70+253500.5   136.573770833 25.5834888889
J090618.38+253448.4   136.576620833 25.5801166667

















In [160]:



    


#manually deleted these from hydra_targets/Beowulf_take2.cat to yield
h1.sdss_image_cutout(targets=wiyn.parse_master('hydra_targets/{h1.name}_take2.cat'.format(**locals()),'target'),
                     scale=30*u.arcmin, imagesize=(1024, 1024),raoffset=0*u.arcmin)



    


















    
Out[160]:
























In [126]:



    


#repeat this as much as necessary to get all configs, running do_whydra on turtle each time
#NOTE: use 20 fibers, not default: do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h1.name}_take2.cat'.format(**locals()), 
                       lst=9,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs [] for removing from master
Writing to hydra_targets/Beowulf_take2_1.ast
SCP commands:
scp hydra_targets/Beowulf_take2_1.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/Beowulf_take2_1.hydra" hydra_targets

















In [127]:



    


#repeat this as much as necessary to get all configs, running do_whydra on turtle each time
#NOTE: use 20 fibers, not default: do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h1.name}_take2.cat'.format(**locals()), 
                       lst=9,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs ['hydra_targets/Beowulf_take2_1.hydra'] for removing from master
Writing to hydra_targets/Beowulf_take2_2.ast
SCP commands:
scp hydra_targets/Beowulf_take2_2.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/Beowulf_take2_2.hydra" hydra_targets

















In [128]:



    


#repeat this as much as necessary to get all configs, running do_whydra on turtle each time
#NOTE: use 20 fibers, not default: do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h1.name}_take2.cat'.format(**locals()), 
                       lst=9,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs ['hydra_targets/Beowulf_take2_1.hydra', 'hydra_targets/Beowulf_take2_2.hydra'] for removing from master
Writing to hydra_targets/Beowulf_take2_3.ast
SCP commands:
scp hydra_targets/Beowulf_take2_3.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/Beowulf_take2_3.hydra" hydra_targets

















In [130]:



    


#repeat this as much as necessary to get all configs, running do_whydra on turtle each time
#NOTE: use 20 fibers, not default: do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h1.name}_take2.cat'.format(**locals()), 
                       lst=9,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs ['hydra_targets/Beowulf_take2_1.hydra', 'hydra_targets/Beowulf_take2_2.hydra', 'hydra_targets/Beowulf_take2_3.hydra'] for removing from master
Writing to hydra_targets/Beowulf_take2_4.ast
SCP commands:
scp hydra_targets/Beowulf_take2_4.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/Beowulf_take2_4.hydra" hydra_targets

















In [131]:



    


#repeat this as much as necessary to get all configs, running do_whydra on turtle each time
#NOTE: use 20 fibers, not default: do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h1.name}_take2.cat'.format(**locals()), 
                       lst=9,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs ['hydra_targets/Beowulf_take2_1.hydra', 'hydra_targets/Beowulf_take2_2.hydra', 'hydra_targets/Beowulf_take2_3.hydra', 'hydra_targets/Beowulf_take2_4.hydra'] for removing from master
Writing to hydra_targets/Beowulf_take2_5.ast
SCP commands:
scp hydra_targets/Beowulf_take2_5.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/Beowulf_take2_5.hydra" hydra_targets

















In [138]:



    


wiyn.imagelist_fibers('hydra_targets/Beowulf_take2_1.hydra','target')
wiyn.imagelist_fibers('hydra_targets/Beowulf_take2_2.hydra','target')
wiyn.imagelist_fibers('hydra_targets/Beowulf_take2_3.hydra','target')
wiyn.imagelist_fibers('hydra_targets/Beowulf_take2_4.hydra','target')
wiyn.imagelist_fibers('hydra_targets/Beowulf_take2_5.hydra','target')



    


















    
Out[138]:








'name ra dec\nJ090659.86+255105.6   136.749425 25.8515583333\nJ090613.17+254349.6   136.554883333 25.7304444444\nJ090711.01+251316.9   136.7959125 25.2213638889\nJ090622.09+253724.4   136.592058333 25.6234555556\nJ090453.08+255133.7   136.221195833 25.8593777778\nJ090341.24+252427.4   135.921833333 25.4076222222\nJ090605.64+251058.3   136.523520833 25.1828666667\nJ090655.04+253501.2   136.729333333 25.5836777778\nJ090428.06+251800.4   136.116954167 25.3001166667\nJ090716.68+251335.8   136.8195125 25.2266277778\nJ090437.66+250352.8   136.15695 25.0646694444\nJ090640.67+253415.5   136.669466667 25.5709722222\nJ090620.46+253310.7   136.5852875 25.5529888889\nJ090600.12+251341.5   136.5005 25.2281944444\nJ090700.08+254523.4   136.750341667 25.7565\nJ090701.78+254321.5   136.7574375 25.7226638889\nJ090504.92+250700.1   136.270504167 25.1167027778\nJ090454.14+250922.9   136.2256 25.1563861111\nJ090357.99+251137.1   135.991654167 25.1936583333\nJ090449.75+252153.8   136.2073125 25.3649694444\nJ090554.75+253124.3   136.478154167 25.5234388889\nJ090348.74+251838.6   135.9530875 25.3107444444\nJ090722.60+252249.2   136.844195833 25.3803444444'

Star Wars





In [132]:



    


h2.altnames



    


















    
Out[132]:








['LukeSkywalker', 'NGC5485', 'NSA53145']

















In [154]:



    


# a master catalog already exists from Feb, with bad fops removed.
# need to remove these already observed, or those overly-clustered
tomatch = ascii.read('hydra_targets/StarWars_successful_premar2014.dat')
gotzcoos = coords.ICRS(tomatch['RA']*u.deg, tomatch['DEC']*u.deg)
removed, reason = wiyn.update_master_catalog('hydra_targets/StarWars.cat'.format(**locals()),
                           'hydra_targets/{h2.name}_take2.cat'.format(**locals()),
                           consolidategrps=None, removecoords=gotzcoos)



    


















    






161 objects in master catalog matched to be removed of 161 possible

















In [148]:



    


#check for duplicate/shredded targets in *original* master cat
print h2.coords
h2.sdss_image_cutout(targets=wiyn.parse_master('hydra_targets/NGC5485.cat'.format(**locals()),'target'),
                     scale=30*u.arcmin, imagesize=(1024, 1024),raoffset=0*u.arcmin)
#only one obvious case to NE



    


















    






<ICRS RA=211.79724 deg, Dec=55.00172 deg>










    
Out[148]:
























In [155]:



    


#and cross with post-observe 
h2.sdss_image_cutout(targets=wiyn.parse_master('hydra_targets/{h2.name}_take2.cat'.format(**locals()),'target'),
                     scale=30*u.arcmin, imagesize=(1024, 1024),raoffset=0*u.arcmin)
#only two obvious case to NE and NWW



    


















    
Out[155]:
























In [157]:



    


coo1 = coords.ICRS(211.85401*u.deg, 55.10302*u.deg)
ang1 = 1*u.arcmin
print wiyn.imagelist_from_master('hydra_targets/{h2.name}_take2.cat'.format(**locals()),'target',nearloc=(coo1, ang1))

coo2 = coords.ICRS(211.37293*u.deg, 55.09409*u.deg)
ang2 = 50*u.arcsec
print wiyn.imagelist_from_master('hydra_targets/{h2.name}_take2.cat'.format(**locals()),'target',nearloc=(coo2, ang2))
#use these with the previous two to figure out what to manually remove because it's the same object twice



    


















    






name ra dec
J140723.93+550627.0   211.849708333 55.1075194444
J140724.64+550537.6   211.852683333 55.0937805556
J140720.91+550617.0   211.837141667 55.1047472222
name ra dec
J140531.10+550516.7   211.379591667 55.0879944444
J140529.91+550529.4   211.374625 55.0915138889

















In [159]:



    


#manually delete the above from the StarWars_take2 catalog to give:
h2.sdss_image_cutout(targets=wiyn.parse_master('hydra_targets/{h2.name}_take2.cat'.format(**locals()),'target'),
                     scale=60*u.arcmin, imagesize=(1024, 1024),raoffset=0*u.arcmin)



    


















    
Out[159]:
























In [162]:



    


#make these with do_hydra, but use *20* sky fibers per target:
# do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h2.name}_take2.cat'.format(**locals()), 
                       lst=11,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs [] for removing from master
Writing to hydra_targets/StarWars_take2_1.ast
SCP commands:
scp hydra_targets/StarWars_take2_1.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/StarWars_take2_1.hydra" hydra_targets

















In [163]:



    


#last one required 
#make these with do_hydra, but use *20* sky fibers per target:
# do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h2.name}_take2.cat'.format(**locals()), 
                       lst=11,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs ['hydra_targets/StarWars_take2_1.hydra'] for removing from master
Writing to hydra_targets/StarWars_take2_2.ast
SCP commands:
scp hydra_targets/StarWars_take2_2.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/StarWars_take2_2.hydra" hydra_targets

















In [164]:



    


#make these with do_hydra, but use *20* sky fibers per target:
# do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h2.name}_take2.cat'.format(**locals()), 
                       lst=11,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs ['hydra_targets/StarWars_take2_1.hydra', 'hydra_targets/StarWars_take2_2.hydra'] for removing from master
Writing to hydra_targets/StarWars_take2_3.ast
SCP commands:
scp hydra_targets/StarWars_take2_3.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/StarWars_take2_3.hydra" hydra_targets

















In [165]:



    


#make these with do_hydra, but use *20* sky fibers per target:
# do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h2.name}_take2.cat'.format(**locals()), 
                       lst=11,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs ['hydra_targets/StarWars_take2_1.hydra', 'hydra_targets/StarWars_take2_2.hydra', 'hydra_targets/StarWars_take2_3.hydra'] for removing from master
Writing to hydra_targets/StarWars_take2_4.ast
SCP commands:
scp hydra_targets/StarWars_take2_4.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/StarWars_take2_4.hydra" hydra_targets

















In [166]:



    


#make these with do_hydra, but use *20* sky fibers per target:
# do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h2.name}_take2.cat'.format(**locals()), 
                       lst=11,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs ['hydra_targets/StarWars_take2_1.hydra', 'hydra_targets/StarWars_take2_2.hydra', 'hydra_targets/StarWars_take2_3.hydra', 'hydra_targets/StarWars_take2_4.hydra'] for removing from master
Writing to hydra_targets/StarWars_take2_5.ast
SCP commands:
scp hydra_targets/StarWars_take2_5.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/StarWars_take2_5.hydra" hydra_targets

















In [167]:



    


#make these with do_hydra, but use *20* sky fibers per target:
# do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h2.name}_take2.cat'.format(**locals()), 
                       lst=11,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs ['hydra_targets/StarWars_take2_1.hydra', 'hydra_targets/StarWars_take2_2.hydra', 'hydra_targets/StarWars_take2_3.hydra', 'hydra_targets/StarWars_take2_4.hydra', 'hydra_targets/StarWars_take2_5.hydra'] for removing from master
Writing to hydra_targets/StarWars_take2_6.ast
SCP commands:
scp hydra_targets/StarWars_take2_6.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/StarWars_take2_6.hydra" hydra_targets

















In [168]:



    


#make these with do_hydra, but use *20* sky fibers per target:
# do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h2.name}_take2.cat'.format(**locals()), 
                       lst=11,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs ['hydra_targets/StarWars_take2_1.hydra', 'hydra_targets/StarWars_take2_2.hydra', 'hydra_targets/StarWars_take2_3.hydra', 'hydra_targets/StarWars_take2_4.hydra', 'hydra_targets/StarWars_take2_5.hydra', 'hydra_targets/StarWars_take2_6.hydra'] for removing from master
Writing to hydra_targets/StarWars_take2_7.ast
SCP commands:
scp hydra_targets/StarWars_take2_7.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/StarWars_take2_7.hydra" hydra_targets

















In [169]:



    


#make these with do_hydra, but use *20* sky fibers per target:
# do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h2.name}_take2.cat'.format(**locals()), 
                       lst=11,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs ['hydra_targets/StarWars_take2_1.hydra', 'hydra_targets/StarWars_take2_2.hydra', 'hydra_targets/StarWars_take2_3.hydra', 'hydra_targets/StarWars_take2_4.hydra', 'hydra_targets/StarWars_take2_5.hydra', 'hydra_targets/StarWars_take2_6.hydra', 'hydra_targets/StarWars_take2_7.hydra'] for removing from master
Writing to hydra_targets/StarWars_take2_8.ast
SCP commands:
scp hydra_targets/StarWars_take2_8.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/StarWars_take2_8.hydra" hydra_targets

















In [170]:



    


#make these with do_hydra, but use *20* sky fibers per target:
# do_whydra -s 20 <fieldname>
wiyn.generate_ast_file('hydra_targets/{h2.name}_take2.cat'.format(**locals()), 
                       lst=11,texp=2, faintmagcut=20, scpname='turtle')



    


















    






Using existing catalogs ['hydra_targets/StarWars_take2_1.hydra', 'hydra_targets/StarWars_take2_2.hydra', 'hydra_targets/StarWars_take2_3.hydra', 'hydra_targets/StarWars_take2_4.hydra', 'hydra_targets/StarWars_take2_5.hydra', 'hydra_targets/StarWars_take2_6.hydra', 'hydra_targets/StarWars_take2_7.hydra', 'hydra_targets/StarWars_take2_8.hydra'] for removing from master
Writing to hydra_targets/StarWars_take2_9.ast
SCP commands:
scp hydra_targets/StarWars_take2_9.ast turtle:/home/ejt26/hydra_simulator/whydra
scp "turtle:/home/ejt26/hydra_simulator/whydra/StarWars_take2_9.hydra" hydra_targets

Make the coordinate cache list





In [172]:



    


wiyn.generate_wiyn_cache('hydra_targets/mar2014.cachelist', 
                         ['hydra_targets/Beowulf_take2*.hydra', 'hydra_targets/StarWars_take2*.hydra'])

Content source: saga-survey/erik

Similar notebooks:

notebook.community | gallery | about