#!/bin/csh
alias echo "echo > /dev/null"
kappa
figaro
unalias echo
gdclear
lutheat
convert
#
echo " "
echo " *** FOCUS ON ARC SPECTRUM *** "
echo " "
echo " Process a group of sequential frames taken with different "
echo " instrument/array focus setting. "
echo " Extract central row from each spectrum and fit gaussian to "
echo " the peak in each spectrum: note FWHM. "
echo " Can then plot the widths of the profiles against UIST focus "
echo " to get best focus for grism/slit combo. "
echo " "
echo " Use QUICK_LOOK _raw frames"
echo " "
set date=20060202
echo "Date set to" $date
echo -n "Enter the number of the first frame in sequence: "
set f1="$<"
echo -n "Enter the number of the last frame in sequence: "
set fn="$<"
set n = $f1
while ($n <= $fn)
if ($n < 10) set numcal=0000${n}
if ($n > 9 && $n < 100) set numcal=000${n}
if ($n > 99 && $n < 1000) set numcal=00${n}
if ($n > 999) set numcal=0${n}
# gaiadisp u${date}_${numcal}_raw
echo " Extract centre row, middle"
## LONG
extract image=u${date}_${numcal}_raw ystart=545 yend=550 spectrum=middle
echo " "
echo " Assuming there is a line in range mask1 to mask2 below !"
echo " "
echo " "
echo " ******************************************** "
echo " ******************************************** "
echo " GAUSSIAN FITTING... "
echo " "
fitgauss in=middle device=xw mask1=360 mask2=400
rm middle.sdf
@ n = $n + 1
end