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
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
|
# SETIMPARS -- Initialize the DAOPHOT task and pset parameters.
procedure setimpars (image, restore, update)
string image {prompt="Image name"}
bool restore {yes, prompt="Restore the last saved parameter set ?"}
bool update {yes, prompt="Update the last saved parameter set ?"}
bool review {no,
prompt="Review the parameters before saving ?"}
file parfile {"", prompt="Input algorithm parameters file"}
file datapars {"", prompt="The input data dependent parameters file"}
file findpars {"",
prompt="The input object detection parameters file"}
file centerpars {"", prompt="The input centering parameters file"}
file fitskypars {"", prompt="The input sky fitting parameters file"}
file photpars {"",
prompt="The input aperture photometry parameters file"}
file daopars {"", prompt="The input psf fitting parameters file"}
bool unlearn {no,
prompt="Unlearn the current algorithm parameters ?"}
begin
# Define some temporary variables.
bool trestore, tupdate
string timage, tinparfile, toutparfile, tpars, tfile1, tfile2
# Read in the image name.
timage = image
trestore = restore
tupdate = update
print ("Setting parameters for image ", timage, " ...")
# Set the input image name.
addstar.image = timage
allstar.image = timage
daoedit.image = timage
daofind.image = timage
group.image = timage
nstar.image = timage
peak.image = timage
pexamine.image = timage
phot.image = timage
psf.image = timage
pstselect.image = timage
substar.image = timage
# Set the input coordinate / sky files back to their defaults.
phot.skyfile = ""
# Set the input photometry files back to their defaults.
addstar.photfile = ""
addstar.simple_text = no
allstar.photfile = "default"
group.photfile = "default"
grpselect.ingroupfile = ""
nstar.groupfile = "default"
peak.photfile = "default"
phot.coords = "default"
psf.photfile = "default"
psf.pstfile = ""
pstselect.photfile = "default"
substar.photfile = "default"
# Set the psfimage back to the default.
addstar.psfimage = "default"
allstar.psfimage = "default"
group.psfimage = "default"
nstar.psfimage = "default"
peak.psfimage = "default"
seepsf.psfimage = ""
substar.psfimage = "default"
# Set the output photometry file names to the default.
allstar.allstarfile = "default"
allstar.rejfile = "default"
daofind.output = "default"
group.groupfile = "default"
grpselect.outgroupfile = ""
nstar.nstarfile = "default"
nstar.rejfile = "default"
peak.peakfile = "default"
peak.rejfile = "default"
phot.output = "default"
psf.groupfile = "default"
psf.opstfile = "default"
pstselect.pstfile = "default"
# Set the output images back to the default.
addstar.addimage = "default"
allstar.subimage = "default"
daofind.starmap = ""
daofind.skymap = ""
psf.psfimage = "default"
seepsf.image = ""
substar.subimage = "default"
# Set any output plot files back to the default.
phot.plotfile = ""
pstselect.plotfile = ""
psf.plotfile = ""
# Get the input parameter file name.
tinparfile = parfile
if (tinparfile == "") {
if (access (timage // ".pars") && trestore)
tinparfile = timage // ".pars"
} else if (! access (tinparfile)) {
print ("File ", tinparfile, " does not exist ...")
return
}
# Read in the input parameters.
if (tinparfile != "") {
print ("Reading algorithm parameters from file ",
tinparfile, " ...")
cl (< tinparfile)
} else {
tpars = datapars
if (access (tpars)) {
print ("Reading datapars parameters from file ", tpars, " ...")
tfile1 = mktemp ("tmp$pars")
tfile2 = mktemp ("tmp$pars")
dparam (tpars, > tfile1)
list = tfile1
while (fscan (list, line) != EOF) {
if (substr (line, 1, 5) != "# EOF")
print ("datapars.", line, >> tfile2)
}
cl (< tfile2)
delete (tfile1 // "," // tfile2, verify-, default_action+,
allversions+, subfiles+, go_ahead+, >& "dev$null")
} else if (unlearn) {
print ("Reading default datapars parameters from disk ...")
unlearn ("daophot.datapars")
} else {
print ("Reading current datapars parameters from disk ...")
}
tpars = findpars
if (access (tpars)) {
print ("Reading findpars parameters from file ", tpars, " ...")
tfile1 = mktemp ("tmp$pars")
tfile2 = mktemp ("tmp$pars")
dparam (tpars, > tfile1)
list = tfile1
while (fscan (list, line) != EOF) {
if (substr (line, 1, 5) != "# EOF")
print ("findpars.", line, >> tfile2)
}
cl (< tfile2)
delete (tfile1 // "," // tfile2, verify-, default_action+,
allversions+, subfiles+, go_ahead+, >& "dev$null")
} else if (unlearn) {
print ("Reading default findpars parameters from disk ...")
unlearn ("daophot.findpars")
} else {
print ("Reading current findpars parameters from disk ...")
}
tpars = centerpars
if (access (tpars)) {
print ("Reading centerpars parameters from file ", tpars,
" ...")
tfile1 = mktemp ("tmp$pars")
tfile2 = mktemp ("tmp$pars")
dparam (tpars, > tfile1)
list = tfile1
while (fscan (list, line) != EOF) {
if (substr (line, 1, 5) != "# EOF")
print ("centerpars.", line, >> tfile2)
}
cl (< tfile2)
delete (tfile1 // "," // tfile2, verify-, default_action+,
allversions+, subfiles+, go_ahead+, >& "dev$null")
} else if (unlearn) {
unlearn ("daophot.centerpars")
print ("Reading default centerpars parameters from disk ...")
} else {
print ("Reading current centerpars parameters from disk ...")
}
tpars = fitskypars
if (access (tpars)) {
print ("Reading fitskypars parameters from file ", tpars,
" ...")
tfile1 = mktemp ("tmp$pars")
tfile2 = mktemp ("tmp$pars")
dparam (tpars, > tfile1)
list = tfile1
while (fscan (list, line) != EOF) {
if (substr (line, 1, 5) != "# EOF")
print ("fitskypars.", line, >> tfile2)
}
cl (< tfile2)
delete (tfile1 // "," // tfile2, verify-, default_action+,
allversions+, subfiles+, go_ahead+, >& "dev$null")
} else if (unlearn) {
unlearn ("daophot.fitskypars")
print ("Reading default fitskypars parameters from disk ...")
} else {
print ("Reading current fitskypars parameters from disk ...")
}
tpars = photpars
if (access (tpars)) {
print ("Reading photpars parameters from file ", tpars,
" ...")
tfile1 = mktemp ("tmp$pars")
tfile2 = mktemp ("tmp$pars")
dparam (tpars, > tfile1)
list = tfile1
while (fscan (list, line) != EOF) {
if (substr (line, 1, 5) != "# EOF")
print ("photpars.", line, >> tfile2)
}
cl (< tfile2)
delete (tfile1 // "," // tfile2, verify-, default_action+,
allversions+, subfiles+, go_ahead+, >& "dev$null")
} else if (unlearn) {
unlearn ("daophot.photpars")
print ("Reading default photpars parameters from disk ...")
} else {
print ("Reading current photpars parameters from disk ...")
}
tpars = daopars
if (access (tpars)) {
print ("Reading psf fitting parameters from file ", tpars,
" ...")
tfile1 = mktemp ("tmp$pars")
tfile2 = mktemp ("tmp$pars")
dparam (tpars, > tfile1)
list = tfile1
while (fscan (list, line) != EOF) {
if (substr (line, 1, 5) != "# EOF")
print ("daopars.", line, >> tfile2)
}
cl (< tfile2)
delete (tfile1 // "," // tfile2, verify-, default_action+,
allversions+, subfiles+, go_ahead+, >& "dev$null")
} else if (unlearn) {
unlearn ("daophot.daopars")
print ("Reading default daopars parameters from disk ...")
} else {
print ("Reading current daopars parameters from disk ...")
}
}
# Review the current values for the algorithm parameters.
if (review) {
eparam ("datapars")
eparam ("findpars")
eparam ("centerpars")
eparam ("fitskypars")
eparam ("photpars")
eparam ("daopars")
}
# Update the output parameter file.
toutparfile = timage // ".pars"
if (tupdate) {
if (access (toutparfile)) {
print ("Updating image parameter file ", toutparfile, " ...")
delete (toutparfile, verify-, default_action+, allversions+,
subfiles+, go_ahead+, >& "dev$null")
} else {
print ("Creating image parameter file ", toutparfile, " ...")
}
dparam ("datapars", "findpars", "centerpars", "fitskypars",
"photpars", "daopars", > toutparfile)
}
end
|
