calib_grism.prg

calib_grism.prg

#!/bin/csh 
alias echo "echo > /dev/null"
kappa
figaro
unalias echo

echo " "
echo " *** SLOPE OF CONTINNUM ***"
echo " "
echo " Extract columns near left and right edges of array "
echo " Fits gaussians to two profiles to check position (and FWHM)"
echo " "
echo " Use QUICK_LOOK _raw frames"
echo " "

set date=20030313
echo "Date set to" $date
echo -n "Number of calibration frame: "
  set in="$<"
    if ($in < 10) set numcal=0000${in}
    if ($in > 9 && $in < 100) set numcal=000${in}
    if ($in > 99 && $in < 1000) set numcal=00${in}
    if ($in > 999) set numcal=0${in}

echo " "

gaiadisp u${date}_${numcal}_raw

echo " Extract two columns, 20-pixels wide, centred on X=100 (left)"
echo " and X=900 (right)"
ystract image=u${date}_${numcal}_raw xstart=90 xend=110 spectrum=left
ystract image=u${date}_${numcal}_raw xstart=890 xend=910 spectrum=right

echo " "
echo " Assuming continuum/pinhole row in range y1=500 to y2=600..."
echo " "

echo " "
echo " ***************************************************  "
echo " ***************************************************  "
echo "             GAUSSIAN FITTING...              "
echo " "

fitgauss in=left device=xw mask1=500 mask2=600
echo " "
echo " ***************************************************  "
echo " ***************************************************  "
echo " *** Make a note of Y-posn. of peak at LEFT (above)"
echo " "

fitgauss in=right device=xw  mask1=500 mask2=600
echo " "
echo " ***************************************************  "
echo " *************************************************** "
echo " *** Make a note of Y-posn. of peak at RIGHT (above) "
echo " "

echo " "
echo " ***************************************************  "
echo " ***************************************************  "
echo " "
echo " Finally..."
echo " "

echo -n " Enter y-value at  LEFT  of spectrum: "
set ya = ( $< )
echo -n " Enter y-value at  RIGHT of spectrum: "
set yb = ( $< )

echo " "
echo " x1 - x2 = 800 (columns 100 and 900 extracted) "
set angle = `calc exp="atand(($ya-$yb)/800)"`

echo " "
echo " >>>> ANGLE OF CONTINUUM =   " $angle
echo " "
echo " "