From 51fb8a6523fe6ec52b393a00bab8c40ebef8af1d Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Mon, 6 Feb 2017 22:28:07 +0800 Subject: Update HOWTO w.r.t. updated xspec scripts And other small updates --- HOWTO_chandra_acis_analysis.txt | 42 ++++++++++++++++++++--------------------- 1 file 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 '_spec.conf' like: +(17) create a config file '_spec.conf' looks like (for batch process): basedir ../.. reg rspec.reg bkgd bkgcorr_blanksky_lbkg.pi | lbkg.pi nh z (18) $ ciao_deproj_spectra.sh reg="rspec.reg" bkgd= basedir="../.." nh= 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> @/xspec_tprofile.tcl - (input 1) + (calculate average temperature and abundance) + xspec> @/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 '_expcorr.conf' with contents: +(20) fix 'NULL' values in 'tprofile.qdp', 'tprofile.txt' & 'tz_average.txt' +(21) $ cd ../../img; +(22) create config '_expcorr.conf' (for batch process): basedir .. reg sbprofile.reg 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 [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 ' to 'global.cfg' (35) $ fitmass global.cfg c # calculate the center values (36) $ fitmass global.cfg 2>&1 | tee global_.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; -- cgit v1.2.2