IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Jan 19, 2006, 2:55:33 PM (21 years ago)
Author:
eugene
Message:

adapted everything to work with pmReadout *readout instead of eamReadout *imdata

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/psphot/src/psphotSourceFits.c

    r5993 r6056  
    44// given a source with an existing modelPSF, attempt a full PSF fit, subtract if successful
    55
    6 bool psphotFitPSF (eamReadout *imdata, pmSource *source) {
     6bool psphotFitPSF (pmReadout *readout, pmSource *source) {
    77
    88    float x, y;
     
    1212
    1313    // extend source radius as needed
    14     psphotCheckRadiusPSF (imdata, source, PSF);
     14    psphotCheckRadiusPSF (readout, source, PSF);
    1515
    1616    x = PSF->params->data.F32[2];
     
    6161}
    6262
    63 bool psphotFitBlob (eamReadout *imdata, pmSource *source, psArray *sources) {
     63bool psphotFitBlob (pmReadout *readout, pmSource *source, psArray *sources) {
    6464
    6565    bool okEXT, okDBL;
     
    7979    pmSource *tmpSrc = pmSourceAlloc ();
    8080
    81     pmModel *EXT = psphotFitEXT (imdata, source);
     81    pmModel *EXT = psphotFitEXT (readout, source);
    8282    okEXT = psphotEvalEXT (tmpSrc, EXT);
    8383    chiEXT = EXT->chisq / EXT->nDOF;
    8484
    85     psArray *DBL = psphotFitDBL (imdata, source);
     85    psArray *DBL = psphotFitDBL (readout, source);
    8686    okDBL  = psphotEvalDBL (tmpSrc, DBL->data[0]);
    8787    okDBL &= psphotEvalDBL (tmpSrc, DBL->data[1]);
     
    129129
    130130// fit a double PSF source to an extended blob
    131 psArray *psphotFitDBL (eamReadout *imdata, pmSource *source) {
     131psArray *psphotFitDBL (pmReadout *readout, pmSource *source) {
    132132
    133133    float x, y, dx, dy;
     
    149149    // save the PSF model from the Ensemble fit
    150150    PSF = source->modelPSF;
    151     psphotCheckRadiusPSF (imdata, source, PSF);
     151    psphotCheckRadiusPSF (readout, source, PSF);
    152152
    153153    modelSet = psArrayAlloc (2);
     
    176176}
    177177
    178 pmModel *psphotFitEXT (eamReadout *imdata, pmSource *source) {
     178pmModel *psphotFitEXT (pmReadout *readout, pmSource *source) {
    179179
    180180    float x, y;
     
    183183    pmModel *EXT = pmSourceModelGuess (source, modelTypeEXT);
    184184
    185     psphotCheckRadiusEXT (imdata, source, EXT);
     185    psphotCheckRadiusEXT (readout, source, EXT);
    186186
    187187    x = EXT->params->data.F32[2];
     
    196196}
    197197
    198 bool psphotFitBlend (eamReadout *imdata, pmSource *source) {
     198bool psphotFitBlend (pmReadout *readout, pmSource *source) {
    199199
    200200    float x, y, dR;
     
    202202    // if this source is not a possible blend, just fit as PSF
    203203    if ((source->blends == NULL) || (source->mode & PM_SOURCE_SATSTAR)) {
    204         bool status = psphotFitPSF (imdata, source);
     204        bool status = psphotFitPSF (readout, source);
    205205        return status;
    206206    }
     
    212212
    213213    // extend source radius as needed
    214     psphotCheckRadiusPSF (imdata, source, PSF);
     214    psphotCheckRadiusPSF (readout, source, PSF);
    215215
    216216    x = PSF->params->data.F32[2];
Note: See TracChangeset for help on using the changeset viewer.