ciao_sbp.sh: Various cleanups

1 files changed, 105 insertions, 151 deletions

diff --git a/scripts/ciao_sbp.sh b/scripts/ciao_sbp.sh
index 5703467..b1831a6 100755
--- a/scripts/ciao_sbp.sh
+++ b/scripts/ciao_sbp.sh
@@ -1,17 +1,15 @@
 #!/bin/sh
 ##
-## Extract `surface brighness profile' after exposure correction finished.
-##
-## NOTES:
-## * Only ACIS-I (chip: 0-3) and ACIS-S (chip: 7) supported
+## Extract `surface brighness profile' specified by the SBP regions and
+## make use of the exposure map.
 ##
 ## Weitian LI <liweitianux@gmail.com>
 ## Created: 2012/08/16
 ##
-VERSION="v4.1"
-UPDATED="2017-02-06"
-##
-## ChangeLogs:
+## Change logs:
+## 2017-02-08, Weitian LI
+##   * Do not generate 'radius_sbp.txt' and 'flux_sbp.txt'
+##   * Small cleanups
 ## v4.1, 2017-02-06, Weitian LI
 ##   * Specify regions format and system for ds9
 ## v4.0, 2015/06/03, Aaron LI
@@ -30,9 +28,6 @@ UPDATED="2017-02-06"
 ##   * provide asol file to correct offset
 ##
 
-unalias -a
-export LC_COLLATE=C
-
 SCRIPT_PATH=`readlink -f $0`
 SCRIPT_DIR=`dirname ${SCRIPT_PATH}`
 CCDGAP_SCRIPT="chandra_ccdgap_rect.py"
@@ -60,8 +55,6 @@ case "$1" in
     -[hH]*|--[hH]*)
         printf "usage:\n"
         printf "    `basename $0` evt_e=<evt_e_file> reg=<sbp_reg> expmap=<exp_map> cellreg=<cell_reg> aspect=<asol_file> [bkg=<bkg> log=<logfile> ]\n"
-        printf "\nversion:\n"
-        printf "    ${VERSION}, ${UPDATED}\n"
         exit ${ERR_USG}
         ;;
 esac
@@ -73,21 +66,17 @@ E_RANGE="700:7000"
 # default ccd edge cut pixel (20pixel)
 DFT_CCD_EDGECUT=25
 # default `event file' which used to match `blanksky' files
-#DFT_EVT="_NOT_EXIST_"
 DFT_EVT_E="`\ls evt2_*_e*.fits 2> /dev/null`"
 # default expmap
 DFT_EXPMAP="`\ls expmap*.fits 2> /dev/null | head -n 1`"
 # default `radial region file' to extract surface brightness
-#DFT_SBP_REG="_NOT_EXIST_"
 DFT_SBP_REG="sbprofile.reg"
 # defalut pointsource region file
 DFT_CELL_REG="`\ls celld*.reg 2> /dev/null`"
 # defalut asol file
 DFT_ASOL_FILE="`\ls -1 pcad*asol*fits 2> /dev/null`"
-
 # default `log file'
 DFT_LOGFILE="sbp_`date '+%Y%m%d'`.log"
-
 ## default parameters }}}
 
 ## functions {{{
@@ -285,33 +274,29 @@ if [ "${ACIS_TYPE}" = "S" ]; then
     punlearn dmlist
     dmlist infile="${SKYFOV}[ccd_id=${CCD}][cols POS]" opt="data,clean" | awk '{for (i=1;i<=NF;i++) print $i }' |sed -e ':a;N;s/\n/,/;ta' | awk -F"]," '{print "polygon("$2}' | awk -F"NaN" '{print $1}' >${TMP_LIST}
     python ${SCRIPT_DIR}/${CCDGAP_SCRIPT} ${TMP_LIST} >${TMP_REC}
