# Reduction of Standard Stars: red spectra
# - steve; Nov. 17, 2004

# grab and preprocess red side data:
ls -1 lred*fits > list.red
vi list.red
:.,$s/.fits//
:x

# verify you got what you want
hselect lred* $I,i_title,LAMPS,SLITNAME,DICHNAME,GRANGLE yes

# twoamp & trim:
twoamp @list.red 
imcopy @list.red//cb[2:2047,101:900] @list.red//cbt

# flatten images & make arc lamp:
imrename lred0175cbt lamp_r

ls -1 lred*cbt.fits > list.flint
vi list.flint
:.,$s/.fits//g

imcombine @list.flint flat_r
response flat_r flat_r flat_rres

ls -1 lred*cbt.fits > list.rsci
vi list.rsci
:.,$s/.fits//g

imarith @list.rsci / flat_rres @list.rsci//f

# extract & wavelength calibrate:
apall lred0120cbtf lower=-15 upper=15 res-
apall lred0170cbtf lower=-15 upper=15 res-
apall lred0315cbtf lower=-15 upper=15 res-
apall lred0368cbtf lower=-15 upper=15 res-

apall lamp_r refer=lred0120cbtf int- fit- tra- rec- back-
identify lamp_r.ms

refspec lred0120cbtf.ms refer=lamp_r.ms
refspec lred0170cbtf.ms refer=lamp_r.ms
refspec lred0315cbtf.ms refer=lamp_r.ms
refspec lred0368cbtf.ms refer=lamp_r.ms
dispcor lred0120cbtf.ms lred0120cbtf.w
dispcor lred0170cbtf.ms lred0170cbtf.w
dispcor lred0315cbtf.ms lred0315cbtf.w
dispcor lred0368cbtf.ms lred0368cbtf.w

# create sensitivity function:
standard lred0120cbtf.w stnd_r560 feige110
standard lred0315cbtf.w stnd_r560 feige110
standard lred0170cbtf.w stnd_r560 hilt600
standard lred0368cbtf.w stnd_r560 hilt600

sensfunc stnd_r560 sens_r560

# try again with 680 dichroic obs
sensfunc stnd_rALL sens_rALL

imrename sens_r560.0001 sens_r560

# calibrate:
calibrate lred0120cbtf.w lred0120cbtf.f sens=sens_r560
calibrate lred0170cbtf.w lred0170cbtf.f sens=sens_r560
calibrate lred0315cbtf.w lred0315cbtf.f sens=sens_r560
calibrate lred0368cbtf.w lred0368cbtf.f sens=sens_r560