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
|
include <gset.h>
include "../lib/daophotdef.h"
include "../lib/psfdef.h"
define FRACTION 0.10
# DP_RADPSF -- Draw a radial profile plot of a data subraster containing a
# candidate psf star.
procedure dp_radpsf (dao, subras, ncols, nlines, x1, y1, title, gp)
pointer dao # pointer to DAOPHOT structure
real subras[ncols,nlines] # data subraster
int ncols, nlines # dimesnions of subraster
int x1, y1 # coordinates of left hand corner
char title[ARB] # title string
pointer gp # pointer to graphics descriptor
int npts
pointer psf, sp, radius, intensity, str
real ymin, ymax, r1, r2, i1, i2
int dp_rivectors()
begin
# Get the pointer to the DAOPHOT PSF fitting substructure.
psf = DP_PSF (dao)
# Allocate temporary space.
call smark (sp)
call salloc (radius, ncols * nlines, TY_REAL)
call salloc (intensity, ncols * nlines, TY_REAL)
call salloc (str, SZ_LINE, TY_CHAR)
# Compute the radial profile.
npts = dp_rivectors (subras, ncols, nlines, x1, y1, DP_CUR_PSFX(psf),
DP_CUR_PSFY(psf), DP_PSFRAD(dao), Memr[radius], Memr[intensity])
call alimr (Memr[intensity], npts, ymin, ymax)
# Make the plot.
call gclear (gp)
r1 = -FRACTION * DP_PSFRAD(dao)
r2 = DP_PSFRAD(dao) + FRACTION * DP_PSFRAD(dao)
i1 = ymin - FRACTION * (ymax - ymin)
i2 = ymax + FRACTION * (ymax - ymin)
call gswind (gp, r1, r2, i1, i2)
call glabax (gp, title, "Radius (pixels)", "Intensity (counts)")
call gpmark (gp, Memr[radius], Memr[intensity], npts, GM_PLUS, 1.0,
1.0)
# Mark the zero radius line.
call gamove (gp, 0.0, i1)
call gadraw (gp, 0.0, i2)
# Mark the sky level.
call gamove (gp, r1, DP_CUR_PSFSKY(psf))
call gadraw (gp, r2, DP_CUR_PSFSKY(psf))
# Mark the half-width at half-maximum.
call gamove (gp, DP_FWHMPSF(dao) / 2.0, i1)
call gadraw (gp, DP_FWHMPSF(dao) / 2.0, i2)
call sprintf (Memc[str], SZ_LINE, "Half-width half-maximum = %0.2f")
call pargr (DP_FWHMPSF(dao) / 2.0)
call gtext (gp, DP_FWHMPSF(dao) / 2.0, i2, Memc[str],
"q=h;u=180;v=t;p=r")
# Mark the fitting radius.
call gamove (gp, DP_FITRAD(dao), i1)
call gadraw (gp, DP_FITRAD(dao), i2)
call sprintf (Memc[str], SZ_LINE, "Fitting radius = %0.2f")
call pargr (DP_FITRAD(dao))
call gtext (gp, DP_FITRAD(dao), i2, Memc[str], "q=h;u=180;v=t;p=r")
call gflush (gp)
call sfree (sp)
end
|