-    XC=` cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $1}'`
-    YC=` cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $2}'`
-    ADD_L=` cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $3/2}'`
-    ANG=` cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $5}'`
-    while [ 1 -eq 1 ]; do
-        if [ `echo "${ANG} < 0" |bc -l ` -eq 1 ] ; then
-            ANG=` echo " ${ANG} + 90 " | bc -l `
-        elif [ `echo "${ANG} >=90" |bc -l ` -eq 1 ] ; then
-            ANG=` echo " ${ANG} - 90 " | bc -l`
-        else
-            break
-        fi
-    done
-    ANG=`echo "${ANG}/180*3.1415926" |bc -l`
-    CCD_1_X_RAW=` echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)-$2*sin($3)}' `
-    CCD_2_X_RAW=` echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)-$2*sin($3)}' `
-    CCD_3_X_RAW=` echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)+$2*sin($3)}' `
-    CCD_4_X_RAW=` echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)+$2*sin($3)}' `
-    CCD_1_Y_RAW=` echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)-$2*sin($3)}' `
-    CCD_2_Y_RAW=` echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)+$2*sin($3)}' `
-    CCD_3_Y_RAW=` echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)-$2*sin($3)}' `
-    CCD_4_Y_RAW=` echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)+$2*sin($3)}' `
-    CCD_1_RAW=` echo "${CCD_1_X_RAW},${CCD_1_Y_RAW}"`
-    CCD_2_RAW=` echo "${CCD_2_X_RAW},${CCD_2_Y_RAW}"`
-    CCD_3_RAW=` echo "${CCD_3_X_RAW},${CCD_3_Y_RAW}"`
-    CCD_4_RAW=` echo "${CCD_4_X_RAW},${CCD_4_Y_RAW}"`
-    REG_CCD_RAW="`echo "polygon(${CCD_1_RAW}, ${CCD_2_RAW}, ${CCD_4_RAW}, ${CCD_3_RAW}) " `"
+    XC=`cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $1}'`
+    YC=`cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $2}'`
+    ADD_L=`cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $3/2}'`
+    ANG=`cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $5}'`
+    if [ `echo "${ANG} < 0" | bc -l` -eq 1 ] ; then
+        ANG=`echo "${ANG} + 90" | bc -l`
+    elif [ `echo "${ANG} >=90" | bc -l` -eq 1 ] ; then
+        ANG=`echo "${ANG} - 90" | bc -l`
+    fi
+    ANG=`echo "${ANG}/180*3.1415926" | bc -l`
+    CCD_1_X_RAW=`echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)-$2*sin($3)}' `
+    CCD_2_X_RAW=`echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)-$2*sin($3)}' `
+    CCD_3_X_RAW=`echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)+$2*sin($3)}' `
+    CCD_4_X_RAW=`echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)+$2*sin($3)}' `
+    CCD_1_Y_RAW=`echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)-$2*sin($3)}' `
+    CCD_2_Y_RAW=`echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)+$2*sin($3)}' `
+    CCD_3_Y_RAW=`echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)-$2*sin($3)}' `
+    CCD_4_Y_RAW=`echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)+$2*sin($3)}' `
+    CCD_1_RAW=`echo "${CCD_1_X_RAW},${CCD_1_Y_RAW}"`
+    CCD_2_RAW=`echo "${CCD_2_X_RAW},${CCD_2_Y_RAW}"`
+    CCD_3_RAW=`echo "${CCD_3_X_RAW},${CCD_3_Y_RAW}"`
+    CCD_4_RAW=`echo "${CCD_4_X_RAW},${CCD_4_Y_RAW}"`
+    REG_CCD_RAW="`echo "polygon(${CCD_1_RAW},${CCD_2_RAW},${CCD_4_RAW},${CCD_3_RAW}) "`"
     DX_2T1=$(echo "`echo ${CCD_2_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_1_RAW} |awk -F\, '{print $1}'`" |bc -l)
     DY_2T1=$(echo "`echo ${CCD_2_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_1_RAW} |awk -F\, '{print $2}'`" |bc -l)
     CCD_WIDTH=`echo "sqrt(${DX_2T1}*${DX_2T1}+${DY_2T1}*${DY_2T1})" | bc -l`
@@ -362,81 +347,77 @@ elif [ "${ACIS_TYPE}" = "I" ]; then
     # ACIS-I
     TMP_REG_FILE_CCD="_ccd_tmp.reg"
     [ -e "${TMP_REG_FILE_CCD}" ] && mv -f ${TMP_REG_FILE_CCD} ${TMP_REG_FILE_CCD}_bak
