diff options
| -rwxr-xr-x | mass_profile/csb_calc_lwt.sh | 160 |
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 - |