aboutsummaryrefslogtreecommitdiff
path: root/noao/digiphot/daophot/lib/allstardef.h
blob: 3f71aa280b67b8dbedb9fcdaa03564bfee54d914 (plain) (blame)
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
# ALLSTAR Structure

define	LEN_ALLSTARSTRUCT (80)

define	DP_ISCACHE	Memi[$1]	# is data cached (default no)
define	DP_ALLCACHE	Memi[$1+2]	# is all the data cached ? (yes)
define	DP_CACHE	Memi[$1+3+$2-1] # which data is cached ?
define	DP_SZCACHE	Memi[$1+6]	# current working set size
define	DP_SZOLDCACHE	Memi[$1+7]	# old working set size
define	DP_DATA		Memi[$1+8]	# pointer to data
define	DP_SUBT		Memi[$1+9]	# pointer to subt
define	DP_WEIGHTS	Memi[$1+10]	# pointer to weights

define	DP_SBUF		Memi[$1+13]	# local subt buffer
define	DP_SNX		Memi[$1+14]	# subt x dimension
define	DP_SNY		Memi[$1+15]	# subt y dimension
define	DP_SLX		Memi[$1+16]	# subt lower left x coord
define	DP_SMX		Memi[$1+17]	# subt lower right x coord
define	DP_SLY		Memi[$1+18]	# subt lower left y coord
define	DP_SMY		Memi[$1+19]	# subt lower right y coord
define	DP_SXOFF	Memi[$1+20]	# subt lower left x offset
define	DP_SYOFF	Memi[$1+21]	# subt lower left y offset

define	DP_WBUF		Memi[$1+23]	# local weight buffer
define	DP_WNX		Memi[$1+24]	# weight x dimension
define	DP_WNY		Memi[$1+25]	# weight y dimension
define	DP_WLX		Memi[$1+26]	# weight lower left x coord
define	DP_WMX		Memi[$1+27]	# weight lower right x coord
define	DP_WLY		Memi[$1+28]	# weight lower left y coord
define	DP_WMY		Memi[$1+29]	# weight lower right y coord
define	DP_WXOFF	Memi[$1+30]	# weight lower left x offset
define	DP_WYOFF	Memi[$1+31]	# weight lower left y offset

define	DP_DBUF		Memi[$1+33]	# local weight buffer
define	DP_DNX		Memi[$1+34]	# weight x dimension
define	DP_DNY		Memi[$1+35]	# weight y dimension
define	DP_DLX		Memi[$1+36]	# weight lower left x coord
define	DP_DMX		Memi[$1+37]	# weight lower right x coord
define	DP_DLY		Memi[$1+38]	# weight lower left y coord
define	DP_DMY		Memi[$1+39]	# weight lower right y coord
define	DP_DXOFF	Memi[$1+40]	# weight lower left x offset
define	DP_DYOFF	Memi[$1+41]	# weight lower left y offset

define	DP_ANUMER	Memi[$1+55]	# pointer to the anumer1 directory
define	DP_ADENOM	Memi[$1+57]	# pointer to the adenom1 directory
define	DP_ARPIXSQ	Memi[$1+59]	# pointer to the rpixsq directory
define	DP_ASUMWT	Memi[$1+60]	# pointer to the sumwt directory
define	DP_ASKIP	Memi[$1+61]	# pointer to the skip directory
define	DP_ALAST	Memi[$1+62]	# pointer to the last directory
define	DP_AXOLD	Memi[$1+63]	# pointer to the xold array
define	DP_AYOLD	Memi[$1+64]	# pointer to the yold array
define	DP_AXCLAMP	Memi[$1+65]	# pointer to the xclamp array
define	DP_AYCLAMP	Memi[$1+66]	# pointer to the yclamp array
define	DP_AX		Memi[$1+67]	# pointer to the x array
define	DP_AV		Memi[$1+68]	# pointer to the v array
define	DP_AC		Memi[$1+69]	# pointer to the c array
define	DP_ANPIX	Memi[$1+70]	# pointer to the npix array
define	DP_AIER		Memi[$1+71]	# pointer to the error array

define	A_SUBT		1
define	A_WEIGHT	2
define	A_DCOPY		3

# Separation criterion

define	FRACTION_MINSEP	0.14
define	CLAMP_FRACTION	0.25

# Limit on N/S ** 2 for two stars to merge.

define	WCRIT0		0.0
define	WCRIT5		1.0
define	WCRIT10		1.5
define	WCRIT15		2.0

# Minimum and maximum sharpness limits.

define	MIN_SHARP	-99.99
define	MAX_SHARP	99.99

# Number of output files columns

define	ALL_NOUTCOLUMN	11

# Define the fitting contstants and constraints

define	RADIUS_FRACTION	     0.95      # % decrease in radius for regrouping
define	DENSE_RADIUS1	     1.2       # 1st limiting radius for regrouping
define	DENSE_RADIUS2	     0.8       # 2nd limiting radius for regrouping
define	MIN_ITER	     4	       # minimum number of iterations
define	MAX_RELERR	     100.0     # maximum relative error	
define	MAX_RSQ		     0.999999  # maximum values of ratio of radii ** 2
define	MAX_RHOSQ	     36.0      # limit on core radius ** 2 for sharp	
define	CHI_NORM	     1.2533141 # sqrt (PI / 2.0)
define	MIN_SUMWT	     3.0       # min value for radial weight sum
define	MIN_NPIX	     4         # min pixels per star for a fit
define	MIN_XYCLAMP	     0.001     # min value of x-y clamp
define	MIN_XYCLAMP_FRACTION 0.5       # min change in value of x-y clamp
define	MAX_XYCLAMP_FRACTION 1.2       # max change in value of x-y clamp
define	MAX_DELTA_FAINTER    0.84      # max permitted brightness decrease
define	MAX_DELTA_BRIGHTER   5.25      # max permitted brightness increase
define	MAX_NEW_ERRMAG	     0.10      # 1st convergence check on mag error
define	MAX_NEW_RELBRIGHT2   0.0005    # 2nd convergenge check on mag
define	MAX_PIXERR2	     4.0e-6    # 2nd convergence check on x/y 
define	MIN_REL_BRIGHT	     1.0e-5    # min relative brightness

# List of ALLSTAR error codes

define	ALLERR_OK		0
define	ALLERR_BIGGROUP		1
define	ALLERR_INDEFSKY		2
define	ALLERR_NOPIX		3
define	ALLERR_SINGULAR		4
define	ALLERR_FAINT		5
define	ALLERR_MERGE		6
define	ALLERR_OFFIMAGE		7