-    for i in `seq 0 3` ; do
-    punlearn dmlist
-    dmlist infile="${SKYFOV}[ccd_id=${i}][cols POS]" opt="data,clean" | awk '{for (i=1;i<=NF;i++) print $i }' |sed -e ':a;N;s/\n/,/;ta' | awk -F"]," '{print "polygon("$2}' | awk -F"NaN" '{print $1}' >${TMP_LIST}
-    python ${SCRIPT_DIR}/${CCDGAP_SCRIPT} ${TMP_LIST} >${TMP_REC}
-    XC=` cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $1}'`
-    YC=` cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $2}'`
-    ADD_L=` cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $3/2}'`
-    ANG=` cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $5}'`
-    while [ 1 -eq 1 ]; do
-        if [ `echo "${ANG} < 0" |bc -l ` -eq 1 ] ; then
-            ANG=` echo " ${ANG} + 90 " | bc -l `
-        elif [ `echo "${ANG} >=90" |bc -l ` -eq 1 ] ; then
-            ANG=` echo " ${ANG} - 90 " | bc -l`
-        else
-            break
+    for i in `seq 0 3`; do
+        punlearn dmlist
+        dmlist infile="${SKYFOV}[ccd_id=${i}][cols POS]" opt="data,clean" | awk '{for (i=1;i<=NF;i++) print $i }' |sed -e ':a;N;s/\n/,/;ta' | awk -F"]," '{print "polygon("$2}' | awk -F"NaN" '{print $1}' >${TMP_LIST}
+        python ${SCRIPT_DIR}/${CCDGAP_SCRIPT} ${TMP_LIST} >${TMP_REC}
+        XC=`cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $1}'`
+        YC=`cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $2}'`
+        ADD_L=`cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $3/2}'`
+        ANG=`cat ${TMP_REC} | awk -F\( '{print $2}' |awk -F\) '{print $1}' |awk -F\, '{print $5}'`
+        if [ `echo "${ANG} < 0" | bc -l` -eq 1 ] ; then
+            ANG=`echo "${ANG} + 90" | bc -l`
+        elif [ `echo "${ANG} >=90" | bc -l` -eq 1 ] ; then
+            ANG=`echo "${ANG} - 90" | bc -l`
         fi
