aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
-rwxr-xr-xscripts/ciao_procevt.sh206
1 files changed, 0 insertions, 206 deletions
diff --git a/scripts/ciao_procevt.sh b/scripts/ciao_procevt.sh
deleted file mode 100755
index 5ee8abd..0000000
--- a/scripts/ciao_procevt.sh
+++ /dev/null
@@ -1,206 +0,0 @@
-#!/bin/sh
-#
-unalias -a
-export LC_COLLATE=C
-###########################################################
-## process `evt2 file' generated by `chandra_repro' ##
-## to produce a `clean evt2 file' ##
-## (remove point source and flares) ##
-## for later image and spectral analysis ##
-## ##
-## NOTES: ##
-## based on previous `ciao_process_evt2.sh' ##
-## support ACIS-I(chip: 0-3) and ACIS-S(chip: 7) ##
-## determine by check `DETNAM' for chip number ##
-## if `DETNAM' has `0123', then `ACIS-I' ##
-## if `DETNAM' has `7', then `ACIS-S' ##
-## ##
-## Weitian LI <liweitianux@gmail.com># ##
-## 2012/08/16 ##
-###########################################################
-##
-VERSION="v3.1"
-UPDATED="2017-02-06"
-##
-## ChangeLogs:
-## v3.1, 2017-02-06, Weitian LI
-## * Specify regions format and system for ds9
-## v3.0, 2015/06/02, Aaron LI
-## * Added 'unalias -a' and 'export LC_COLLATE=C'
-## * Replaced 'grep' with '\grep', 'ls' with '\ls'
-## * Copy needed pfiles to current working directory, and
-## set environment variable $PFILES to use these first.
-## v2.2, 2014/10/30, Weitian LI
-## small fix to the generation of chips script:
-## changed '>>' to '>'
-## v2.1, 2012/08/16, Weitian LI
-## improve invoke `chips', run it in a separate terminal
-##
-
-
-## error code {{{
-ERR_USG=1
-ERR_DIR=11
-ERR_EVT=12
-ERR_BKG=13
-ERR_REG=14
-ERR_ASOL=21
-ERR_BPIX=22
-ERR_PBK=23
-ERR_MSK=24
-ERR_BKGTY=31
-ERR_SPEC=32
-ERR_DET=41
-## error code }}}
-
-## usage, help {{{
-case "$1" in
- -[hH]*|--[hH]*)
- printf "usage:\n"
- printf " `basename $0` evt=<raw_evt_file>\n"
- printf "\nversion:\n"
- printf " ${VERSION}, ${UPDATED}\n"
- exit ${ERR_USG}
- ;;
-esac
-## usage, help }}}
-
-## default parameters {{{
-# default `event file' which used to match `blanksky' files
-#DFT_EVT="_NOT_EXIST_"
-DFT_EVT="`\ls acisf?????*_repro_evt2.fits 2> /dev/null`"
-## default parameters }}}
-
-## functions {{{
-# process commandline arguments
-# cmdline arg format: `KEY=VALUE'
-getopt_keyval() {
- until [ -z "$1" ]
- do
- key=${1%%=*} # extract key
- val=${1#*=} # extract value
- keyval="${key}=\"${val}\""
- echo "## getopt: eval '${keyval}'"
- eval ${keyval}
- shift # shift, process next one
- done
-}
-## functions }}}
-
-## parameters {{{
-# process cmdline args using `getopt_keyval'
-getopt_keyval "$@"
-
-# check given parameters
-# check evt file
-if [ -r "${evt}" ]; then
- EVT=${evt}
-elif [ -r "${DFT_EVT}" ]; then
- EVT=${DFT_EVT}
-else
- read -p "evt2 file: " EVT
- if ! [ -r "${EVT}" ]; then
- printf "ERROR: cannot access given \`${EVT}' evt file\n"
- exit ${ERR_EVT}
- fi
-fi
-printf "## use evt file: \`${EVT}'\n"
-## parameters }}}
-
-## prepare parameter files (pfiles) {{{
-CIAO_TOOLS="dmkeypar dmcopy celldetect dmextract"
-
-# 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 }}}
-
-## determine ACIS type {{{
-punlearn dmkeypar
-DETNAM=`dmkeypar ${EVT} DETNAM echo=yes`
-if echo ${DETNAM} | \grep -q 'ACIS-0123'; then
- printf "## \`DETNAM' (${DETNAM}) has chips 0123\n"
- printf "## ACIS-I\n"
- ACIS_TYPE="ACIS-I"
- CCD="0:3"
- ROOTNAME="evt2_c`echo ${CCD} | tr ':' '-'`"
-elif echo ${DETNAM} | \grep -q 'ACIS-[0-6]*7'; then
- printf "## \`DETNAM' (${DETNAM}) has chip 7\n"
- printf "## ACIS-S\n"
- ACIS_TYPE="ACIS-S"
- CCD="7"
- ROOTNAME="evt2_c${CCD}"
-else
- printf "ERROR: unknown detector type: ${DETNAM}\n"
- exit ${ERR_DET}
-fi
-## ACIS type }}}
-
-## main process {{{
-printf "filter raw evt2 file by ccd_id ...\n"
-EVT2_ORIG="${ROOTNAME}_orig.fits"
-punlearn dmcopy
-dmcopy infile="${EVT}[ccd_id=${CCD}]" outfile=${EVT2_ORIG} clobber=yes
-
-# detect point sources
-printf "using \`celldetect' to find the point sources ...\n"
-CELLD="celld_${ROOTNAME}"
-[ -e "${CELLD}.reg" ] && mv -fv ${CELLD}.reg ${CELLD}.reg_bak
-punlearn celldetect
-celldetect infile=${EVT2_ORIG} outfile="${CELLD}.fits" \
- regfile="${CELLD}.reg" clobber=yes
-
-printf "check the result of \`celldetect' ...\n"
-printf "modify if necessary and save as the same name, \`${CELLD}.reg'\n"
-cp -fv ${CELLD}.reg ${CELLD}_orig.reg
-ds9 ${EVT2_ORIG} -regions format ciao \
- -regions system physical \
- -regions ${CELLD}.reg \
- -cmap he -bin factor 2
-
-EVT2_RMSRCS="${ROOTNAME}_rmsrcs.fits"
-punlearn dmcopy
-dmcopy infile="${EVT2_ORIG}[exclude sky=region(${CELLD}.reg)]" \
- outfile=${EVT2_RMSRCS} clobber=yes
-
-LC_REG="ex_bkg.reg"
-printf "filter flares ...\n"
-printf "select a big source region and save as \`${LC_REG}'\n"
-touch ${LC_REG}
-ds9 ${EVT2_RMSRCS} -cmap he -bin factor 2
-
-printf "create the lightcurve ...\n"
-LC="${LC_REG%.reg}.lc"
-[ -e "${LC}" ] && mv -fv ${LC} ${LC}_bak
-punlearn dmextract
-dmextract infile="${EVT2_RMSRCS}[exclude sky=region(${LC_REG})][bin time=::200]" \
- outfile="${LC}" opt=ltc1 clobber=yes
-
-# generate a script for `chips'
-# for convenience and backup
-LC_SCRIPT="_${LC%.lc}.chips"
-LC_SCALE=1.2
-GTI="${LC%.lc}.gti"
-cat > ${LC_SCRIPT} << _EOF_
-from lightcurves import *
-lc_clean("${LC}")
-lc_clean("${LC}", scale=${LC_SCALE}, outfile="${GTI}")
-print_window("${LC%.lc}_lc.jpg", ["format", "jpg", "clobber", "True"])
-_EOF_
-
-printf "generate GTI ...\n"
-# [-x] let chips run in a separate terminal
-chips -x ${LC_SCRIPT}
-
-printf "remove the background flare ...\n"
-EVT2_CLEAN="${ROOTNAME}_clean.fits"
-punlearn dmcopy
-dmcopy infile="${EVT2_RMSRCS}[@${GTI}]" outfile=${EVT2_CLEAN} clobber=yes
-## main process }}}
-
-printf "FINISHED\n"