Superseded by 'scripts/ciao_calc_csb.sh'

1 files changed, 0 insertions, 160 deletions

diff --git a/mass_profile/csb_calc_lwt.sh b/mass_profile/csb_calc_lwt.sh
deleted file mode 100755
index fa6e654..0000000
--- a/mass_profile/csb_calc_lwt.sh
+++ /dev/null
@@ -1,160 +0,0 @@
-#!/bin/sh
-#
-# Calculate the surface brightness concentration (C_SB) for objects
-# meet the conditions: 'z>0.3' or 'counts_in_0.048R500<500'.
-#
-# Aaron LI
-# Updated: 2016-05-28
-#
-
-ERR_CALC=1
-ERR_DIR=2
-ERR_JSON=3
-ERR_Z=4
-ERR_CNT=5
-
-# default basedir relative to 'spc/profile'
-DFT_BASEDIR="../.."
-# default imgdir relative to 'basedir'
-DFT_IMGDIR="img"
-DFT_JSON_PAT="*_INFO.json"
-RSPEC_REG="rspec.reg"
-CSB_RES="csb_results_`date '+%Y%m%d'`.txt"
-
-## functions {{{
-# process commandline arguments
-# cmdline arg format: `KEY=VALUE'
-getopt_keyval() {
-    until [ -z "$1" ]
-    do
-        key=${1%%=*}                    # extract key
-        val=${1#*=}                     # extract value
-        keyval="${key}=\"${val}\""
-        echo "## getopt: eval '${keyval}'"
-        eval ${keyval}
-        shift                           # shift, process next one
-    done
-}
-## functions }}}
-
-## parameters {{{
-# process cmdline args using `getopt_keyval'
-getopt_keyval "$@"
-
-# basedir
-if [ -d "${basedir}" ] && ls ${basedir}/*repro_evt2.fits > /dev/null 2>&1; then
-    BASEDIR=${basedir}
-elif [ -d "${DFT_BASEDIR}" ] && ls ${DFT_BASEDIR}/*repro_evt2.fits > /dev/null 2>&1; then
-    BASEDIR=${DFT_BASEDIR}
-else
-    read -p "> basedir (contains info json): " BASEDIR
-    if [ ! -d "${BASEDIR}" ] || ! ls ${BASEDIR}/*repro_evt2.fits >/dev/null 2>&1; then
-        printf "ERROR: given \`${BASEDIR}' invalid!\n"
-        exit ${ERR_DIR}
-    fi
-fi
-BASEDIR=`( cd ${BASEDIR} && pwd -P )`
-printf "## use basedir: \`${BASEDIR}'\n"
-# mass dir
-if [ ! -z "${imgdir}" ] && [ -d "${BASEDIR}/${imgdir}" ]; then
-    IMG_DIR=`( cd ${BASEDIR}/${imgdir} && pwd -P )`
-elif [ -d "${BASEDIR}/${DFT_IMGDIR}" ]; then
-    IMG_DIR=`( cd ${BASEDIR}/${DFT_IMGDIR} && pwd -P )`
-else
-    read -p "> img dir (relative to basedir): " IMG_DIR
-    if [ ! -d "${BASEDIR}/${IMG_DIR}" ]; then
-        printf "ERROR: given \`${IMG_DIR}' invalid\n"
-        exit ${ERR_DIR}
-    else
-        IMG_DIR="${BASEDIR}/${IMG_DIR}"
-    fi
-fi
-printf "## use imgdir: \`${IMG_DIR}'\n"
-# info json
-if [ ! -z "${json}" ] && [ -r "${BASEDIR}/${json}" ]; then
-    JSON_FILE="${BASEDIR}/${json}"
-elif [ `ls -1 ${BASEDIR}/${DFT_JSON_PAT} 2>/dev/null | wc -l` -eq 1 ]; then
-    JSON_FILE="`ls ${BASEDIR}/${DFT_JSON_PAT} 2>/dev/null`"
-else
-    read -p "> info json: " JSON_FILE
-    if [ ! -r "${BASEDIR}/${JSON_FILE}" ]; then
-        printf "ERROR: given \`${JSON_FILE}' not exist!\n"
-        exit ${ERR_JSON}
-    fi
-fi
-printf "## use json_file: \`${JSON_FILE}'\n"
-## }}}
-
-## main {{{
-# in 'spc/profile'
-X=`grep -iE '(pie|annulus)' ${RSPEC_REG} | head -n 1 | awk -F'(' '{ print $2 }' | awk -F',' '{ print $1 }'`
-Y=`grep -iE '(pie|annulus)' ${RSPEC_REG} | head -n 1 | awk -F'(' '{ print $2 }' | awk -F',' '{ print $2 }'`
-# json file
-Z=`grep -i '"redshift"' ${JSON_FILE} | awk -F':' '{ print $2 }' | tr -d ' ,'`
-R500=`grep '"R500.*kpc' ${JSON_FILE} | awk -F':' '{ print $2 }' | tr -d ' ,'`
-OBS_ID=`grep '"Obs.*ID' ${JSON_FILE} | awk -F':' '{ print $2 }' | tr -d ' ,'`
-OBJ_NAME=`grep '"Source\ Name' ${JSON_FILE} | awk -F':' '{ print $2 }' | sed -e 's/\ *"//' -e 's/"\ *,$//'`
-CT=`grep '"Cooling_time' ${JSON_FILE} | awk -F':' '{ print $2 }' | tr -d ' ,'`
-
-cd ${IMG_DIR}
-printf "entered img directory\n"
-EXPMAP=`ls expmap*.fits 2> /dev/null`
-EVT_E=`ls evt*e700-7000*.fits 2> /dev/null` 
-
-### test Z>0.3?
-if [ `echo "${Z} < 0.3" | bc -l` -eq 1 ]; then
-    F_WZ=true
-    WZ="WZ"
-    printf "*** WARNING: redshift z=${Z} < 0.3 ***\n"
-#    exit ${ERR_Z}
-fi
-
-KPC_PER_PIXEL=`cosmo_calc ${Z} | grep 'kpc/pixel' | awk '{ print $3 }'`
-RC_PIX=`echo "scale=2; 0.048 * ${R500} / ${KPC_PER_PIXEL}" | bc -l`
-# test counts_in_0.048R500<500?
-RC_REG="pie(${X},${Y},0,${RC_PIX},0,360)"
-punlearn dmlist
-CNT_RC=`dmlist infile="${EVT_E}[sky=${RC_REG}]" opt=block | grep 'EVENTS' | awk '{ print $8 }'`
-printf "R500=${R500}, 0.048R500_pix=${RC_PIX}, counts_in_0.048R500=${CNT_RC}\n"
-if [ ${CNT_RC} -gt 500 ]; then
-    F_WC=true
-    WC="WC"
-    printf "*** WARNING: counts_in_0.048R500=${CNT_RC} > 500 ***\n"
-#    exit ${ERR_CNT}
-fi
-
-TMP_REG="_tmp_csb.reg"
-TMP_S="_tmp_csb.fits"
-
-R1=`echo "scale=2;  40 / ${KPC_PER_PIXEL}" | bc -l`
-R2=`echo "scale=2; 400 / ${KPC_PER_PIXEL}" | bc -l` 
-cat > ${TMP_REG} << _EOF_
-pie(${X},${Y},0,${R1},0,360)
-pie(${X},${Y},0,${R2},0,360)
-_EOF_
-
-printf "CHECK the regions (R1=${R1}, R2=${R2}) ...\n"
-ds9 ${EVT_E} -regions ${TMP_REG} -cmap sls -bin factor 4
-
-punlearn dmextract
-dmextract infile="${EVT_E}[bin sky=@${TMP_REG}]" outfile="${TMP_S}" exp=${EXPMAP} opt=generic clobber=yes
-punlearn dmlist
-S1=`dmlist "${TMP_S}[cols SUR_FLUX]" opt="data,clean" | grep -v '#' | sed -n -e 's/\ *//' -e '1p'`
-S2=`dmlist "${TMP_S}[cols SUR_FLUX]" opt="data,clean" | grep -v '#' | sed -n -e 's/\ *//' -e '2p'`
-CSB=`echo "${S1} ${S2}" | awk '{ print $1/$2/100 }'` 
-
-[ -e ${CSB_RES} ] && mv -f ${CSB_RES} ${CSB_RES}_bak
-printf "\n==============================\n"
-printf "z=${Z}, R500=${R500} (kpc)\n" | tee -a ${CSB_RES}
-printf "0.048R500=${RC_PIX}, counts=${CNT_RC}\n" | tee -a ${CSB_RES}
-printf "R1=${R1}, R2=${R2} (pixel)\n" | tee -a ${CSB_RES}
-printf "S1=${S1}, S2=${S2} (sur_flux)\n" | tee -a ${CSB_RES}
-printf "C_sb: ${CSB}\n" | tee -a ${CSB_RES}
-[ "x${F_WZ}" = "xtrue" ] && printf "${WZ}\n" | tee -a ${CSB_RES}
-[ "x${F_WC}" = "xtrue" ] && printf "${WC}\n" | tee -a ${CSB_RES}
-printf "# OBS_ID,OBJ_NAME,Z,R500,RC_PIX,CNT_RC,CT,R1_PIX,R2_PIX,S1,S2,CSB,WRDSHFT,WCNT\n" | tee -a ${CSB_RES}
-printf "# $OBS_ID,$OBJ_NAME,$Z,$R500,$RC_PIX,$CNT_RC,$CT,$R1,$R2,$S1,$S2,$CSB,$WZ,$WC\n\n" | tee -a ${CSB_RES}
-## main }}}
-
-exit 0
-