aboutsummaryrefslogtreecommitdiffstats
path: root/mass_profile/calc_lx_dbeta.sh
blob: 104e4b39ec390d7e1b79f90b38fcd56804b7c0b3 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
#!/bin/sh

echo "### \$#: $#"
echo "### `pwd -P`"
if [ $# -gt 1 ]
then
    :
else
    echo "Usage:$0 <cfg file> <rout in kpc> [c]"
    echo "If central value only, append a \"c\""
    exit
fi
export PATH="/usr/local/bin:/usr/bin:/bin:$PATH"
export PGPLOT_FONT="${HEADAS}/lib/grfont.dat"

# blist file
# (energy bands to calc cooling function data)
BLIST="blist.txt"
[ -e "${BLIST}" ] && mv -f ${BLIST} ${BLIST}_bak
cat > ${BLIST} << _EOF_
bolo
0.7 7
0.1 2.4
_EOF_

cfg_file=$1
rout=$2
base_path=`dirname $0`
echo $base_path
#initialize sbp config file
sbp_cfg=`grep '^sbp_cfg' $cfg_file|awk '{print $2}'`
#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 data file
sbp_data_file=`grep '^sbp_file' $sbp_cfg | awk '{ print $2 }'`
#initialize the temperature profile file
T_file=`grep '^T_file' $sbp_cfg|awk '{print $2}'`
#initialize the rmin_kpc for nfw mass profile fitting
nfw_rmin_kpc=`grep '^nfw_rmin_kpc' $cfg_file|awk '{print $2}'`
#echo $t_profile_type
cm_per_pixel=`grep '^cm_per_pixel' $sbp_cfg|awk '{print $2}'`
da=`python -c "print($cm_per_pixel/(.492/3600/180*3.1415926))"`

#determine which temperature profile to be used, and fit the T profile
if [ "$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 2> /dev/null
    mv -f wang2012_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}'`
dl=`python -c "print($da*(1+$z)**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
#$base_path/coolfunc_calc_bolo.sh ${T_file} $abund $nh $z cfunc_bolo.dat
#$base_path/coolfunc_calc_0.7-7.sh ${T_file} $abund $nh $z cfunc_0.7-7.dat
#$base_path/coolfunc_calc_0.1-2.4.sh ${T_file} $abund $nh $z cfunc_0.1-2.4.dat
$base_path/coolfunc_calc_erg.sh ${T_file} $abund $nh $z "cfunc_" ${BLIST}
mv flux_cnt_ratio.txt flux_cnt_ratio_center.txt
#fit sbp
prog="calc_lx_dbeta"
lx_res="lx_dbeta_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
echo $cfunc_file
#exit

#store central value
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

#rm -f summary_lx.dat
#calculate cooling time
#echo $dl

## calculate center values
if [ $# -eq 3 ]; then
    $base_path/analyze_lx.py
    $base_path/analyze_fx.py
    exit 0
fi


###########################################################
#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
    #t_data_file=temp_shuffled_t.dat
    #exit

    if [ "$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 2> /dev/null
	mv -f wang2012_dump.qdp ${T_file}
    else
	echo temperature profile name invalid!
	exit
    fi

    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

    #exit

    echo "### `pwd -P`"
    echo "### $i ###"

    $base_path/coolfunc_calc.sh ${T_file} $abund $nh $z $cfunc_file
    #$base_path/coolfunc_calc_bolo.sh ${T_file} $abund $nh $z cfunc_bolo.dat
    #$base_path/coolfunc_calc_0.7-7.sh ${T_file} $abund $nh $z cfunc_0.7-7.dat
    #$base_path/coolfunc_calc_0.1-2.4.sh ${T_file} $abund $nh $z cfunc_0.1-2.4.dat
    $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
    #grep Lx $lx_res |awk '{print $2}' >>summary_lx.dat
    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