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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <mach.h>
include <imhdr.h>
include <imio.h>
include <fio.h>
include <qpioset.h>
include "qpf.h"
# ZFIOQP -- QPF virtual file driver. This driver presents to the caller a
# virtual file space containing a two dimensional array of type short or int
# pixels, wherein each "pixel" is a count of the number of events from a
# QPOE event list which map into that pixel. An i/o request results in
# runtime filtering and sampling of the event list, mapping each event which
# passes the filter into the corresponding output pixel, and incrementing the
# value of that pixel to count the event.
# QPFZOP -- Open the file driver for i/o on the QPIO descriptor opened at
# qpf_open time.
procedure qpfzop (pkfn, mode, status)
char pkfn[ARB] #I packed virtual filename from FIO
int mode #I file access mode (ignored)
int status #O output status - i/o channel if successful
int ip
pointer sp, fn, qpf
int ctoi()
begin
call smark (sp)
call salloc (fn, SZ_FNAME, TY_CHAR)
# The QPF descriptor is passed encoded in the pseudo filename as
# "QPFxxxx" (decimal). Extract this and return it as the i/o
# channel for the driver.
ip = 4
call strupk (pkfn, Memc[fn], SZ_FNAME)
if (ctoi (Memc[fn], ip, qpf) <= 0)
status = ERR
else
status = qpf
QPF_IOSTAT(qpf) = 0
call sfree (sp)
end
# QPFZCL -- Close the QPF binary file driver.
procedure qpfzcl (chan, status)
int chan #I QPF i/o channel
int status #O output status
begin
status = OK
end
# QPFZRD -- Read a segment of the virtual pixel array into the output buffer,
# i.e., zero the output buffer and sample the event list, accumulating counts
# in the output array.
procedure qpfzrd (chan, obuf, nbytes, boffset)
int chan #I QPF i/o channel
char obuf[ARB] #O output buffer
int nbytes #I nbytes to be read
int boffset #I file offset at which read commences
pointer qpf, im, io
int vs[2], ve[2]
real xblock, yblock
int szb_pixel, ncols, pixel, nev, xoff, yoff
int qpio_readpixs(), qpio_readpixi()
include <szpixtype.inc>
begin
qpf = chan
im = QPF_IM(qpf)
io = QPF_IO(qpf)
xblock = QPF_XBLOCK(qpf)
yblock = QPF_YBLOCK(qpf)
ncols = IM_PHYSLEN(im,1)
xoff = QPF_VS(qpf,1)
yoff = QPF_VS(qpf,2)
szb_pixel = pix_size[IM_PIXTYPE(im)] * SZB_CHAR
# Convert boffset, nbytes to vs, ve.
pixel = (boffset - 1) / szb_pixel
vs[1] = (mod (pixel, ncols)) * xblock + xoff
vs[2] = (pixel / ncols) * yblock + yoff
pixel = (boffset-1 + nbytes - szb_pixel) / szb_pixel
ve[1] = (mod (pixel, ncols)) * xblock + (xblock-1) + xoff
ve[2] = (pixel / ncols) * yblock + (yblock-1) + yoff
# Call readpix to sample image into the output buffer. Zero the buffer
# first since the read is additive.
call aclrc (obuf, nbytes / SZB_CHAR)
iferr {
switch (IM_PIXTYPE(im)) {
case TY_SHORT:
nev = qpio_readpixs (io, obuf, vs, ve, 2, xblock, yblock)
case TY_INT:
nev = qpio_readpixi (io, obuf, vs, ve, 2, xblock, yblock)
}
} then {
QPF_IOSTAT(qpf) = ERR
} else
QPF_IOSTAT(qpf) = nbytes
end
# QPFZWR -- Write to the virtual pixel array. QPF permits only read-only
# access, but we ignore write requests, so return OK and do nothing if this
# routine is called.
procedure qpfzwr (chan, ibuf, nbytes, boffset)
int chan #I QPF i/o channel
char ibuf[ARB] #O datg buffer
int nbytes #I nbytes to be written
int boffset #I file offset to write at
pointer qpf
begin
qpf = chan
QPF_IOSTAT(qpf) = nbytes
end
# QPFZWT -- Return the number of virtual bytes transferred in the last i/o
# request.
procedure qpfzwt (chan, status)
int chan #I QPF i/o channel
int status #O i/o channel status
pointer qpf
begin
qpf = chan
status = QPF_IOSTAT(qpf)
end
# QPFZST -- Query device/file parameters.
procedure qpfzst (chan, param, value)
int chan #I QPF i/o channel
int param #I parameter to be returned
int value #O parameter value
pointer qpf, im, io
int szb_pixel, npix
int qpio_stati()
include <szpixtype.inc>
begin
qpf = chan
im = QPF_IM(qpf)
io = QPF_IO(qpf)
npix = IM_PHYSLEN(im,1) * IM_PHYSLEN(im,2)
szb_pixel = pix_size[IM_PIXTYPE(im)] * SZB_CHAR
switch (param) {
case FSTT_BLKSIZE:
value = 1
case FSTT_FILSIZE:
value = npix * szb_pixel
case FSTT_OPTBUFSIZE:
value = min (npix*szb_pixel, qpio_stati(io,QPIO_OPTBUFSIZE))
case FSTT_MAXBUFSIZE:
value = npix * szb_pixel
default:
value = ERR
}
end
|