+        ANG=`echo "${ANG}/180*3.1415926" | bc -l`
+        CCD_1_X_RAW=`echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)-$2*sin($3)}' `
+        CCD_2_X_RAW=`echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)-$2*sin($3)}' `
+        CCD_3_X_RAW=`echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)+$2*sin($3)}' `
+        CCD_4_X_RAW=`echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)+$2*sin($3)}' `
+        CCD_1_Y_RAW=`echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)-$2*sin($3)}' `
+        CCD_2_Y_RAW=`echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)+$2*sin($3)}' `
+        CCD_3_Y_RAW=`echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)-$2*sin($3)}' `
+        CCD_4_Y_RAW=`echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)+$2*sin($3)}' `
+        CCD_1_RAW=`echo "${CCD_1_X_RAW},${CCD_1_Y_RAW}"`
+        CCD_2_RAW=`echo "${CCD_2_X_RAW},${CCD_2_Y_RAW}"`
+        CCD_3_RAW=`echo "${CCD_3_X_RAW},${CCD_3_Y_RAW}"`
+        CCD_4_RAW=`echo "${CCD_4_X_RAW},${CCD_4_Y_RAW}"`
+        REG_CCD_RAW=`echo "polygon(${CCD_1_RAW},${CCD_2_RAW},${CCD_4_RAW},${CCD_3_RAW}) "`
+        DX_2T1=$(echo "`echo ${CCD_2_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_1_RAW} |awk -F\, '{print $1}'`" |bc -l)
+        DY_2T1=$(echo "`echo ${CCD_2_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_1_RAW} |awk -F\, '{print $2}'`" |bc -l)
+        CCD_WIDTH=`echo "sqrt(${DX_2T1}*${DX_2T1}+${DY_2T1}*${DY_2T1})" | bc -l`
+        CCD_2T1_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_2T1}/${CCD_WIDTH}" | bc -l`
+        CCD_2T1_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_2T1}/${CCD_WIDTH}" | bc -l`
+        DX_3T1=$(echo "`echo ${CCD_3_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_1_RAW} |awk -F\, '{print $1}'`" |bc -l)
+        DY_3T1=$(echo "`echo ${CCD_3_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_1_RAW} |awk -F\, '{print $2}'`" |bc -l)
+        CCD_3T1_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_3T1}/${CCD_WIDTH}" | bc -l`
+        CCD_3T1_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_3T1}/${CCD_WIDTH}" | bc -l`
+        CCD_1_X=$(echo "` echo ${CCD_1_RAW} |awk -F\, '{print $1}' ` + `echo ${CCD_2T1_MOV_X}` +`echo ${CCD_3T1_MOV_X}` "| bc -l)
+        CCD_1_Y=$(echo "` echo ${CCD_1_RAW} |awk -F\, '{print $2}' ` + `echo ${CCD_2T1_MOV_Y}` +`echo ${CCD_3T1_MOV_Y}` "| bc -l)
+        DX_1T2=$(echo "`echo ${CCD_1_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_2_RAW} |awk -F\, '{print $1}'`" |bc -l)
+        DY_1T2=$(echo "`echo ${CCD_1_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_2_RAW} |awk -F\, '{print $2}'`" |bc -l)
+        CCD_1T2_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_1T2}/${CCD_WIDTH}" | bc -l`
+        CCD_1T2_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_1T2}/${CCD_WIDTH}" | bc -l`
+        DX_4T2=$(echo "`echo ${CCD_4_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_2_RAW} |awk -F\, '{print $1}'`" |bc -l)
+        DY_4T2=$(echo "`echo ${CCD_4_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_2_RAW} |awk -F\, '{print $2}'`" |bc -l)
+        CCD_4T2_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_4T2}/${CCD_WIDTH}" | bc -l`
+        CCD_4T2_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_4T2}/${CCD_WIDTH}" | bc -l`
+        CCD_2_X=$(echo "` echo ${CCD_2_RAW} |awk -F\, '{print $1}' ` + `echo ${CCD_1T2_MOV_X}` +`echo ${CCD_4T2_MOV_X}` "| bc -l)
+        CCD_2_Y=$(echo "` echo ${CCD_2_RAW} |awk -F\, '{print $2}' ` + `echo ${CCD_1T2_MOV_Y}` +`echo ${CCD_4T2_MOV_Y}` "| bc -l)
+        DX_1T3=$(echo "`echo ${CCD_1_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_3_RAW} |awk -F\, '{print $1}'`" |bc -l)
+        DY_1T3=$(echo "`echo ${CCD_1_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_3_RAW} |awk -F\, '{print $2}'`" |bc -l)
+        CCD_1T3_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_1T3}/${CCD_WIDTH}" | bc -l`
+        CCD_1T3_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_1T3}/${CCD_WIDTH}" | bc -l`
+        DX_4T3=$(echo "`echo ${CCD_4_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_3_RAW} |awk -F\, '{print $1}'`" |bc -l)
+        DY_4T3=$(echo "`echo ${CCD_4_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_3_RAW} |awk -F\, '{print $2}'`" |bc -l)
+        CCD_4T3_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_4T3}/${CCD_WIDTH}" | bc -l`
+        CCD_4T3_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_4T3}/${CCD_WIDTH}" | bc -l`
+        CCD_3_X=$(echo "` echo ${CCD_3_RAW} |awk -F\, '{print $1}' ` + `echo ${CCD_1T3_MOV_X}` +`echo ${CCD_4T3_MOV_X}` "| bc -l)
+        CCD_3_Y=$(echo "` echo ${CCD_3_RAW} |awk -F\, '{print $2}' ` + `echo ${CCD_1T3_MOV_Y}` +`echo ${CCD_4T3_MOV_Y}` "| bc -l)
+        DX_2T4=$(echo "`echo ${CCD_2_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_4_RAW} |awk -F\, '{print $1}'`" |bc -l)
+        DY_2T4=$(echo "`echo ${CCD_2_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_4_RAW} |awk -F\, '{print $2}'`" |bc -l)
+        CCD_2T4_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_2T4}/${CCD_WIDTH}" | bc -l`
+        CCD_2T4_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_2T4}/${CCD_WIDTH}" | bc -l`
+        DX_3T4=$(echo "`echo ${CCD_3_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_4_RAW} |awk -F\, '{print $1}'`" |bc -l)
+        DY_3T4=$(echo "`echo ${CCD_3_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_4_RAW} |awk -F\, '{print $2}'`" |bc -l)
+        CCD_3T4_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_3T4}/${CCD_WIDTH}" | bc -l`
+        CCD_3T4_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_3T4}/${CCD_WIDTH}" | bc -l`
+        CCD_4_X=$(echo "` echo ${CCD_4_RAW} |awk -F\, '{print $1}' ` + `echo ${CCD_2T4_MOV_X}` +`echo ${CCD_3T4_MOV_X}` "| bc -l)
+        CCD_4_Y=$(echo "` echo ${CCD_4_RAW} |awk -F\, '{print $2}' ` + `echo ${CCD_2T4_MOV_Y}` +`echo ${CCD_3T4_MOV_Y}` "| bc -l)
+        REG_CCD_CUT=`echo "polygon(${CCD_1_X},${CCD_1_Y},${CCD_2_X},${CCD_2_Y},${CCD_4_X},${CCD_4_Y},${CCD_3_X},${CCD_3_Y})"`
+        echo ${REG_CCD_CUT} >>${TMP_REG_FILE_CCD}
     done
