star_flags) # Write flags file prihdu = fits.PrimaryHDU(zpick.zwarning) thdulist = fits.HDUList([prihdu]) thdulist.writeto( '/uufs/astro.utah.edu/common/home/u0814744/scratch/screens/flags-%s.fits' % plate, overwrite=True) output = io.WriteRedmonster( zpick, dest='/uufs/astro.utah.edu/common/home/u0814744/scratch/screens', overwrite=True) # Read yanny file x = y.yanny(filename='spInspect_alltest_bolton.par.txt', np=True) # Get fibers, zpipe, zperson for each plate args = n.where(x['BOSSOBJECT']['plate'] == 3686)[0] fibers3686 = [] zpipe3686 = [] zperson3686 = [] for i in args: fibers3686.append(x['BOSSOBJECT'][i][2] - 1) zpipe3686.append(x['BOSSOBJECT'][i][5]) zperson3686.append(x['BOSSOBJECT'][i][6]) args = n.where(x['BOSSOBJECT']['plate'] == 3687)[0] fibers3687 = [] zpipe3687 = [] zperson3687 = []
import matplotlib.pyplot as p p.interactive(True) from redmonster.sandbox import yanny as y from astropy.io import fits from redmonster.datamgr import spec, io from redmonster.physics import zfinder, zfitter, zpicker from redmonster.physics import misc from time import gmtime, strftime from os.path import join from os import environ from scipy.optimize import curve_fit from os.path import join from os import environ # Read yanny file x = y.yanny(filename='spInspect_alltest_bolton.par.txt', np=True) # Get fibers, zpipe, zperson for each plate args = n.where(x['BOSSOBJECT']['plate'] == 3686)[0] fibers3686 = [] zpipe3686 = [] zperson3686 = [] classperson3686 = [] comments3686 = [] for i in args: fibers3686.append( x['BOSSOBJECT'][i][2]-1) zpipe3686.append( x['BOSSOBJECT'][i][5]) zperson3686.append( x['BOSSOBJECT'][i][6]) classperson3686.append( x['BOSSOBJECT'][i][7]) comments3686.append( x['BOSSOBJECT'][i][8])
import numpy as n from redmonster.sandbox import yanny as y # Read yanny file x = y.yanny(filename='/uufs/astro.utah.edu/common/home/u0814744/boss/\ spInspect_alltest_bolton.par.txt', np=True) # Get fibers, zpipe, zperson for each plate args = n.where(x['BOSSOBJECT']['plate'] == 3686)[0] fibers3686 = [] zpipe3686 = [] zperson3686 = [] for i in args: fibers3686.append( x['BOSSOBJECT'][i][2]) zpipe3686.append( x['BOSSOBJECT'][i][5]) zperson3686.append( x['BOSSOBJECT'][i][6]) args = n.where(x['BOSSOBJECT']['plate'] == 3687)[0] fibers3687 = [] zpipe3687 = [] zperson3687 = [] for i in args: fibers3687.append( x['BOSSOBJECT'][i][2]) zpipe3687.append( x['BOSSOBJECT'][i][5]) zperson3687.append( x['BOSSOBJECT'][i][6]) args = n.where(x['BOSSOBJECT']['plate'] == 3804)[0] fibers3804 = [] zpipe3804 = [] zperson3804 = []
import matplotlib.pyplot as p p.interactive(True) from redmonster.sandbox import yanny as y from astropy.io import fits from redmonster.datamgr import spec, io from redmonster.physics import zfinder, zfitter, zpicker from redmonster.physics import misc from time import gmtime, strftime import multiprocessing as mp startfib = 100 nfib = 1 # Read yanny file x = y.yanny(filename='/Users/boltonlab3/Downloads/spInspect_alltest_bolton.par.txt', np=True) # Get fibers, zpipe, zperson for each plate args = n.where(x['BOSSOBJECT']['plate'] == 3686)[0] fibers3686 = [] zpipe3686 = [] zperson3686 = [] for i in args: fibers3686.append( x['BOSSOBJECT'][i][2]) zpipe3686.append( x['BOSSOBJECT'][i][5]) zperson3686.append( x['BOSSOBJECT'][i][6]) # Only keep 50 fibers fiberid = fibers3686[startfib:startfib+nfib] zpipe = zpipe3686[startfib:startfib+nfib] zperson = zperson3686[startfib:startfib+nfib]
import numpy as n from redmonster.sandbox import yanny as y # Read yanny file x = y.yanny(filename='/uufs/astro.utah.edu/common/home/u0814744/boss/\ spInspect_alltest_bolton.par.txt', np=True) # Get fibers, zpipe, zperson for each plate args = n.where(x['BOSSOBJECT']['plate'] == 3686)[0] fibers3686 = [] zpipe3686 = [] zperson3686 = [] for i in args: fibers3686.append(x['BOSSOBJECT'][i][2]) zpipe3686.append(x['BOSSOBJECT'][i][5]) zperson3686.append(x['BOSSOBJECT'][i][6]) args = n.where(x['BOSSOBJECT']['plate'] == 3687)[0] fibers3687 = [] zpipe3687 = [] zperson3687 = [] for i in args: fibers3687.append(x['BOSSOBJECT'][i][2]) zpipe3687.append(x['BOSSOBJECT'][i][5]) zperson3687.append(x['BOSSOBJECT'][i][6]) args = n.where(x['BOSSOBJECT']['plate'] == 3804)[0] fibers3804 = [] zpipe3804 = []