Changeset 36108 for trunk/psphot/src
- Timestamp:
- Sep 11, 2013, 11:00:36 AM (13 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotExtendedSourceFits.c (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotExtendedSourceFits.c
r36086 r36108 475 475 } 476 476 477 bool testObject = false; 478 // testObject |= ((fabs(source->peak->xf - 179) < 5) && (fabs(source->peak->yf - 1138) < 5)); 479 // testObject |= ((fabs(source->peak->xf - 5047) < 5) && (fabs(source->peak->yf - 151) < 5)); 480 // testObject |= ((fabs(source->peak->xf - 3929) < 5) && (fabs(source->peak->yf - 4109) < 5)); 481 // testObject |= ((fabs(source->peak->xf - 915) < 5) && (fabs(source->peak->yf - 5998) < 5)); 482 // testObject |= ((fabs(source->peak->xf - 5406) < 5) && (fabs(source->peak->yf - 326) < 5)); 483 if (testObject) { 484 fprintf (stderr, "test object @ %f, %f\n", source->peak->xf, source->peak->yf); 485 psTraceSetLevel ("psModules.objects.pmPCM_MinimizeChisq", 5); 486 psTraceSetLevel ("psphot.psphotExtendedSourceFits_Threaded", 5); 487 } 488 477 489 // loop here over the models chosen for each source (exclude by S/N) 478 490 // Reset the iterator … … 531 543 psFree (source->modelFlux); 532 544 source->modelFlux = NULL; 533 // if ((fabs(source->peak->xf - 2572) < 20) && (fabs(source->peak->yf - 5874) < 20)) {534 // fprintf (stderr, "test object\n");535 // }536 537 545 modelFit = psphotFitEXT (modelFit, readout, source, fitOptions, modelType, maskVal, markVal); 538 546 if (!modelFit) { … … 646 654 psTrace ("psphot", 4, "best ext model for %f, %f : %s chisq = %f\n", source->moments->Mx, source->moments->My, pmModelClassGetName (source->modelEXT->type), source->modelEXT->chisq); 647 655 psTrace ("psphot", 5, "extended source model for source at %7.1f, %7.1f", source->moments->Mx, source->moments->My); 656 657 if (testObject) { 658 psTraceSetLevel ("psModules.objects.pmPCM_MinimizeChisq", 0); 659 psTraceSetLevel ("psphot.psphotExtendedSourceFits_Threaded", 0); 660 } 648 661 } 649 662
Note:
See TracChangeset
for help on using the changeset viewer.
