aboutsummaryrefslogtreecommitdiffstats
path: root/scripts
diff options
context:
space:
mode:
authorAaron LI <aaronly.me@outlook.com>2017-02-26 14:26:26 +0800
committerAaron LI <aaronly.me@outlook.com>2017-02-26 14:26:26 +0800
commitecbf372e3e2239516d7d29a35640e077f9915ba9 (patch)
treeb3cc8beb63db63603982a220836457eed2bd93dc /scripts
parentdf176d449931fcb4c96c082c2385623be0cb11c2 (diff)
downloadchandra-acis-analysis-ecbf372e3e2239516d7d29a35640e077f9915ba9.tar.bz2
ciao_bkg_spectra.sh: Update to use 'manifest.py' and 'results.py'
Diffstat (limited to 'scripts')
-rwxr-xr-xscripts/ciao_bkg_spectra.sh147
1 files changed, 29 insertions, 118 deletions
diff --git a/scripts/ciao_bkg_spectra.sh b/scripts/ciao_bkg_spectra.sh
index 96560e6..dc12199 100755
--- a/scripts/ciao_bkg_spectra.sh
+++ b/scripts/ciao_bkg_spectra.sh
@@ -1,26 +1,24 @@
#!/bin/sh
#
-unalias -a
-export LC_COLLATE=C
-###########################################################
-## extract background spectra from src and blanksky ##
-## renormalization the blank spectrum ##
-## ##
-## Ref: Chandra spectrum analysis ##
-## http://cxc.harvard.edu/ciao/threads/extended/ ##
-## Ref: specextract ##
-## http://cxc.harvard.edu/ciao/ahelp/specextract.html ##
-## Ref: CIAO v4.4 region bugs ##
-## http://cxc.harvard.edu/ciao/bugs/regions.html#bug-12187
-## ##
-## Weitian LI ##
-## 2012/07/24 ##
-###########################################################
-
-VERSION="v5.0"
-UPDATED="2015/06/02"
-# ChangeLogs:
-# v5.0, 2015/06/02, Aaron LI
+# extract background spectra from src and blanksky
+# renormalization the blank spectrum
+#
+# References:
+# * Chandra spectrum analysis
+# http://cxc.harvard.edu/ciao/threads/extended/
+# * specextract
+# http://cxc.harvard.edu/ciao/ahelp/specextract.html
+# * CIAO v4.4 region bugs
+# http://cxc.harvard.edu/ciao/bugs/regions.html#bug-12187
+#
+# Weitian LI
+# 2012-07-24
+
+# Change logs:
+# 2017-02-26, Weitian LI
+# * Use 'manifest.py' and 'results.py'
+# * Simplify 'specextract' parameters handling
+# v5.0, 2015/06/02, Weitian LI
# * Removed 'GREP_OPTIONS' and replace 'grep' with '\grep'
# * Removed 'trap INT date'
# * Copy needed pfiles to current working directory, and
@@ -47,12 +45,10 @@ UPDATED="2015/06/02"
case "$1" in
-[hH]*|--[hH]*)
printf "usage:\n"
- printf " `basename $0` evt=<evt2_clean> reg=<reglist> blank=<blanksky_evt> basedir=<base_dir> nh=<nH> z=<redshift> [ grouptype=<NUM_CTS|BIN> grouptypeval=<number> binspec=<binspec> log=<log_file> ]\n"
+ printf " `basename $0` evt=<evt2_clean> reg=<reglist> blank=<blanksky_evt> nh=<nH> z=<redshift> [ grouptype=<NUM_CTS|BIN> grouptypeval=<number> binspec=<binspec> log=<log_file> ]\n"
printf "\nNotes:\n"
printf " If grouptype=NUM_CTS, then grouptypeval required.\n"
printf " If grouptype=BIN, then binspec required.\n"
- printf "\nversion:\n"
- printf " ${VERSION}, ${UPDATED}\n"
exit ${ERR_USG}
;;
esac
@@ -60,14 +56,9 @@ esac
## default parameters {{{
# default `event file' which used to match `blanksky' files
-#DFT_EVT="_NOT_EXIST_"
-DFT_EVT="`\ls evt2*_clean.fits`"
+DFT_EVT=$(manifest.py -b getpath -r evt2_clean)
# default `blanksky file'
-#DFT_BLANK="_NOT_EXIST_"
-DFT_BLANK="`\ls blanksky*.fits`"
-# default dir which contains `asols, asol.lis, ...' files
-#DFT_BASEDIR="_NOT_EXIST_"
-DFT_BASEDIR=".."
+DFT_BLANK=$(manifest.py -b getpath -r bkg_blank)
# default parameters for 'dmgroup'
DFT_GROUPTYPE="NUM_CTS"
DFT_GROUPTYPEVAL="20"
@@ -76,15 +67,9 @@ DFT_BINSPEC="1:128:2,129:256:4,257:512:8,513:1024:16"
# default `log file'
DFT_LOGFILE="bkg_spectra_`date '+%Y%m%d'`.log"
-## howto find files in `basedir'
-# default `asol.lis pattern'
-DFT_ASOLIS_PAT="acis*asol?.lis"
-# default `bad pixel filename pattern'
-DFT_BPIX_PAT="acis*repro*bpix?.fits"
-# default `pbk file pattern'
-DFT_PBK_PAT="acis*pbk?.fits"
-# default `msk file pattern'
-DFT_MSK_PAT="acis*msk?.fits"
+ASOL=$(manifest.py -b -s "," getpath -r asol)
+BPIX=$(manifest.py -b getpath -r bpix)
+MSK=$(manifest.py -b getpath -r msk)
## default parameters }}}
## error code {{{
@@ -203,33 +188,18 @@ else
fi
# check given nH
if [ -z "${nh}" ]; then
- read -p "> value of nH: " N_H
+ N_H=$(results.py get nh)
else
N_H=${nh}
fi
printf "## use nH: ${N_H}\n" | ${TOLOG}
# check given redshift
if [ -z "${z}" ]; then
- read -p "> value of redshift: " REDSHIFT
+ REDSHIFT=$(results.py get z)
else
REDSHIFT=${z}
fi
printf "## use redshift: ${REDSHIFT}\n" | ${TOLOG}
-# check given dir
-if [ -d "${basedir}" ]; then
- BASEDIR=${basedir}
-elif [ -d "${DFT_BASEDIR}" ]; then
- BASEDIR=${DFT_BASEDIR}
-else
- read -p "> basedir (contains asol files): " BASEDIR
- if [ ! -d "${BASEDIR}" ]; then
- printf "ERROR: given \`${BASEDIR}' NOT a directory\n"
- exit ${ERR_DIR}
- fi
-fi
-# remove the trailing '/'
-BASEDIR=`echo ${BASEDIR} | sed 's/\/*$//'`
-printf "## use basedir: \`${BASEDIR}'\n" | ${TOLOG}
# check given dmgroup parameters: grouptype, grouptypeval, binspec
if [ -z "${grouptype}" ]; then
GROUPTYPE="${DFT_GROUPTYPE}"
@@ -270,37 +240,6 @@ if [ "x${INVALID}" != "x" ]; then
printf "WARNING: some pie region's END_ANGLE > 360\n" | ${TOLOG}
printf " CIAO v4.4 tools may run into trouble\n"
fi
-
-# check files in `basedir'
-printf "check needed files in basedir \`${BASEDIR}' ...\n"
-# check asolis files
-ASOLIS=`\ls -1 ${BASEDIR}/${DFT_ASOLIS_PAT} | head -n 1`
-if [ -z "${ASOLIS}" ]; then
- printf "ERROR: cannot find \"${DFT_ASOLIS_PAT}\" in dir \`${BASEDIR}'\n"
- exit ${ERR_ASOL}
-fi
-printf "## use asolis: \`${ASOLIS}'\n" | ${TOLOG}
-# check badpixel file
-BPIX=`\ls -1 ${BASEDIR}/${DFT_BPIX_PAT} | head -n 1`
-if [ -z "${BPIX}" ]; then
- printf "ERROR: cannot find \"${DFT_BPIX_PAT}\" in dir \`${BASEDIR}'\n"
- exit ${ERR_BPIX}
-fi
-printf "## use badpixel: \`${BPIX}'\n" | ${TOLOG}
-# check pbk file
-PBK=`\ls -1 ${BASEDIR}/${DFT_PBK_PAT} | head -n 1`
-if [ -z "${PBK}" ]; then
- printf "ERROR: cannot find \"${DFT_PBK_PAT}\" in dir \`${BASEDIR}'\n"
- exit ${ERR_PBK}
-fi
-printf "## use pbk: \`${PBK}'\n" | ${TOLOG}
-# check msk file
-MSK=`\ls -1 ${BASEDIR}/${DFT_MSK_PAT} | head -n 1`
-if [ -z "${MSK}" ]; then
- printf "ERROR: cannot find \"${DFT_MSK_PAT}\" in dir \`${BASEDIR}'\n"
- exit ${ERR_MSK}
-fi
-printf "## use msk: \`${MSK}'\n" | ${TOLOG}
## check files }}}
## prepare parameter files (pfiles) {{{
@@ -347,39 +286,11 @@ for reg_i in ${REGLIST}; do
# NO background response files
# NO background spectrum (generate by self)
# NO spectrum grouping (group by self using `dmgroup')
- # Determine parameters for different versions of specextract {{{
- # 'pbkfile' parameter deprecated in CIAO-4.6
- if `pget specextract pbkfile >/dev/null 2>&1`; then
- P_PBKFILE="pbkfile=${PBK}"
- else
- P_PBKFILE=""
- fi
- # specextract: revision 2013-06:
- # 'correct' parameter renamed to 'correctpsf'
- if `pget specextract correct >/dev/null 2>&1`; then
- P_CORRECT="correct=no"
- else
- P_CORRECT="correctpsf=no"
- fi
- # specextract: revision 2013-12:
- # 'weight' parameter controls whether ONLY ARFs are weighted.
- # 'weight_rmf' added to control whether RMFs are weighted.
- # NOTE:
- # (1) only when 'weight=yes' will the 'weight_rmf' parameter be used.
- # (2) no longer distingush between unweighted and weighted reponses
- # in file extension; only .arf & .rmf are now used.
- if `pget specextract weight_rmf >/dev/null 2>&1`; then
- P_WEIGHTRMF="weight_rmf=yes"
- else
- P_WEIGHTRMF=""
- fi
- # }}}
punlearn specextract
specextract infile="${EVT}[sky=region(${REG_TMP})]" \
- outroot=${LBKG_PI%.pi} bkgfile="" asp="@${ASOLIS}" \
+ outroot=${LBKG_PI%.pi} bkgfile="" asp="${ASOL}" \
mskfile="${MSK}" badpixfile="${BPIX}" \
- ${P_PBKFILE} ${P_CORRECT} \
- weight=yes ${P_WEIGHTRMF} \
+ correctpsf=no weight=yes weight_rmf=yes \
energy="0.3:11.0:0.01" channel="1:1024:1" \
bkgresp=no combine=no binarfwmap=2 \
grouptype=NONE binspec=NONE \