imcopy fibre_flat.fits[0] fibre_flat_img.imh
imcopy fibre_calibg.fits[0] fibre_calib_img.imh

imred
crutil
cosmicrays obj_none_30_15.fits[0] obj_none_cr.fits thresho=500 fluxrat=4 npasses=5 window=7
cosmicrays obj_QSOz14_H19_30_15.fits[0] obj_QSOz14_H19_cr.fits thresho=500 fluxrat=4 npasses=5 window=7
cosmicrays obj_Em_15_16_30_15.fits[0] obj_Em_15_16_cr.fits thresho=500 fluxrat=4 npasses=5 window=7
cosmicrays obj_Em_16_16_30_15.fits[0] obj_Em_16_16_cr.fits thresho=500 fluxrat=4 npasses=5 window=7
cosmicrays obj_flat_AB20_30_15.fits[0] obj_flat_AB20_cr.fits thresho=500 fluxrat=4 npasses=5 window=7
cosmicrays obj_flat_AB21_30_15.fits[0] obj_flat_AB21_cr.fits thresho=500 fluxrat=4 npasses=5 window=7

imarith obj_QSOz14_H19_crg.fits - obj_none_cr.fits obj_QSOz14_H19_ss.imh
imarith obj_Em_15_16_crg.fits - obj_none_cr.fits obj_Em_15_16_ss.imh
imarith obj_Em_16_16_crg.fits - obj_none_cr.fits obj_Em_16_16_ss.imh
imarith obj_flat_AB20_crg.fits - obj_none_cr.fits obj_flat_AB20_ss.imh
imarith obj_flat_AB21_crg.fits - obj_none_cr.fits obj_flat_AB21_ss.imh

imdel *_crg.fits

hydra
unlearn hydra

params.coordli = "arcXe.txt"
params.t_funct = "cheby"
params.t_order = 4
params.weights = "variance"
params.i_funct = "cheby"
params.sort = "MJD"
params.group = "MJD"
reid.search = 100.0
reid.cradius = 150.0

awk '{print NR,1}' obj_QSOz14_H19.list > list200

dohydra obj_QSOz14_H19_ss.imh apref=fibre_flat_img.imh flat=fibre_flat_img.imh \
arcs1=fibre_calib_img.imh readnoi=RON gain=GAIN datamax=INDEF fibers=200 width=6 \
minsep=8 maxsep=12 apidtab=list200 objbeam=0,1 skybeam=0 \
scatter- fitflat+ clean- dispcor+ savearc+ skysubt- skyedit- savesky- \
splot+ redo+ update+ batch- listonl-

dohydra obj_Em_15_16_ss.imh apref=fibre_flat_img.imh flat=fibre_flat_img.imh \
arcs1=fibre_calib_img.imh readnoi=RON gain=GAIN datamax=INDEF fibers=200 width=6 \
minsep=8 maxsep=12 apidtab=list200 objbeam=0,1 skybeam=0 \
scatter- fitflat+ clean- dispcor+ savearc+ skysubt- skyedit- savesky- \
splot+ redo+ update+ batch- listonl-

dohydra obj_Em_16_16_ss.imh apref=fibre_flat_img.imh flat=fibre_flat_img.imh \
arcs1=fibre_calib_img.imh readnoi=RON gain=GAIN datamax=INDEF fibers=200 width=6 \
minsep=8 maxsep=12 apidtab=list200 objbeam=0,1 skybeam=0 \
scatter- fitflat+ clean- dispcor+ savearc+ skysubt- skyedit- savesky- \
splot+ redo+ update+ batch- listonl-

dohydra obj_flat_AB20_ss.imh apref=fibre_flat_img.imh flat=fibre_flat_img.imh \
arcs1=fibre_calib_img.imh readnoi=RON gain=GAIN datamax=INDEF fibers=200 width=6 \
minsep=8 maxsep=12 apidtab=list200 objbeam=0,1 skybeam=0 \
scatter- fitflat+ clean- dispcor+ savearc+ skysubt- skyedit- savesky- \
splot+ redo+ update+ batch- listonl-

dohydra obj_flat_AB21_ss.imh apref=fibre_flat_img.imh flat=fibre_flat_img.imh \
arcs1=fibre_calib_img.imh readnoi=RON gain=GAIN datamax=INDEF fibers=200 width=6 \
minsep=8 maxsep=12 apidtab=list200 objbeam=0,1 skybeam=0 \
scatter- fitflat+ clean- dispcor+ savearc+ skysubt- skyedit- savesky- \
splot+ redo+ update+ batch- listonl-

low_rej = 0.5 niterate=3 order=5

imarith obj_QSOz14_H19_ss.msg.imh / obj_QSOz14_H19_ss.ms.imh[*,50] obj_QSOz14_H19_nr.fits
imarith obj_Em_15_16_ss.msg.imh / obj_Em_15_16_ss.ms.imh[*,50] obj_Em_15_16_nr.fits
imarith obj_Em_16_16_ss.msg.imh / obj_Em_16_16_ss.ms.imh[*,50] obj_Em_16_16_nr.fits
imarith obj_flat_AB20_ss.msg.imh / obj_flat_AB20_ss.ms.imh[*,50] obj_flat_AB20_nr.fits
imarith obj_flat_AB21_ss.msg.imh / obj_flat_AB21_ss.ms.imh[*,50] obj_flat_AB21_nr.fits

oned
rspectext ../template/obj_F0V_J15 obj_F0V_J15.fits dtype=linear crval1=8000.0 cdelt1=1.0
dispcor obj_F0V_J15.fits obj_F0V_J15_cor.fits w1 =9984.014 w2=18194.830 dw=3.9955 nw=2056 flux-

imarith obj_QSOz14_H19_nrg.fits * obj_F0V_J15_cor3.fits obj_QSOz14_H19_cal.fits
imarith obj_Em_15_16_nrg.fits * obj_F0V_J15_cor.fits obj_Em_15_16_cal.fits
imarith obj_Em_16_16_nrg.fits * obj_F0V_J15_cor.fits obj_Em_16_16_cal.fits
imarith obj_flat_AB20_nrg.fits * obj_F0V_J15_cor.fits obj_flat_AB20_cal.fits
imarith obj_flat_AB21_nrg.fits * obj_F0V_J15_cor.fits obj_flat_AB21_cal.fits

/* Simulate 5 hours exposure : 30min * 10 */

blkavg obj_Em_16_16_calg.fits[*,2:11] obj_Em_16_16_10h.fits 1 10 option=average

/* Fit emission line width and strength */

fitprofs obj_Em_15_16_calg.fits region="15800 16200" \
lines="!1,50,100,150,200" logfile=obj_flat_Em_15_16_log \
positio="Em_16" profile="gaussian" gfwhm=20.0 fitback=yes fitposi=all fitgfwh=all

awk '{if($1!="#") print $1,$2,$3,$6}' obj_Em_15_16_log
