diff options
Diffstat (limited to 'bin/fit_mass.sh')
-rwxr-xr-x | bin/fit_mass.sh | 13 |
1 files changed, 10 insertions, 3 deletions
diff --git a/bin/fit_mass.sh b/bin/fit_mass.sh index 8df5966..e4e0614 100755 --- a/bin/fit_mass.sh +++ b/bin/fit_mass.sh @@ -74,6 +74,8 @@ z=`grep '^z' ${sbp_cfg} | awk '{ print $2 }'` cm_per_pixel=`cosmo_calc.py -b --cm-per-pixel ${z}` sed -i'' "s/^cm_per_pixel.*$/cm_per_pixel ${cm_per_pixel}/" ${sbp_cfg} +cfunc_table="coolfunc_table_photon.txt" + if grep -q '^beta2' $sbp_cfg; then MODEL="dbeta" MODEL_NAME="double-beta" @@ -94,8 +96,12 @@ cp -fv ${tprofile_dump} ${tprofile} mv -fv ${tprofile_dump} ${tprofile_center} mv -fv fit_result.qdp ${tprofile_fit_center} -${base_path}/calc_coolfunc.sh ${tprofile_center} \ - ${abund} ${nh} ${z} ${cfunc_profile} +if [ ! -f ${cfunc_table} ]; then + ${base_path}/calc_coolfunc_table.py -Z ${abund} -n ${nh} -z ${z} \ + -u photon -o ${cfunc_table} +fi +${base_path}/calc_coolfunc_profile.py -C -t ${cfunc_table} \ + -T ${tprofile_center} -o ${cfunc_profile} cfunc_profile_center="coolfunc_profile_center.txt" cp -f ${cfunc_profile} ${cfunc_profile_center} @@ -163,7 +169,8 @@ for i in `seq 1 ${MC_TIMES}`; do printf "## ${i} / ${MC_TIMES} ##\n" printf "## `pwd -P` ##\n" - ${base_path}/calc_coolfunc.sh ${tprofile} ${abund} ${nh} ${z} ${cfunc_profile} + ${base_path}/calc_coolfunc_profile.py -C -t ${cfunc_table} \ + -T ${tprofile} -o ${cfunc_profile} ${base_path}/${PROG_SBPFIT} ${TMP_SBP_CFG} 2> /dev/null cat ${RES_SBPFIT} ${base_path}/fit_nfw_mass mass_int.dat ${z} ${nfw_rmin_kpc} 2> /dev/null |