blob: 3e2d8817bb312e0fa72adfcbf39820b4af4cc81e (
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
|
#!/bin/sh
#
# Make an image from evt2 file by binning the 'TDET' coordinates.
# If the 'TDET' coordinate information is not available (e.g., blanksky),
# the 'DET' coordinate is used for binning.
#
# Aaron LI
# 2015/08/21
#
if [ $# -ne 2 ]; then
printf "Usage:\n"
printf " `basename $0` <input:evt2> <output:img>\n"
exit 1
fi
EVT2="$1"
OUTIMG="$2"
# Test whether 'TEDT' coordinate exists?
if dmlist ${EVT2} cols | grep -q tdet; then
COORDX="tdetx"
COORDY="tdety"
else
printf "WARNING: tdet coordinate NOT exist! use det instead.\n"
COORDX="detx"
COORDY="dety"
fi
# Get TDET/DET coordinate min & max values
# COORDX:
punlearn dmstat
dmstat "${EVT2}[cols ${COORDX}]" >/dev/null 2>&1
XMIN=`pget dmstat out_min`
XMAX=`pget dmstat out_max`
dmstat "${EVT2}[cols ${COORDY}]" >/dev/null 2>&1
YMIN=`pget dmstat out_min`
YMAX=`pget dmstat out_max`
BINSPEC="[bin ${COORDX}=${XMIN}:${XMAX}:1,${COORDY}=${YMIN}:${YMAX}:1]"
punlearn dmcopy
CMD="dmcopy \"${EVT2}${BINSPEC}\" ${OUTIMG}"
printf "CMD: ${CMD}\n"
eval ${CMD}
|
