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
|
# SDOARCS -- Determine dispersion relation for spectrum based on reference arcs.
procedure sdoarcs (spec, arcref, reextract, arcap, logfile, batch)
file spec
file arcref
bool reextract
bool arcap
file logfile
bool batch
struct *fd
begin
string imtype, ectype
int i, j, k, n
file temp, arc1, arc2, str1, str2, arctype, apref, arc, arcec, logs
file specec, specarc
bool verbose1
imtype = "." // envget ("imtype")
i = stridx (",", imtype)
if (i > 0)
imtype = substr (imtype, 1, i-1)
ectype = ".ec" // imtype
n = strlen (imtype)
temp = mktemp ("tmp$iraf")
if (batch)
verbose1 = no
else
verbose1 = verbose
if (verbose1)
logs = logfile//",STDOUT"
else
logs = logfile
for (j=1; j<=2; j+=1) {
# The reference spectra refer initially to the 2D image. At the
# end we will reset them to refer to the 1D spectra.
hselect (spec, "refspec"//j, yes, > temp)
fd = temp
k = fscan (fd, arc1, str1)
fd = ""; delete (temp, verify=no)
if (k < 1)
break
# Strip possible image extension.
i = strlen (arc1)
if (i > n && substr (arc1, i-n+1, i) == imtype)
arc1 = substr (arc1, 1, i-n)
# Set extraction output and aperture reference depending on whether
# the arcs are to be rextracted using recentered or retraced object
# apertures.
if (arcap) {
arc2 = spec // arc1
apref = spec
if (access (arc2//ectype))
imdelete (arc2//ectype, verify=no)
delete (database//"/ec"//arc2//".ec*", verify = no)
} else {
arc2 = arc1
apref = apslitproc.references
if (reextract && access (arc2//ectype)) {
if (arc2 != arcref)
imdelete (arc2//ectype, verify=no)
}
}
# Extract and determine dispersion function if necessary.
if (!access (arc2//ectype)) {
delete (database//"/ec"//arc2//".ec*", verify = no)
if (!batch)
print ("Extract and reidentify arc spectrum ", arc1)
print ("Extract and reidentify arc spectrum ", arc1, >> logfile)
apslitproc (arc1, output=arc2//".ec", references=apref,
background="none", clean=no, weights="none",
verbose=verbose1)
ecreidentify (arc2//".ec", arcref//".ec", shift=0.,
cradius=sparams.cradius, threshold=sparams.threshold,
refit=sparams.refit, database=database, logfiles=logs)
# If not reextracting arcs based on object apertures
# then save the extracted arc to avoid doing it again.
if (arc1 != arc2)
imdelete (arc2//".ec", verify=no)
}
# Set the REFSPEC parameters for echelle spectrum.
if (k == 1)
hedit (spec//".ec", "refspec"//j, arc2//".ec", add=yes,
verify=no, show=no, update=yes)
else
hedit (spec//".ec", "refspec"//j, arc2//".ec "//str1, add=yes,
verify=no, show=no, update=yes)
}
end
|