-    ANG=`echo "${ANG}/180*3.1415926" |bc -l`
-    CCD_1_X_RAW=` echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)-$2*sin($3)}' `
-    CCD_2_X_RAW=` echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)-$2*sin($3)}' `
-    CCD_3_X_RAW=` echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)+$2*sin($3)}' `
-    CCD_4_X_RAW=` echo " ${XC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)+$2*sin($3)}' `
-    CCD_1_Y_RAW=` echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)-$2*sin($3)}' `
-    CCD_2_Y_RAW=` echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1+$2*cos($3)+$2*sin($3)}' `
-    CCD_3_Y_RAW=` echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)-$2*sin($3)}' `
-    CCD_4_Y_RAW=` echo " ${YC}  ${ADD_L} ${ANG} "| awk '{print $1-$2*cos($3)+$2*sin($3)}' `
-    CCD_1_RAW=` echo "${CCD_1_X_RAW},${CCD_1_Y_RAW}"`
-    CCD_2_RAW=` echo "${CCD_2_X_RAW},${CCD_2_Y_RAW}"`
-    CCD_3_RAW=` echo "${CCD_3_X_RAW},${CCD_3_Y_RAW}"`
-    CCD_4_RAW=` echo "${CCD_4_X_RAW},${CCD_4_Y_RAW}"`
-    REG_CCD_RAW=`echo "polygon(${CCD_1_RAW}, ${CCD_2_RAW}, ${CCD_4_RAW}, ${CCD_3_RAW}) " `
-    DX_2T1=$(echo "`echo ${CCD_2_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_1_RAW} |awk -F\, '{print $1}'`" |bc -l)
-    DY_2T1=$(echo "`echo ${CCD_2_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_1_RAW} |awk -F\, '{print $2}'`" |bc -l)
-    CCD_WIDTH=`echo "sqrt(${DX_2T1}*${DX_2T1}+${DY_2T1}*${DY_2T1})" | bc -l`
-    CCD_2T1_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_2T1}/${CCD_WIDTH}" | bc -l`
-    CCD_2T1_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_2T1}/${CCD_WIDTH}" | bc -l`
-    DX_3T1=$(echo "`echo ${CCD_3_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_1_RAW} |awk -F\, '{print $1}'`" |bc -l)
-    DY_3T1=$(echo "`echo ${CCD_3_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_1_RAW} |awk -F\, '{print $2}'`" |bc -l)
-    CCD_3T1_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_3T1}/${CCD_WIDTH}" | bc -l`
-    CCD_3T1_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_3T1}/${CCD_WIDTH}" | bc -l`
-    CCD_1_X=$(echo "` echo ${CCD_1_RAW} |awk -F\, '{print $1}' ` + `echo ${CCD_2T1_MOV_X}` +`echo ${CCD_3T1_MOV_X}` "| bc -l)
-    CCD_1_Y=$(echo "` echo ${CCD_1_RAW} |awk -F\, '{print $2}' ` + `echo ${CCD_2T1_MOV_Y}` +`echo ${CCD_3T1_MOV_Y}` "| bc -l)
-    DX_1T2=$(echo "`echo ${CCD_1_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_2_RAW} |awk -F\, '{print $1}'`" |bc -l)
-    DY_1T2=$(echo "`echo ${CCD_1_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_2_RAW} |awk -F\, '{print $2}'`" |bc -l)
-    CCD_1T2_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_1T2}/${CCD_WIDTH}" | bc -l`
-    CCD_1T2_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_1T2}/${CCD_WIDTH}" | bc -l`
-    DX_4T2=$(echo "`echo ${CCD_4_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_2_RAW} |awk -F\, '{print $1}'`" |bc -l)
-    DY_4T2=$(echo "`echo ${CCD_4_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_2_RAW} |awk -F\, '{print $2}'`" |bc -l)
-    CCD_4T2_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_4T2}/${CCD_WIDTH}" | bc -l`
-    CCD_4T2_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_4T2}/${CCD_WIDTH}" | bc -l`
-    CCD_2_X=$(echo "` echo ${CCD_2_RAW} |awk -F\, '{print $1}' ` + `echo ${CCD_1T2_MOV_X}` +`echo ${CCD_4T2_MOV_X}` "| bc -l)
-    CCD_2_Y=$(echo "` echo ${CCD_2_RAW} |awk -F\, '{print $2}' ` + `echo ${CCD_1T2_MOV_Y}` +`echo ${CCD_4T2_MOV_Y}` "| bc -l)
-    DX_1T3=$(echo "`echo ${CCD_1_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_3_RAW} |awk -F\, '{print $1}'`" |bc -l)
-    DY_1T3=$(echo "`echo ${CCD_1_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_3_RAW} |awk -F\, '{print $2}'`" |bc -l)
-    CCD_1T3_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_1T3}/${CCD_WIDTH}" | bc -l`
-    CCD_1T3_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_1T3}/${CCD_WIDTH}" | bc -l`
-    DX_4T3=$(echo "`echo ${CCD_4_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_3_RAW} |awk -F\, '{print $1}'`" |bc -l)
-    DY_4T3=$(echo "`echo ${CCD_4_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_3_RAW} |awk -F\, '{print $2}'`" |bc -l)
-    CCD_4T3_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_4T3}/${CCD_WIDTH}" | bc -l`
-    CCD_4T3_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_4T3}/${CCD_WIDTH}" | bc -l`
-    CCD_3_X=$(echo "` echo ${CCD_3_RAW} |awk -F\, '{print $1}' ` + `echo ${CCD_1T3_MOV_X}` +`echo ${CCD_4T3_MOV_X}` "| bc -l)
-    CCD_3_Y=$(echo "` echo ${CCD_3_RAW} |awk -F\, '{print $2}' ` + `echo ${CCD_1T3_MOV_Y}` +`echo ${CCD_4T3_MOV_Y}` "| bc -l)
-    DX_2T4=$(echo "`echo ${CCD_2_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_4_RAW} |awk -F\, '{print $1}'`" |bc -l)
-    DY_2T4=$(echo "`echo ${CCD_2_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_4_RAW} |awk -F\, '{print $2}'`" |bc -l)
-    CCD_2T4_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_2T4}/${CCD_WIDTH}" | bc -l`
-    CCD_2T4_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_2T4}/${CCD_WIDTH}" | bc -l`
-    DX_3T4=$(echo "`echo ${CCD_3_RAW} | awk -F\, '{print $1}'`-`echo  ${CCD_4_RAW} |awk -F\, '{print $1}'`" |bc -l)
-    DY_3T4=$(echo "`echo ${CCD_3_RAW} | awk -F\, '{print $2}'`-`echo  ${CCD_4_RAW} |awk -F\, '{print $2}'`" |bc -l)
-    CCD_3T4_MOV_X=`echo "${DFT_CCD_EDGECUT}*${DX_3T4}/${CCD_WIDTH}" | bc -l`
-    CCD_3T4_MOV_Y=`echo "${DFT_CCD_EDGECUT}*${DY_3T4}/${CCD_WIDTH}" | bc -l`
-    CCD_4_X=$(echo "` echo ${CCD_4_RAW} |awk -F\, '{print $1}' ` + `echo ${CCD_2T4_MOV_X}` +`echo ${CCD_3T4_MOV_X}` "| bc -l)
-    CCD_4_Y=$(echo "` echo ${CCD_4_RAW} |awk -F\, '{print $2}' ` + `echo ${CCD_2T4_MOV_Y}` +`echo ${CCD_3T4_MOV_Y}` "| bc -l)
-    REG_CCD_CUT=`echo "polygon(${CCD_1_X},${CCD_1_Y},${CCD_2_X},${CCD_2_Y},${CCD_4_X},${CCD_4_Y},${CCD_3_X},${CCD_3_Y})"`
-    echo ${REG_CCD_CUT} >>${TMP_REG_FILE_CCD}
-done
     REG_FILE_CCD="_ccd.reg"
     [ -e "${REG_FILE_CCD}" ] && mv -fv ${REG_FILE_CCD} ${REG_FILE_CCD}_bak
  #   echo "` cat ${TMP_REG_FILE_CCD} | head -n 1 | tail -n 1` + ` cat ${TMP_REG_FILE_CCD} | head -n 2 | tail -n 1` +`cat ${TMP_REG_FILE_CCD} | head -n 3 | tail -n 1`+`cat ${TMP_REG_FILE_CCD} | head -n 4 | tail -n 1 `" >${REG_FILE_CCD}
