diff options
| -rw-r--r-- | files/empty_v3.json | 1 | ||||
| -rwxr-xr-x | scripts/calc_lx_simple.sh (renamed from scripts/calc_lx_simple_v2.sh) | 3 | ||||
| -rwxr-xr-x | scripts/chandra_gensbpreg.sh | 40 |
3 files changed, 36 insertions, 8 deletions
diff --git a/files/empty_v3.json b/files/empty_v3.json index 1095c1e..c150611 100644 --- a/files/empty_v3.json +++ b/files/empty_v3.json @@ -1,6 +1,7 @@ { "Obs. ID": , "Source Name": "", + "Unified Name": "", "Obs. Date": "", "Detector": "", "Exposure (ks)": , diff --git a/scripts/calc_lx_simple_v2.sh b/scripts/calc_lx_simple.sh index 576a43e..a13c67b 100755 --- a/scripts/calc_lx_simple_v2.sh +++ b/scripts/calc_lx_simple.sh @@ -14,6 +14,8 @@ export PATH="/usr/bin:$PATH" ## 2014/06/18: use env variable 'MASS_PROFILE_DIR' ## ########################################################### +UPDATED="2014/06/18" + ## usage, `path_conffile' is the configuration file ## which contains the `path' to each `repro/mass' directory if [ $# -ne 1 ]; then @@ -21,6 +23,7 @@ if [ $# -ne 1 ]; then printf " `basename $0` <mass_dir>\n" printf "\nNOTE:\n" printf " script cannot handle \`~' in path\n" + printf "\nVersion: ${UPDATED}\n" exit 1 fi diff --git a/scripts/chandra_gensbpreg.sh b/scripts/chandra_gensbpreg.sh index df56501..ecc1fe9 100755 --- a/scripts/chandra_gensbpreg.sh +++ b/scripts/chandra_gensbpreg.sh @@ -1,19 +1,41 @@ #!/bin/sh -# -# v2.0, 2013/03/06, LIweitiaNux -# add param `-t' to test STN -# -########################################################### +## +## This script generate a series of regions for the extraction of +## radial surface brightness profile (SBP). +## +## Regions geneartion algorithm: +## (1) innermost 10 regions, we require a mininal of 5 pixel as well +## as 50 counts within the 0.7-7.0 keV range. +## (2) following regions: R_out = R_in * 1.2, and require STN > 1.5. +## +## Reference: +## [1] region generate algorithm ??? (TODO) +## +## Author: Zhenghao ZHU +## Created: ??? (TODO) +## +## ChangeLogs: +## v2.1, 2015/02/13, Weitian LI +## * added '1' to denominators when calculate STN to avoid division by zero +## * added script description +## v2.0, 2013/03/06, Weitian LI +## * added the parameter `-t' to print STN results for testing +## + +UPDATED="2015/02/13" # minimal counts CNT_MIN=50 + # energy: 700-7000eV -- channel 49:479 CH_LOW=49 CH_HI=479 + # energy 9.5-12keV -- channel 651:822 CH_BKG_LOW=651 CH_BKG_HI=822 + if [ $# -lt 6 ]; then printf "usage:\n" printf " `basename $0` <evt> <evt_e> <x> <y> <bkg_pi> <reg_out>\n" @@ -57,6 +79,7 @@ X=$3 Y=$4 BKGSPC=$5 REG_OUT=$6 + [ -f "${REG_OUT}" ] && mv -fv ${REG_OUT} ${REG_OUT}_bak echo "EVT: ${EVT}" echo "EVT_E: ${EVT_E}" @@ -117,10 +140,11 @@ while [ `echo "${STN} > 1.5" | bc -l` -eq 1 ]; do COUNT_SRC=`dmstat "${TMP_SPC}[channel=${CH_LOW}:${CH_HI}][cols counts]" | grep "sum:" | awk '{print $2}'` COUNT_BKG=`dmstat "${BKGSPC}[channel=${CH_LOW}:${CH_HI}][cols counts]" | grep "sum:" | awk '{print $2}'` - # echo "CNT_SRC: ${COUNT_SRC}, IDX_SRC: ${INDEX_SRC}, CNT_BKG: ${COUNT_BKG}, IDX_BKG: ${INDEX_BKG}" - # exit + #echo "CNT_SRC: ${COUNT_SRC}, IDX_SRC: ${INDEX_SRC}, CNT_BKG: ${COUNT_BKG}, IDX_BKG: ${INDEX_BKG}" + #exit 99 - STN=`echo ${COUNT_SRC} ${INDEX_SRC} ${COUNT_BKG} ${INDEX_BKG} | awk '{ printf("%f",$1/$2/$3*$4) }'` + # Add '1' to the denominators to avoid division by zero. + STN=`echo ${COUNT_SRC} ${INDEX_SRC} ${COUNT_BKG} ${INDEX_BKG} | awk '{ printf("%f", ($1 / ($2 + 1)) / ($3 / ($4 + 1))) }'` CNT=`dmlist "${EVT_E}[sky=${TMP_REG}]" blocks | grep 'EVENTS' | awk '{ print $8 }'` echo "CNT: ${CNT}" echo "CNT_MIN: ${CNT_MIN}" |