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
|
# DO3FIBERS -- Process Coude fiber spectra from 2D to wavelength calibrated 1D.
#
# The task PROC does all of the interactive work and BATCH does the
# background work. This procedure is organized this way to minimize the
# dictionary space when the background task is submitted.
procedure do3fiber (objects)
string objects = "" {prompt="List of object spectra"}
file apref = "" {prompt="Aperture reference spectrum"}
file flat = "" {prompt="Flat field spectrum"}
string arcs = "" {prompt="List of arc spectra"}
file arctable = "" {prompt="Arc assignment table (optional)\n"}
string readnoise = "RDNOISE" {prompt="Read out noise sigma (photons)"}
string gain = "GAIN" {prompt="Photon gain (photons/data number)"}
real datamax = INDEF {prompt="Max data value / cosmic ray threshold"}
int fibers = 3 {prompt="Number of fibers"}
real width = 6. {prompt="Width of profiles (pixels)"}
string crval = "INDEF" {prompt="Approximate central wavelength"}
string cdelt = "INDEF" {prompt="Approximate dispersion"}
string objaps = "2" {prompt="Object apertures"}
string arcaps = "1,3" {prompt="Arc apertures\n"}
bool scattered = no {prompt="Subtract scattered light?"}
bool fitflat = yes {prompt="Fit and ratio flat field spectrum?"}
bool recenter = yes {prompt="Recenter object apertures?"}
bool edit = no {prompt="Edit/review object apertures?"}
bool clean = no {prompt="Detect and replace bad pixels?"}
bool dispcor = yes {prompt="Dispersion correct spectra?"}
bool splot = yes {prompt="Plot the final spectrum?"}
bool redo = no {prompt="Redo operations if previously done?"}
bool update = yes {prompt="Update spectra if cal data changes?"}
bool batch = no {prompt="Extract objects in batch?"}
bool listonly = no {prompt="List steps but don't process?\n"}
pset params = "" {prompt="Algorithm parameters"}
begin
apscript.readnoise = readnoise
apscript.gain = gain
apscript.nfind = fibers
apscript.width = width
apscript.t_width = width
apscript.radius = width
apscript.clean = clean
apscript.order = "increasing"
proc.datamax = datamax
proc (objects, apref, flat, "", arcs, "", "",
arctable, fibers, "", crval, cdelt, objaps, "", arcaps, "",
"", "", scattered, fitflat, recenter, edit, no, no, clean,
dispcor, no, no, no, no, no, splot, redo, update, batch, listonly)
if (proc.dobatch) {
print ("-- Do remaining spectra as a batch job --")
print ("batch&batch") | cl
}
end
|
