Changeset 34382
- Timestamp:
- Sep 4, 2012, 6:31:26 AM (14 years ago)
- Location:
- branches/eam_branches/ipp-20120805/Ohana/src
- Files:
-
- 53 edited
-
addstar/src/ReadStarsFITS.c (modified) (9 diffs)
-
addstar/src/SEDfit.c (modified) (2 diffs)
-
addstar/src/load2mass_catalog.c (modified) (1 diff)
-
addstar/src/mkcmf.c (modified) (5 diffs)
-
addstar/src/update_coords.c (modified) (1 diff)
-
dvomerge/src/merge_catalogs_old.c (modified) (1 diff)
-
libautocode/Makefile.Targets (modified) (2 diffs)
-
libautocode/def/Stars.d (modified) (1 diff)
-
libautocode/def/average-ps1-v4.d (modified) (2 diffs)
-
libautocode/def/average.d (modified) (2 diffs)
-
libautocode/def/cmf-ps1-dev-0.d (modified) (1 diff)
-
libautocode/def/cmf-ps1-dev-1.d (modified) (1 diff)
-
libautocode/def/cmf-ps1-sv1.d (modified) (2 diffs)
-
libautocode/def/cmf-ps1-v1.d (modified) (1 diff)
-
libautocode/def/cmf-ps1-v2.d (modified) (1 diff)
-
libautocode/def/cmf-ps1-v3.d (modified) (2 diffs)
-
libautocode/def/measure-panstarrs-dev-0.d (modified) (1 diff)
-
libautocode/def/measure-panstarrs-dev-1.d (modified) (1 diff)
-
libautocode/def/measure-ps1-dev-1.d (modified) (1 diff)
-
libautocode/def/measure-ps1-dev-2.d (modified) (1 diff)
-
libautocode/def/measure-ps1-v1.d (modified) (1 diff)
-
libautocode/def/measure-ps1-v2.d (modified) (1 diff)
-
libautocode/def/measure-ps1-v3.d (modified) (1 diff)
-
libautocode/def/measure-ps1-v4.d (modified) (2 diffs)
-
libautocode/def/measure.d (modified) (2 diffs)
-
libdvo/src/dbExtractAverages.c (modified) (1 diff)
-
libdvo/src/dbExtractMeasures.c (modified) (1 diff)
-
libdvo/src/dbFields.c (modified) (1 diff)
-
libdvo/src/dvo_catalog.c (modified) (2 diffs)
-
libdvo/src/dvo_convert_PS1_DEV_1.c (modified) (6 diffs)
-
libdvo/src/dvo_convert_PS1_DEV_2.c (modified) (6 diffs)
-
libdvo/src/dvo_convert_PS1_REF.c (modified) (4 diffs)
-
libdvo/src/dvo_convert_PS1_V1.c (modified) (6 diffs)
-
libdvo/src/dvo_convert_PS1_V2.c (modified) (6 diffs)
-
libdvo/src/dvo_convert_PS1_V3.c (modified) (6 diffs)
-
libdvo/src/dvo_convert_PS1_V4.c (modified) (4 diffs)
-
libdvo/src/dvo_convert_elixir.c (modified) (5 diffs)
-
libdvo/src/dvo_convert_loneos.c (modified) (5 diffs)
-
libdvo/src/dvo_convert_panstarrs_DEV_0.c (modified) (6 diffs)
-
libdvo/src/dvo_convert_panstarrs_DEV_1.c (modified) (6 diffs)
-
opihi/dvo/avextract.c (modified) (1 diff)
-
opihi/dvo/avmatch.c (modified) (1 diff)
-
opihi/dvo/gstar.c (modified) (3 diffs)
-
opihi/dvo/photometry.c (modified) (3 diffs)
-
relastro/src/BrightCatalog.c (modified) (7 diffs)
-
relastro/src/UpdateObjectOffsets.c (modified) (1 diff)
-
relastro/src/UpdateObjects.c (modified) (1 diff)
-
relastro/src/load_catalogs.c (modified) (1 diff)
-
relastro/src/relastro_objects.c (modified) (1 diff)
-
relphot/src/bcatalog.c (modified) (1 diff)
-
relphot/src/load_catalogs.c (modified) (1 diff)
-
relphot/src/relphot_objects.c (modified) (1 diff)
-
relphot/src/setMrelFinal.c (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/ipp-20120805/Ohana/src/addstar/src/ReadStarsFITS.c
r34260 r34382 120 120 121 121 stars[i].measure.psfChisq = NAN; // not provided by SMPDATA: 122 stars[i].measure.psfQ ual= NAN; // not provided by SMPDATA:122 stars[i].measure.psfQF = NAN; // not provided by SMPDATA: 123 123 stars[i].measure.psfNdof = NAN_S_INT; // not provided by SMPDATA: 124 124 stars[i].measure.psfNpix = NAN_S_INT; // not provided by SMPDATA: … … 189 189 190 190 stars[i].measure.psfChisq = ps1data[i].psfChisq; 191 stars[i].measure.psfQ ual = ps1data[i].psfQual;191 stars[i].measure.psfQF = ps1data[i].psfQF; 192 192 stars[i].measure.psfNdof = NAN_S_INT; // not provided by PS1_DEV_0: 193 193 stars[i].measure.psfNpix = NAN_S_INT; // not provided by PS1_DEV_0: … … 261 261 262 262 stars[i].measure.psfChisq = ps1data[i].psfChisq; 263 stars[i].measure.psfQ ual = ps1data[i].psfQual;263 stars[i].measure.psfQF = ps1data[i].psfQF; 264 264 stars[i].measure.psfNdof = NAN_S_INT; // not provided by PS1_DEV_1: 265 265 stars[i].measure.psfNpix = NAN_S_INT; // not provided by PS1_DEV_1: … … 345 345 346 346 stars[i].measure.psfChisq = ps1data[i].psfChisq; 347 stars[i].measure.psfQ ual = ps1data[i].psfQual;347 stars[i].measure.psfQF = ps1data[i].psfQF; 348 348 stars[i].measure.psfNdof = ps1data[i].psfNdof; 349 349 stars[i].measure.psfNpix = ps1data[i].psfNpix; … … 430 430 431 431 stars[i].measure.psfChisq = ps1data[i].psfChisq; 432 stars[i].measure.psfQ ual = ps1data[i].psfQual;432 stars[i].measure.psfQF = ps1data[i].psfQF; 433 433 stars[i].measure.psfNdof = ps1data[i].psfNdof; 434 434 stars[i].measure.psfNpix = ps1data[i].psfNpix; … … 507 507 508 508 stars[i].measure.psfChisq = ps1data[i].psfChisq; 509 stars[i].measure.psfQ ual = ps1data[i].psfQual;509 stars[i].measure.psfQF = ps1data[i].psfQF; 510 510 stars[i].measure.psfNdof = ps1data[i].psfNdof; 511 511 stars[i].measure.psfNpix = ps1data[i].psfNpix; … … 584 584 585 585 stars[i].measure.psfChisq = ps1data[i].psfChisq; 586 stars[i].measure.psfQ ual = ps1data[i].psfQual;586 stars[i].measure.psfQF = ps1data[i].psfQF; 587 587 stars[i].measure.psfNdof = ps1data[i].psfNdof; 588 588 stars[i].measure.psfNpix = ps1data[i].psfNpix; … … 669 669 670 670 stars[i].measure.psfChisq = ps1data[i].psfChisq; 671 stars[i].measure.psfQ ual = ps1data[i].psfQual;671 stars[i].measure.psfQF = ps1data[i].psfQF; 672 672 stars[i].measure.psfNdof = ps1data[i].psfNdof; 673 673 stars[i].measure.psfNpix = ps1data[i].psfNpix; … … 748 748 749 749 stars[i].measure.psfChisq = ps1data[i].psfChisq; 750 stars[i].measure.psfQ ual = ps1data[i].psfQual;750 stars[i].measure.psfQF = ps1data[i].psfQF; 751 751 stars[i].measure.psfNdof = ps1data[i].psfNdof; 752 752 stars[i].measure.psfNpix = ps1data[i].psfNpix; -
branches/eam_branches/ipp-20120805/Ohana/src/addstar/src/SEDfit.c
r21508 r34382 147 147 // construct an average object for this object 148 148 // XXX for now, the output objects will have limited astrometric interpretation... 149 dvo_average_init (&outcat[0].average[Nave]); 149 150 outcat[0].average[Nave].R = incat[0].average[i].R; 150 151 outcat[0].average[Nave].D = incat[0].average[i].D; 151 outcat[0].average[Nave].dR = 0;152 outcat[0].average[Nave].dD = 0;153 outcat[0].average[Nave].uR = 0;154 outcat[0].average[Nave].uD = 0;155 outcat[0].average[Nave].duR = 0;156 outcat[0].average[Nave].duD = 0;157 outcat[0].average[Nave].P = 0;158 outcat[0].average[Nave].dP = 0;159 160 // XXX for now, set the average mag data to NULL161 outcat[0].average[Nave].Nmeasure = 0;162 outcat[0].average[Nave].Nmissing = 0;163 outcat[0].average[Nave].Xp = NAN_S_SHORT;164 outcat[0].average[Nave].measureOffset = Nmeas;165 outcat[0].average[Nave].missingOffset = -1;166 outcat[0].average[Nave].flags = 0;167 152 168 153 for (j = 0; j < Nsec; j++) { 169 outcat[0].secfilt[Nave*Nsec+j].M = NAN; 170 outcat[0].secfilt[Nave*Nsec+j].dM = NAN; 171 outcat[0].secfilt[Nave*Nsec+j].Xm = NAN_S_SHORT; 154 dvo_secfilt_init (&outcat[0].secfilt[Nave*Nsec+j]); 172 155 } 173 156 … … 181 164 for (j = 0; valid && (j < Nmodel); j++) { 182 165 n = modelRow[j]; 166 dvo_measure_init (&outcat[0].measure[Nmeas]); 183 167 outcat[0].measure[Nmeas].dR = 0.0; 184 168 outcat[0].measure[Nmeas].dD = 0.0; -
branches/eam_branches/ipp-20120805/Ohana/src/addstar/src/load2mass_catalog.c
r21508 r34382 34 34 catalog[0].average[Nave].Nmeasure = 0; 35 35 catalog[0].average[Nave].Nmissing = 0; 36 catalog[0].average[Nave]. Xp= NAN_S_SHORT;36 catalog[0].average[Nave].ChiSqAve = NAN_S_SHORT; 37 37 catalog[0].average[Nave].measureOffset = Nmeas; 38 38 catalog[0].average[Nave].missingOffset = -1; -
branches/eam_branches/ipp-20120805/Ohana/src/addstar/src/mkcmf.c
r34260 r34382 419 419 stars[i].fy = FY; 420 420 stars[i].df = DF; 421 stars[i].psfQ ual= PSFQUAL;421 stars[i].psfQF = PSFQUAL; 422 422 stars[i].nFrames = 1; 423 423 } … … 466 466 stars[i].fy = FY; 467 467 stars[i].df = DF; 468 stars[i].psfQ ual= PSFQUAL;468 stars[i].psfQF = PSFQUAL; 469 469 stars[i].nFrames = 1; 470 470 stars[i].flags = FLAGS; … … 514 514 stars[i].fy = FY; 515 515 stars[i].df = DF; 516 stars[i].psfQ ual= PSFQUAL;516 stars[i].psfQF = PSFQUAL; 517 517 stars[i].nFrames = 1; 518 518 stars[i].flags = FLAGS; … … 551 551 // randomly give poor PSFQF values 552 552 if ((BAD_PSFQF_FRAC > 0.0) && (drand48() < BAD_PSFQF_FRAC)) { 553 stars[i].psfQ ual= 0.25;553 stars[i].psfQF = 0.25; 554 554 } else { 555 stars[i].psfQ ual= PSFQUAL;555 stars[i].psfQF = PSFQUAL; 556 556 } 557 557 … … 606 606 // randomly give poor PSFQF values 607 607 if ((BAD_PSFQF_FRAC > 0.0) && (drand48() < BAD_PSFQF_FRAC)) { 608 stars[i].psfQ ual= 0.25;608 stars[i].psfQF = 0.25; 609 609 } else { 610 stars[i].psfQ ual= PSFQUAL;610 stars[i].psfQF = PSFQUAL; 611 611 } 612 612 -
branches/eam_branches/ipp-20120805/Ohana/src/addstar/src/update_coords.c
r27435 r34382 47 47 dR2 = r2 / Npt - r*r; 48 48 dD2 = d2 / Npt - d*d; 49 average[0]. Xp= sqrt (dD2 + dR2 / SQ(cos(d*RAD_DEG)));50 /* Xp is scatter in position in hundredths of arcsec*/49 average[0].ChiSqAve = sqrt (dD2 + dR2 / SQ(cos(d*RAD_DEG))); 50 /* ChiSqAve is supposed to be a chisq */ 51 51 52 52 return; -
branches/eam_branches/ipp-20120805/Ohana/src/dvomerge/src/merge_catalogs_old.c
r33963 r34382 320 320 output[0].average[Nave].dP = 0; 321 321 322 output[0].average[Nave]. Xp= 0;322 output[0].average[Nave].stargal = 0; 323 323 output[0].average[Nave].ChiSqAve = 0.0; 324 324 output[0].average[Nave].ChiSqPM = 0.0; -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/Makefile.Targets
r34260 r34382 70 70 $(ASRC)/cmf-ps1-v2.$(ARCH).o \ 71 71 $(ASRC)/cmf-ps1-v3.$(ARCH).o \ 72 $(ASRC)/cmf-ps1-v4.$(ARCH).o \ 72 73 $(ASRC)/cmf-ps1-sv1.$(ARCH).o \ 74 $(ASRC)/cmf-ps1-sv2.$(ARCH).o \ 73 75 $(ASRC)/cmf-smpdata.$(ARCH).o \ 74 76 $(ASRC)/getstar-ps1-dev-0.$(ARCH).o \ … … 154 156 $(AINC)/cmf-ps1-v2.h \ 155 157 $(AINC)/cmf-ps1-v3.h \ 158 $(AINC)/cmf-ps1-v4.h \ 156 159 $(AINC)/cmf-ps1-sv1.h \ 160 $(AINC)/cmf-ps1-sv2.h \ 157 161 $(AINC)/cmf-smpdata.h \ 158 162 $(AINC)/getstar-ps1-dev-0.h \ -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/Stars.d
r21508 r34382 1 *** note that this file / structure has been deprecated 2 3 1 4 STRUCT Stars 2 5 EXTNAME STARS -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/average-ps1-v4.d
r34260 r34382 2 2 EXTNAME DVO_AVERAGE_PS1_V4 3 3 TYPE BINTABLE 4 SIZE 12 04 SIZE 128 5 5 DESCRIPTION DVO Average Object Table 6 6 … … 25 25 FIELD Trange, TIME_RANGE, int, mean epoch (PM,PAR ref), unix time seconds 26 26 27 FIELD Xp, SIGMA_POS, float, position scatter, 1/100 arcsec 27 FIELD psfQF, PSF_QF, float, psf coverage (bad masks) 28 FIELD psfQFperf, PSF_QF_PERF, float, psf coverage (all masks) 29 30 FIELD stargal, STARGAL_SEP, float, star / galaxy separator, 1/100 arcsec 28 31 FIELD Npos, NUMBER_POS, unsigned short, number of detections used for astrometry 29 32 -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/average.d
r33647 r34382 2 2 EXTNAME DVO_AVERAGE 3 3 TYPE BINTABLE 4 SIZE 12 04 SIZE 128 5 5 DESCRIPTION DVO Average Object Table 6 6 … … 25 25 FIELD Trange, TIME_RANGE, int, mean epoch (PM,PAR ref), unix time seconds 26 26 27 FIELD Xp, SIGMA_POS, float, position scatter, 1/100 arcsec 27 FIELD psfQF, PSF_QF, float, psf coverage (bad masks) 28 FIELD psfQFperf, PSF_QF_PERF, float, psf coverage (all masks) 29 30 FIELD stargal, STARGAL_SEP, float, star / galaxy separator, 1/100 arcsec 28 31 FIELD Npos, NUMBER_POS, unsigned short, number of detections used for astrometry 29 32 -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/cmf-ps1-dev-0.d
r21508 r34382 20 20 FIELD fy, PSF_WIDTH_Y, float, semi-minor, pixels 21 21 FIELD df, PSF_THETA, float, ellipse angle, degrees 22 FIELD psfQ ual, PSF_QF, float, quality factor22 FIELD psfQF, PSF_QF, float, quality factor 23 23 FIELD nFrames, N_FRAMES, short, images overlapping peak 24 24 FIELD dummy, DUMMY, short, padding -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/cmf-ps1-dev-1.d
r21508 r34382 22 22 FIELD fy, PSF_WIDTH_Y, float, semi-minor, pixels 23 23 FIELD df, PSF_THETA, float, ellipse angle, degrees 24 FIELD psfQ ual, PSF_QF, float, quality factor24 FIELD psfQF, PSF_QF, float, quality factor 25 25 FIELD nFrames, N_FRAMES, short, images overlapping peak 26 26 FIELD flags, FLAGS, short, padding -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/cmf-ps1-sv1.d
r34260 r34382 51 51 52 52 FIELD df, PSF_THETA, float, ellipse angle, degrees 53 FIELD psfQ ual,PSF_QF, float, quality factor54 FIELD psfQ ualPerfect,PSF_QF_PERFECT, float, quality factor perfect53 FIELD psfQF, PSF_QF, float, quality factor 54 FIELD psfQFperf, PSF_QF_PERFECT, float, quality factor perfect 55 55 FIELD psfNdof, PSF_NDOF, int, psf degrees of freedom 56 56 … … 112 112 # FIELD fy, PSF_MINOR, PSF_MINOR 1E float, psf fit minor axis, pixels 113 113 # FIELD df, PSF_THETA, PSF_THETA 1E float, ellipse angle, degrees 114 # FIELD psfQ ual,PSF_QF, PSF_QF 1E float, quality factor115 # FIELD psfQ ualPerfect,PSF_QF_PERFECT, PSF_QF_PERFECT 1E float, quality factor perfect114 # FIELD psfQF, PSF_QF, PSF_QF 1E float, quality factor 115 # FIELD psfQFperf, PSF_QF_PERFECT, PSF_QF_PERFECT 1E float, quality factor perfect 116 116 # FIELD psfNdof, PSF_NDOF, PSF_NDOF 1J int, psf degrees of freedom 117 117 # FIELD psfNpix, PSF_NPIX, PSF_NPIX 1J int, psf number of pixels -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/cmf-ps1-v1.d
r34260 r34382 30 30 FIELD fy, PSF_MINOR, float, psf fit minor axis, pixels 31 31 FIELD df, PSF_THETA, float, ellipse angle, degrees 32 FIELD psfQ ual, PSF_QF, float, quality factor32 FIELD psfQF, PSF_QF, float, quality factor 33 33 FIELD psfNdof, PSF_NDOF, int, psf degrees of freedom 34 34 FIELD psfNpix, PSF_NPIX, int, psf number of pixels -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/cmf-ps1-v2.d
r34260 r34382 30 30 FIELD fy, PSF_MINOR, float, psf fit minor axis, pixels 31 31 FIELD df, PSF_THETA, float, ellipse angle, degrees 32 FIELD psfQ ual, PSF_QF, float, quality factor32 FIELD psfQF, PSF_QF, float, quality factor 33 33 FIELD psfNdof, PSF_NDOF, int, psf degrees of freedom 34 34 FIELD psfNpix, PSF_NPIX, int, psf number of pixels -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/cmf-ps1-v3.d
r34260 r34382 36 36 FIELD fy, PSF_MINOR, float, psf fit minor axis, pixels 37 37 FIELD df, PSF_THETA, float, ellipse angle, degrees 38 FIELD psfQ ual,PSF_QF, float, quality factor39 FIELD psfQ ualPerfect,PSF_QF_PERFECT, float, quality factor perfect38 FIELD psfQF, PSF_QF, float, quality factor 39 FIELD psfQFperf, PSF_QF_PERFECT, float, quality factor perfect 40 40 FIELD psfNdof, PSF_NDOF, int, psf degrees of freedom 41 41 FIELD psfNpix, PSF_NPIX, int, psf number of pixels … … 68 68 # FIELD dFlux, PSF_INST_FLUX_SIG, float, psf flux error, counts 69 69 # FIELD MapRaw, AP_MAG_RAW, float, raw aperture mag, mags 70 # FIELD psfQ ualPerfect,PSF_QF_PERFECT, float, quality factor perfect70 # FIELD psfQFperf, PSF_QF_PERFECT, float, quality factor perfect 71 71 # FIELD M3c, MOMENTS_M3C, float, third moment cos(t), pixels^3 72 72 # FIELD M3s, MOMENTS_M3S, float, third moment sin(t), pixels^3 -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/measure-panstarrs-dev-0.d
r21508 r34382 46 46 47 47 # do we need more resolution than a short? should this be a log? 48 FIELD psfQ ual,PSF_QF, short, psf coverage/quality factor48 FIELD psfQF, PSF_QF, short, psf coverage/quality factor 49 49 50 50 FIELD dophot, DOPHOT, char, dophot type -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/measure-panstarrs-dev-1.d
r21508 r34382 46 46 47 47 # do we need more resolution than a short? should this be a log? 48 FIELD psfQ ual,PSF_QF, short, psf coverage/quality factor48 FIELD psfQF, PSF_QF, short, psf coverage/quality factor 49 49 50 50 FIELD dophot, DOPHOT, char, dophot type -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/measure-ps1-dev-1.d
r21508 r34382 33 33 34 34 # do we need more resolution than a short? should this be a log? 35 FIELD psfQ ual,PSF_QF, float, psf coverage/quality factor35 FIELD psfQF, PSF_QF, float, psf coverage/quality factor 36 36 FIELD psfChisq, PSF_CHISQ, float, psf coverage/quality factor 37 37 FIELD crNsigma, CR_NSIGMA, float, psf coverage/quality factor -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/measure-ps1-dev-2.d
r21508 r34382 34 34 35 35 # do we need more resolution than a short? should this be a log? 36 FIELD psfQ ual,PSF_QF, float, psf coverage/quality factor36 FIELD psfQF, PSF_QF, float, psf coverage/quality factor 37 37 FIELD psfChisq, PSF_CHISQ, float, psf fit chisq 38 38 FIELD crNsigma, CR_NSIGMA, float, Nsigma deviation towards CR -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/measure-ps1-v1.d
r21508 r34382 39 39 40 40 # do we need more resolution than a short? should this be a log? 41 FIELD psfQ ual,PSF_QF, float, psf coverage/quality factor41 FIELD psfQF, PSF_QF, float, psf coverage/quality factor 42 42 FIELD psfChisq, PSF_CHISQ, float, psf fit chisq 43 43 FIELD psfNdof, PSF_NDOF, int, psf degrees of freedom -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/measure-ps1-v2.d
r30604 r34382 39 39 40 40 # do we need more resolution than a short? should this be a log? 41 FIELD psfQ ual,PSF_QF, float, psf coverage/quality factor41 FIELD psfQF, PSF_QF, float, psf coverage/quality factor 42 42 FIELD psfChisq, PSF_CHISQ, float, psf fit chisq 43 43 FIELD psfNdof, PSF_NDOF, int, psf degrees of freedom -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/measure-ps1-v3.d
r34260 r34382 39 39 40 40 # do we need more resolution than a short? should this be a log? 41 FIELD psfQ ual,PSF_QF, float, psf coverage/quality factor41 FIELD psfQF, PSF_QF, float, psf coverage/quality factor 42 42 FIELD psfChisq, PSF_CHISQ, float, psf fit chisq 43 43 FIELD psfNdof, PSF_NDOF, int, psf degrees of freedom -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/measure-ps1-v4.d
r34260 r34382 47 47 48 48 # do we need more resolution than a short? should this be a log? 49 FIELD psfQual, PSF_QF, float, psf coverage/quality factor 50 FIELD psfChisq, PSF_CHISQ, float, psf fit chisq 51 FIELD psfNdof, PSF_NDOF, int, psf degrees of freedom 52 FIELD psfNpix, PSF_NPIX, int, psf number of pixels 53 FIELD crNsigma, CR_NSIGMA, float, Nsigma deviation towards CR 54 FIELD extNsigma, EXT_NSIGMA, float, Nsigma deviation towards EXT 49 FIELD psfQF, PSF_QF, float, psf coverage/quality factor 50 FIELD psfQFPerf, PSF_QF_PEFECT, float, psf coverage / quality factor (all mask bits) 51 FIELD psfChisq, PSF_CHISQ, float, psf fit chisq 52 53 FIELD psfNdof, PSF_NDOF, int, psf degrees of freedom 54 FIELD psfNpix, PSF_NPIX, int, psf number of pixels 55 FIELD crNsigma, CR_NSIGMA, float, Nsigma deviation towards CR 56 FIELD extNsigma, EXT_NSIGMA, float, Nsigma deviation towards EXT 55 57 56 58 # model shape parameters … … 73 75 FIELD dRsys, POS_SYS_ERR, short, systematic error from astrom, 1/100 of pixels 74 76 75 FIELD pad, PAD, char[4], padding76 77 77 # local astrometry scales 78 78 FIELD posangle, POSANGLE, short, position angle sky to chip, (0xffff/360) deg -
branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/measure.d
r34260 r34382 2 2 EXTNAME DVO_MEASURE 3 3 TYPE BINTABLE 4 SIZE 17 24 SIZE 176 5 5 DESCRIPTION DVO Detection Measurement Table 6 6 … … 47 47 48 48 # do we need more resolution than a short? should this be a log? 49 FIELD psfQual, PSF_QF, float, psf coverage/quality factor 50 FIELD psfChisq, PSF_CHISQ, float, psf fit chisq 51 FIELD psfNdof, PSF_NDOF, int, psf degrees of freedom 52 FIELD psfNpix, PSF_NPIX, int, psf number of pixels 53 FIELD crNsigma, CR_NSIGMA, float, Nsigma deviation towards CR 54 FIELD extNsigma, EXT_NSIGMA, float, Nsigma deviation towards EXT 49 FIELD psfQF, PSF_QF, float, psf coverage/quality factor 50 FIELD psfQFPerf, PSF_QF_PEFECT, float, psf coverage / quality factor (all mask bits) 51 FIELD psfChisq, PSF_CHISQ, float, psf fit chisq 52 53 FIELD psfNdof, PSF_NDOF, int, psf degrees of freedom 54 FIELD psfNpix, PSF_NPIX, int, psf number of pixels 55 FIELD crNsigma, CR_NSIGMA, float, Nsigma deviation towards CR 56 FIELD extNsigma, EXT_NSIGMA, float, Nsigma deviation towards EXT 55 57 56 58 # model shape parameters -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dbExtractAverages.c
r34260 r34382 147 147 case AVE_TRANGE: 148 148 value.Flt = TimeValue (average[0].Trange, 0, TimeFormat); 149 break;150 151 case AVE_Xp:152 value.Flt = 0.01*average[0].Xp;153 149 break; 154 150 -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dbExtractMeasures.c
r34346 r34382 480 480 break; 481 481 case MEAS_PSF_QF: /* OK */ 482 value.Flt = measure[0].psfQ ual;482 value.Flt = measure[0].psfQF; 483 483 break; 484 484 case MEAS_PSF_QF_PERFECT: /* OK */ -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dbFields.c
r34260 r34382 335 335 if (!strcasecmp (fieldName, "EXTID_LO")) ESCAPE (AVE_EXTID_LO, MAG_NONE, OPIHI_INT); 336 336 337 // if (!strcasecmp (fieldName, "Xp")) ESCAPE (AVE_Xp, MAG_NONE, OPIHI_FLT);338 // if (!strcasecmp (fieldName, "TYPE")) ESCAPE (AVE_TYPE, MAG_NONE, OPIHI_INT);339 // for words that don't parse, try a photcode340 341 337 // check for code:mode in photcode name 342 338 code = ParsePhotcodeField (fieldName, &mode, MAG_AVE); -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dvo_catalog.c
r34260 r34382 89 89 average->Trange = 0; 90 90 91 average-> Xp =0;91 average->stargal = 0.0; 92 92 average->Npos = 0; 93 93 … … 174 174 measure->extID = 0; 175 175 176 measure->psfQ ual= NAN;176 measure->psfQF = NAN; 177 177 measure->psfChisq = NAN; 178 178 measure->psfNdof = 0; -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dvo_convert_PS1_DEV_1.c
r34260 r34382 27 27 out[i].detID = in[i].detID; 28 28 out[i].imageID = in[i].imageID; 29 out[i].psfQ ual = in[i].psfQual;29 out[i].psfQF = in[i].psfQF; 30 30 out[i].psfChisq = in[i].psfChisq; 31 31 out[i].crNsigma = in[i].crNsigma; … … 65 65 out[i].FluxKron = NAN; 66 66 out[i].dFluxKron = NAN; 67 out[i].psfQFperf = NAN; 67 68 } 68 69 return (out); … … 93 94 out[i].detID = in[i].detID; 94 95 out[i].imageID = in[i].imageID; 95 out[i].psfQ ual = in[i].psfQual;96 out[i].psfQF = in[i].psfQF; 96 97 out[i].psfChisq = in[i].psfChisq; 97 98 out[i].crNsigma = in[i].crNsigma; … … 126 127 out[i].R = in[i].R; 127 128 out[i].D = in[i].D; 128 out[i].Xp = in[i].Xp;129 129 out[i].dR = in[i].dR; 130 130 out[i].dD = in[i].dD; … … 159 159 out[i].photFlagsUpper = 0; 160 160 out[i].photFlagsLower = 0; 161 162 // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway) 163 out[i].psfQF = NAN; 164 out[i].psfQFperf = NAN; 165 out[i].stargal = NAN; 161 166 } 162 167 return (out); … … 174 179 out[i].R = in[i].R; 175 180 out[i].D = in[i].D; 176 out[i].Xp = in[i].Xp;177 181 out[i].dR = in[i].dR; 178 182 out[i].dD = in[i].dD; -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dvo_convert_PS1_DEV_2.c
r34260 r34382 29 29 out[i].detID = in[i].detID; 30 30 out[i].imageID = in[i].imageID; 31 out[i].psfQ ual = in[i].psfQual;31 out[i].psfQF = in[i].psfQF; 32 32 out[i].psfChisq = in[i].psfChisq; 33 33 out[i].crNsigma = in[i].crNsigma; … … 62 62 out[i].FluxKron = NAN; 63 63 out[i].dFluxKron = NAN; 64 } 64 out[i].psfQFperf = NAN; 65 } 65 66 return (out); 66 67 } … … 92 93 out[i].detID = in[i].detID; 93 94 out[i].imageID = in[i].imageID; 94 out[i].psfQ ual = in[i].psfQual;95 out[i].psfQF = in[i].psfQF; 95 96 out[i].psfChisq = in[i].psfChisq; 96 97 out[i].crNsigma = in[i].crNsigma; … … 130 131 out[i].P = in[i].P; 131 132 out[i].dP = in[i].dP; 132 out[i].Xp = in[i].Xp;133 133 out[i].Nmeasure = in[i].Nmeasure; 134 134 out[i].Nmissing = in[i].Nmissing; … … 153 153 out[i].photFlagsUpper = 0; 154 154 out[i].photFlagsLower = 0; 155 156 // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway) 157 out[i].psfQF = NAN; 158 out[i].psfQFperf = NAN; 159 out[i].stargal = NAN; 155 160 } 156 161 return (out); … … 176 181 out[i].P = in[i].P; 177 182 out[i].dP = in[i].dP; 178 out[i].Xp = in[i].Xp;179 183 out[i].Nmeasure = in[i].Nmeasure; 180 184 out[i].Nmissing = in[i].Nmissing; -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dvo_convert_PS1_REF.c
r34260 r34382 33 33 out[i].catID = in[i].catID; 34 34 out[i].extID = 0; 35 out[i].psfQ ual= 1.0; // we do not save a value, but some operations (relphot/average) filter on this value..35 out[i].psfQF = 1.0; // we do not save a value, but some operations (relphot/average) filter on this value.. 36 36 out[i].psfChisq = NAN; 37 37 out[i].psfNdof = 0; … … 60 60 out[i].FluxKron = NAN; 61 61 out[i].dFluxKron = NAN; 62 out[i].psfQFperf = NAN; 62 63 } 63 64 return (out); … … 105 106 out[i].P = NAN; 106 107 out[i].dP = NAN; 107 out[i].Xp = NAN_S_SHORT;108 108 out[i].ChiSqAve = NAN; 109 109 out[i].ChiSqPM = NAN; … … 126 126 out[i].photFlagsUpper = 0; 127 127 out[i].photFlagsLower = 0; 128 129 // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway) 130 out[i].psfQF = NAN; 131 out[i].psfQFperf = NAN; 132 out[i].stargal = NAN; 128 133 } 129 134 return (out); -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dvo_convert_PS1_V1.c
r34260 r34382 33 33 out[i].catID = in[i].catID; 34 34 out[i].extID = in[i].extID; 35 out[i].psfQ ual = in[i].psfQual;35 out[i].psfQF = in[i].psfQF; 36 36 out[i].psfChisq = in[i].psfChisq; 37 37 out[i].psfNdof = in[i].psfNdof; … … 60 60 out[i].FluxKron = NAN; 61 61 out[i].dFluxKron = NAN; 62 out[i].psfQFperf = NAN; 62 63 } 63 64 return (out); … … 94 95 out[i].catID = in[i].catID; 95 96 out[i].extID = in[i].extID; 96 out[i].psfQ ual = in[i].psfQual;97 out[i].psfQF = in[i].psfQF; 97 98 out[i].psfChisq = in[i].psfChisq; 98 99 out[i].psfNdof = in[i].psfNdof; … … 137 138 out[i].P = in[i].P; 138 139 out[i].dP = in[i].dP; 139 out[i].Xp = in[i].Xp;140 140 out[i].ChiSqAve = in[i].ChiSq; 141 141 out[i].ChiSqPM = NAN; … … 158 158 out[i].photFlagsUpper = 0; 159 159 out[i].photFlagsLower = 0; 160 161 // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway) 162 out[i].psfQF = NAN; 163 out[i].psfQFperf = NAN; 164 out[i].stargal = NAN; 160 165 } 161 166 return (out); … … 181 186 out[i].P = in[i].P; 182 187 out[i].dP = in[i].dP; 183 out[i].Xp = in[i].Xp;184 188 out[i].ChiSq = in[i].ChiSqAve; 185 189 out[i].Npos = in[i].Npos; -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dvo_convert_PS1_V2.c
r34260 r34382 33 33 out[i].catID = in[i].catID; 34 34 out[i].extID = in[i].extID; 35 out[i].psfQ ual = in[i].psfQual;35 out[i].psfQF = in[i].psfQF; 36 36 out[i].psfChisq = in[i].psfChisq; 37 37 out[i].psfNdof = in[i].psfNdof; … … 61 61 out[i].FluxKron = NAN; 62 62 out[i].dFluxKron = NAN; 63 out[i].psfQFperf = NAN; 63 64 } 64 65 return (out); … … 95 96 out[i].catID = in[i].catID; 96 97 out[i].extID = in[i].extID; 97 out[i].psfQ ual = in[i].psfQual;98 out[i].psfQF = in[i].psfQF; 98 99 out[i].psfChisq = in[i].psfChisq; 99 100 out[i].psfNdof = in[i].psfNdof; … … 138 139 out[i].P = in[i].P; 139 140 out[i].dP = in[i].dP; 140 out[i].Xp = in[i].Xp;141 141 out[i].ChiSqAve = in[i].ChiSqAve; 142 142 out[i].ChiSqPM = in[i].ChiSqPM; … … 159 159 out[i].photFlagsUpper = 0; 160 160 out[i].photFlagsLower = 0; 161 162 // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway) 163 out[i].psfQF = NAN; 164 out[i].psfQFperf = NAN; 165 out[i].stargal = NAN; 161 166 } 162 167 return (out); … … 182 187 out[i].P = in[i].P; 183 188 out[i].dP = in[i].dP; 184 out[i].Xp = in[i].Xp;185 189 out[i].ChiSqAve = in[i].ChiSqAve; 186 190 out[i].ChiSqPM = in[i].ChiSqPM; -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dvo_convert_PS1_V3.c
r34260 r34382 33 33 out[i].catID = in[i].catID; 34 34 out[i].extID = in[i].extID; 35 out[i].psfQ ual = in[i].psfQual;35 out[i].psfQF = in[i].psfQF; 36 36 out[i].psfChisq = in[i].psfChisq; 37 37 out[i].psfNdof = in[i].psfNdof; … … 61 61 out[i].FluxKron = NAN; 62 62 out[i].dFluxKron = NAN; 63 out[i].psfQFperf = NAN; 63 64 } 64 65 return (out); … … 95 96 out[i].catID = in[i].catID; 96 97 out[i].extID = in[i].extID; 97 out[i].psfQ ual = in[i].psfQual;98 out[i].psfQF = in[i].psfQF; 98 99 out[i].psfChisq = in[i].psfChisq; 99 100 out[i].psfNdof = in[i].psfNdof; … … 138 139 out[i].P = in[i].P; 139 140 out[i].dP = in[i].dP; 140 out[i].Xp = in[i].Xp;141 141 out[i].ChiSqAve = in[i].ChiSqAve; 142 142 out[i].ChiSqPM = in[i].ChiSqPM; … … 157 157 out[i].catID = in[i].catID; 158 158 out[i].extID = in[i].extID; 159 160 // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway) 161 out[i].psfQF = NAN; 162 out[i].psfQFperf = NAN; 163 out[i].stargal = NAN; 159 164 } 160 165 return (out); … … 180 185 out[i].P = in[i].P; 181 186 out[i].dP = in[i].dP; 182 out[i].Xp = in[i].Xp;183 187 out[i].ChiSqAve = in[i].ChiSqAve; 184 188 out[i].ChiSqPM = in[i].ChiSqPM; -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dvo_convert_PS1_V4.c
r34260 r34382 39 39 out[i].catID = in[i].catID; 40 40 out[i].extID = in[i].extID; 41 out[i].psfQual = in[i].psfQual; 41 out[i].psfQF = in[i].psfQF; 42 out[i].psfQFperf = in[i].psfQFperf; 42 43 out[i].psfChisq = in[i].psfChisq; 43 44 out[i].psfNdof = in[i].psfNdof; … … 99 100 out[i].catID = in[i].catID; 100 101 out[i].extID = in[i].extID; 101 out[i].psfQual = in[i].psfQual; 102 out[i].psfQF = in[i].psfQF; 103 out[i].psfQFperf = in[i].psfQFperf; 102 104 out[i].psfChisq = in[i].psfChisq; 103 105 out[i].psfNdof = in[i].psfNdof; … … 142 144 out[i].P = in[i].P; 143 145 out[i].dP = in[i].dP; 144 out[i].Xp = in[i].Xp; 146 out[i].psfQF = in[i].psfQF; 147 out[i].psfQFperf = in[i].psfQFperf; 148 out[i].stargal = in[i].stargal; 145 149 out[i].ChiSqAve = in[i].ChiSqAve; 146 150 out[i].ChiSqPM = in[i].ChiSqPM; … … 184 188 out[i].P = in[i].P; 185 189 out[i].dP = in[i].dP; 186 out[i].Xp = in[i].Xp; 190 out[i].psfQF = in[i].psfQF; 191 out[i].psfQFperf = in[i].psfQFperf; 192 out[i].stargal = in[i].stargal; 187 193 out[i].ChiSqAve = in[i].ChiSqAve; 188 194 out[i].ChiSqPM = in[i].ChiSqPM; -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dvo_convert_elixir.c
r34260 r34382 35 35 out[i].Sky = 0; 36 36 out[i].dSky = 0; 37 out[i].psfQ ual= 0;37 out[i].psfQF = 0; 38 38 out[i].psfChisq = 0; 39 39 out[i].crNsigma = 0; … … 72 72 out[i].FluxKron = NAN; 73 73 out[i].dFluxKron = NAN; 74 out[i].psfQFperf = NAN; 74 75 } 75 76 return (out); … … 127 128 out[i].R = in[i].R; 128 129 out[i].D = in[i].D; 129 out[i].Xp = in[i].Xp;130 130 131 131 // changed for PANSTARRS_DEV_0 (moved from Average to Measure) … … 170 170 out[i].photFlagsUpper = 0; 171 171 out[i].photFlagsLower = 0; 172 173 // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway) 174 out[i].psfQF = NAN; 175 out[i].psfQFperf = NAN; 176 out[i].stargal = NAN; 172 177 } 173 178 return (out); … … 185 190 out[i].R = in[i].R; 186 191 out[i].D = in[i].D; 187 out[i].Xp = in[i].Xp;188 192 189 193 // changed for PANSTARRS_DEV_0 (moved from Average to Measure) -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dvo_convert_loneos.c
r34260 r34382 39 39 out[i].Sky = 0; 40 40 out[i].dSky = 0; 41 out[i].psfQ ual= 0;41 out[i].psfQF = 0; 42 42 out[i].psfChisq = 0; 43 43 out[i].crNsigma = 0; … … 74 74 out[i].FluxKron = NAN; 75 75 out[i].dFluxKron = NAN; 76 out[i].psfQFperf = NAN; 76 77 } 77 78 return (out); … … 118 119 out[i].R = in[i].R; 119 120 out[i].D = in[i].D; 120 out[i].Xp = in[i].Xp;121 121 122 122 // added for ELIXIR … … 162 162 out[i].photFlagsUpper = 0; 163 163 out[i].photFlagsLower = 0; 164 165 // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway) 166 out[i].psfQF = NAN; 167 out[i].psfQFperf = NAN; 168 out[i].stargal = NAN; 164 169 } 165 170 return (out); … … 177 182 out[i].R = in[i].R; 178 183 out[i].D = in[i].D; 179 out[i].Xp = in[i].Xp;180 184 181 185 // changed for PANSTARRS_DEV_0 (moved from Average to Measure) -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dvo_convert_panstarrs_DEV_0.c
r34260 r34382 36 36 37 37 // changed or added for PS1_DEV_1 (2008.02.26) 38 out[i].psfQ ual = (in[i].psfQual == NAN_S_SHORT) ? NAN : in[i].psfQual;38 out[i].psfQF = (in[i].psfQF == NAN_S_SHORT) ? NAN : in[i].psfQF; 39 39 out[i].dbFlags = in[i].flags; 40 40 out[i].detID = in[i].detID_lo; … … 65 65 out[i].Mkron = NAN; 66 66 out[i].dMkron = NAN; 67 out[i].FluxPSF = NAN; 68 out[i].dFluxPSF = NAN; 69 out[i].FluxKron = NAN; 70 out[i].dFluxKron = NAN; 67 out[i].FluxPSF = NAN; 68 out[i].dFluxPSF = NAN; 69 out[i].FluxKron = NAN; 70 out[i].dFluxKron = NAN; 71 out[i].psfQFperf = NAN; 71 72 } 72 73 return (out); … … 105 106 // changed or added for PS1_DEV_1 (2008.02.26) 106 107 out[i].flags = in[i].dbFlags; 107 out[i].psfQ ual = isnan(in[i].psfQual) ? NAN_S_SHORT : in[i].psfQual;108 out[i].psfQF = isnan(in[i].psfQF) ? NAN_S_SHORT : in[i].psfQF; 108 109 out[i].detID_hi = 0; 109 110 out[i].detID_lo = in[i].detID; … … 133 134 out[i].R = in[i].R; 134 135 out[i].D = in[i].D; 135 out[i].Xp = in[i].Xp;136 136 out[i].dR = in[i].dR; 137 137 out[i].dD = in[i].dD; … … 166 166 out[i].photFlagsUpper = 0; 167 167 out[i].photFlagsLower = 0; 168 169 // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway) 170 out[i].psfQF = NAN; 171 out[i].psfQFperf = NAN; 172 out[i].stargal = NAN; 168 173 } 169 174 return (out); … … 182 187 out[i].R = in[i].R; 183 188 out[i].D = in[i].D; 184 out[i].Xp = in[i].Xp;185 189 out[i].dR = in[i].dR; 186 190 out[i].dD = in[i].dD; -
branches/eam_branches/ipp-20120805/Ohana/src/libdvo/src/dvo_convert_panstarrs_DEV_1.c
r34260 r34382 36 36 37 37 // changed or added for PS1_DEV_1 (2008.02.26) 38 out[i].psfQ ual = (in[i].psfQual == NAN_S_SHORT) ? NAN : in[i].psfQual;38 out[i].psfQF = (in[i].psfQF == NAN_S_SHORT) ? NAN : in[i].psfQF; 39 39 out[i].dbFlags = in[i].flags; 40 40 out[i].detID = in[i].detID_lo; … … 65 65 out[i].Mkron = NAN; 66 66 out[i].dMkron = NAN; 67 out[i].FluxPSF = NAN; 68 out[i].dFluxPSF = NAN; 69 out[i].FluxKron = NAN; 70 out[i].dFluxKron = NAN; 67 out[i].FluxPSF = NAN; 68 out[i].dFluxPSF = NAN; 69 out[i].FluxKron = NAN; 70 out[i].dFluxKron = NAN; 71 out[i].psfQFperf = NAN; 71 72 } 72 73 return (out); … … 105 106 // changed or added for PS1_DEV_1 (2008.02.26) 106 107 out[i].flags = in[i].dbFlags; 107 out[i].psfQ ual = isnan(in[i].psfQual) ? NAN_S_SHORT : in[i].psfQual;108 out[i].psfQF = isnan(in[i].psfQF) ? NAN_S_SHORT : in[i].psfQF; 108 109 out[i].detID_hi = 0; 109 110 out[i].detID_lo = in[i].detID; … … 133 134 out[i].R = in[i].R; 134 135 out[i].D = in[i].D; 135 out[i].Xp = in[i].Xp;136 136 out[i].dR = in[i].dR; 137 137 out[i].dD = in[i].dD; … … 166 166 out[i].photFlagsUpper = 0; 167 167 out[i].photFlagsLower = 0; 168 169 // added for PS1_V4 (Xp dropped in V4 onward, was not really used anyway) 170 out[i].psfQF = NAN; 171 out[i].psfQFperf = NAN; 172 out[i].stargal = NAN; 168 173 } 169 174 return (out); … … 182 187 out[i].R = in[i].R; 183 188 out[i].D = in[i].D; 184 out[i].Xp = in[i].Xp;185 189 out[i].dR = in[i].dR; 186 190 out[i].dD = in[i].dD; -
branches/eam_branches/ipp-20120805/Ohana/src/opihi/dvo/avextract.c
r34088 r34382 325 325 gprint (GP_ERR, " <photcode>:nphot : number of measurements used for average magnitude in this photcode\n"); 326 326 327 // gprint (GP_ERR, " Xp : NOT VALID\n");328 327 // gprint (GP_ERR, " type : dophot type (unused)\n"); 329 328 // gprint (GP_ERR, " typefrac : dophot type fraction (unused)\n"); -
branches/eam_branches/ipp-20120805/Ohana/src/opihi/dvo/avmatch.c
r33963 r34382 351 351 gprint (GP_ERR, " <photcode>:nphot : number of measurements used for average magnitude in this photcode\n"); 352 352 353 // gprint (GP_ERR, " Xp : NOT VALID\n");354 353 // gprint (GP_ERR, " type : dophot type (unused)\n"); 355 354 // gprint (GP_ERR, " typefrac : dophot type fraction (unused)\n"); -
branches/eam_branches/ipp-20120805/Ohana/src/opihi/dvo/gstar.c
r33963 r34382 190 190 gprint (GP_LOG, "%5.2f ", 3600.0*sqrt(r)); 191 191 gprint (GP_LOG, "%3d ", catalog.average[k].Nmeasure); 192 gprint (GP_LOG, "%4.1f ", 0.01*catalog.average[k].Xp);192 gprint (GP_LOG, "%4.1f ", catalog.average[k].ChiSqAve); 193 193 gprint (GP_LOG, "%5x ", catalog.average[k].flags); 194 194 gprint (GP_LOG, "%x ", catalog.average[k].objID); … … 205 205 gprint (GP_LOG, "%f ", catalog.average[k].dP); 206 206 207 gprint (GP_LOG, "%f ", catalog.average[k].ChiSqAve);208 207 gprint (GP_LOG, "%f ", catalog.average[k].ChiSqPM); 209 208 gprint (GP_LOG, "%f ", catalog.average[k].ChiSqPar); … … 292 291 gprint (GP_LOG, "0x%08x ", catalog.measure[m].detID); 293 292 gprint (GP_LOG, "0x%08x ", catalog.measure[m].imageID); 294 gprint (GP_LOG, "%.3f ", catalog.measure[m].psfQ ual);293 gprint (GP_LOG, "%.3f ", catalog.measure[m].psfQF); 295 294 gprint (GP_LOG, "%.1f ", catalog.measure[m].psfChisq); 296 295 gprint (GP_LOG, "%.1f ", catalog.measure[m].crNsigma); -
branches/eam_branches/ipp-20120805/Ohana/src/opihi/dvo/photometry.c
r31635 r34382 226 226 if (!strcasecmp (parname, "Nmeas")) param = AVE_NMEAS; 227 227 if (!strcasecmp (parname, "Nmiss")) param = AVE_NMISS; 228 if (!strcasecmp (parname, "Xp")) param = AVE_Xp;229 228 if (!strcasecmp (parname, "Xm")) param = AVE_Xm; 230 229 if (!strcasecmp (parname, "flag")) param = AVE_OBJ_FLAGS; … … 294 293 void GetAverageParamHelp () { 295 294 gprint (GP_ERR, "value may be one of the following:\n"); 296 gprint (GP_ERR, " ra dec dmag Nmeas Nmiss Xm XpNphot Ncode flag type typefrac\n\n");295 gprint (GP_ERR, " ra dec dmag Nmeas Nmiss Xm Nphot Ncode flag type typefrac\n\n"); 297 296 gprint (GP_ERR, "value may also be a valid photcode\n"); 298 297 gprint (GP_ERR, "photcodes or 'mag' may have optional magnitude mode: mag,[Mave, Mref]\n"); … … 637 636 case AVE_NMISS: 638 637 value = average[0].Nmissing; 639 break;640 case AVE_Xp:641 value = 0.01*average[0].Xp;642 638 break; 643 639 case AVE_OBJ_FLAGS: -
branches/eam_branches/ipp-20120805/Ohana/src/relastro/src/BrightCatalog.c
r33652 r34382 154 154 GET_COLUMN(Tmean, "MEAN_EPOCH", int); 155 155 GET_COLUMN(Trange, "TIME_RANGE", int); 156 GET_COLUMN( Xp, "SIGMA_POS",float);156 GET_COLUMN(stargal, "STARGAL_SEP", float); 157 157 GET_COLUMN(Npos, "NUMBER_POS", short); 158 158 GET_COLUMN(Nmeasure, "NMEASURE", short); … … 189 189 average[i].Tmean = Tmean[i] ; 190 190 average[i].Trange = Trange[i] ; 191 average[i]. Xp = Xp[i];191 average[i].stargal = stargal[i] ; 192 192 average[i].Npos = Npos[i] ; 193 193 average[i].Nmeasure = Nmeasure[i] ; … … 221 221 free (Tmean); 222 222 free (Trange); 223 free ( Xp);223 free (stargal); 224 224 free (Npos); 225 225 free (Nmeasure); … … 487 487 int *Tmean ; ALLOCATE (Tmean , int , catalog->Naverage); 488 488 int *Trange ; ALLOCATE (Trange , int , catalog->Naverage); 489 float * Xp ; ALLOCATE (Xp, float , catalog->Naverage);489 float *stargal ; ALLOCATE (stargal , float , catalog->Naverage); 490 490 short *Npos ; ALLOCATE (Npos , short , catalog->Naverage); 491 491 short *Nmeasure ; ALLOCATE (Nmeasure , short , catalog->Naverage); … … 520 520 Tmean[i] = average[i].Tmean ; 521 521 Trange[i] = average[i].Trange ; 522 Xp[i] = average[i].Xp;522 stargal[i] = average[i].stargal ; 523 523 Npos[i] = average[i].Npos ; 524 524 Nmeasure[i] = average[i].Nmeasure ; … … 552 552 gfits_set_bintable_column (&theader, &ftable, "MEAN_EPOCH", Tmean, catalog->Naverage); 553 553 gfits_set_bintable_column (&theader, &ftable, "TIME_RANGE", Trange, catalog->Naverage); 554 gfits_set_bintable_column (&theader, &ftable, "S IGMA_POS", Xp,catalog->Naverage);554 gfits_set_bintable_column (&theader, &ftable, "STARGAL_SEP", stargal, catalog->Naverage); 555 555 gfits_set_bintable_column (&theader, &ftable, "NUMBER_POS", Npos, catalog->Naverage); 556 556 gfits_set_bintable_column (&theader, &ftable, "NMEASURE", Nmeasure, catalog->Naverage); … … 582 582 free (Tmean); 583 583 free (Trange); 584 free ( Xp);584 free (stargal); 585 585 free (Npos); 586 586 free (Nmeasure); -
branches/eam_branches/ipp-20120805/Ohana/src/relastro/src/UpdateObjectOffsets.c
r33963 r34382 95 95 // RESET (-reset) 96 96 // TimeSelect -time 97 // (note that psfQ ualis applied rigidly at 0.85, as is the galaxy test)97 // (note that psfQF is applied rigidly at 0.85, as is the galaxy test) 98 98 // ImagSelect, ImagMin, ImagMax 99 99 // MaxDensityUse, MaxDensityValue -
branches/eam_branches/ipp-20120805/Ohana/src/relastro/src/UpdateObjects.c
r33652 r34382 352 352 catalog[i].average[j].dP = fit.dp; // parallax error in arcsec 353 353 354 // Xp is supposed to be the position scatter, not the chisq : fix this:355 // catalog[i].average[j].Xp = (fit.Nfit > 1) ? 100.0*log10(fit.chisq) : NAN_S_SHORT;356 354 catalog[i].average[j].ChiSqAve = fitAve.chisq; 357 355 catalog[i].average[j].ChiSqPM = fitPM.chisq; 358 356 catalog[i].average[j].ChiSqPar = fitPAR.chisq; 359 catalog[i].average[j].Xp = 0.0;360 357 catalog[i].average[j].Tmean = (Tmean * 86400 * 365.26) + T2000; 361 358 catalog[i].average[j].Trange = (Trange * 86400 * 365.26); -
branches/eam_branches/ipp-20120805/Ohana/src/relastro/src/load_catalogs.c
r33963 r34382 137 137 // RESET (-reset) 138 138 // TimeSelect -time 139 // (note that psfQ ualis applied rigidly at 0.85, as is the galaxy test)139 // (note that psfQF is applied rigidly at 0.85, as is the galaxy test) 140 140 // ImagSelect, ImagMin, ImagMax 141 141 // MaxDensityUse, MaxDensityValue -
branches/eam_branches/ipp-20120805/Ohana/src/relastro/src/relastro_objects.c
r34260 r34382 123 123 // RESET (-reset) 124 124 // TimeSelect -time 125 // (note that psfQ ualis applied rigidly at 0.85, as is the galaxy test)125 // (note that psfQF is applied rigidly at 0.85, as is the galaxy test) 126 126 // ImagSelect, ImagMin, ImagMax 127 127 // MaxDensityUse, MaxDensityValue -
branches/eam_branches/ipp-20120805/Ohana/src/relphot/src/bcatalog.c
r33963 r34382 117 117 118 118 // skip garbage measurements 119 if (catalog[0].measure[offset].psfQ ual< 0.85) { Npsfqf ++; continue; }119 if (catalog[0].measure[offset].psfQF < 0.85) { Npsfqf ++; continue; } 120 120 if (isnan(catalog[0].measure[offset].M)) { Nnan ++; continue; } 121 121 -
branches/eam_branches/ipp-20120805/Ohana/src/relphot/src/load_catalogs.c
r33963 r34382 136 136 // TimeSelect -time 137 137 // DophotSelect 138 // (note that psfQ ualis applied rigidly at 0.85, as is the galaxy test)138 // (note that psfQF is applied rigidly at 0.85, as is the galaxy test) 139 139 // MAG_LIM 140 140 // SIGMA_LIM -
branches/eam_branches/ipp-20120805/Ohana/src/relphot/src/relphot_objects.c
r33963 r34382 143 143 // TimeSelect -time 144 144 // DophotSelect 145 // (note that psfQ ualis applied rigidly at 0.85, as is the galaxy test)145 // (note that psfQF is applied rigidly at 0.85, as is the galaxy test) 146 146 // MAG_LIM 147 147 // SIGMA_LIM -
branches/eam_branches/ipp-20120805/Ohana/src/relphot/src/setMrelFinal.c
r34088 r34382 237 237 if (catalog[0].measure[m].photFlags & code->photomBadMask) goto skip; 238 238 if ((catalog[0].measure[m].photcode > 10000) && (catalog[0].measure[m].photcode < 10500)) { 239 if (catalog[0].measure[m].psfQ ual< 0.85) goto skip;239 if (catalog[0].measure[m].psfQF < 0.85) goto skip; 240 240 } 241 241 }
Note:
See TracChangeset
for help on using the changeset viewer.
