-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathextract-spec-esas-pn.sh
executable file
·100 lines (80 loc) · 2.88 KB
/
extract-spec-esas-pn.sh
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
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
######################################################################
# runs spectral extraction in the esas bands
#
# Snowden & Kuntz, 2011:
#
# pn-spectra and pn-spectra processes the filtered event files to
# produce background spectra for the entire energy range and selected
# region and background images for the selected region and selected
# band for the individual ccds. The region selection expression is in
# an input file and should be in detector coordinates. If the input
# file does not exist, reg.txt in this case, the default is to process
# the entire FOV. The input energies are in eV. Change the caldb
# directory string to whatever is appropriate. In this step the
# processing is to create images in two bands for the instruments.
dir=$1
here=`pwd`
cd $dir
######################################################################
# extract pn spectra - soft band
prefix=$PN_EV_PREFIX_LIST
elow='400' # detection bands minima [eV]
ehigh='1250' # detection bands maxima [eV]
regfile=$PN_SRC_REGFILE
pattern=4
quad1=$PN_QUAD1
quad2=$PN_QUAD2
quad3=$PN_QUAD3
quad4=$PN_QUAD4
pn-spectra prefix="$prefix" caldb=$esas_caldb region=$regfile mask=1 elow=$elow ehigh=$ehigh pattern=4 quad1=$quad1 quad2=$quad2 quad3=$quad3 quad4=$quad4
if [[ $? -ne 0 ]]
then
echo -e "\n** error: pn spectral extraction in soft band failed!"
echo -e "*** error in script: $0\n"
cd $startdir
exit 1
fi
######################################################################
# extract pn spectra - hard band
prefix=$PN_EV_PREFIX_LIST
elow='2000' # detection bands minima [eV]
ehigh='7200' # detection bands maxima [eV]
regfile=$PN_SRC_REGFILE
pattern=4
quad1=$PN_QUAD1
quad2=$PN_QUAD2
quad3=$PN_QUAD3
quad4=$PN_QUAD4
pn-spectra prefix="$prefix" caldb=$esas_caldb region=$regfile mask=1 elow=$elow ehigh=$ehigh pattern=4 quad1=$quad1 quad2=$quad2 quad3=$quad3 quad4=$quad4
if [[ $? -ne 0 ]]
then
echo -e "\n** error: pn spectral extraction in hard band failed!"
echo -e "*** error in script: $0\n"
cd $startdir
exit 1
fi
# # part refactored to its own script in extract-spec-band
# ######################################################################
# # extract pn spectra - standard band
# prefix=$PN_EV_PREFIX_LIST
# elow='500' # detection bands minima [eV]
# ehigh='2000' # detection bands maxima [eV]
# regfile=$PN_SRC_REGFILE
# pattern=4
# quad1=$PN_QUAD1
# quad2=$PN_QUAD2
# quad3=$PN_QUAD3
# quad4=$PN_QUAD4
# pn-spectra prefix="$prefix" caldb=$esas_caldb region=$regfile mask=1 elow=$elow ehigh=$ehigh pattern=4 quad1=$quad1 quad2=$quad2 quad3=$quad3 quad4=$quad4
# if [[ $? -ne 0 ]]
# then
# echo -e "\n** error: pn spectral extraction in standard band failed!"
# echo -e "*** error in script: $0\n"
# cd $startdir
# exit 1
# fi
######################################################################
# exit
cd $here
echo -e "\n$0 in $obsid done!"
exit 0