irs1_sample1_lryp.fits

display range: -20,20 Divide irs1_sample1_lrxpb.fits into 9 1D images Fit and subtract the residual backbround using object mask Merge 9 1D images

#Extract the parts of the spectra using 9 masks
!echo "1 0 0 0 0 0 0 0 0" > extract.dat
imdel irs1_sample1_lry1.fits
!../extract2 irs1_sample1_lrxpb.fits irs1_sample1_lry1.fits
!echo "0 1 0 0 0 0 0 0 0" > extract.dat
imdel irs1_sample1_lry2.fits
!../extract2 irs1_sample1_lrxpb.fits irs1_sample1_lry2.fits
!echo "0 0 1 0 0 0 0 0 0" > extract.dat
imdel irs1_sample1_lry3.fits
!../extract2 irs1_sample1_lrxpb.fits irs1_sample1_lry3.fits
!echo "0 0 0 1 0 0 0 0 0" > extract.dat
imdel irs1_sample1_lry4.fits
!../extract2 irs1_sample1_lrxpb.fits irs1_sample1_lry4.fits
!echo "0 0 0 0 1 0 0 0 0" > extract.dat
imdel irs1_sample1_lry5.fits
!../extract2 irs1_sample1_lrxpb.fits irs1_sample1_lry5.fits
!echo "0 0 0 0 0 1 0 0 0" > extract.dat
imdel irs1_sample1_lry6.fits
!../extract2 irs1_sample1_lrxpb.fits irs1_sample1_lry6.fits
!echo "0 0 0 0 0 0 1 0 0" > extract.dat
imdel irs1_sample1_lry7.fits
!../extract2 irs1_sample1_lrxpb.fits irs1_sample1_lry7.fits
!echo "0 0 0 0 0 0 0 1 0" > extract.dat
imdel irs1_sample1_lry8.fits
!../extract2 irs1_sample1_lrxpb.fits irs1_sample1_lry8.fits
!echo "0 0 0 0 0 0 0 0 1" > extract.dat
imdel irs1_sample1_lry9.fits
!../extract2 irs1_sample1_lrxpb.fits irs1_sample1_lry9.fits
!echo "0.5 1 1 1 1 1 1 1 0.5" > extract.dat
#Prepare object mask
imdel tmp1,tmp2
imcomb irs1_sample1_lry4,irs1_sample1_lry5,irs1_sample1_lry6 tmp1
imcomb irs1_sample1_lry1,irs1_sample1_lry2,irs1_sample1_lry3,irs1_sample1_lry7,irs1_sample1_lry8,irs1_sample1_lry9 tmp2
imarith tmp1 - tmp2 tmp1
imfunc tmp1 tmp1 abs
fit1d tmp1 tmp1 fit axis=1 inter- function=legendre order=1
imrepl tmp1 -1. upper=10
imrepl tmp1  0. lower=10
imarith tmp1 * -1 tmp1
imdel mask,mask2
imcopy tmp1 mask.fits
imarith mask.fits * irs1_sample1_lr_mask.fits mask.fits
#Fit and subtract the background using object mask
imdel tmp1
imarith irs1_sample1_lry1 / irs1_cont_lra irs1_sample1_lry1
!../fixbadpix irs1_sample1_lry1.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=20 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1_lry1 - tmp1 irs1_sample1_lry1
imarith irs1_sample1_lry1 * irs1_cont_lra irs1_sample1_lry1
imdel tmp1
imarith irs1_sample1_lry2 / irs1_cont_lra irs1_sample1_lry2
!../fixbadpix irs1_sample1_lry2.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=20 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1_lry2 - tmp1 irs1_sample1_lry2
imarith irs1_sample1_lry2 * irs1_cont_lra irs1_sample1_lry2
imdel tmp1
imarith irs1_sample1_lry3 / irs1_cont_lra irs1_sample1_lry3
!../fixbadpix irs1_sample1_lry3.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=20 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1_lry3 - tmp1 irs1_sample1_lry3
imarith irs1_sample1_lry3 * irs1_cont_lra irs1_sample1_lry3
imdel tmp1
imarith irs1_sample1_lry4 / irs1_cont_lra irs1_sample1_lry4
!../fixbadpix irs1_sample1_lry4.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=20 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1_lry4 - tmp1 irs1_sample1_lry4
imarith irs1_sample1_lry4 * irs1_cont_lra irs1_sample1_lry4
imdel tmp1
imarith irs1_sample1_lry5 / irs1_cont_lra irs1_sample1_lry5
!../fixbadpix irs1_sample1_lry5.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=20 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1_lry5 - tmp1 irs1_sample1_lry5
imarith irs1_sample1_lry5 * irs1_cont_lra irs1_sample1_lry5
imdel tmp1
imarith irs1_sample1_lry6 / irs1_cont_lra irs1_sample1_lry6
!../fixbadpix irs1_sample1_lry6.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=20 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1_lry6 - tmp1 irs1_sample1_lry6
imarith irs1_sample1_lry6 * irs1_cont_lra irs1_sample1_lry6
imdel tmp1
imarith irs1_sample1_lry7 / irs1_cont_lra irs1_sample1_lry7
!../fixbadpix irs1_sample1_lry7.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=20 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1_lry7 - tmp1 irs1_sample1_lry7
imarith irs1_sample1_lry7 * irs1_cont_lra irs1_sample1_lry7
imdel tmp1
imarith irs1_sample1_lry8 / irs1_cont_lra irs1_sample1_lry8
!../fixbadpix irs1_sample1_lry8.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=20 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1_lry8 - tmp1 irs1_sample1_lry8
imarith irs1_sample1_lry8 * irs1_cont_lra irs1_sample1_lry8
imdel tmp1
imarith irs1_sample1_lry9 / irs1_cont_lra irs1_sample1_lry9
!../fixbadpix irs1_sample1_lry9.fits tmp1.fits y
fit1d tmp1 tmp1 fit axis=2 inter- function=spline3 order=20 high_r=3 low_r=3 nite=10 grow=0
imarith irs1_sample1_lry9 - tmp1 irs1_sample1_lry9
imarith irs1_sample1_lry9 * irs1_cont_lra irs1_sample1_lry9
#Merge 9 images
imdel irs1_sample1_lryp.fits
imrename mask.fits irs1_sample1_lr_mask.fits
!../pileup irs1_sample1_lry
displ irs1_sample1_lryp 2