aboutsummaryrefslogtreecommitdiffstats
path: root/scripts
diff options
context:
space:
mode:
authorAaron LI <aaronly.me@outlook.com>2017-02-18 16:28:35 +0800
committerAaron LI <aaronly.me@outlook.com>2017-02-18 16:28:35 +0800
commit89adbc474890c50963537e8e32bdeff8aa3e10d1 (patch)
treed3c0b63ca7467af8a6530acff55cbf735eb05030 /scripts
parent7dfd1127a5a63863de5bbfbab26bd29d482dec8d (diff)
downloadchandra-acis-analysis-89adbc474890c50963537e8e32bdeff8aa3e10d1.tar.bz2
Remove unused ciao_check_offset.sh
Diffstat (limited to 'scripts')
-rwxr-xr-xscripts/ciao_check_offset.sh269
1 files changed, 0 insertions, 269 deletions
diff --git a/scripts/ciao_check_offset.sh b/scripts/ciao_check_offset.sh
deleted file mode 100755
index c57bbaf..0000000
--- a/scripts/ciao_check_offset.sh
+++ /dev/null
@@ -1,269 +0,0 @@
-#!/bin/sh
-#
-unalias -a
-export LC_COLLATE=C
-###########################################################
-## based on `ciao_update_xcentroid.sh' ##
-## get `xcentroid' from region `sbprofile.reg' ##
-## convert from physical coords to WCS corrds ##
-## calculate offset with chandra pointing coords ##
-## ##
-## Weitian LI <liweitianux@gmail.com> ##
-## 2014/06/25 ##
-###########################################################
-##
-VERSION="v2.0"
-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.
-## * Replace 'grep' with '\grep', 'ls' with '\ls'
-##
-
-## error code {{{
-ERR_USG=1
-ERR_DIR=11
-ERR_EVT=12
-ERR_BKG=13
-ERR_REG=14
-ERR_INFO=15
-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, help {{{
-case "$1" in
- -[hH]*|--[hH]*)
- printf "usage:\n"
- printf " `basename $0` evt=<evt_file> reg=<sbp_reg> basedir=<base_dir>\n"
- printf "\nversion:\n"
- printf " ${VERSION}, ${UPDATED}\n"
- exit ${ERR_USG}
- ;;
-esac
-## usage, help }}}
-
-## default parameters {{{
-# critical offset values (unit: dd:mm:ss)
-OFF_CRIT_S="0:3:0" # 0deg 3min 0sec
-OFF_CRIT_I="0:7:0" # 0deg 7min 0sec
-
-# default `event file' which used to match `blanksky' files
-#DFT_EVT="_NOT_EXIST_"
-DFT_EVT="`\ls evt2*_clean.fits 2> /dev/null`"
-# default dir which contains `asols, asol.lis, ...' files
-# DFT_BASEDIR="_NOT_EXIST_"
-DFT_BASEDIR=".."
-# default `radial region file' to extract surface brightness
-#DFT_SBP_REG="_NOT_EXIST_"
-DFT_SBP_REG="sbprofile.reg"
-
-## howto find files in `basedir'
-# default `asol.lis pattern'
-DFT_ASOLIS_PAT="acis*asol?.lis"
-## 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
-}
-
-## unit/format conversion
-ddmmss2deg() {
- # convert 'dd:mm:ss' to 'deg'
- echo "$1" | awk -F':' '
- function abs(x) { return ((x<0.0) ? (-x) : x) }
- function sign(x) { return ((x<0.0) ? (-1.0) : 1.0) }
- {
- value = abs($1) + ($2)/60.0 + ($3)/3600.0;
- printf("%.8f", sign($1)*value);
- }'
-}
-
-deg2ddmmss() {
- # convert 'deg' to 'dd:mm:ss'
- echo "$1" | awk '
- function abs(x) { return ((x<0.0) ? (-x) : x) }
- {
- deg = $1;
- dd = int(deg);
- mm = int(abs(deg-dd)*60.0);
- ss = (abs(deg-dd)*60.0 - mm)*60.0;
- printf("%d:%d:%.2f", dd, mm, ss);
- }'
-}
-## functions }}}
-
-## check CIAO init {{{
-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}"
-## check CIAO }}}
-
-## 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 "clean 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" | ${TOLOG}
-# check given region file(s)
-if [ -r "${reg}" ]; then
- SBP_REG="${reg}"
-elif [ -r "${DFT_SBP_REG}" ]; then
- SBP_REG=${DFT_SBP_REG}
-else
- read -p "> surface brighness radial region file: " SBP_REG
- if [ ! -r "${SBP_REG}" ]; then
- printf "ERROR: cannot access given \`${SBP_REG}' region file\n"
- exit ${ERR_REG}
- fi
-fi
-printf "## use reg file(s): \`${SBP_REG}'\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}
-## parameters }}}
-
-## prepare parameter files (pfiles) {{{
-CIAO_TOOLS="dmkeypar 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 }}}
-
-## main process {{{
-# asolis
-ASOLIS=`( cd ${BASEDIR} && \ls ${DFT_ASOLIS_PAT} 2> /dev/null )`
-
-# get (x,y) from sbp region
-printf "get (x,y) from ${SBP_REG}\n"
-X=`\grep -iE '(pie|annulus)' ${SBP_REG} | head -n 1 | awk -F',' '{ print $1 }' | tr -d 'a-zA-Z() '`
-Y=`\grep -iE '(pie|annulus)' ${SBP_REG} | head -n 1 | awk -F',' '{ print $2 }' | tr -d 'a-zA-Z() '`
-
-# dmcoords to convert (x,y) to (ra,dec)
-printf "\`dmcoords' to convert (x,y) to (ra,dec) ...\n"
-punlearn dmcoords
-dmcoords infile="${EVT}" asolfile="@${BASEDIR}/${ASOLIS}" option=sky celfmt=deg x=${X} y=${Y}
-RA=`pget dmcoords ra`
-DEC=`pget dmcoords dec`
-
-# get observation pointing coordinates
-punlearn dmkeypar
-RA_PNT=`dmkeypar infile="${EVT}" keyword="RA_PNT" echo=yes`
-DEC_PNT=`dmkeypar infile="${EVT}" keyword="DEC_PNT" echo=yes`
-
-## 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"
-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"
-else
- printf "ERROR: unknown detector type: ${DETNAM}\n"
- exit ${ERR_DET}
-fi
-## ACIS type }}}
-
-# calculate offset
-OFF_DEG=`echo "${RA} ${RA_PNT} ${DEC} ${DEC_PNT}" | awk '{ printf("%.8f", sqrt(($1-$2)^2 + ($3-$4)^2)) }'`
-OFF_DDMMSS=`deg2ddmmss ${OFF_DEG}`
-
-# compare offset with the given critical values
-if [ "${ACIS_TYPE}" = "ACIS-S" ]; then
- DEG_OFF_CRIT=`ddmmss2deg ${OFF_CRIT_S}`
-else
- DEG_OFF_CRIT=`ddmmss2deg ${OFF_CRIT_I}`
-fi
-if [ `echo "${OFF_DEG} > ${DEG_OFF_CRIT}" | bc -l` -eq 1 ]; then
- LARGE_OFFSET="YES"
-else
- LARGE_OFFSET="NO"
-fi
-
-# output results
-printf "###################################################\n"
-printf "## Detector type: ${ACIS_TYPE}\n"
-printf "## Critical offset values: (format: dd:mm:ss)\n"
-printf "## OFF_CRIT_S: ${OFF_CRIT_S}\n"
-printf "## OFF_CRIT_I: ${OFF_CRIT_I}\n"
-printf "##\n"
-printf "## Our Results:\n"
-printf "## Chandra sky coordinates:\n"
-printf "## (x,y): ($X,$Y)\n"
-printf "## WCS coordinates in 'deg' unit:\n"
-printf "## (ra,dec): (${RA},${DEC})\n"
-printf "##\n"
-printf "## Observation pointing coordinates:\n"
-printf "## (ra,dec): (${RA_PNT},${DEC_PNT})\n"
-printf "##\n"
-printf "## Offset (format: dd:mm:ss):\n"
-printf "## offset: ${OFF_DDMMSS}\n"
-printf "## offset_deg: ${OFF_DEG}\n"
-if [ "${LARGE_OFFSET}" = "YES" ]; then
-printf "## *** WARNING: LARGE OFFSET ***\n"
-fi
-printf "###################################################\n"
-
-## main }}}
-
-exit 0
-