diff options
Diffstat (limited to 'scripts')
| -rwxr-xr-x | scripts/chandra_genspcreg.sh | 138 |
1 files changed, 72 insertions, 66 deletions
diff --git a/scripts/chandra_genspcreg.sh b/scripts/chandra_genspcreg.sh index 2c0a8d3..22754df 100755 --- a/scripts/chandra_genspcreg.sh +++ b/scripts/chandra_genspcreg.sh @@ -1,80 +1,83 @@ #!/bin/sh -## -## This script generate a series of regions for the extraction of -## radial surface brightness profile (SBP). -## -## Regions geneartion algorithm: -## (TODO) -## -## Author: Zhenghao ZHU -## Created: ??? -## -UPDATED="2015/06/03" -## ChangeLogs: -## v2.0, 2015/06/03, Aaron LI -## * Copy needed pfiles to current working directory, and -## set environment variable $PFILES to use these first. -## * Added missing punlearn -## * Removed the section of dmlist.par & dmextract.par deletion -## +# +# This script generate a series of regions for the extraction of +# radial surface brightness profile (SBP). +# +# Regions geneartion algorithm: +# (TODO) +# +# Author: Zhenghao ZHU +# Created: ??? +# +# Change logs: +# 2017-02-26, Weitian LI +# * Further simplify arguments handling +# * Remove test support (-t) for simplification +# * Rename 'stn' to 'SNR' +# * Add ds9 view +# v2.0, 2015/06/03, Weitian LI +# * Copy needed pfiles to current working directory, and +# set environment variable $PFILES to use these first. +# * Added missing punlearn +# * Removed the section of dmlist.par & dmextract.par deletion -## prepare parameter files (pfiles) {{{ -CIAO_TOOLS="dmstat dmlist dmextract" +# minimal counts +CNT_MIN=2500 -# Copy necessary pfiles for localized usage -for tool in ${CIAO_TOOLS}; do - pfile=`paccess ${tool}` - [ -n "${pfile}" ] && punlearn ${tool} && cp -Lvf ${pfile} . -done +# energy: 700-7000eV -- channel 49:479 +ERANGE=700:7000 +CH_LOW=49 +CH_HI=479 -# Modify environment variable 'PFILES' to use local pfiles first -export PFILES="./:${PFILES}" -## pfiles }}} +# energy 9.5-12keV -- channel 651:822 +CH_BKG_LOW=651 +CH_BKG_HI=822 -if [ $# -ne 6 ] ; then +if [ $# -ne 4 ] ; then printf "usage:\n" - printf " `basename $0` <evt> <evt_e> <bkg_pi> <x> <y> <reg_out>\n" - printf "updated:\n" - printf " ${UPDATED}\n" + printf " `basename $0` <evt> <bkg_pi> <reg_in> <reg_out>\n" exit 1 fi EVT=$1 -EVT_E=$2 -BKGSPC=$3 -X=$4 -Y=$5 -REG_OUT=$6 -[ -f "${REG_OUT}" ] && mv -fv ${REG_OUT} ${REG_OUT}_bak +BKGSPC=$2 +REG_IN=$3 +REG_OUT=$4 + +X=`\grep -i 'point' ${REG_IN} | head -n 1 | tr -d 'a-zA-Z() ' | awk -F',' '{ print $1 }'` +Y=`\grep -i 'point' ${REG_IN} | head -n 1 | tr -d 'a-zA-Z() ' | awk -F',' '{ print $2 }'` echo "EVT: ${EVT}" -echo "EVT_E: ${EVT_E}" +echo "ERANGE: ${ERANGE}" +echo "Center: (${X},${Y})" echo "BKGSPC: ${BKGSPC}" -echo "X: ${X}" -echo "Y: ${Y}" echo "" -#min counts -CNT_MIN=2500 -#singal to noise -STN=10 +## prepare parameter files (pfiles) {{{ +CIAO_TOOLS="dmstat dmlist dmextract" -#energy700:7000 -- channel 49:479 -CH_LOW=49 -CH_HI=479 -#energy 9,5kev-12kev -- channel 651:822 -CH_BKG_LOW=651 -CH_BKG_HI=822 +# Copy necessary pfiles for localized usage +for tool in ${CIAO_TOOLS}; do + pfile=`paccess ${tool}` + [ -n "${pfile}" ] && punlearn ${tool} && cp -Lvf ${pfile} . +done + +# Modify environment variable 'PFILES' to use local pfiles first +export PFILES="./:${PFILES}" +## pfiles }}} + +[ -f "${REG_OUT}" ] && mv -fv ${REG_OUT} ${REG_OUT}_bak RIN=0 ROUT=0 CNTS=0 ROUT_MAX=1500 +SNR=10 -STN_FILE="spc_stn.dat" -[ -e ${STN_FILE} ] && mv -fv ${STN_FILE} ${STN_FILE}_bak +SNR_FILE="spc_snr.dat" +[ -e ${SNR_FILE} ] && mv -fv ${SNR_FILE} ${SNR_FILE}_bak i=0 -while [ `echo "$STN > 2 "| bc -l` -eq 1 ] ; do +while [ `echo "$SNR > 2 "| bc -l` -eq 1 ] ; do if [ `echo "$ROUT > $ROUT_MAX" | bc -l` -eq 1 ]; then break fi @@ -84,13 +87,13 @@ while [ `echo "$STN > 2 "| bc -l` -eq 1 ] ; do echo "${TMP_REG}" >> ${REG_OUT} fi i=`expr $i + 1` - printf "gen reg#$i ...\n" - if [ ${ROUT} -eq 0 ] ; then + printf "Generate region #$i ...\n" + if [ ${ROUT} -eq 0 ] ; then ROUT=5 fi TMP_REG="pie($X,$Y,$RIN,$ROUT,0,360)" punlearn dmlist - CNTS=`dmlist "${EVT_E}[sky=${TMP_REG}]" blocks | \grep 'EVENTS' | awk '{print $8}'` + CNTS=`dmlist "${EVT}[energy=${ERANGE}][sky=${TMP_REG}]" blocks | \grep 'EVENTS' | awk '{print $8}'` while [ ${CNTS} -lt ${CNT_MIN} ]; do ROUT=`expr $ROUT + 1 ` if [ `echo "$ROUT > $ROUT_MAX" | bc -l` -eq 1 ]; then @@ -98,7 +101,7 @@ while [ `echo "$STN > 2 "| bc -l` -eq 1 ] ; do fi TMP_REG="pie($X,$Y,$RIN,$ROUT,0,360)" punlearn dmlist - CNTS=`dmlist "${EVT_E}[sky=${TMP_REG}]" blocks | \grep 'EVENTS' | awk '{print $8}'` + CNTS=`dmlist "${EVT}[energy=${ERANGE}][sky=${TMP_REG}]" blocks | \grep 'EVENTS' | awk '{print $8}'` done TMP_SPC=_tmpspc.pi punlearn dmextract @@ -109,13 +112,13 @@ while [ `echo "$STN > 2 "| 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}' ` - if [ ${INDEX_SRC} -eq 0 ] ;then - STN=10000 + if [ ${INDEX_SRC} -eq 0 ] ;then + SNR=10000 else - STN=`echo ${COUNT_SRC} ${INDEX_SRC} ${COUNT_BKG} ${INDEX_BKG} | awk '{ printf("%f",$1/$2/$3*$4) }' ` + SNR=`echo ${COUNT_SRC} ${INDEX_SRC} ${COUNT_BKG} ${INDEX_BKG} | awk '{ printf("%f",$1/$2/$3*$4) }' ` fi - echo " STN: ${STN}" - echo "${STN}" >> "${STN_FILE}" + echo " SNR: ${SNR}" + echo "${SNR}" >> "${SNR_FILE}" done ## fix 'i', to consistent with the actual annuluses @@ -133,7 +136,7 @@ elif [ $i -gt 6 ]; then mv -fv ${REG_OUT} ${REG_OUT}_2500bak CNTS=0 punlearn dmlist - CNTS_TOTAL=`dmlist "${EVT_E}[sky=pie($X,$Y,0,$RIN,0,360)]" blocks | \grep 'EVENTS' | awk '{print $8}'` + CNTS_TOTAL=`dmlist "${EVT}[energy=${ERANGE}][sky=pie($X,$Y,0,$RIN,0,360)]" blocks | \grep 'EVENTS' | awk '{print $8}'` CNTS_USE=`echo "${CNTS_TOTAL} 6" | awk '{printf("%d", $1/$2)}'` echo "CNT_USE: ${CNT_USE}" printf "*** too many annulus ***\n" @@ -144,12 +147,12 @@ elif [ $i -gt 6 ]; then while [ $j -le 6 ] ; do while [ ${CNTS} -lt ${CNTS_USE} ] ; do ROUT=`expr ${ROUT} + 1 ` - if [ ${ROUT} -gt ${ROUT_MAX} ]; then + if [ ${ROUT} -gt ${ROUT_MAX} ]; then break fi TMP_REG="pie($X,$Y,$RIN,$ROUT,0,360)" punlearn dmlist - CNTS=`dmlist "${EVT_E}[sky=${TMP_REG}]" blocks | \grep 'EVENTS' | awk '{print $8}'` + CNTS=`dmlist "${EVT}[energy=${ERANGE}][sky=${TMP_REG}]" blocks | \grep 'EVENTS' | awk '{print $8}'` done j=`expr $j + 1 ` echo "${TMP_REG}" >> ${REG_OUT} @@ -158,3 +161,6 @@ elif [ $i -gt 6 ]; then done fi +printf "check SBP regions ...\n" +ds9 ${EVT} -regions format ciao -regions system physical \ + -regions ${REG_OUT} -cmap he -bin factor 4 |