Changeset 6715 for trunk/psphot/src/psphotEvalPSF.c
- Timestamp:
- Mar 27, 2006, 4:31:17 PM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotEvalPSF.c (modified) (16 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotEvalPSF.c
r5993 r6715 54 54 // do we actually have a valid PSF model? 55 55 if (model == NULL) { 56 source->mode &= ~PM_SOURCE_ FITTED;56 source->mode &= ~PM_SOURCE_MODE_FITTED; 57 57 return false; 58 58 } … … 63 63 break; 64 64 case PM_MODEL_UNTRIED: 65 source->mode &= ~PM_SOURCE_ FITTED;65 source->mode &= ~PM_SOURCE_MODE_FITTED; 66 66 return false; 67 67 case PM_MODEL_BADARGS: … … 69 69 case PM_MODEL_OFFIMAGE: 70 70 default: 71 source->mode |= PM_SOURCE_ FAIL;71 source->mode |= PM_SOURCE_MODE_FAIL; 72 72 return false; 73 73 } 74 74 75 75 // unless we prove otherwise, this object is a star. 76 source->type = PM_SOURCE_ STAR;76 source->type = PM_SOURCE_TYPE_STAR; 77 77 78 78 // the following source->mode information pertains to modelPSF: 79 source->mode |= PM_SOURCE_ PSFMODEL;79 source->mode |= PM_SOURCE_MODE_PSFMODEL; 80 80 81 81 // if the object has fitted peak above saturation, label as SATSTAR … … 84 84 // XXX no extended object can saturate and stay extended... 85 85 if (model->params->data.F32[1] >= SATURATION) { 86 if (source->mode & PM_SOURCE_ PSFSTAR) {86 if (source->mode & PM_SOURCE_MODE_PSFSTAR) { 87 87 psLogMsg ("psphot", 5, "PSFSTAR marked SATSTAR\n"); 88 88 } 89 source->mode |= PM_SOURCE_ SATSTAR;89 source->mode |= PM_SOURCE_MODE_SATSTAR; 90 90 return true; 91 91 } … … 93 93 // if the object has a fitted peak below 0, the fit did not converge cleanly 94 94 if (model->params->data.F32[1] <= 0) { 95 source->mode |= PM_SOURCE_ FAIL;95 source->mode |= PM_SOURCE_MODE_FAIL; 96 96 return false; 97 97 } … … 99 99 // if the source was predicted to be a SATSTAR, but it fitted below saturation, 100 100 // make a note to the user 101 if (source->mode & PM_SOURCE_ SATSTAR) {101 if (source->mode & PM_SOURCE_MODE_SATSTAR) { 102 102 psLogMsg ("psphot", 5, "SATSTAR marked normal (fitted peak below saturation)\n"); 103 source->mode &= ~PM_SOURCE_ SATSTAR;103 source->mode &= ~PM_SOURCE_MODE_SATSTAR; 104 104 } 105 105 … … 115 115 nSy = dSY / hypot (MIN_DS, 1 / (SY * SN)); 116 116 117 // assign PM_SOURCE_ GOODSTAR to bright objects within PSF region of dparams[]117 // assign PM_SOURCE_MODE_GOODSTAR to bright objects within PSF region of dparams[] 118 118 keep = TRUE; 119 119 keep &= (fabs(nSx) < PSF_SHAPE_NSIGMA); … … 124 124 125 125 // this source is not a star, warn if it was a PSFSTAR 126 if (source->mode & PM_SOURCE_ PSFSTAR) {126 if (source->mode & PM_SOURCE_MODE_PSFSTAR) { 127 127 psLogMsg ("psphot", 5, "PSFSTAR demoted based on fit quality (%f, %f : %f %f %f %f)\n", 128 128 model->params->data.F32[2], model->params->data.F32[3], nSx, nSy, SN, Chi); … … 134 134 // object appears to be small, suspected defect 135 135 if ((nSx <= -PSF_SHAPE_NSIGMA) || (nSy <= -PSF_SHAPE_NSIGMA)) { 136 source->type = PM_SOURCE_ DEFECT;136 source->type = PM_SOURCE_TYPE_DEFECT; 137 137 return false; 138 138 } … … 140 140 // object appears to be large, suspected extended source 141 141 if ((nSx >= PSF_SHAPE_NSIGMA) || (nSy >= PSF_SHAPE_NSIGMA)) { 142 source->type = PM_SOURCE_ EXTENDED;142 source->type = PM_SOURCE_TYPE_EXTENDED; 143 143 return false; 144 144 } 145 145 146 146 // poor-quality fit; only keep if nothing else works... 147 source->mode |= PM_SOURCE_ POOR;147 source->mode |= PM_SOURCE_MODE_POOR; 148 148 return false; 149 149 } … … 155 155 // do we actually have a valid PSF model? 156 156 if (model == NULL) { 157 source->mode &= ~PM_SOURCE_ FITTED;157 source->mode &= ~PM_SOURCE_MODE_FITTED; 158 158 return false; 159 159 } … … 164 164 break; 165 165 case PM_MODEL_UNTRIED: 166 source->mode &= ~PM_SOURCE_ FITTED;166 source->mode &= ~PM_SOURCE_MODE_FITTED; 167 167 return false; 168 168 case PM_MODEL_BADARGS: … … 170 170 case PM_MODEL_OFFIMAGE: 171 171 default: 172 source->mode |= PM_SOURCE_ FAIL;172 source->mode |= PM_SOURCE_MODE_FAIL; 173 173 return false; 174 174 } 175 175 176 176 // unless we prove otherwise, this object is a star. 177 source->type = PM_SOURCE_ STAR;177 source->type = PM_SOURCE_TYPE_STAR; 178 178 179 179 // the following source->mode information pertains to modelPSF: 180 source->mode |= PM_SOURCE_ PSFMODEL;180 source->mode |= PM_SOURCE_MODE_PSFMODEL; 181 181 182 182 // if the object has fitted peak above saturation, label as SATSTAR … … 185 185 // XXX no extended object can saturate and stay extended... 186 186 if (model->params->data.F32[1] >= SATURATION) { 187 if (source->mode & PM_SOURCE_ PSFSTAR) {187 if (source->mode & PM_SOURCE_MODE_PSFSTAR) { 188 188 psLogMsg ("psphot", 5, "PSFSTAR marked SATSTAR\n"); 189 189 } 190 source->mode |= PM_SOURCE_ SATSTAR;190 source->mode |= PM_SOURCE_MODE_SATSTAR; 191 191 return true; 192 192 } … … 194 194 // if the object has a fitted peak below 0, the fit did not converge cleanly 195 195 if (model->params->data.F32[1] <= 0) { 196 source->mode |= PM_SOURCE_ FAIL;196 source->mode |= PM_SOURCE_MODE_FAIL; 197 197 return false; 198 198 } … … 200 200 // if the source was predicted to be a SATSTAR, but it fitted below saturation, 201 201 // make a note to the user 202 if (source->mode & PM_SOURCE_ SATSTAR) {202 if (source->mode & PM_SOURCE_MODE_SATSTAR) { 203 203 psLogMsg ("psphot", 5, "SATSTAR marked normal (fitted peak below saturation)\n"); 204 source->mode &= ~PM_SOURCE_ SATSTAR;204 source->mode &= ~PM_SOURCE_MODE_SATSTAR; 205 205 } 206 206 return true;
Note:
See TracChangeset
for help on using the changeset viewer.