@@ -515,41 +496,14 @@ punlearn dmlist
 dmlist infile="${SBP_RMID}[cols RMID,R_ERR,SUR_FLUX,SUR_FLUX_ERR]" \
     outfile="${SBP_TXT}" opt="data,clean"
 
-## QDP for sbp {{{
 printf "generate a handy QDP file for sbp ...\n"
-cp -fv ${SBP_TXT} ${SBP_QDP}
-# change comment sign
-sed -i'' 's/#/!/g' ${SBP_QDP}
-# add QDP commands
-sed -i'' '1 i\
-SKIP SINGLE' ${SBP_QDP}
-sed -i'' '1 i\
-READ SERR 1 2' ${SBP_QDP}
-sed -i'' '2 i\
-LABEL Y "Surface Flux (photons/cm\\u2\\d/pixel\\u2\\d/s)"' ${SBP_QDP}
-sed -i'' '2 i\
-LABEL X "Radius (pixel)"' ${SBP_QDP}
-sed -i'' '2 i\
-LABEL T "Surface Brightness Profile"' ${SBP_QDP}
-## QDP }}}
-
-printf "generate sbp fitting needed files ...\n"
-SBP_RADIUS="radius_sbp.txt"
-SBP_FLUX="flux_sbp.txt"
-[ -e "${SBP_RADIUS}" ] && mv -fv ${SBP_RADIUS} ${SBP_RADIUS}_bak
-[ -e "${SBP_FLUX}" ] && mv -fv ${SBP_FLUX} ${SBP_FLUX}_bak
-punlearn dmlist
-dmlist infile="${SBP_RMID}[cols R]" \
-    opt="data,clean" | awk '{ print $2 }' > ${SBP_RADIUS}
-# change the first line `R[2]' to `0.0'
-sed -i'' 's/R.*/0\.0/' ${SBP_RADIUS}
-dmlist infile="${SBP_RMID}[cols SUR_FLUX,SUR_FLUX_ERR]" \
-    opt="data,clean" > ${SBP_FLUX}
-# remove the first comment line
-sed -i'' '/#.*/d' ${SBP_FLUX}
+echo "\
+READ SERR 1 2
+SKIP SINGLE
+LABEL T \"Surface Brightness Profile\"
+LABEL X \"Radius (pixel)\"
+LABEL Y \"Surface Flux (photons/cm\\u2\\d/pixel\\u2\\d/s)\"
+!$(cat ${SBP_TXT})" > ${SBP_QDP}  # NOTE the '!' which comments the header
 
 ## sbp data }}}
-
 ## main }}}
-
-exit 0