Remove obsolete 'fit_{d,}beta_mass_profile.sh'

These two scripts are superseded by 'fit_nfwmass_{d,}beta.sh', which
are already merged together as 'fit_nfwmass.sh'.

2 files changed, 0 insertions, 398 deletions

diff --git a/mass_profile/fit_beta_mass_profile.sh b/mass_profile/fit_beta_mass_profile.sh
deleted file mode 100755
index 3795189..0000000
--- a/mass_profile/fit_beta_mass_profile.sh
+++ /dev/null
@@ -1,200 +0,0 @@
-#!/bin/bash
-
-echo $#
-if [ $# -eq 1 ]
-then
-    :
-else
-    echo "Usage:$0 <cfg file>"
-    exit
-fi
-export PGPLOT_FONT=`locate grfont.dat|head -1`
-
-cfg_file=$1
-base_path=`dirname $0`
-echo $base_path
-#initialize profile type name
-t_profile_type=`grep t_profile $cfg_file|awk '{print $2}'`
-#initialize data file name
-t_data_file=`grep t_data_file $cfg_file|awk '{print $2}'`
-#initialize sbp config file
-sbp_cfg=`grep sbp_cfg $cfg_file|awk '{print $2}'`
-#initialize the temperature profile file
-T_file=`grep '^T_file' $sbp_cfg|awk '{print $2}'`
-#echo $t_profile_type
-cm_per_pixel=`grep '^cm_per_pixel' $sbp_cfg|awk '{print $2}'`
-
-#determine which temperature profile to be used, and fit the T profile
-if [ $t_profile_type == zyy ]
-then
-    t_param_file=`grep t_param_file $cfg_file|awk '{print $2}'`
-    $base_path/fit_zyy_model $t_data_file $t_param_file $cm_per_pixel
-    mv -f zyy_dump.qdp ${T_file}
-elif [ $t_profile_type == m0603246 ]
-then
-    $base_path/fit_m0603246 $t_data_file $cm_per_pixel
-    mv -f m0603246_dump.qdp ${T_file}
-elif [ $t_profile_type == wang2012 ]
-then
-    t_param_file=`grep t_param_file $cfg_file|awk '{print $2}'`
-    $base_path/fit_wang2012_model $t_data_file $t_param_file $cm_per_pixel
-    mv -f wang2012_dump.qdp ${T_file}
-elif [ $t_profile_type == allen ]
-then
-    $base_path/fit_allen_model $t_data_file $cm_per_pixel
-    mv -f allen_dump.qdp ${T_file}
-elif [ $t_profile_type == zzl ]
-then
-    t_param_file=`grep t_param_file $cfg_file|awk '{print $2}'`
-    $base_path/fit_zzl_model $t_data_file $t_param_file
-    mv -f zzl_dump.qdp ${T_file}
-else
-    echo temperature profile name invalid!
-    exit
-fi
-
-cfunc_file=`grep '^cfunc_file' ${sbp_cfg} |awk '{print $2}'`
-z=`grep  '^z' ${sbp_cfg}|awk '{print $2}'`
-abund=`grep '^abund' ${cfg_file} |awk '{print $2}'`
-nh=`grep '^nh' ${cfg_file} |awk '{print $2}'`
-$base_path/coolfunc_calc.sh ${T_file} $abund $nh $z $cfunc_file
-mv flux_cnt_ratio.txt flux_cnt_ratio_center.txt
-#fit sbp
-$base_path/fit_beta_sbp $sbp_cfg
-echo $cfunc_file
-#exit
-
-#store central valu
-mv sbp_fit.qdp sbp_fit_center.qdp
-mv mass_int.qdp mass_int_center.qdp
-mv overdensity.qdp overdensity_center.qdp
-mv gas_mass_int.qdp gas_mass_int_center.qdp
-sbp_data_file=`grep sbp_file $sbp_cfg|awk '{print $2}'`
-radius_sbp_file=`grep radius_sbp_file ${cfg_file}|awk '{print $2}'`
-
-if [ x"$radius_sbp_file" == x ]
-then
-    echo "Error, must have radius_sbp_file assigned, this file should be a 4-column file, which contains the radius, radius err, sbp, and sbp err"
-    exit
-fi
-
-cat ${radius_sbp_file} | sed 's/#.*$//' | grep -Ev '^\s*$' > .tmp.txt
-mv .tmp.txt ${radius_sbp_file}
-
-rm -f summary_shuffle_mass_profile.qdp
-rm -f summary_overdensity.qdp
-rm -f summary_mass_profile.qdp
-rm -f summary_gas_mass_profile.qdp
-
-#100 times of Monte-carlo simulation to determine error
-#just repeat above steps
-for i in `seq 1 100`
-do
-    echo $t_data_file
-    $base_path/shuffle_T.py $t_data_file temp_shuffled_t.dat
-    $base_path/shuffle_sbp.py $sbp_data_file temp_shuffled_sbp.dat
-    
-#exit
-    if [ $t_profile_type == zyy ]
-    then
-	t_param_file=`grep t_param_file $cfg_file|awk '{print $2}'`
-	$base_path/fit_zyy_model temp_shuffled_t.dat $t_param_file $cm_per_pixel
-	mv -f zyy_dump.qdp ${T_file}
-    elif [ $t_profile_type == m0603246 ]
-    then
-	$base_path/fit_m0603246 temp_shuffled_t.dat $cm_per_pixel
-	mv -f m0603246_dump.qdp ${T_file}
-    elif [ $t_profile_type == wang2012 ]
-    then
-	t_param_file=`grep t_param_file $cfg_file|awk '{print $2}'`
-	$base_path/fit_wang2012_model temp_shuffled_t.dat $t_param_file $cm_per_pixel
-    mv -f wang2012_dump.qdp ${T_file}
-    elif [ $t_profile_type == allen ]
-    then
-	$base_path/fit_allen_model temp_shuffled_t.dat $cm_per_pixel
-	mv -f allen_dump.qdp ${T_file}
-    elif [ $t_profile_type == zzl ]
-    then
-	t_param_file=`grep t_param_file $cfg_file|awk '{print $2}'`
-    $base_path/fit_zzl_model temp_shuffled_t.dat $t_param_file
-    mv -f zzl_dump.qdp ${T_file}
-    else
-	echo temperature profile name invalid!
-	exit
-    fi
-    
-#exit
-    echo >temp_sbp.cfg
-    
-    cat $sbp_cfg|while read l
-do
-    if echo $l|grep sbp_file >/dev/null
-    then
-	echo sbp_file temp_shuffled_sbp.dat >>temp_sbp.cfg
-    elif echo $l|grep T_file >/dev/null
-    then
-	echo T_file ${T_file} >>temp_sbp.cfg
-    else
-	echo $l >>temp_sbp.cfg
-    fi
-    
-done
-
-$base_path/coolfunc_calc.sh ${T_file} $abund $nh $z $cfunc_file
-
-$base_path/fit_beta_sbp temp_sbp.cfg
-cat mass_int.qdp >>summary_mass_profile.qdp
-echo no no no >>summary_mass_profile.qdp
-
-cat overdensity.qdp >>summary_overdensity.qdp
-echo no no no >>summary_overdensity.qdp
-
-cat gas_mass_int.qdp >>summary_gas_mass_profile.qdp
-echo no no no >>summary_gas_mass_profile.qdp
-done
-#analys the errors
-$base_path/analyze_mass_profile.py 200
-$base_path/analyze_mass_profile.py 500
-$base_path/analyze_mass_profile.py 1500
-#$base_path/analyze_mass_profile.py 2500
-
-r500=`$base_path/analyze_mass_profile.py 500|grep r500|awk '{print $2}'`
-#$base_path/calc_lx $radius_sbp_file flux_cnt_ratio_center.txt $z $r500 $t_data_file
-r200=`$base_path/analyze_mass_profile.py 200|grep r200|awk '{print $2}'`
-#$base_path/calc_lx $radius_sbp_file flux_cnt_ratio_center.txt $z $r200 $t_data_file
-
-r500e=`$base_path/analyze_mass_profile.py 500|grep '^r500' 2>/dev/null|awk '{print $2,$3}'`
-m500e=`$base_path/analyze_mass_profile.py 500|grep '^m500' 2>/dev/null|awk '{print $2,$3}'`
-L500=`$base_path/calc_lx $radius_sbp_file flux_cnt_ratio_center.txt $z $r500 Tprofile.dat 2>/dev/null|awk '{print $2,$3,$4}'`
-mg500e=`$base_path/analyze_mass_profile.py 500|grep '^gas_m' 2>/dev/null|awk '{print $2,$3}'`
-fg500e=`$base_path/analyze_mass_profile.py 500|grep '^gas_fraction' 2>/dev/null|awk '{print $2,$3}'`
-
-
-r200e=`$base_path/analyze_mass_profile.py 200|grep '^r200' 2>/dev/null|awk '{print $2,$3}'`
-m200e=`$base_path/analyze_mass_profile.py 200|grep '^m200' 2>/dev/null|awk '{print $2,$3}'`
-L200=`$base_path/calc_lx $radius_sbp_file flux_cnt_ratio_center.txt $z $r200 Tprofile.dat 2>/dev/null|awk '{print $2,$3,$4}'`
-mg200e=`$base_path/analyze_mass_profile.py 200|grep '^gas_m' 2>/dev/null|awk '{print $2,$3}'`
-fg200e=`$base_path/analyze_mass_profile.py 200|grep '^gas_fraction' 2>/dev/null|awk '{print $2,$3}'`
-
-
-
-echo "******************"
-echo "Final results:"
-echo "******************"
-echo 
-echo 
-
-echo r500= $r500e  kpc
-echo m500= $m500e M_sun
-echo L500= $L500 erg/s
-echo gas mass 500= $mg500e M_sun
-echo gas fractho 500= $fg500e x100%
-
-echo r200= $r200e kpc
-echo m200= $m200e M_sun
-echo L200= $L200 erg/s
-echo gas mass 200= $mg200e M_sun
-echo gas fractho 200= $fg200e x100%
-
-
-
diff --git a/mass_profile/fit_dbeta_mass_profile.sh b/mass_profile/fit_dbeta_mass_profile.sh
deleted file mode 100755
index a57c3a8..0000000
--- a/mass_profile/fit_dbeta_mass_profile.sh
+++ /dev/null
@@ -1,198 +0,0 @@
-#!/bin/bash
-
-echo $#
-if [ $# -eq 1 ]
-then
-    :
-else
-    echo "Usage:$0 <cfg file>"
-    exit
-fi
-export PGPLOT_FONT=`locate grfont.dat|head -1`
-
-cfg_file=$1
-base_path=`dirname $0`
-echo $base_path
-#initialize profile type name
-t_profile_type=`grep t_profile $cfg_file|awk '{print $2}'`
-#initialize data file name
-t_data_file=`grep t_data_file $cfg_file|awk '{print $2}'`
-#initialize sbp config file
-sbp_cfg=`grep sbp_cfg $cfg_file|awk '{print $2}'`
-#initialize the temperature profile file
-T_file=`grep '^T_file' $sbp_cfg|awk '{print $2}'`
-#echo $t_profile_type
-cm_per_pixel=`grep '^cm_per_pixel' $sbp_cfg|awk '{print $2}'`
-#determine which temperature profile to be used, and fit the T profile
-if [ $t_profile_type == zyy ]
-then
-    t_param_file=`grep t_param_file $cfg_file|awk '{print $2}'`
-    $base_path/fit_zyy_model $t_data_file $t_param_file $cm_per_pixel
-    mv -f zyy_dump.qdp ${T_file}
-elif [ $t_profile_type == m0603246 ]
-then
-    $base_path/fit_m0603246 $t_data_file $cm_per_pixel
-    mv -f m0603246_dump.qdp ${T_file}
-elif [ $t_profile_type == wang2012 ]
-then
-    t_param_file=`grep t_param_file $cfg_file|awk '{print $2}'`
-    $base_path/fit_wang2012_model $t_data_file $t_param_file $cm_per_pixel
-    mv -f wang2012_dump.qdp ${T_file}
-elif [ $t_profile_type == allen ]
-then
-    $base_path/fit_allen_model $t_data_file $cm_per_pixel
-    mv -f allen_dump.qdp ${T_file}
-elif [ $t_profile_type == zzl ]
-then
-    t_param_file=`grep t_param_file $cfg_file|awk '{print $2}'`
-    $base_path/fit_zzl_model $t_data_file $t_param_file
-    mv -f zzl_dump.qdp ${T_file}
-else
-    echo temperature profile name invalid!
-    exit
-fi
-
-cfunc_file=`grep '^cfunc_file' ${sbp_cfg} |awk '{print $2}'`
-z=`grep  '^z' ${sbp_cfg}|awk '{print $2}'`
-abund=`grep '^abund' ${cfg_file} |awk '{print $2}'`
-nh=`grep '^nh' ${cfg_file} |awk '{print $2}'`
-$base_path/coolfunc_calc.sh ${T_file} $abund $nh $z $cfunc_file
-mv flux_cnt_ratio.txt flux_cnt_ratio_center.txt
-#fit sbp
-$base_path/fit_dbeta_sbp $sbp_cfg
-echo $cfunc_file
-#exit
-
-#store central valu
-mv sbp_fit.qdp sbp_fit_center.qdp
-mv mass_int.qdp mass_int_center.qdp
-mv overdensity.qdp overdensity_center.qdp
-mv gas_mass_int.qdp gas_mass_int_center.qdp
-sbp_data_file=`grep sbp_file $sbp_cfg|awk '{print $2}'`
-radius_sbp_file=`grep radius_sbp_file ${cfg_file}|awk '{print $2}'`
-
-if [ x"$radius_sbp_file" == x ]
-then
-    echo "Error, must have radius_sbp_file assigned, this file should be a 4-column file, which contains the radius, radius err, sbp, and sbp err"
-    exit
-fi
-
-cat ${radius_sbp_file} | sed 's/#.*$//' | grep -Ev '^\s*$' > .tmp.txt
-mv .tmp.txt ${radius_sbp_file}
-
-rm -f summary_shuffle_mass_profile.qdp
-rm -f summary_overdensity.qdp
-rm -f summary_mass_profile.qdp
-rm -f summary_gas_mass_profile.qdp
-
-#100 times of Monte-carlo simulation to determine error
-#just repeat above steps
-for i in `seq 1 100`
-do
-    echo $t_data_file
-    $base_path/shuffle_T.py $t_data_file temp_shuffled_t.dat
-    $base_path/shuffle_sbp.py $sbp_data_file temp_shuffled_sbp.dat
-
-#exit
-    if [ $t_profile_type == zyy ]
-    then
-	t_param_file=`grep t_param_file $cfg_file|awk '{print $2}'`
-	$base_path/fit_zyy_model temp_shuffled_t.dat $t_param_file $cm_per_pixel
-	mv -f zyy_dump.qdp ${T_file}
-    elif [ $t_profile_type == m0603246 ]
-    then
-	$base_path/fit_m0603246 temp_shuffled_t.dat $cm_per_pixel
-	mv -f m0603246_dump.qdp ${T_file}
-    elif [ $t_profile_type == wang2012 ]
-    then
-	t_param_file=`grep t_param_file $cfg_file|awk '{print $2}'`
-	$base_path/fit_wang2012_model temp_shuffled_t.dat $t_param_file $cm_per_pixel
-	mv -f wang2012_dump.qdp ${T_file}
-    elif [ $t_profile_type == allen ]
-    then
-	$base_path/fit_allen_model temp_shuffled_t.dat $cm_per_pixel
-	mv -f allen_dump.qdp ${T_file}
-    elif [ $t_profile_type == zzl ]
-    then
-	t_param_file=`grep t_param_file $cfg_file|awk '{print $2}'`
-	$base_path/fit_zzl_model temp_shuffled_t.dat $t_param_file
-	mv -f zzl_dump.qdp ${T_file}
-    else
-	echo temperature profile name invalid!
-	exit
-    fi
-    
-#exit
-    echo >temp_sbp.cfg
-    
-    cat $sbp_cfg|while read l
-do
-    if echo $l|grep sbp_file >/dev/null
-    then
-	echo sbp_file temp_shuffled_sbp.dat >>temp_sbp.cfg
-    elif echo $l|grep T_file >/dev/null
-    then
-	echo T_file ${T_file} >>temp_sbp.cfg
-    else
-	echo $l >>temp_sbp.cfg
-    fi
-    
-done
-
-$base_path/coolfunc_calc.sh ${T_file} $abund $nh $z $cfunc_file
-
-$base_path/fit_dbeta_sbp temp_sbp.cfg
-cat mass_int.qdp >>summary_mass_profile.qdp
-echo no no no >>summary_mass_profile.qdp
-
-cat overdensity.qdp >>summary_overdensity.qdp
-echo no no no >>summary_overdensity.qdp
-
-cat gas_mass_int.qdp >>summary_gas_mass_profile.qdp
-echo no no no >>summary_gas_mass_profile.qdp
-done
-#analys the errors
-$base_path/analyze_mass_profile.py 200
-$base_path/analyze_mass_profile.py 500
-$base_path/analyze_mass_profile.py 1500
-#$base_path/analyze_mass_profile.py 2500
-
-r500=`$base_path/analyze_mass_profile.py 500|grep r500|awk '{print $2}'`
-#$base_path/calc_lx $radius_sbp_file flux_cnt_ratio_center.txt $z $r500 $t_data_file
-r200=`$base_path/analyze_mass_profile.py 200|grep r200|awk '{print $2}'`
-#$base_path/calc_lx $radius_sbp_file flux_cnt_ratio_center.txt $z $r200 $t_data_file
-
-r500e=`$base_path/analyze_mass_profile.py 500|grep '^r500' 2>/dev/null|awk '{print $2,$3}'`
-m500e=`$base_path/analyze_mass_profile.py 500|grep '^m500' 2>/dev/null|awk '{print $2,$3}'`
-L500=`$base_path/calc_lx $radius_sbp_file flux_cnt_ratio_center.txt $z $r500 Tprofile.dat 2>/dev/null|awk '{print $2,$3,$4}'`
-mg500e=`$base_path/analyze_mass_profile.py 500|grep '^gas_m' 2>/dev/null|awk '{print $2,$3}'`
-fg500e=`$base_path/analyze_mass_profile.py 500|grep '^gas_fraction' 2>/dev/null|awk '{print $2,$3}'`
-
-
-r200e=`$base_path/analyze_mass_profile.py 200|grep '^r200' 2>/dev/null|awk '{print $2,$3}'`
-m200e=`$base_path/analyze_mass_profile.py 200|grep '^m200' 2>/dev/null|awk '{print $2,$3}'`
-L200=`$base_path/calc_lx $radius_sbp_file flux_cnt_ratio_center.txt $z $r200 Tprofile.dat 2>/dev/null|awk '{print $2,$3,$4}'`
-mg200e=`$base_path/analyze_mass_profile.py 200|grep '^gas_m' 2>/dev/null|awk '{print $2,$3}'`
-fg200e=`$base_path/analyze_mass_profile.py 200|grep '^gas_fraction' 2>/dev/null|awk '{print $2,$3}'`
-
-
-
-echo "******************"
-echo "Final results:"
-echo "******************"
-echo 
-echo 
-
-echo r500= $r500e  kpc
-echo m500= $m500e M_sun
-echo L500= $L500 erg/s
-echo gas mass 500= $mg500e M_sun
-echo gas fractho 500= $fg500e x100%
-
-echo r200= $r200e kpc
-echo m200= $m200e M_sun
-echo L200= $L200 erg/s
-echo gas mass 200= $mg200e M_sun
-echo gas fractho 200= $fg200e x100%
-
-