diff options
| -rw-r--r-- | HOWTO_chandra_acis_analysis.txt | 41 | ||||
| -rw-r--r-- | README.md | 2 |
2 files changed, 24 insertions, 19 deletions
diff --git a/HOWTO_chandra_acis_analysis.txt b/HOWTO_chandra_acis_analysis.txt index 06b08af..4c91682 100644 --- a/HOWTO_chandra_acis_analysis.txt +++ b/HOWTO_chandra_acis_analysis.txt @@ -66,17 +66,18 @@ Step-by-step guide to analyze ACIS data: temp <avg_temp> abund <avg_abund> (23) $ ciao_expcorr_sbp.sh basedir=.. nh=<nh> z=<redshift> temp=<avg_temp> abund=<avg_abund> -## -------------------------------------------------------- -(24) $ cd ../mass; clean_massdir.sh; -(25) check needed linked files: - * sbprofile.txt - * tprofile.txt -(26) check needed config files: +## -------------------------------------------------------- +(24) $ cd ../mass + $ ln -s ../img/sbprofile.txt . + $ ln -s ../spc/profile/tprofile.txt . +(25) Copy the sample config files located at 'files' directory: * mass.conf * wang2012_param.txt * sbp_sbeta.conf * sbp_dbeta.conf +(26) Fill 'nH', 'abund' in 'mass.conf'; + and 'z' in 'sbp_sbeta.conf' and 'sbp_dbeta.conf' (27) $ fittp tprofile.txt wang2012_param.txt (28) $ qdp fit_result.qdp (check fitted temperature profile, and adjust parameter accordingly) @@ -87,32 +88,36 @@ Step-by-step guide to analyze ACIS data: $ ln -s sbp_dbeta.conf sbp.cfg # use double-beta (32) $ fitnfw <z> [rmin_kpc] (33) $ qdp nfw_fit_result.qdp # check fitted nfw profile, and ajust 'rmin_kpc' -(34) add 'nfw_rmin_kpc <rmin_kpc>' to 'mass.conf -(35) $ fitmass mass.conf c # calculate the center values +(34) Update 'nfw_rmin_kpc' in 'mass.conf +(35) $ fitmass mass.conf c # calculate the central values (36) $ fitmass mass.conf 2>&1 | tee mass_<date>.log # calculate mass data (37) $ collect_infodata.sh -(38) $ cd ../img; +(38) $ cd ../img $ chandra_update_xcentroid.sh -(39) $ cd ../spc/profile; +(39) $ cd ../spc/profile (40) $ ciao_r500avgt.sh inner=0.1 outer=0.5 # check 0.5R500 range $ ciao_r500avgt.sh inner=0.2 outer=0.5 (41) $ xspec xspec> @xspec_r500avgt_0.1-0.5.xcm - xspec> fit; fit; - xspec> cpd /xs; pl l del; - xspec> error 1. 2 3 # calc errors (68%) of kT and Abundanc -(42) update the data of 'T(0.1-0.5 R500)', 'Z(0.1-0.5 R500)' ... + xspec> fit; cpd /xs; pl l del; + xspec> error 1.0 2 3 + (calculate the 1 sigma errors for temperature and abundance) +(42) update the following values in the INFO.json file: + * 'T(0.1-0.5 R500)' + * 'Z(0.1-0.5 R500)' +(43) repeat the above two steps for region "0.2-0.5 R500", and update + the values of 'T(0.2-0.5 R500)' and 'Z(0.2-0.5 R500)'. ## -------------------------------------------------------- -(43) $ cd ../..; # in 'repro' dir +(44) $ cd ../..; # in 'repro' dir $ cp -a mass lxfx; cd lxfx; -(44) $ calclxfx global.cfg c 500 200 # 'c' for center values +(45) $ calclxfx global.cfg c 500 200 # 'c' for center values $ calclxfx global.cfg 500 200 # calc 'errors' -(45) $ getlxfx . c 500 200 # for center results +(46) $ getlxfx . c 500 200 # for center results $ getlxfx . 500 200 # for all results (with errors) ## -------------------------------------------------------- -(46) $ cd ..; cd spc/profile; +(47) $ cd ..; cd spc/profile; calculate cooling time and Csb value: $ ciao_calc_ct_csb.sh check the region whether exceeds ccd edge, and answer y/n @@ -94,7 +94,7 @@ alias fitmass="${MASS_PROFILE_DIR}/fit_mass.sh" alias fitnfw="${MASS_PROFILE_DIR}/fit_nfw_mass mass_int.dat" alias fitsbp="${MASS_PROFILE_DIR}/fit_sbp.sh" alias fittp="${MASS_PROFILE_DIR}/fit_wang2012_model" -alias calclxfx="${MASS_PROFILE_DIR}/calc_lxfx_simple.sh" +alias calclxfx="${MASS_PROFILE_DIR}/calc_lxfx_wrapper.sh" alias getlxfx="${MASS_PROFILE_DIR}/get_lxfx_data.sh" ``` |