calc_lxfx.sh: Update to use the new style config files

1 files changed, 47 insertions, 56 deletions

diff --git a/mass_profile/calc_lxfx.sh b/mass_profile/calc_lxfx.sh
index 922c8fc..e240bbe 100755
--- a/mass_profile/calc_lxfx.sh
+++ b/mass_profile/calc_lxfx.sh
@@ -17,24 +17,27 @@
 # Weitian LI
 # 2016-06-07
 #
+# Change logs:
+# 2017-02-07, Weitian LI
+#   * Update to use the new style configuration files
+#
 
 if [ $# -eq 2 ] || [ $# -eq 3 ]; then
     :
 else
     echo "usage:"
-    echo "    `basename $0` <global.cfg> <rout_kpc> [c]"
+    echo "    `basename $0` <mass.conf> <rout_kpc> [c]"
     echo ""
     echo "arguments:"
-    echo "    <global.cfg>: main config file used for mass calculation"
+    echo "    <mass.conf>: config file for mass profile calculation"
     echo "    <rout_kpc>: outer/cut radius within which to calculate Lx & Fx"
     echo "                e.g., r500, r200 (unit: kpc)"
     echo "    [c]: optional; if specified, do not calculate the errors"
     exit 1
 fi
 export PATH="/usr/local/bin:/usr/bin:/bin:$PATH"
-export PGPLOT_FONT="${HEADAS}/lib/grfont.dat"
 
-cfg_file="$1"
+mass_cfg="$1"
 rout="$2"
 case "$3" in
     [cC])
@@ -48,24 +51,17 @@ esac
 base_path=$(dirname $(realpath $0))
 
 ## Extract settings/values from the config file
-# model name of the temperature profile
-t_profile_type=`grep '^t_profile' $cfg_file | awk '{ print $2 }'`
-# parameter file for temperature profile fitting
-t_param_file=`grep '^t_param_file' $cfg_file | awk '{ print $2 }'`
-# data points of the temperature profile
-t_data_file=`grep '^t_data_file' $cfg_file | awk '{ print $2 }'`
-# config file for SBP fitting
-sbp_cfg=`grep '^sbp_cfg' $cfg_file | awk '{ print $2 }'`
-# SBP flux data
-sbp_data_file=`grep '^sbp_file' $sbp_cfg | awk '{ print $2 }'`
-# fitted temperature profile
-T_file=`grep '^T_file' $sbp_cfg | awk '{ print $2 }'`
-# filename of the cooling function data
-cfunc_file=`grep '^cfunc_file' $sbp_cfg |awk '{ print $2 }'`
-cm_per_pixel=`grep '^cm_per_pixel' $sbp_cfg | awk '{ print $2 }'`
-z=`grep '^z' $sbp_cfg | awk '{ print $2 }'`
-nh=`grep '^nh' $cfg_file |awk '{ print $2 }'`
-abund=`grep '^abund' $cfg_file |awk '{ print $2 }'`
+nh=`grep            '^nh'            ${mass_cfg} | awk '{ print $2 }'`
+abund=`grep         '^abund'         ${mass_cfg} | awk '{ print $2 }'`
+tprofile_data=`grep '^tprofile_data' ${mass_cfg} | awk '{ print $2 }'`
+tprofile_cfg=`grep  '^tprofile_cfg'  ${mass_cfg} | awk '{ print $2 }'`
+sbp_cfg=`grep       '^sbp_cfg'       ${mass_cfg} | awk '{ print $2 }'`
+
+sbp_data=`grep      '^sbp_data'      ${sbp_cfg} | awk '{ print $2 }'`
+tprofile=`grep      '^tprofile'      ${sbp_cfg} | awk '{ print $2 }'`
+cfunc_profile=`grep '^cfunc_profile' ${sbp_cfg} | awk '{ print $2 }'`
+z=`grep             '^z'             ${sbp_cfg} | awk '{ print $2 }'`
+cm_per_pixel=`grep  '^cm_per_pixel'  ${sbp_cfg} | awk '{ print $2 }'`
 
 if grep -q '^beta2' $sbp_cfg; then
     MODEL="dbeta"
@@ -73,13 +69,11 @@ else
     MODEL="beta"
 fi
 
-# only 'wang2012' model supported
-if [ "X${t_profile_type}" != "Xwang2012" ]; then
-    echo "ERROR: unsupported temperature profile: ${t_profile_type}"
-    exit 11
-fi
-${base_path}/fit_wang2012_model $t_data_file $t_param_file $cm_per_pixel 2> /dev/null
-mv -f wang2012_dump.qdp ${T_file}
+PROG_TPROFILE="fit_wang2012_model"
+tprofile_dump="wang2012_dump.qdp"
+${base_path}/${PROG_TPROFILE} ${tprofile_data} ${tprofile_cfg} \
+            ${cm_per_pixel} 2> /dev/null
+mv -fv ${tprofile_dump} ${tprofile}
 
 # energy bands for which the cooling function data will be calculated
 BLIST="blist.txt"
