Remove 'scripts/ciao_genregs.sh'

1 files changed, 0 insertions, 307 deletions

diff --git a/scripts/ciao_genregs.sh b/scripts/ciao_genregs.sh
deleted file mode 100755
index 7ab3f41..0000000
--- a/scripts/ciao_genregs.sh
+++ /dev/null
@@ -1,307 +0,0 @@
-#!/bin/sh
-##
-## Generate regions for 'radial spectra analysis' (rspec.reg),
-## and regions for 'surface brightness profile' extraction.
-##
-## Author: Weitian LI
-## Created: 2013/10/12
-##
-VERSION="v2.1"
-UPDATED="2017-02-06"
-##
-## Changelogs:
-## v2.1, 2017-02-06, Weitian LI
-##   * Specify regions format and system for ds9
-## v2.0, 2015/06/03, Aaron LI
-##   * Updated script description
-##   * Replaced 'grep' with '\grep', 'ls' with '\ls'
-##   * ds9 colormap changed from 'sls' to 'he'
-##   * Copy needed pfiles to current working directory, and
-##     set environment variable $PFILES to use these first.
-## v1.0, 2013/10/12, Weitian LI
-##   split from 'ciao_expcorr_sbp_v3.sh'
-##
-
-unalias -a
-export LC_COLLATE=C
-
-SCRIPT_PATH=`readlink -f $0`
-SCRIPT_DIR=`dirname ${SCRIPT_PATH}`
-XCENTROID_SCRIPT="chandra_xcentroid.sh"
-GEN_SPCREG_SCRIPT="chandra_genspcreg.sh"
-GEN_SBPREG_SCRIPT="chandra_gensbpreg.sh"
-
-## err 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
-ERR_ENG=42
-ERR_CIAO=100
-## error code }}}
-
-## usage {{{
-case "$1" in
-    -[hH]*|--[hH]*)
-        printf "usage:\n"
-        printf "    `basename $0` evt=<evt2_clean> reg_in=<reg_in> bkgd=<bkgd_spec> ds9=<Yes|no>\n"
-        printf "\nversion:\n"
-        printf "    ${VERSION}, ${UPDATED}\n"
-        exit ${ERR_USG}
-        ;;
-esac
-## usage }}}
-
-## link needed files {{{
-BKGD_FILE=`\ls ../bkg/bkgcorr_bl*.pi 2> /dev/null | head -n 1`
-if [ -r "${BKGD_FILE}" ]; then
-    ln -svf ${BKGD_FILE} .
-fi
-ASOL_FILE=`\ls ../pcad*_asol?.fits 2> /dev/null`
-if [ -r "${ASOL_FILE}" ]; then
-    ln -svf ${ASOL_FILE} .
-fi
-CELL_REG_FILE=`\ls ../evt/celld*.reg 2> /dev/null | \grep -v 'orig'`
-if [ -r "${CELL_REG_FILE}" ]; then
-    ln -svf ${CELL_REG_FILE} .
-fi
-# }}}
-
-## default parameters {{{
-## clean evt2 file
-DFT_EVT=`\ls evt2*_clean.fits 2> /dev/null`
-## the repro dir
-DFT_BASEDIR=".."
-# default `asol file'
-ASOL="`\ls pcadf*_asol1.fits 2> /dev/null | head -n 1`"
-
-## energy range
-# format: `E_START:E_END:E_WIDTH'
-DFT_ENERGY=700:7000:100
-E_START=`echo ${DFT_ENERGY} | awk -F':' '{ print $1 }'`
-E_END=`echo ${DFT_ENERGY} | awk -F':' '{ print $2 }'`
-
-## log file
-DFT_LOGFILE="genreg_`date '+%Y%m%d'`.log"
-
-## background spectra
-DFT_BKGD=`\ls bkgcorr_bl*.pi | head -n 1`
-## default parameters }}}
-
-## functions {{{
-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 }}}
-
-## check ciao init and set path to solve conflit with heasoft {{{
-if [ -z "${ASCDS_INSTALL}" ]; then
-    printf "ERROR: CIAO NOT initialized\n"
-    exit ${ERR_CIAO}
-fi
-
-## XXX: heasoft's `pget' etc. tools conflict with some CIAO tools
-printf "set \$PATH to avoid conflicts between HEAsoft and CIAO\n"
-export PATH="${ASCDS_BIN}:${ASCDS_CONTRIB}:${PATH}"
-printf "## PATH: ${PATH}\n"
-## check ciao&heasoft}}}
-
-## parameters {{{
-getopt_keyval "$@"
-
-##log file
-LOGFILE="not_exist"
-
-##check evt file
-
-if [ -r "${evt}" ]; then
-    EVT=${evt}
-elif [ -r "${DFT_EVT}" ]; then
-    EVT=${DFT_EVT}
-else
-    read -p "clean evt2 file: " EVT
-    if [ ! -r "${EVT}" ]; then
-        printf "ERROR: cannot access given \`${EVT}' evt file\n"
-        exit ${ERR_EVT}
-    fi
-fi
-
-##check ori sbp region
-if [ -r "${reg_in}" ]; then
-    REG_IN="${reg_in}"
-else
-    read -p "input region file: " REG_IN
-    if [ ! -r "${REG_IN}" ]; then
-        printf "ERROR: cannot access given \`${REG_IN}' evt file\n"
-        exit ${ERR_REG}
-    fi
-fi
-
-# check give basedir
-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/\/*$//'`
-
-#background spectrum
-if [ -r "${bkgd}" ] ;then
-    BKGD=${bkgd}
-elif [ -r "${DFT_BKGD}" ] ; then
-    BKGD="${DFT_BKGD}"
-    #ln -svf ${DFT_BKGD} .
-else
-    read -p ">background spectrum file: " BKGD
-    if [ ! -d ${BKGD} ] ; then
-        printf "ERROR on background spectrum file"
-        exit ${ERR_BKG}
-    fi
-fi
-
-## ds9: flag to determine whether to use ds9 check centroid and regions
-if [ ! -z "${ds9}" ]; then
-    case "${ds9}" in
-        [nN][oO]|[fF]*)
-            F_DS9="NO"
-            ;;
-        *)
-            F_DS9="YES"
-            ;;
-    esac
-else
-    F_DS9="YES"
-fi
-## parameters }}}
-
-## prepare parameter files (pfiles) {{{
-CIAO_TOOLS="dmkeypar dmcopy dmcoords"
-
-# 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 {{{
-# consistent with `ciao_procevt'
-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"
-    NEW_DETNAM="ACIS-0123"
-    ROOTNAME="c0-3_e${E_START}-${E_END}"
-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"
-    NEW_DETNAM="ACIS-7"
-    ROOTNAME="c7_e${E_START}-${E_END}"
-else
-    printf "ERROR: unknown detector type: ${DETNAM}\n"
-    exit ${ERR_DET}
-fi
-## ACIS type }}}
-
-## filter by energy band
-printf "filter out events in energy band: \`${E_START}:${E_END}' ...\n"
-EVT_E="evt2_${ROOTNAME}.fits"
-if [ ! -r "${EVT_E}" ]; then
-    punlearn dmcopy
-    dmcopy infile="${EVT}[energy=${E_START}:${E_END}]" outfile="${EVT_E}" clobber=yes
-fi
-
-printf "======== X-RAY CENTROID =======\n"
-CMD="${SCRIPT_DIR}/${XCENTROID_SCRIPT} evt=${EVT} reg=${REG_IN} conv=yes 2>&1 | tee xcentroid.dat"
-printf "CMD: $CMD\n"
-${SCRIPT_DIR}/${XCENTROID_SCRIPT} evt=${EVT} reg=${REG_IN} conv=yes 2>&1 | tee xcentroid.dat
-X=`\grep '(X,Y)' xcentroid.dat | tr -d ' XY():' | awk -F',' '{ print $2 }'`
-Y=`\grep '(X,Y)' xcentroid.dat | tr -d ' XY():' | awk -F',' '{ print $3 }'`
-CNTRD_WCS_REG="centroid_wcs.reg"
-CNTRD_PHY_REG="centroid_phy.reg"
-printf "## X centroid: ($X,$Y)\n"
-if [ "${F_DS9}" = "YES" ]; then
-    printf "check the X centroid ...\n"
-    ds9 ${EVT_E} -regions format ciao \
-        -regions system physical \
-        -regions ${CNTRD_PHY_REG} \
-        -cmap he -bin factor 4
-fi
-X0=$X
-Y0=$Y
-X=`\grep -i 'point' ${CNTRD_PHY_REG} | head -n 1 | tr -d 'a-zA-Z() ' | awk -F',' '{ print $1 }'`
-Y=`\grep -i 'point' ${CNTRD_PHY_REG} | head -n 1 | tr -d 'a-zA-Z() ' | awk -F',' '{ print $2 }'`
-if [ "x${X}" != "x${X0}" ] || [ "x${Y}" != "x${Y0}" ]; then
-    printf "## X CENTROID CHANGED -> ($X,$Y)\n"
-    # update ${CNTRD_WCS_REG}
-    printf "update ${CNTRD_WCS_REG} ...\n"
-    rm -f ${CNTRD_WCS_REG}
-    punlearn dmcoords
-    dmcoords infile="${EVT}" asolfile="${ASOL}" option=sky x=${X} y=${Y}
-    RA=`pget dmcoords ra`
-    DEC=`pget dmcoords dec`
-    echo "point(${RA},${DEC})" > ${CNTRD_WCS_REG}
-fi
-printf "======== X-RAY CENTROID FINISHED =======\n\n"
-
-SPC_REG=rspec.reg
-SBP_REG=sbprofile.reg
-
-printf "======== GENERATE SBPROFILE REGIONS =======\n"
-CMD="${SCRIPT_DIR}/${GEN_SBPREG_SCRIPT} ${EVT} ${EVT_E} ${X} ${Y} ${BKGD} ${SBP_REG}"
-printf "CMD: ${CMD}\n"
-${SCRIPT_DIR}/${GEN_SBPREG_SCRIPT} ${EVT} ${EVT_E} ${X} ${Y} ${BKGD} ${SBP_REG}
-if [ "${F_DS9}" = "YES" ]; then
-    printf "check SBP regions ...\n"
-    ds9 ${EVT_E} -regions format ciao \
-        -regions system physical \
-        -regions ${SBP_REG} \
-        -cmap he -bin factor 4
-fi
-printf "======== GENERATE SBPROFILE REGIONS FINISHED =======\n\n"
-
-printf "======== GENERATE SPECTRUM REGIONS =======\n"
-CMD="${SCRIPT_DIR}/${GEN_SPCREG_SCRIPT} ${EVT} ${EVT_E} ${BKGD} ${X} ${Y} ${SPC_REG}"
-printf "CMD: ${CMD}\n"
-${SCRIPT_DIR}/${GEN_SPCREG_SCRIPT} ${EVT} ${EVT_E} ${BKGD} ${X} ${Y} ${SPC_REG}
-if [ "${F_DS9}" = "YES" ]; then
-    printf "check SPC regions ...\n"
-    ds9 ${EVT_E} -regions format ciao \
-        -regions system physical \
-        -regions ${SPC_REG} \
-        -cmap he -bin factor 4
-fi
-printf "======== GENERATE SPECTRUM REGIONS FINISHED =======\n\n"
-
-exit 0