diff options
| -rw-r--r-- | HOWTO_chandra_acis_analysis.txt | 42 |
1 files changed, 21 insertions, 21 deletions
diff --git a/HOWTO_chandra_acis_analysis.txt b/HOWTO_chandra_acis_analysis.txt index 81cc959..85ce39b 100644 --- a/HOWTO_chandra_acis_analysis.txt +++ b/HOWTO_chandra_acis_analysis.txt @@ -62,30 +62,26 @@ nH tool: https://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl $ ln -s ../../img/rspec.reg img_rspec.reg (16) ds9 open 'evt2*_clean.fits' with regs 'img_rspec.reg'; adjust the regions and save as 'rspec.reg' -(17) create conf file '<NAME>_spec.conf' like: +(17) create a config file '<NAME>_spec.conf' looks like (for batch process): basedir ../.. reg rspec.reg bkgd bkgcorr_blanksky_lbkg.pi | lbkg.pi nh <nh> z <redshift> (18) $ ciao_deproj_spectra.sh reg="rspec.reg" bkgd=<bkgd> basedir="../.." nh=<nh> z=<z> -(19) Fit the radial spectra to derive the radial temperature profile: +(19) Fit the radial spectra to derive the radial temperature profile, + as well as the average temperature and abundance: $ xspec xspec> @xspec_deproj.xcm xspec> fit - (more tweaks when necessary) + (tweaks parameters when necessary) xspec> @<path>/xspec_tprofile.tcl - (input 1) + (calculate average temperature and abundance) + xspec> @<path>/xspec_avg_tz.tcl xspec> exit -(20) fix 'NULL' values in 'tcl_temp_profile.qdp' & 'tcl_temp_profile.txt' -(21) Determine the average temperature and abundance: - $ xspec - xspec> newpar <..>=<..> (bind all kT's and Abundanc's) - xspec> fit - (get average temperature and abundance) - xspec> exit -(22) $ cd ../../img; - create '<NAME>_expcorr.conf' with contents: +(20) fix 'NULL' values in 'tprofile.qdp', 'tprofile.txt' & 'tz_average.txt' +(21) $ cd ../../img; +(22) create config '<NAME>_expcorr.conf' (for batch process): basedir .. reg sbprofile.reg nh <nh> @@ -97,12 +93,16 @@ nH tool: https://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl (24) $ cd ../mass; clean_massdir.sh; (25) check needed linked files: - sbprofile.txt, flux_sbp.txt, radius_sbp.txt, tcl_temp_profile.txt + * sbprofile.txt + * flux_sbp.txt + * radius_sbp.txt + * tprofile.txt (26) check needed config files: - wang2012_param.txt, - fitting_mass.conf, fitting_sbp.conf, - fitting_dbeta_mass.conf, fitting_dbeta_sbp.conf, - global.cfg + * wang2012_param.txt + * fitting_mass.conf + * fitting_sbp.conf + * fitting_dbeta_mass.conf + * fitting_dbeta_sbp.conf (27) $ fitwang wang2012_param.txt (28) $ qdp fit_result.qdp # check fitted temperature profile (29) $ fitsbp fitting_mass.conf # single-beta sbp @@ -112,15 +112,15 @@ nH tool: https://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl $ ln -s fitting_dbeta_mass.conf global.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' to 'global.cfg' +(34) add 'nfw_rmin_kpc <rmin_kpc>' to 'global.cfg' (35) $ fitmass global.cfg c # calculate the center values (36) $ fitmass global.cfg 2>&1 | tee global_<date>.log # calculate mass data (37) $ collect_infodata.sh (38) $ cd ../img; $ chandra_update_xcentroid.sh (39) $ cd ../spc/profile; -(40) $ chr500 inner=0.1 outer=0.5 # check 0.5R500 range - $ chr500 inner=0.2 outer=0.5 +(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; |