blob: e5230991dfd848fdaa9c7ba7536c355ce20daf18 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
|
#!/bin/sh
#
# Copyright (c) 2016 Aaron LI
# MIT license
#
# Generate the exposure map and apply exposure correction
# using CIAO `fluximage`.
#
# NOTE:
# The existing "instmap_weights.txt" is ued for exposure map generation.
#
# Created: 2016-04-28
#
case "$1" in
-[hH]*)
echo "Usage:"
echo " `basename $0` <evt_e> <img> [ instmap_weights.txt ]"
exit 1
;;
esac
EVT_E="$1"
IMG="$2"
SPEC_WGT="${3:-instmap_weights.txt}"
REPRO_DIR=".."
ASOLIS=`\ls ${REPRO_DIR}/acisf*asol?.lis`
BPIX=`\ls ${REPRO_DIR}/acisf*repro_bpix?.fits`
MSK=`\ls ${REPRO_DIR}/acisf*msk?.fits`
ROOTNAME=`echo "${EVT_E%.fits}" | sed -e 's/^evt2_//'`
## get `xygrid' for image
punlearn get_sky_limits
get_sky_limits image="${IMG}" verbose=0
XYGRID=`pget get_sky_limits xygrid`
echo "xygrid: ${XYGRID}"
punlearn ardlib
echo "invoking fluximage to generate expmap and apply correction ..."
punlearn fluximage
fluximage infile="${EVT_E}" outroot="${ROOTNAME}" \
binsize=1 bands="${SPEC_WGT}" xygrid="${XYGRID}" \
asol="@${ASOLIS}" badpixfile="${BPIX}" \
maskfile="${MSK}" clobber=yes
# make symbolic links
# clipped counts image
ln -svf ${ROOTNAME}*band*thresh.img img_${ROOTNAME}_thresh.fits
# clipped exposure map
ln -svf ${ROOTNAME}*band*thresh.expmap expmap_${ROOTNAME}.fits
# exposure-corrected image
ln -svf ${ROOTNAME}*band*flux.img img_expcorr_${ROOTNAME}.fits
|