@@ -90,9 +84,11 @@ bolo
 0.1 2.4
 _EOF_
 
-${base_path}/coolfunc_calc.sh ${T_file} $abund $nh $z $cfunc_file
-${base_path}/coolfunc_calc_erg.sh ${T_file} $abund $nh $z "cfunc_" ${BLIST}
-mv flux_cnt_ratio.txt flux_cnt_ratio_center.txt
+${base_path}/coolfunc_calc.sh ${tprofile} ${abund} \
+            ${nh} ${z} ${cfunc_profile}
+${base_path}/coolfunc_calc_erg.sh ${tprofile} ${abund} \
+            ${nh} ${z} "cfunc_" ${BLIST}
+mv -fv flux_cnt_ratio.txt flux_cnt_ratio_center.txt
 
 PROG="calc_lx_${MODEL}"
 LXF_RES="lx_${MODEL}_param.txt"
@@ -128,36 +124,33 @@ fi
 # Estimate the errors of Lx and Fx by Monte Carlo simulation
 MC_TIMES=100
 for i in `seq 1 ${MC_TIMES}`; do
-    $base_path/shuffle_T.py $t_data_file temp_shuffled_t.dat
-    $base_path/shuffle_sbp.py $sbp_data_file temp_shuffled_sbp.dat
-
-    if [ "$t_profile_type" = "wang2012" ]; then
-        ${base_path}/fit_wang2012_model temp_shuffled_t.dat $t_param_file $cm_per_pixel 2> /dev/null
-        mv -f wang2012_dump.qdp ${T_file}
-    else
-        echo "ERROR: unsupported temperature profile: ${t_profile_type}"
-        exit 11
-    fi
-
-    # clear ${TMP_SBP_CFG}
-    TMP_SBP_CFG="temp_sbp.cfg"
-    # : > ${TMP_SBP_CFG}
+    ${base_path}/shuffle_profile.py ${tprofile_data} tmp_tprofile.txt
+    ${base_path}/shuffle_profile.py ${sbp_data} tmp_sbprofile.txt
+
+    # temperature profile
+    ${base_path}/${PROG_TPROFILE} tmp_tprofile.txt ${tprofile_cfg} \
+                ${cm_per_pixel} 2> /dev/null
+    mv -f ${tprofile_dump} ${tprofile}
+
+    TMP_SBP_CFG="tmp_sbp.cfg"
     [ -e "${TMP_SBP_CFG}" ] && rm -f ${TMP_SBP_CFG}
     cat ${sbp_cfg} | while read l; do
-        if echo "${l}" | grep -q '^sbp_file' >/dev/null; then
-            echo "sbp_file  temp_shuffled_sbp.dat" >> ${TMP_SBP_CFG}
-        elif echo "${l}" | grep -q '^T_file' >/dev/null; then
-            echo "T_file  ${T_file}" >> ${TMP_SBP_CFG}
+        if echo "${l}" | grep -q '^sbp_data' >/dev/null; then
+            echo "sbp_data  tmp_sbprofile.txt" >> ${TMP_SBP_CFG}
+        elif echo "${l}" | grep -q '^tprofile' >/dev/null; then
+            echo "tprofile  ${tprofile}" >> ${TMP_SBP_CFG}
         else
             echo "${l}" >> ${TMP_SBP_CFG}
         fi
     done
 
     echo "### `pwd -P`"
-    echo "### $i / ${MC_TIMES} ###"
-    ${base_path}/coolfunc_calc.sh ${T_file} $abund $nh $z $cfunc_file
-    ${base_path}/coolfunc_calc_erg.sh ${T_file} $abund $nh $z "cfunc_" ${BLIST}
-    ${base_path}/${PROG} temp_sbp.cfg ${rout} \
+    echo "### ${i} / ${MC_TIMES} ###"
+    ${base_path}/coolfunc_calc.sh ${tprofile} ${abund} \
+                ${nh} ${z} ${cfunc_profile}
+    ${base_path}/coolfunc_calc_erg.sh ${tprofile} ${abund} \
+                ${nh} ${z} "cfunc_" ${BLIST}
+    ${base_path}/${PROG} ${TMP_SBP_CFG} ${rout} \
                 cfunc_bolo.dat \
                 cfunc_0.7-7.dat \
                 cfunc_0.1-2.4.dat 2> /dev/null
@@ -175,5 +168,3 @@ done # end of 'for'
 # analyze Lx & Fx Monte Carlo results
 ${base_path}/analyze_lxfx.py "Lx" summary_lx.dat lx_result.txt ${BLIST}
 ${base_path}/analyze_lxfx.py "Fx" summary_fx.dat fx_result.txt ${BLIST}
-
-exit 0