Update 'chandra_gensbpreg.sh'.

* 'chandra_gensbpreg.sh':
  - Add '1' to denominators when calculate STN to avoid division by zero
  - Add script description
* 'empty_v3.json':
  - Add 'Unified Name' field
* Rename 'calc_lx_simple_v2.sh' -> 'calc_lx_simple.sh'.

2 files changed, 35 insertions, 8 deletions

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}"