Rename 3 scripts; Clean up HOWTO a bit

1 files changed, 630 insertions, 0 deletions

diff --git a/scripts/collect_infodata.sh b/scripts/collect_infodata.sh
new file mode 100755
index 0000000..d78d1da
--- /dev/null
+++ b/scripts/collect_infodata.sh
@@ -0,0 +1,630 @@
+#!/bin/sh
+##
+unalias -a
+export LC_COLLATE=C
+# fix path for python
+export PATH=/usr/bin:/usr/local/bin:$PATH
+#########################################################################
+## Collect the calculated mass data and write to the JSON file.
+##
+## Weitian LI <liweitianux@gmail.com>
+## August 31, 2012
+##
+## ChangeLogs:
+##   check JSON syntax, modify output to agree with the syntax
+##   Ref: http://json.parser.online.fr/
+## v1.1, 2012/09/05, Weitian LI
+##   add `T_avg(0.2-0.5 R500)' and `T_err'
+## v2.0, 2012/10/14, Weitian LI
+##   add parameters
+## v2.1, 2012/11/07, Weitian LI
+##   account for `fitting_dbeta'
+## v2.2, 2012/12/18, Weitian LI
+##   add `beta' and `cooling time' parameters
+## v3.0, 2013/02/09, Weitian LI
+##   modified for `new sample info format'
+## v3.1, 2013/05/18, Weitian LI
+##   add key `Feature'
+## v3.2, 2013/05/29, Weitian LI
+##   add key `XCNTRD_RA, XCNTRD_DEC'
+## v3.3, 2013/10/14, Weitian LI
+##   add key `Unified Name'
+## v3.4, 2015/06/03, Aaron LI
+##   * Copy needed pfiles to current working directory, and
+##     set environment variable $PFILES to use these first.
+##   * Replaced 'grep' with '\grep', 'ls' with '\ls'
+##
+VERSION="v3.4"
+UPDATED="2015/06/03"
+
+## error code {{{
+ERR_USG=1
+ERR_DIR=11
+ERR_CFG=12
+ERR_RES=13
+ERR_COSC=14
+ERR_BETA=101
+ERR_BETA2=102
+ERR_JSON=201
+ERR_MLOG=202
+ERR_BLOG=203
+ERR_CLOG=204
+## error code }}}
+
+## usage, help {{{
+case "$1" in
+    -[hH]*|--[hH]*)
+        printf "usage:\n"
+        printf "    `basename $0` json=<info_json> cfg=<main_cfg> res=<final_result> basedir=<repro_dir> massdir=<mass_dir>\n"
+        printf "\nversion:\n"
+        printf "${VERSION}, ${UPDATED}\n"
+        exit ${ERR_USG}
+        ;;
+esac
+## usage, help }}}
+
+## default parameters {{{
+# default basedir
+DFT_BASEDIR=".."
+# default mass dir
+DFT_MASSDIR="mass"
+# default pattern for json info file
+DFT_JSON_PAT="*_INFO.json"
+# main config file
+DFT_CFG_PAT="global.cfg"
+# default result file
+DFT_RES_FINAL_PAT="final_result.txt"
+# default sbprofile region file
+DFT_SBPROFILE_REG="sbprofile.reg"
+# default radial spectra region file
+DFT_RSPEC_REG="rspec.reg"
+# default CMD to determine the `spc/profile' dir
+DFT_SPC_DIR_CMD='dirname `readlink ${T_DATA_FILE}`'
+# default CMD to determine the `img' dir
+DFT_IMG_DIR_CMD='dirname `readlink ${RADIUS_SBP_FILE}`'
+# default config file pattern for `expcorr'
+DFT_EXPCORR_CONF_PAT="*_expcorr.conf"
+## default parameters }}}
+
+## 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 "$@"
+
+# init directory
+INIT_DIR=`pwd -P`
+
+# check given parameters
+# check given dir
+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}" ]; then
+        printf "ERROR: given \`${BASEDIR}' NOT a directory\n"
+        exit ${ERR_DIR}
+    elif ! \ls ${BASEDIR}/*repro_evt2.fits > /dev/null 2>&1; then
+        printf "ERROR: given \`${BASEDIR}' NOT contains needed evt files\n"
+        exit ${ERR_DIR}
+    fi
+fi
+# remove the trailing '/'
+BASEDIR=`( cd ${INIT_DIR}/${BASEDIR} && pwd -P )`
+printf "## use basedir: \`${BASEDIR}'\n"
+# mass dir
+if [ ! -z "${massdir}" ] && [ -d "${BASEDIR}/${massdir}" ]; then
+    MASS_DIR=`( cd ${BASEDIR}/${massdir} && pwd -P )`
+elif [ -d "${BASEDIR}/${DFT_MASSDIR}" ]; then
+    MASS_DIR=`( cd ${BASEDIR}/${DFT_MASSDIR} && pwd -P )`
+else
+    read -p "> mass dir (relative to basedir): " MASS_DIR
+    if [ ! -d "${BASEDIR}/${MASS_DIR}" ]; then
+        printf "ERROR: given \`${BASEDIR}/${MASS_DIR}' NOT a directory\n"
+        exit ${ERR_DIR}
+    fi
+fi
+# remove the trailing '/'
+MASS_DIR=`echo ${MASS_DIR} | sed 's/\/*$//'`
+printf "## use massdir: \`${MASS_DIR}'\n"
+# check json file
+if [ ! -z "${json}" ] && [ -r "${BASEDIR}/${json}" ]; then
+    JSON_FILE="${json}"
+elif [ "`\ls ${BASEDIR}/${DFT_JSON_PAT} | wc -l`" -eq 1 ]; then
+    JSON_FILE=`( cd ${BASEDIR} && \ls ${DFT_JSON_PAT} )`
+else
+    read -p "> info json file: " JSON_FILE
+    if ! [ -r "${BASEDIR}/${JSON_FILE}" ]; then
+        printf "ERROR: cannot access given \`${BASEDIR}/${INFO_JSON}' file\n"
+        exit ${ERR_JSON}
+    fi
+fi
+printf "## use info json file: \`${JSON_FILE}'\n"
+# check main config file
+if [ ! -z "${cfg}" ] && [ -r "${cfg}" ]; then
+    CFG_FILE="${cfg}"
+elif [ -r "${DFT_CFG_PAT}" ]; then
+    CFG_FILE="${DFT_CFG_PAT}"
+else
+    read -p "> main config file: " CFG_FILE
+    if [ ! -r "${CFG_FILE}" ]; then
+        printf "ERROR: cannot access given \`${CFG_JSON}' file\n"
+        exit ${ERR_CFG}
+    fi
+fi
+printf "## use main config file: \`${CFG_FILE}'\n"
+SBP_CFG=`\grep '^sbp_cfg' ${CFG_FILE} | awk '{ print $2 }'`
+printf "## sbp config file: \`${SBP_CFG}'\n"
+# check final result file
+if [ ! -z "${res}" ] && [ -r "${res}" ]; then
+    RES_FINAL="${res}"
+elif [ -r "${DFT_RES_FINAL_PAT}" ]; then
+    RES_FINAL="${DFT_RES_FINAL_PAT}"
+else
+    read -p "> final result file: " RES_FINAL
+    if [ ! -r "${RES_FINAL}" ]; then
+        printf "ERROR: cannot access given \`${RES_FINAL}' file\n"
+        exit ${ERR_RES}
+    fi
+fi
+printf "## use final result file: \`${RES_FINAL}'\n"
+## parameters }}}
+
+## directory & file {{{
+BASE_PATH=`dirname $0`
+printf "## BASE_PATH: ${BASE_PATH}\n"
+COSMO_CALC="${BASE_PATH}/cosmo_calc"
+if [ ! -x "${COSMO_CALC}" ]; then
+    printf "*** ERROR: \`${COSMO_CALC}' not executable\n"
+    exit ${ERR_COSC}
+fi
+EVT_DIR="${BASEDIR}/evt"
+IMG_DIR="${BASEDIR}/img"
+SPEC_DIR="${BASEDIR}/spc/profile"
+## dir & file }}}
+
+## prepare parameter files (pfiles) {{{
+CIAO_TOOLS="dmkeypar"
+
+# Copy necessary pfiles for localized usage
+for tool in ${CIAO_TOOLS}; do
+    pfile=`paccess ${tool}`
+    [ -n "${pfile}" ] && punlearn ${tool} && cp -Lvf ${pfile} .
+done
+
+# Modify environment variable 'PFILES' to use local pfiles first
+export PFILES="./:${PFILES}"
+## pfiles }}}
+
+cd ${BASEDIR}
+printf "## enter directory: `pwd -P`\n"
+
+## in dir `repro' {{{
+punlearn dmkeypar
+EVT_RAW=`\ls *repro_evt2.fits`
+OBS_ID=`dmkeypar ${EVT_RAW} OBS_ID echo=yes`
+DATE_OBS=`dmkeypar ${EVT_RAW} DATE-OBS echo=yes`
+EXPOSURE_RAW=`dmkeypar ${EVT_RAW} EXPOSURE echo=yes | awk '{ print $1/1000 }'`
+## ACIS_TYPE
+DETNAM=`dmkeypar ${EVT_RAW} DETNAM echo=yes`
+if echo ${DETNAM} | \grep -q 'ACIS-0123'; then
+    ACIS_TYPE="ACIS-I"
+elif echo ${DETNAM} | \grep -q 'ACIS-[0-6]*7'; then
+    ACIS_TYPE="ACIS-S"
+else
+    printf "*** ERROR: unknown detector type: ${DETNAM}\n"
+    ACIS_TYPE="UNKNOWN"
+fi
+## dir `repro' }}}
+
+## in dir `repro/evt' {{{
+cd ${EVT_DIR}
+EVT_CLEAN=`\ls evt2_c*_clean.fits`
+EXPOSURE_CLEAN=`dmkeypar ${EVT_CLEAN} EXPOSURE echo=yes | awk '{ print $1/1000 }'`
+## dir `repro/evt' }}}
+
+## in dir `repro/mass' {{{
+cd ${MASS_DIR}
+
+# misc {{{
+N_H=`\grep '^nh' ${CFG_FILE} | awk '{ print $2 }'`
+Z=`\grep '^z' ${SBP_CFG} | awk '{ print $2 }'`
+T_DATA_FILE=`\grep '^t_data_file' ${CFG_FILE} | awk '{ print $2 }'`
+NFW_RMIN_KPC=`\grep '^nfw_rmin_kpc' ${CFG_FILE} | awk '{ print $2 }'`
+E_Z=`${COSMO_CALC} ${Z} | \grep -i 'Hubble_parameter' | awk '{ print $3 }'`
+KPC_PER_PIXEL=`${COSMO_CALC} ${Z} | \grep 'kpc/pixel' | awk '{ print $3 }'`
+RADIUS_SBP_FILE=`\grep '^radius_sbp_file' ${CFG_FILE} | awk '{ print $2 }'`
+RMAX_SBP_PIX=`tail -n 1 ${RADIUS_SBP_FILE} | awk '{ print $1+$2 }'`
+RMAX_SBP_KPC=`echo "${RMAX_SBP_PIX} ${KPC_PER_PIXEL}" | awk '{ printf("%.2f", $1*$2) }'`
+SPC_DIR=`eval ${DFT_SPC_DIR_CMD}`
+if [ -f "${SPC_DIR}/${DFT_RSPEC_REG}" ]; then
+    RMAX_TPRO_PIX=`\grep -iE '(pie|annulus)' ${SPC_DIR}/${DFT_RSPEC_REG} | tail -n 1 | awk -F',' '{ print $4 }'`
+    RMAX_TPRO_KPC=`echo "${RMAX_TPRO_PIX} ${KPC_PER_PIXEL}" | awk '{ printf("%.2f", $1*$2) }'`
+fi
+IMG_DIR=`eval ${DFT_IMG_DIR_CMD}`
+EXPCORR_CONF=`\ls ${IMG_DIR}/${DFT_EXPCORR_CONF_PAT} 2> /dev/null`
+echo "EXPCORR_CONF: ${EXPCORR_CONF}"
+if [ -f "${EXPCORR_CONF}" ]; then
+    T_REF=`\grep '^temp' ${EXPCORR_CONF} | awk '{ print $2 }'`
+    Z_REF=`\grep '^abund' ${EXPCORR_CONF} | awk '{ print $2 }'`
+fi
+[ -z "${NFW_RMIN_KPC}" ] && NFW_RMIN_KPC="null"
+[ -z "${RMAX_SBP_PIX}" ] && RMAX_SBP_PIX="null"
+[ -z "${RMAX_SBP_KPC}" ] && RMAX_SBP_KPC="null"
+[ -z "${T_REF}" ]        && T_REF="null"
+[ -z "${Z_REF}" ]        && Z_REF="null"
+# misc }}}
+
+## determine single/double beta {{{
+if \grep -q '^beta2' ${SBP_CFG}; then
+    MODEL_SBP="double-beta"
+    n01=`\grep '^n01' ${RES_FINAL} | awk '{ print $2 }'`
+    beta1=`\grep '^beta1' ${RES_FINAL} | awk '{ print $2 }'`
+    rc1=`\grep -E '^rc1\s' ${RES_FINAL} | awk '{ print $2 }'`
+    rc1_kpc=`\grep '^rc1_kpc' ${RES_FINAL} | awk '{ print $2 }'`
+    n02=`\grep '^n02' ${RES_FINAL} | awk '{ print $2 }'`
+    beta2=`\grep '^beta2' ${RES_FINAL} | awk '{ print $2 }'`
+    rc2=`\grep -E '^rc2\s' ${RES_FINAL} | awk '{ print $2 }'`
+    rc2_kpc=`\grep '^rc2_kpc' ${RES_FINAL} | awk '{ print $2 }'`
+    BKG=`\grep '^bkg' ${RES_FINAL} | awk '{ print $2 }'`
+    # beta1 -> smaller rc; beta2 -> bigger rc
+    if [ `echo "${rc1} < ${rc2}" | bc -l` -eq 1 ]; then
+        N01=${n01}
+        BETA1=${beta1}
+        RC1=${rc1}
+        RC1_KPC=${rc1_kpc}
+        N02=${n02}
+        BETA2=${beta2}
+        RC2=${rc2}
+        RC2_KPC=${rc2_kpc}
+    else
+        # beta1 <-> beta2 (swap)
+        N01=${n02}
+        BETA1=${beta2}
+        RC1=${rc2}
+        RC1_KPC=${rc2_kpc}
+        N02=${n01}
+        BETA2=${beta1}
+        RC2=${rc1}
+        RC2_KPC=${rc1_kpc}
+    fi
+else
+    MODEL_SBP="single-beta"
+    N01="null"
+    BETA1="null"
+    RC1="null"
+    RC1_KPC="null"
+    N02=`\grep '^n0' ${RES_FINAL} | awk '{ print $2 }'`
+    BETA2=`\grep '^beta' ${RES_FINAL} | awk '{ print $2 }'`
+    RC2=`\grep -E '^rc\s' ${RES_FINAL} | awk '{ print $2 }'`
+    RC2_KPC=`\grep '^rc_kpc' ${RES_FINAL} | awk '{ print $2 }'`
+    BKG=`\grep '^bkg' ${RES_FINAL} | awk '{ print $2 }'`
+fi
+## single/double beta }}}
+
+## get `mass/virial_radius/luminosity' {{{
+# 200 data {{{
+R200_VAL=`\grep '^r200' ${RES_FINAL} | awk '{ print $2 }'`
+R200_ERR_L=`\grep '^r200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+R200_ERR_U=`\grep '^r200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+M200_VAL=`\grep '^m200' ${RES_FINAL} | awk '{ print $2 }'`
+M200_ERR_L=`\grep '^m200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+M200_ERR_U=`\grep '^m200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+L200_VAL=`\grep '^L200' ${RES_FINAL} | awk '{ print $2 }'`
+L200_ERR=`\grep '^L200' ${RES_FINAL} | awk '{ print $4 }'`
+MGAS200_VAL=`\grep '^gas_m200' ${RES_FINAL} | awk '{ print $2 }'`
+MGAS200_ERR_L=`\grep '^gas_m200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+MGAS200_ERR_U=`\grep '^gas_m200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+FGAS200_VAL=`\grep '^gas_fraction200' ${RES_FINAL} | awk '{ print $2 }'`
+FGAS200_ERR_L=`\grep '^gas_fraction200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+FGAS200_ERR_U=`\grep '^gas_fraction200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+[ -z "${R200_VAL}" ]       && R200_VAL="null"
+[ -z "${R200_ERR_L}" ]     && R200_ERR_L="null"
+[ -z "${R200_ERR_U}" ]     && R200_ERR_U="null"
+[ -z "${M200_VAL}" ]       && M200_VAL="null"
+[ -z "${M200_ERR_L}" ]     && M200_ERR_L="null"
+[ -z "${M200_ERR_U}" ]     && M200_ERR_U="null"
+[ -z "${L200_VAL}" ]       && L200_VAL="null"
+[ -z "${L200_ERR}" ]       && L200_ERR="null"
+[ -z "${MGAS200_VAL}" ]    && MGAS200_VAL="null"
+[ -z "${MGAS200_ERR_L}" ]  && MGAS200_ERR_L="null"
+[ -z "${MGAS200_ERR_U}" ]  && MGAS200_ERR_U="null"
+[ -z "${FGAS200_VAL}" ]    && FGAS200_VAL="null"
+[ -z "${FGAS200_ERR_L}" ]  && FGAS200_ERR_L="null"
+[ -z "${FGAS200_ERR_U}" ]  && FGAS200_ERR_U="null"
+# 200 }}}
+# 500 data {{{
+R500_VAL=`\grep '^r500' ${RES_FINAL} | awk '{ print $2 }'`
+R500_ERR_L=`\grep '^r500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+R500_ERR_U=`\grep '^r500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+M500_VAL=`\grep '^m500' ${RES_FINAL} | awk '{ print $2 }'`
+M500_ERR_L=`\grep '^m500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+M500_ERR_U=`\grep '^m500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+L500_VAL=`\grep '^L500' ${RES_FINAL} | awk '{ print $2 }'`
+L500_ERR=`\grep '^L500' ${RES_FINAL} | awk '{ print $4 }'`
+MGAS500_VAL=`\grep '^gas_m500' ${RES_FINAL} | awk '{ print $2 }'`
+MGAS500_ERR_L=`\grep '^gas_m500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+MGAS500_ERR_U=`\grep '^gas_m500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+FGAS500_VAL=`\grep '^gas_fraction500' ${RES_FINAL} | awk '{ print $2 }'`
+FGAS500_ERR_L=`\grep '^gas_fraction500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+FGAS500_ERR_U=`\grep '^gas_fraction500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+[ -z "${R500_VAL}" ]       && R500_VAL="null"
+[ -z "${R500_ERR_L}" ]     && R500_ERR_L="null"
+[ -z "${R500_ERR_U}" ]     && R500_ERR_U="null"
+[ -z "${M500_VAL}" ]       && M500_VAL="null"
+[ -z "${M500_ERR_L}" ]     && M500_ERR_L="null"
+[ -z "${M500_ERR_U}" ]     && M500_ERR_U="null"
+[ -z "${L500_VAL}" ]       && L500_VAL="null"
+[ -z "${L500_ERR}" ]       && L500_ERR="null"
+[ -z "${MGAS500_VAL}" ]    && MGAS500_VAL="null"
+[ -z "${MGAS500_ERR_L}" ]  && MGAS500_ERR_L="null"
+[ -z "${MGAS500_ERR_U}" ]  && MGAS500_ERR_U="null"
+[ -z "${FGAS500_VAL}" ]    && FGAS500_VAL="null"
+[ -z "${FGAS500_ERR_L}" ]  && FGAS500_ERR_L="null"
+[ -z "${FGAS500_ERR_U}" ]  && FGAS500_ERR_U="null"
+# 500 }}}
+# 1500 data {{{
+R1500_VAL=`\grep '^r1500' ${RES_FINAL} | awk '{ print $2 }'`
+R1500_ERR_L=`\grep '^r1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+R1500_ERR_U=`\grep '^r1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+M1500_VAL=`\grep '^m1500' ${RES_FINAL} | awk '{ print $2 }'`
+M1500_ERR_L=`\grep '^m1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+M1500_ERR_U=`\grep '^m1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+L1500_VAL=`\grep '^L1500' ${RES_FINAL} | awk '{ print $2 }'`
+L1500_ERR=`\grep '^L1500' ${RES_FINAL} | awk '{ print $4 }'`
+MGAS1500_VAL=`\grep '^gas_m1500' ${RES_FINAL} | awk '{ print $2 }'`
+MGAS1500_ERR_L=`\grep '^gas_m1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+MGAS1500_ERR_U=`\grep '^gas_m1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+FGAS1500_VAL=`\grep '^gas_fraction1500' ${RES_FINAL} | awk '{ print $2 }'`
+FGAS1500_ERR_L=`\grep '^gas_fraction1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+FGAS1500_ERR_U=`\grep '^gas_fraction1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+[ -z "${R1500_VAL}" ]      && R1500_VAL="null"
+[ -z "${R1500_ERR_L}" ]    && R1500_ERR_L="null"
+[ -z "${R1500_ERR_U}" ]    && R1500_ERR_U="null"
+[ -z "${M1500_VAL}" ]      && M1500_VAL="null"
+[ -z "${M1500_ERR_L}" ]    && M1500_ERR_L="null"
+[ -z "${M1500_ERR_U}" ]    && M1500_ERR_U="null"
+[ -z "${L1500_VAL}" ]      && L1500_VAL="null"
+[ -z "${L1500_ERR}" ]      && L1500_ERR="null"
+[ -z "${MGAS1500_VAL}" ]   && MGAS1500_VAL="null"
+[ -z "${MGAS1500_ERR_L}" ] && MGAS1500_ERR_L="null"
+[ -z "${MGAS1500_ERR_U}" ] && MGAS1500_ERR_U="null"
+[ -z "${FGAS1500_VAL}" ]   && FGAS1500_VAL="null"
+[ -z "${FGAS1500_ERR_L}" ] && FGAS1500_ERR_L="null"
+[ -z "${FGAS1500_ERR_U}" ] && FGAS1500_ERR_U="null"
+# 1500 }}}
+# 2500 data {{{
+R2500_VAL=`\grep '^r2500' ${RES_FINAL} | awk '{ print $2 }'`
+R2500_ERR_L=`\grep '^r2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+R2500_ERR_U=`\grep '^r2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+M2500_VAL=`\grep '^m2500' ${RES_FINAL} | awk '{ print $2 }'`
+M2500_ERR_L=`\grep '^m2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+M2500_ERR_U=`\grep '^m2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+L2500_VAL=`\grep '^L2500' ${RES_FINAL} | awk '{ print $2 }'`
+L2500_ERR=`\grep '^L2500' ${RES_FINAL} | awk '{ print $4 }'`
+MGAS2500_VAL=`\grep '^gas_m2500' ${RES_FINAL} | awk '{ print $2 }'`
+MGAS2500_ERR_L=`\grep '^gas_m2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+MGAS2500_ERR_U=`\grep '^gas_m2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+FGAS2500_VAL=`\grep '^gas_fraction2500' ${RES_FINAL} | awk '{ print $2 }'`
+FGAS2500_ERR_L=`\grep '^gas_fraction2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
+FGAS2500_ERR_U=`\grep '^gas_fraction2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+[ -z "${R2500_VAL}" ]      && R2500_VAL="null"
+[ -z "${R2500_ERR_L}" ]    && R2500_ERR_L="null"
+[ -z "${R2500_ERR_U}" ]    && R2500_ERR_U="null"
+[ -z "${M2500_VAL}" ]      && M2500_VAL="null"
+[ -z "${M2500_ERR_L}" ]    && M2500_ERR_L="null"
+[ -z "${M2500_ERR_U}" ]    && M2500_ERR_U="null"
+[ -z "${L2500_VAL}" ]      && L2500_VAL="null"
+[ -z "${L2500_ERR}" ]      && L2500_ERR="null"
+[ -z "${MGAS2500_VAL}" ]   && MGAS2500_VAL="null"
+[ -z "${MGAS2500_ERR_L}" ] && MGAS2500_ERR_L="null"
+[ -z "${MGAS2500_ERR_U}" ] && MGAS2500_ERR_U="null"
+[ -z "${FGAS2500_VAL}" ]   && FGAS2500_VAL="null"
+[ -z "${FGAS2500_ERR_L}" ] && FGAS2500_ERR_L="null"
+[ -z "${FGAS2500_ERR_U}" ] && FGAS2500_ERR_U="null"
+# 2500 }}}
+FGRR=`\grep '^gas_fraction.*r2500.*r500=' ${RES_FINAL} | sed 's/^.*r500=//' | awk '{ print $1 }'`
+FGRR_ERR_L=`\grep '^gas_fraction.*r2500.*r500=' ${RES_FINAL} | sed 's/^.*r500=//' | awk '{ print $2 }' | awk -F'/' '{ print $1 }'`
+FGRR_ERR_U=`\grep '^gas_fraction.*r2500.*r500=' ${RES_FINAL} | sed 's/^.*r500=//' | awk '{ print $2 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
+[ -z "${FGRR}" ]       && FGRR="null"
+[ -z "${FGRR_ERR_L}" ] && FGRR_ERR_L="null"
+[ -z "${FGRR_ERR_U}" ] && FGRR_ERR_U="null"
+## mrl }}}
+
+## rcool & cooling time {{{
+RCOOL=`\grep '^cooling' ${RES_FINAL} | awk '{ print $4 }'`
+COOLING_TIME=`\grep '^cooling' ${RES_FINAL} | awk -F'=' '{ print $2 }' | tr -d ' Gyr'`
+[ -z "${RCOOL}" ] && RCOOL="null"
+[ -z "${COOLING_TIME}" ] && COOLING_TIME="null"
+## cooling time }}}
+## repro/mass }}}
+
+cd ${BASEDIR}
+## orig json file {{{
+printf "## collect data from original info file ...\n"
+OBJ_NAME=`\grep '"Source\ Name' ${JSON_FILE} | sed 's/.*"Source.*":\ //' | sed 's/^"//' | sed 's/"\ *,$//'`
+OBJ_UNAME=`\grep '"Unified\ Name' ${JSON_FILE} | sed 's/.*"Unified.*":\ //' | sed 's/^"//' | sed 's/"\ *,$//'`
+OBJ_RA=`\grep '"R\.\ A\.' ${JSON_FILE} | sed 's/.*"R\.\ A\.":\ //' | sed 's/^"//' | sed 's/"\ *,$//'`
+OBJ_DEC=`\grep '"Dec\.' ${JSON_FILE} | sed 's/.*"Dec\.":\ //' | sed 's/^"//' | sed 's/"\ *,$//'`
+OBJ_XCRA=`\grep '"XCNTRD_RA' ${JSON_FILE} | sed 's/.*"XCNTRD_RA":\ //' | sed 's/^"//' | sed 's/"\ *,$//'`
+OBJ_XCDEC=`\grep '"XCNTRD_DEC' ${JSON_FILE} | sed 's/.*"XCNTRD_DEC":\ //' | sed 's/^"//' | sed 's/"\ *,$//'`
+REDSHIFT=`\grep '"redshift' ${JSON_FILE} | sed 's/.*"redshift.*":\ //' | sed 's/\ *,$//'`
+COOLCORE=`\grep -i '"cool.*core' ${JSON_FILE} | sed 's/.*"[cC]ool.*":\ //' | sed 's/\ *,$//'`
+[ -z "${COOLCORE}" ] && COOLCORE="null"
+OBJ_FEATURE=`\grep '"Feature' ${JSON_FILE} | sed 's/.*"Feature":\ //' | sed 's/^"//' | sed 's/"\ *,*$//'`
+OBJ_NOTE=`\grep '"NOTE' ${JSON_FILE} | sed 's/.*"NOTE":\ //' | sed 's/^"//' | sed 's/"\ *,*$//'`
+
+# T & Z {{{
+T_1R500=`\grep '"T(0\.1.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
+T_1ERR=`\grep '"T_err(.*0\.1.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
+T_1ERR_L=`\grep '"T_err_l.*0\.1.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
+T_1ERR_U=`\grep '"T_err_u.*0\.1.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
+Z_1R500=`\grep '"Z(0\.1.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
+Z_1ERR=`\grep '"Z_err(.*0\.1.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
+Z_1ERR_L=`\grep '"Z_err_l.*0\.1.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
+Z_1ERR_U=`\grep '"Z_err_u.*0\.1.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
+T_2R500=`\grep '"T(0\.2.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
+T_2ERR=`\grep '"T_err(.*0\.2.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
+T_2ERR_L=`\grep '"T_err_l.*0\.2.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
+T_2ERR_U=`\grep '"T_err_u.*0\.2.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
+Z_2R500=`\grep '"Z(0\.2.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
+Z_2ERR=`\grep '"Z_err(.*0\.2.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
+Z_2ERR_L=`\grep '"Z_err_l.*0\.2.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
+Z_2ERR_U=`\grep '"Z_err_u.*0\.2.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
+[ -z "${T_1R500}" ]  && T_1R500="null"
+[ -z "${T_1ERR}" ]   && T_1ERR="null"
+[ -z "${T_1ERR_L}" ] && T_1ERR_L="null"
+[ -z "${T_1ERR_U}" ] && T_1ERR_U="null"
+[ -z "${Z_1R500}" ]  && Z_1R500="null"
+[ -z "${Z_1ERR}" ]   && Z_1ERR="null"
+[ -z "${Z_1ERR_L}" ] && Z_1ERR_L="null"
+[ -z "${Z_1ERR_U}" ] && Z_1ERR_U="null"
+[ -z "${T_2R500}" ]  && T_2R500="null"
+[ -z "${T_2ERR}" ]   && T_2ERR="null"
+[ -z "${T_2ERR_L}" ] && T_2ERR_L="null"
+[ -z "${T_2ERR_U}" ] && T_2ERR_U="null"
+[ -z "${Z_2R500}" ]  && Z_2R500="null"
+[ -z "${Z_2ERR}" ]   && Z_2ERR="null"
+[ -z "${Z_2ERR_L}" ] && Z_2ERR_L="null"
+[ -z "${Z_2ERR_U}" ] && Z_2ERR_U="null"
+# T & Z }}}
+## json data }}}
+
+## output data to JSON file {{{
+printf "## save collected data into file \`${JSON_FILE}' ...\n"
+mv -fv ${JSON_FILE} ${JSON_FILE}_bak
+cat > ${JSON_FILE} << _EOF_
+{
+    "Obs. ID": ${OBS_ID},
+    "Source Name": "${OBJ_NAME}",
+    "Unified Name": "${OBJ_UNAME}",
+    "Obs. Date": "${DATE_OBS}",
+    "Detector": "${ACIS_TYPE}",
+    "Exposure (ks)": ${EXPOSURE_RAW},
+    "Clean Exposure (ks)": ${EXPOSURE_CLEAN},
+    "R. A.": "${OBJ_RA}",
+    "Dec.": "${OBJ_DEC}",
+    "XCNTRD_RA": "${OBJ_XCRA}",
+    "XCNTRD_DEC": "${OBJ_XCDEC}",
+    "nH (10^22 cm^-2)": ${N_H},
+    "redshift": ${REDSHIFT},
+    "E(z)": ${E_Z},
+    "T_ref (keV)": ${T_REF},
+    "Z_ref (solar)": ${Z_REF},
+    "Rmax_SBP (pixel)": ${RMAX_SBP_PIX},
+    "Rmax_Tpro (pixel)": ${RMAX_TPRO_PIX},
+    "Rmax_SBP (kpc)": ${RMAX_SBP_KPC},
+    "Rmax_Tpro (kpc)": ${RMAX_TPRO_KPC},
+    "NFW_Rmin (kpc)": ${NFW_RMIN_KPC},
+    "Model_SBP": "${MODEL_SBP}",
+    "n01": ${N01},
+    "beta1": ${BETA1},
+    "rc1": ${RC1},
+    "rc1_kpc": ${RC1_KPC},
+    "n02": ${N02},
+    "beta2": ${BETA2},
+    "rc2": ${RC2},
+    "rc2_kpc": ${RC2_KPC},
+    "bkg": ${BKG},
+    "R200 (kpc)": ${R200_VAL},
+    "R200_err_lower (1sigma)": ${R200_ERR_L},
+    "R200_err_upper (1sigma)": ${R200_ERR_U},
+    "M200 (M_sun)": ${M200_VAL},
+    "M200_err_lower (1sigma)": ${M200_ERR_L},
+    "M200_err_upper (1sigma)": ${M200_ERR_U},
+    "L200 (erg/s)": ${L200_VAL},
+    "L200_err (1sigma)": ${L200_ERR},
+    "M_gas200 (M_sun)": ${MGAS200_VAL},
+    "M_gas200_err_lower (1sigma)": ${MGAS200_ERR_L},
+    "M_gas200_err_upper (1sigma)": ${MGAS200_ERR_U},
+    "F_gas200": ${FGAS200_VAL},
+    "F_gas200_err_lower (1sigma)": ${FGAS200_ERR_L},
+    "F_gas200_err_upper (1sigma)": ${FGAS200_ERR_U},
+    "R500 (kpc)": ${R500_VAL},
+    "R500_err_lower (1sigma)": ${R500_ERR_L},
+    "R500_err_upper (1sigma)": ${R500_ERR_U},
+    "M500 (M_sun)": ${M500_VAL},
+    "M500_err_lower (1sigma)": ${M500_ERR_L},
+    "M500_err_upper (1sigma)": ${M500_ERR_U},
+    "L500 (erg/s)": ${L500_VAL},
+    "L500_err (1sigma)": ${L500_ERR},
+    "M_gas500 (M_sun)": ${MGAS500_VAL},
+    "M_gas500_err_lower (1sigma)": ${MGAS500_ERR_L},
+    "M_gas500_err_upper (1sigma)": ${MGAS500_ERR_U},
+    "F_gas500": ${FGAS500_VAL},
+    "F_gas500_err_lower (1sigma)": ${FGAS500_ERR_L},
+    "F_gas500_err_upper (1sigma)": ${FGAS500_ERR_U},
+    "R1500": ${R1500_VAL},
+    "R1500_err_lower": ${R1500_ERR_L},
+    "R1500_err_upper": ${R1500_ERR_U},
+    "M1500": ${M1500_VAL},
+    "M1500_err_lower": ${M1500_ERR_L},
+    "M1500_err_upper": ${M1500_ERR_U},
+    "L1500": ${L1500_VAL},
+    "L1500_err": ${L1500_ERR},
+    "M_gas1500": ${MGAS1500_VAL},
+    "M_gas1500_err_lower": ${MGAS1500_ERR_L},
+    "M_gas1500_err_upper": ${MGAS1500_ERR_U},
+    "F_gas1500": ${FGAS1500_VAL},
+    "F_gas1500_err_lower": ${FGAS1500_ERR_L},
+    "F_gas1500_err_upper": ${FGAS1500_ERR_U},
+    "R2500": ${R2500_VAL},
+    "R2500_err_lower": ${R2500_ERR_L},
+    "R2500_err_upper": ${R2500_ERR_U},
+    "M2500": ${M2500_VAL},
+    "M2500_err_lower": ${M2500_ERR_L},
+    "M2500_err_upper": ${M2500_ERR_U},
+    "L2500": ${L2500_VAL},
+    "L2500_err": ${L2500_ERR},
+    "M_gas2500": ${MGAS2500_VAL},
+    "M_gas2500_err_lower": ${MGAS2500_ERR_L},
+    "M_gas2500_err_upper": ${MGAS2500_ERR_U},
+    "F_gas2500": ${FGAS2500_VAL},
+    "F_gas2500_err_lower": ${FGAS2500_ERR_L},
+    "F_gas2500_err_upper": ${FGAS2500_ERR_U},
+    "T(0.1-0.5 R500)": ${T_1R500},
+    "T_err(0.1-0.5 R500)": ${T_1ERR},
+    "T_err_l(0.1-0.5 R500)": ${T_1ERR_L},
+    "T_err_u(0.1-0.5 R500)": ${T_1ERR_U},
+    "Z(0.1-0.5 R500)": ${Z_1R500},
+    "Z_err(0.1-0.5 R500)": ${Z_1ERR},
+    "Z_err_l(0.1-0.5 R500)": ${Z_1ERR_L},
+    "Z_err_u(0.1-0.5 R500)": ${Z_1ERR_U},
+    "T(0.2-0.5 R500)": ${T_2R500},
+    "T_err(0.2-0.5 R500)": ${T_2ERR},
+    "T_err_l(0.2-0.5 R500)": ${T_2ERR_L},
+    "T_err_u(0.2-0.5 R500)": ${T_2ERR_U},
+    "Z(0.2-0.5 R500)": ${Z_2R500},
+    "Z_err(0.2-0.5 R500)": ${Z_2ERR},
+    "Z_err_l(0.2-0.5 R500)": ${Z_2ERR_L},
+    "Z_err_u(0.2-0.5 R500)": ${Z_2ERR_U},
+    "F_gas(R2500-R500)": ${FGRR},
+    "F_gas_err_l(R2500-R500)": ${FGRR_ERR_L},
+    "F_gas_err_u(R2500-R500)": ${FGRR_ERR_U},
+    "R_cool (kpc)": ${RCOOL},
+    "Cooling_time (Gyr)": ${COOLING_TIME},
+    "Cool_core": ${COOLCORE},
+    "Feature": "${OBJ_FEATURE}",
+    "NOTE": "${OBJ_NOTE}"
+},
+_EOF_
+## output JSON }}}
+
+exit 0
+