Skip to content

Instantly share code, notes, and snippets.

@mattgiguere

mattgiguere/vumps_comparison_code

Created September 9, 2015 19:06
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save mattgiguere/45a91878b77f170eb51e to your computer and use it in GitHub Desktop.
Save mattgiguere/45a91878b77f170eb51e to your computer and use it in GitHub Desktop.
Download ZIP
code for analyzing fringing in the MAOS/VUMPS quartz exposures
Raw
vumps_comparison_code
ft0826 = readfits('/tous/vumps/flats/vumps150826.hgh_flat.fits')
q0826 = readfits('/tous/vumps/fitspec/150826/rvumps150826.1034.fits')
tststr = readfits('/tous/vumps/fitspec/150826/rvumps150826.1031.fits')
vega = readfits('/tous/vumps/fitspec/150901/rvumps150901.1035.fits')
q0901 = readfits('/tous/vumps/fitspec/150901/rvumps150901.1023.fits')
items = ['150901 flat', '150826 flat', '150826 qtz', '150826 teststar', '150901 vega', '150901 qtz']
plot, ft29[*,40,1], col=0, /xsty, /ysty
oplot, q0901[1,*,40], col=120
oplot, ft0826[*,40,1], col=80
oplot, q0826[1,*,40], col=250
oplot, tststr[1,*,40], col=70
oplot, vega[1,*,40], col=210
al_legend, items, color=[0, 80, 250, 70, 210, 120], linestyle=[0,0,0,0, 0, 0]
; Now plot the normalized central region
begv=1500
endv=2500
ps_open, '~/Desktop/vumps_plots', /encaps, /color
plot, ft29[begv:endv,40,1]/max(ft29[begv:endv,40,1]) + 0.5, col=0, /xsty, /ysty, yran=[0.9, 1.5]
oplot, q0901[1,begv:endv,40]/max(q0901[1,begv:endv,40]) + 0.45, col=120
oplot, ft0826[begv:endv,40,1]/max(ft0826[begv:endv,40,1]) + 0.4, col=80
oplot, q0826[1,begv:endv,40]/max(q0826[1,begv:endv,40]) + 0.35, col=250
oplot, tststr[1,begv:endv,40]/max(tststr[1,begv:endv,40]) + 0.3, col=70
oplot, vega[1,begv:endv,40]/max(vega[1,begv:endv,40])+0.15, col=210
al_legend, items, color=[0, 80, 250, 70, 210, 120], linestyle=[0,0,0,0, 0, 0], /bottom
ps_close
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment