diff options
Diffstat (limited to 'mass_profile/calc_lxfx.sh')
| -rwxr-xr-x | mass_profile/calc_lxfx.sh | 176 |
1 files changed, 176 insertions, 0 deletions
diff --git a/mass_profile/calc_lxfx.sh b/mass_profile/calc_lxfx.sh new file mode 100755 index 0000000..a943a42 --- /dev/null +++ b/mass_profile/calc_lxfx.sh @@ -0,0 +1,176 @@ +#!/bin/sh +# +# Wrapper script used to calculate the luminosity (Lx) and flux (Fx) data, +# which invokes the programming 'calc_lx_beta' (single-beta SBP) or +# 'calc_lx_dbeta' (double-beta SBP). +# +# Output: +# * lx_result.txt +# * fx_result.txt +# +# Author: Junhua GU +# Created: 2013-06-24 +# +# Weitian LI +# 2016-06-07 +# + +if [ $# -eq 2 ] || [ $# -eq 3 ]; then + : +else + echo "usage:" + echo " `basename $0` <global.cfg> <rout_kpc> [c]" + echo "" + echo "arguments:" + echo " <global.cfg>: main config file used for mass 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" +rout="$2" +case "$3" in + [cC]) + F_C="YES" + ;; + *) + F_C="NO" + ;; +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 }'` +nfw_rmin_kpc=`grep '^nfw_rmin_kpc' $cfg_file | 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 }'` +da=`python -c "print($cm_per_pixel/(0.492/3600/180*3.1415926))"` +dl=`python -c "print($da*(1+$z)**2)"` + +if grep -q '^beta2' $sbp_cfg; then + MODEL="dbeta" +else + MODEL="beta" +fi + +if [ "$t_profile_type" = "wang2012" ]; then + ${base_path}/fit_wang2012_model $t_data_file $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 + +# energy bands for which the cooling function data will be calculated +BLIST="blist.txt" +[ -e "${BLIST}" ] && mv -f ${BLIST} ${BLIST}_bak +cat > ${BLIST} << _EOF_ +bolo +0.7 7 +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 + +PROG="calc_lx_${MODEL}" +LXF_RES="lx_${MODEL}_param.txt" +${base_path}/${PROG} ${sbp_cfg} ${rout} \ + cfunc_bolo.dat \ + cfunc_0.7-7.dat \ + cfunc_0.1-2.4.dat 2> /dev/null +LX1=`grep '^Lx1' ${LX_RES} | awk '{ print $2 }'` +LX2=`grep '^Lx2' ${LX_RES} | awk '{ print $2 }'` +LX3=`grep '^Lx3' ${LX_RES} | awk '{ print $2 }'` +FX1=`grep '^Fx1' ${LX_RES} | awk '{ print $2 }'` +FX2=`grep '^Fx2' ${LX_RES} | awk '{ print $2 }'` +FX3=`grep '^Fx3' ${LX_RES} | awk '{ print $2 }'` + +echo $LX1 $LX2 $LX3 >summary_lx.dat +echo $FX1 $FX2 $FX3 >summary_fx.dat + +# save the calculated central values +mv ${LX_RES} ${LX_RES%.txt}_center.txt +mv lx_sbp_fit.qdp lx_sbp_fit_center.qdp +mv lx_rho_fit.dat lx_rho_fit_center.dat + +# only calculate the central values +if [ "${F_C}" = "YES" ]; then + echo "Calculate the central values only ..." + ${base_path}/analyze_lx.py + ${base_path}/analyze_fx.py + exit 0 +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 + + echo >temp_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" >>temp_sbp.cfg + elif echo "${l}" | grep -q '^T_file' >/dev/null; then + echo "T_file ${T_file}" >>temp_sbp.cfg + else + echo "${l}" >>temp_sbp.cfg + fi + done + + echo "### `pwd -P`" + echo "### $i ###" + ${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} \ + cfunc_bolo.dat \ + cfunc_0.7-7.dat \ + cfunc_0.1-2.4.dat 2> /dev/null + LX1=`grep '^Lx1' ${LX_RES} | awk '{ print $2 }'` + LX2=`grep '^Lx2' ${LX_RES} | awk '{ print $2 }'` + LX3=`grep '^Lx3' ${LX_RES} | awk '{ print $2 }'` + FX1=`grep '^Fx1' ${LX_RES} | awk '{ print $2 }'` + FX2=`grep '^Fx2' ${LX_RES} | awk '{ print $2 }'` + FX3=`grep '^Fx3' ${LX_RES} | awk '{ print $2 }'` + + echo $LX1 $LX2 $LX3 >>summary_lx.dat + echo $FX1 $FX2 $FX3 >>summary_fx.dat +done # end of 'for' + +# analyze lx & fx +${base_path}/analyze_lx.py +${base_path}/analyze_fx.py + +exit 0 |