IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 34383


Ignore:
Timestamp:
Sep 4, 2012, 1:22:26 PM (14 years ago)
Author:
eugene
Message:

added cmf formats PS1_V4 and PS1_SV2 to both psModules and Ohana/libautocode; checked that they both work; checked that dvo schema PS1_V4 can be loaded and read

Location:
branches/eam_branches/ipp-20120805
Files:
3 added
12 edited

Legend:

Unmodified
Added
Removed

  • branches/eam_branches/ipp-20120805/Ohana/src/addstar/include/addstar.h

    r34260 r34383  
    269269int        NewReflist_Thread      PROTO((int BindSocket));
    270270
    271 Stars     *Convert_SMPDATA        PROTO((FTable *table, unsigned int *nstars));
    272 Stars     *Convert_PS1_DEV_0      PROTO((FTable *table, unsigned int *nstars));
    273 Stars     *Convert_PS1_DEV_1      PROTO((FTable *table, unsigned int *nstars));
    274 Stars     *Convert_PS1_V1         PROTO((FTable *table, unsigned int *nstars));
    275 Stars     *Convert_PS1_V1_Alt     PROTO((FTable *table, unsigned int *nstars));
    276 Stars     *Convert_PS1_V2         PROTO((FTable *table, unsigned int *nstars));
    277 Stars     *Convert_PS1_V3         PROTO((FTable *table, unsigned int *nstars));
    278 Stars     *Convert_PS1_SV1        PROTO((FTable *table, unsigned int *nstars));
    279 Stars     *Convert_PS1_SV1_Alt    PROTO((FTable *table, unsigned int *nstars));
    280 
    281271int        InitStar               PROTO((Stars *star));
    282272

  • branches/eam_branches/ipp-20120805/Ohana/src/addstar/src/MatchHeaders.c

    r31160 r34383  
    5858    if (!strcmp (exttype, "PS1_V2")) goto keep;
    5959    if (!strcmp (exttype, "PS1_V3")) goto keep;
     60    if (!strcmp (exttype, "PS1_V4")) goto keep;
    6061    if (!strcmp (exttype, "PS1_SV1")) goto keep;
     62    if (!strcmp (exttype, "PS1_SV2")) goto keep;
    6163    continue;
    6264

  • branches/eam_branches/ipp-20120805/Ohana/src/addstar/src/ReadStarsFITS.c

    r34382 r34383  
    11# include "addstar.h"
     2
     3// this is the complete list of FITS format input files the addstar knows
     4// (excluding SDSS data and reference database info, such as 2MASS)
     5// NOTE: these must also be listed in MatchHeaders.c (line ~ 62)
     6Stars     *Convert_SMPDATA        PROTO((FTable *table, unsigned int *nstars));
     7Stars     *Convert_PS1_DEV_0      PROTO((FTable *table, unsigned int *nstars));
     8Stars     *Convert_PS1_DEV_1      PROTO((FTable *table, unsigned int *nstars));
     9Stars     *Convert_PS1_V1         PROTO((FTable *table, unsigned int *nstars));
     10Stars     *Convert_PS1_V1_Alt     PROTO((FTable *table, unsigned int *nstars));
     11Stars     *Convert_PS1_V2         PROTO((FTable *table, unsigned int *nstars));
     12Stars     *Convert_PS1_V3         PROTO((FTable *table, unsigned int *nstars));
     13Stars     *Convert_PS1_V4         PROTO((FTable *table, unsigned int *nstars));
     14Stars     *Convert_PS1_SV1        PROTO((FTable *table, unsigned int *nstars));
     15Stars     *Convert_PS1_SV1_Alt    PROTO((FTable *table, unsigned int *nstars));
     16Stars     *Convert_PS1_SV2        PROTO((FTable *table, unsigned int *nstars));
    217
    318// given a file with the pointer at the start of the table block and the
    … …  
    5469    stars = Convert_PS1_V3 (&table, &Nstars);
    5570  }
     71  if (!strcmp (type, "PS1_V4")) {
     72    stars = Convert_PS1_V4 (&table, &Nstars);
     73  }
    5674  if (!strcmp (type, "PS1_SV1")) {
    5775    stars = Convert_PS1_SV1 (&table, &Nstars);
     76  }
     77  if (!strcmp (type, "PS1_SV2")) {
     78    stars = Convert_PS1_SV2 (&table, &Nstars);
    5879  }
    5980  if (stars == NULL) {
    … …  
    616637}
    617638
    618 Stars *Convert_PS1_SV1 (FTable *table, unsigned int *nstars) {
     639Stars *Convert_PS1_V4 (FTable *table, unsigned int *nstars) {
    619640
    620641  off_t Nstars;
    … …  
    622643  double ZeroPt;
    623644  Stars *stars;
    624   CMF_PS1_SV1 *ps1data;
    625 
    626   if (table[0].header[0].Naxis[0] == 196) {
    627     stars = Convert_PS1_SV1_Alt (table, nstars);
    628     return (stars);
    629   }
    630 
    631   ps1data = gfits_table_get_CMF_PS1_SV1 (table, &Nstars, NULL);
     645  CMF_PS1_V4 *ps1data;
     646
     647  ps1data = gfits_table_get_CMF_PS1_V4 (table, &Nstars, NULL);
    632648  if (!ps1data) {
    633649    fprintf (stderr, "skipping inconsistent entry\n");
    … …  
    658674    stars[i].measure.Mkron      = (ps1data[i].kronFlux > 0.0) ? -2.5*log10(ps1data[i].kronFlux) + ZeroPt : NAN;
    659675    stars[i].measure.dMkron     = (ps1data[i].kronFlux > 0.0) ? ps1data[i].kronFluxErr / ps1data[i].kronFlux : NAN;
    660 
     676                       
    661677    // these fluxes are converted from counts to counts/sec in FilterStars.c
    662678    stars[i].measure.FluxPSF    = ps1data[i].Flux;
    … …  
    689705
    690706    // the Average fields and the following Measure fields are set in FilterStars after
    691     // the image metadata is in hand:  dR, dD, Mcal, dt, airmass, az, t, imageID, extID,
    692     // averef is set in find_matches, dbFlags is zero on ingest.
     707    // the image metadata is in hand:  dR, dD, Mcal, dt, airmass, az, t, imageID, extID.
     708
     709    // averef is set in find_matches
     710
     711    // dbFlags is zero on ingest.
    693712
    694713    // the following fields are currently not being set anywhere: t_msec
    … …  
    698717}
    699718
    700 Stars *Convert_PS1_SV1_Alt (FTable *table, unsigned int *nstars) {
     719Stars *Convert_PS1_SV1 (FTable *table, unsigned int *nstars) {
    701720
    702721  off_t Nstars;
    … …  
    706725  CMF_PS1_SV1 *ps1data;
    707726
    708   // some test output files were produced called CMF_PS1_SV1 but with mismatch byte boundaries
    709 
    710   ps1data = gfits_table_get_CMF_PS1_SV1_Alt (table, &Nstars, NULL);
     727  if (table[0].header[0].Naxis[0] == 196) {
     728    stars = Convert_PS1_SV1_Alt (table, nstars);
     729    return (stars);
     730  }
     731
     732  ps1data = gfits_table_get_CMF_PS1_SV1 (table, &Nstars, NULL);
    711733  if (!ps1data) {
    712734    fprintf (stderr, "skipping inconsistent entry\n");
    … …  
    777799}
    778800
     801Stars *Convert_PS1_SV1_Alt (FTable *table, unsigned int *nstars) {
     802
     803  off_t Nstars;
     804  unsigned int i;
     805  double ZeroPt;
     806  Stars *stars;
     807  CMF_PS1_SV1 *ps1data;
     808
     809  // some test output files were produced called CMF_PS1_SV1 but with mismatch byte boundaries
     810
     811  ps1data = gfits_table_get_CMF_PS1_SV1_Alt (table, &Nstars, NULL);
     812  if (!ps1data) {
     813    fprintf (stderr, "skipping inconsistent entry\n");
     814    return (NULL);
     815  }
     816  ZeroPt = GetZeroPoint();
     817
     818  ALLOCATE (stars, Stars, Nstars);
     819  for (i = 0; i < Nstars; i++) {
     820    InitStar (&stars[i]);
     821    stars[i].measure.Xccd       = ps1data[i].X;
     822    stars[i].measure.Yccd       = ps1data[i].Y;
     823    stars[i].measure.dXccd      = ToShortPixels(ps1data[i].dX);
     824    stars[i].measure.dYccd      = ToShortPixels(ps1data[i].dY);
     825
     826    stars[i].measure.posangle   = ToShortDegrees(ps1data[i].posangle);
     827    stars[i].measure.pltscale   = ps1data[i].pltscale;
     828
     829    if ((ps1data[i].M >= 0.0) || isnan(ps1data[i].M)) {
     830      stars[i].measure.M      = NAN;
     831    } else {
     832      stars[i].measure.M      = ps1data[i].M + ZeroPt;
     833    }
     834    stars[i].measure.dM         = ps1data[i].dM;
     835    stars[i].measure.dMcal      = ps1data[i].dMcal;
     836    stars[i].measure.Map        = ps1data[i].Map + ZeroPt;
     837                       
     838    stars[i].measure.Mkron      = (ps1data[i].kronFlux > 0.0) ? -2.5*log10(ps1data[i].kronFlux) + ZeroPt : NAN;
     839    stars[i].measure.dMkron     = (ps1data[i].kronFlux > 0.0) ? ps1data[i].kronFluxErr / ps1data[i].kronFlux : NAN;
     840
     841    // these fluxes are converted from counts to counts/sec in FilterStars.c
     842    stars[i].measure.FluxPSF    = ps1data[i].Flux;
     843    stars[i].measure.dFluxPSF   = ps1data[i].dFlux;
     844    stars[i].measure.FluxKron   = ps1data[i].kronFlux;
     845    stars[i].measure.dFluxKron  = ps1data[i].kronFluxErr;
     846
     847    stars[i].measure.Sky        = ps1data[i].sky;
     848    stars[i].measure.dSky       = ps1data[i].dSky;
     849                       
     850    stars[i].measure.psfChisq   = ps1data[i].psfChisq;
     851    stars[i].measure.psfQF      = ps1data[i].psfQF;
     852    stars[i].measure.psfNdof    = ps1data[i].psfNdof;
     853    stars[i].measure.psfNpix    = ps1data[i].psfNpix;
     854    stars[i].measure.crNsigma   = ps1data[i].crNsigma;
     855    stars[i].measure.extNsigma  = ps1data[i].extNsigma;
     856
     857    stars[i].measure.FWx        = ToShortPixels(ps1data[i].fx);
     858    stars[i].measure.FWy        = ToShortPixels(ps1data[i].fy);
     859    stars[i].measure.theta      = ToShortDegrees(ps1data[i].df);
     860
     861    stars[i].measure.Mxx        = ToShortPixels(ps1data[i].Mxx);
     862    stars[i].measure.Mxy        = ToShortPixels(ps1data[i].Mxy);
     863    stars[i].measure.Myy        = ToShortPixels(ps1data[i].Myy);
     864                       
     865    stars[i].measure.photFlags  = ps1data[i].flags;
     866
     867    // this is may optionally be replaced by the internal sequence (see FilterStars.c)
     868    stars[i].measure.detID      = ps1data[i].detID;
     869
     870    // the Average fields and the following Measure fields are set in FilterStars after
     871    // the image metadata is in hand:  dR, dD, Mcal, dt, airmass, az, t, imageID, extID,
     872    // averef is set in find_matches, dbFlags is zero on ingest.
     873
     874    // the following fields are currently not being set anywhere: t_msec
     875  }   
     876  *nstars = Nstars;
     877  return (stars);
     878}
     879
     880Stars *Convert_PS1_SV2 (FTable *table, unsigned int *nstars) {
     881
     882  off_t Nstars;
     883  unsigned int i;
     884  double ZeroPt;
     885  Stars *stars;
     886  CMF_PS1_SV2 *ps1data;
     887
     888  ps1data = gfits_table_get_CMF_PS1_SV2 (table, &Nstars, NULL);
     889  if (!ps1data) {
     890    fprintf (stderr, "skipping inconsistent entry\n");
     891    return (NULL);
     892  }
     893  ZeroPt = GetZeroPoint();
     894
     895  ALLOCATE (stars, Stars, Nstars);
     896  for (i = 0; i < Nstars; i++) {
     897    InitStar (&stars[i]);
     898    stars[i].measure.Xccd       = ps1data[i].X;
     899    stars[i].measure.Yccd       = ps1data[i].Y;
     900    stars[i].measure.dXccd      = ToShortPixels(ps1data[i].dX);
     901    stars[i].measure.dYccd      = ToShortPixels(ps1data[i].dY);
     902
     903    stars[i].measure.posangle   = ToShortDegrees(ps1data[i].posangle);
     904    stars[i].measure.pltscale   = ps1data[i].pltscale;
     905
     906    if ((ps1data[i].M >= 0.0) || isnan(ps1data[i].M)) {
     907      stars[i].measure.M      = NAN;
     908    } else {
     909      stars[i].measure.M      = ps1data[i].M + ZeroPt;
     910    }
     911    stars[i].measure.dM         = ps1data[i].dM;
     912    stars[i].measure.dMcal      = ps1data[i].dMcal;
     913    stars[i].measure.Map        = ps1data[i].Map + ZeroPt;
     914                       
     915    stars[i].measure.Mkron      = (ps1data[i].kronFlux > 0.0) ? -2.5*log10(ps1data[i].kronFlux) + ZeroPt : NAN;
     916    stars[i].measure.dMkron     = (ps1data[i].kronFlux > 0.0) ? ps1data[i].kronFluxErr / ps1data[i].kronFlux : NAN;
     917
     918    // these fluxes are converted from counts to counts/sec in FilterStars.c
     919    stars[i].measure.FluxPSF    = ps1data[i].Flux;
     920    stars[i].measure.dFluxPSF   = ps1data[i].dFlux;
     921    stars[i].measure.FluxKron   = ps1data[i].kronFlux;
     922    stars[i].measure.dFluxKron  = ps1data[i].kronFluxErr;
     923
     924    stars[i].measure.Sky        = ps1data[i].sky;
     925    stars[i].measure.dSky       = ps1data[i].dSky;
     926                       
     927    stars[i].measure.psfChisq   = ps1data[i].psfChisq;
     928    stars[i].measure.psfQF      = ps1data[i].psfQF;
     929    stars[i].measure.psfNdof    = ps1data[i].psfNdof;
     930    stars[i].measure.psfNpix    = ps1data[i].psfNpix;
     931    stars[i].measure.crNsigma   = ps1data[i].crNsigma;
     932    stars[i].measure.extNsigma  = ps1data[i].extNsigma;
     933
     934    stars[i].measure.FWx        = ToShortPixels(ps1data[i].fx);
     935    stars[i].measure.FWy        = ToShortPixels(ps1data[i].fy);
     936    stars[i].measure.theta      = ToShortDegrees(ps1data[i].df);
     937
     938    stars[i].measure.Mxx        = ToShortPixels(ps1data[i].Mxx);
     939    stars[i].measure.Mxy        = ToShortPixels(ps1data[i].Mxy);
     940    stars[i].measure.Myy        = ToShortPixels(ps1data[i].Myy);
     941                       
     942    stars[i].measure.photFlags  = ps1data[i].flags;
     943
     944    // this is may optionally be replaced by the internal sequence (see FilterStars.c)
     945    stars[i].measure.detID      = ps1data[i].detID;
     946
     947    // the Average fields and the following Measure fields are set in FilterStars after
     948    // the image metadata is in hand:  dR, dD, Mcal, dt, airmass, az, t, imageID, extID,
     949    // averef is set in find_matches, dbFlags is zero on ingest.
     950
     951    // the following fields are currently not being set anywhere: t_msec
     952  }   
     953  *nstars = Nstars;
     954  return (stars);
     955}
     956
     957

  • branches/eam_branches/ipp-20120805/Ohana/src/addstar/src/mkcmf.c

    r34382 r34383  
    1414void gauss_init (int Nbin);
    1515double rnd_gauss (double mean, double sigma);
     16void writeStars_PS1_V4 (FTable *ftable, double *X, double *Y, double *M, unsigned int *Flag, int Nstars);
    1617void writeStars_PS1_V3 (FTable *ftable, double *X, double *Y, double *M, unsigned int *Flag, int Nstars);
    1718void writeStars_PS1_V2 (FTable *ftable, double *X, double *Y, double *M, unsigned int *Flag, int Nstars);
    … …  
    281282  if (!strcmp(type, "PS1_V3")) {
    282283    writeStars_PS1_V3 (&ftable, X, Y, M, Flag, Nstars);
     284    found = TRUE;
     285  }
     286  if (!strcmp(type, "PS1_V4")) {
     287    writeStars_PS1_V4 (&ftable, X, Y, M, Flag, Nstars);
    283288    found = TRUE;
    284289  }
    … …  
    635640}
    636641
     642void writeStars_PS1_V4 (FTable *ftable, double *X, double *Y, double *M, unsigned int *Flag, int Nstars) {
     643
     644  int i;
     645  CMF_PS1_V4 *stars;
     646  float flux, fSN;
     647
     648  // XXX add gaussian-distributed noise based on counts
     649  // this needs to make different output 'stars' entries depending on the desired type
     650  ALLOCATE (stars, CMF_PS1_V4, Nstars);
     651  gauss_init (2048);
     652  for (i = 0; i < Nstars; i++) {
     653    stars[i].detID = i;
     654
     655    flux = pow (10.0, -0.4*M[i]);
     656    fSN = 1.0 / sqrt(flux);
     657
     658    stars[i].X = X[i];
     659    stars[i].Y = Y[i];
     660    stars[i].M = M[i];
     661    stars[i].Map = M[i] - 0.05;
     662
     663    if (ADDNOISE) {
     664      stars[i].X += FX * fSN * rnd_gauss(0.0, 1.0);
     665      stars[i].Y += FY * fSN * rnd_gauss(0.0, 1.0);
     666      stars[i].M += fSN*rnd_gauss(0.0, 1.0);
     667    }
     668
     669    // randomly give poor PSFQF values
     670    if ((BAD_PSFQF_FRAC > 0.0) && (drand48() < BAD_PSFQF_FRAC)) {
     671      stars[i].psfQF     = 0.25;
     672      stars[i].psfQFperf = 0.24;
     673    } else {
     674      stars[i].psfQF     = PSFQUAL;
     675      stars[i].psfQFperf = MAX(PSFQUAL - 0.01, 0.0);
     676    }
     677   
     678    stars[i].dX = FX * fSN;
     679    stars[i].dY = FY * fSN;
     680    stars[i].dM = fSN;
     681
     682    stars[i].Mpeak     = M[i] + 1.0;
     683    stars[i].sky       = SKY;
     684    stars[i].dSky      = DSKY;
     685    stars[i].psfChisq  = PSFCHI;
     686    stars[i].crNsigma  = CRN;
     687    stars[i].extNsigma = EXTN;
     688    stars[i].fx        = FX;
     689    stars[i].fy        = FY;
     690    stars[i].df        = DF;
     691    stars[i].nFrames   = 1;
     692    stars[i].flags     = Flag[i];
     693
     694    stars[i].kronFlux  = flux * 1.25;
     695    stars[i].kronFluxErr = fSN * flux * 1.25;
     696  }
     697
     698  gfits_table_set_CMF_PS1_V4 (ftable, stars, Nstars);
     699  gfits_modify (ftable->header, "EXTTYPE",   "%s", 1, "PS1_V4");
     700}
     701

  • branches/eam_branches/ipp-20120805/Ohana/src/addstar/test/simple.dvo

    r34260 r34383  
    44
    55macro test.all
     6  # test.fields (DVO schema) (CMF format)
    67  test.fields PS1_DEV_0 PS1_V1
    78  test.fields PS1_DEV_1 PS1_V1
    … …  
    910  test.fields PS1_V2    PS1_V1
    1011  test.fields PS1_V3    PS1_V1
     12  test.fields PS1_V4    PS1_V1
    1113
    1214  test.fields PS1_DEV_0 PS1_V2
    … …  
    1517  test.fields PS1_V2    PS1_V2
    1618  test.fields PS1_V3    PS1_V2
     19  test.fields PS1_V4    PS1_V4
    1720
    1821  test.fields PS1_DEV_0 PS1_V3
    … …  
    2124  test.fields PS1_V2    PS1_V3
    2225  test.fields PS1_V3    PS1_V3
     26  test.fields PS1_V4    PS1_V4
    2327
    2428  test.fields PS1_DEV_0 PS1_V4
    … …  
    2731  test.fields PS1_V2    PS1_V4
    2832  test.fields PS1_V3    PS1_V4
     33  test.fields PS1_V4    PS1_V4
    2934end 
    3035
    … …  
    3641  end
    3742
    38   tapPLAN 51
     43  echo ""
     44  echo "testing --- CMF = $1, DVO = $2 ---"
    3945
    4046  exec rm -rf catdir.test
    … …  
    4551  mkinput
    4652  exec mkcmf test.in.txt test.cmf -date 2008/1/1 -time 01:00:00 -radec $RA $DEC -type $1
    47   exec addstar -D CATDIR catdir.test -D CAMERA simtest test.cmf -D CATFORMAT $2 -quick-airmass
     53  if ($TAP_VERBOSE)
     54    exec addstar -D CATDIR catdir.test -D CAMERA simtest test.cmf -D CATFORMAT $2 -quick-airmass
     55  else
     56    exec addstar -D CATDIR catdir.test -D CAMERA simtest test.cmf -D CATFORMAT $2 -quick-airmass >& tmp.log
     57  end
    4858
    4959  list testfields -copy testfields_$1
     60
     61  # set the number of tests based on the testfields list
     62  tapPLAN {3 + 3*($testfields:n - 1)}
    5063
    5164  catdir catdir.test
    … …  
    7790  for i 1 $testfields:n
    7891    list name -split $testfields:$i
    79     if ("$name:0" == "SKIP") continue
    80     if ("$name:2" == "SKIP") continue
     92    if ("$name:0" == "SKIP")
     93      tapSKIP 3
     94      continue
     95    end
     96    if ("$name:2" == "SKIP")
     97      tapSKIP 3
     98      continue
     99    end
    81100    # must be the same length
    82101    set v1 = $name:0
    … …  
    106125    # THETA is stored to only (360/65536) deg accuracy
    107126    if ("$name:0" == "PSF_THETA")
    108       echo $MEAN
     127      # echo $MEAN
    109128      tapOK {abs($MEAN)  < 0.006} "$name:0 vs $name:2 (MEAN)"
    110129      tapOK {abs($SIGMA) < 0.001} "$name:0 vs $name:2 (SIGMA)"
    … …  
    304323  N_FRAMES          : SKIP # not ingested into DVO
    305324end
     325
     326# this list is good for PS1_V4
     327list testfields_PS1_V4
     328  IPP_IDET          : detid     
     329  X_PSF             : xccd     
     330  Y_PSF             : yccd     
     331  X_PSF_SIG         : xccd:err # FAIL
     332  Y_PSF_SIG         : yccd:err # FAIL
     333  POSANGLE          : SKIP # astrometry is not calibrated in the cmf
     334  PLTSCALE          : SKIP # astrometry is not calibrated in the cmf
     335  PSF_INST_MAG      : mag:inst 
     336  PSF_INST_MAG_SIG  : mag:err   
     337  PSF_INST_FLUX     : SKIP # not ingested into DVO
     338  PSF_INST_FLUX_SIG : SKIP # not ingested into DVO
     339  AP_MAG_STANDARD   : mag:aperinst # FAIL
     340  AP_MAG_RAW        : SKIP # not ingested into DVO
     341  AP_MAG_RADIUS     : SKIP # not ingested into DVO
     342  CAL_PSF_MAG       : SKIP # photometry is not calibrated in the cmf
     343  CAL_PSF_MAG_SIG   : SKIP # photometry is not calibrated in the cmf
     344  RA_PSF            : SKIP # astrometry is not calibrated in the cmf
     345  DEC_PSF           : SKIP # astrometry is not calibrated in the cmf
     346  PEAK_FLUX_AS_MAG  : SKIP # not ingested into DVO
     347  SKY               : sky       
     348  SKY_SIG           : sky_err   
     349  PSF_CHISQ         : psf_chisq
     350  CR_NSIGMA         : cr_nsigma
     351  EXT_NSIGMA        : ext_nsigma
     352  PSF_MAJOR         : FWHM_MAJ 
     353  PSF_MINOR         : FWHM_MIN 
     354  PSF_THETA         : THETA # FAIL
     355  PSF_QF            : PSF_QF   
     356  PSF_QF_PERFECT    : SKIP # not ingested into DVO
     357  PSF_NDOF          : PSF_NDOF
     358  PSF_NPIX          : PSF_NPIX
     359  MOMENTS_XX        : MXX
     360  MOMENTS_XY        : MXY
     361  MOMENTS_YY        : MYY
     362  MOMENTS_M3C       : SKIP # not ingested into DVO
     363  MOMENTS_M3S       : SKIP # not ingested into DVO
     364  MOMENTS_M4C       : SKIP # not ingested into DVO
     365  MOMENTS_M4S       : SKIP # not ingested into DVO
     366  MOMENTS_R1        : SKIP # not ingested into DVO
     367  MOMENTS_RH        : SKIP # not ingested into DVO
     368  KRON_FLUX         : mag:kroninst
     369  KRON_FLUX_ERR     : mag:kronerr
     370  KRON_FLUX_INNER   : SKIP # not ingested into DVO
     371  KRON_FLUX_OUTER   : SKIP # not ingested into DVO
     372  FLAGS             : phot_flags
     373  N_FRAMES          : SKIP # not ingested into DVO
     374end

  • branches/eam_branches/ipp-20120805/Ohana/src/addstar/test/tap.dvo

    r33653 r34383  
    22
    33if (not($?TAP_BREAK)) set TAP_BREAK = 0
     4if (not($?TAP_VERBOSE)) set TAP_VERBOSE = 0
    45
    56macro tapOK
    … …  
    1011
    1112  if ($1)
    12     echo "ok : $2"
     13    if ($TAP_VERBOSE)
     14      echo "ok : $2"
     15    end
    1316    $TAP_LAST = 1
    1417  else
    … …  
    4245  end
    4346
    44   $TAP_NSKIP ++
    45   $TAP_NDONE ++
     47  $TAP_NSKIP += $1
     48  $TAP_NDONE += $1
    4649end
    4750

  • branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/measure-ps1-v4.d

    r34382 r34383  
    4848# do we need more resolution than a short? should this be a log?
    4949FIELD psfQF,          PSF_QF,        float,          psf coverage/quality factor
    50 FIELD psfQFPerf,      PSF_QF_PEFECT, float,          psf coverage / quality factor (all mask bits)
     50FIELD psfQFperf,      PSF_QF_PEFECT, float,          psf coverage / quality factor (all mask bits)
    5151FIELD psfChisq,       PSF_CHISQ,     float,          psf fit chisq
    5252

  • branches/eam_branches/ipp-20120805/Ohana/src/libautocode/def/measure.d

    r34382 r34383  
    4848# do we need more resolution than a short? should this be a log?
    4949FIELD psfQF,          PSF_QF,        float,          psf coverage/quality factor
    50 FIELD psfQFPerf,      PSF_QF_PEFECT, float,          psf coverage / quality factor (all mask bits)
     50FIELD psfQFperf,      PSF_QF_PEFECT, float,          psf coverage / quality factor (all mask bits)
    5151FIELD psfChisq,       PSF_CHISQ,     float,          psf fit chisq
    5252

  • branches/eam_branches/ipp-20120805/psModules/src/objects/Makefile.am

    r34381 r34383  
    4848        pmSourceIO_CMF_PS1_V4.c \
    4949        pmSourceIO_CMF_PS1_SV1.c \
     50        pmSourceIO_CMF_PS1_SV2.c \
    5051        pmSourceIO_CMF_PS1_DV1.c \
    5152        pmSourceIO_CMF_PS1_DV2.c \
    … …  
    136137
    137138# pmSourceID_CMF_* functions use a common framework
    138 BUILT_SOURCES = pmSourceIO_CMF_PS1_V1.c pmSourceIO_CMF_PS1_V2.c pmSourceIO_CMF_PS1_V3.c pmSourceIO_CMF_PS1_DV1.c pmSourceIO_CMF_PS1_DV2.c pmSourceIO_CMF_PS1_DV3.c pmSourceIO_CMF_PS1_SV1.c
     139BUILT_SOURCES = \
     140pmSourceIO_CMF_PS1_V1.c \
     141pmSourceIO_CMF_PS1_V2.c \
     142pmSourceIO_CMF_PS1_V3.c \
     143pmSourceIO_CMF_PS1_V4.c \
     144pmSourceIO_CMF_PS1_DV1.c \
     145pmSourceIO_CMF_PS1_DV2.c \
     146pmSourceIO_CMF_PS1_DV3.c \
     147pmSourceIO_CMF_PS1_SV1.c \
     148pmSourceIO_CMF_PS1_SV2.c
    139149
    140150pmSourceIO_CMF_PS1_V1.c : pmSourceIO_CMF.c.in mksource.pl

  • branches/eam_branches/ipp-20120805/psModules/src/objects/pmSourceIO.c

    r34376 r34383  
    575575            PM_SOURCES_WRITE("PS1_V4",    CMF_PS1_V4);
    576576            PM_SOURCES_WRITE("PS1_SV1",   CMF_PS1_SV1);
     577            PM_SOURCES_WRITE("PS1_SV2",   CMF_PS1_SV2);
    577578            PM_SOURCES_WRITE("PS1_DV1",   CMF_PS1_DV1);
    578579            PM_SOURCES_WRITE("PS1_DV2",   CMF_PS1_DV2);
     580            PM_SOURCES_WRITE("PS1_DV3",   CMF_PS1_DV3);
    579581
    580582            psFree (outhead);
    … …  
    10481050            PM_SOURCES_READ_PSF("PS1_V4",    CMF_PS1_V4);
    10491051            PM_SOURCES_READ_PSF("PS1_SV1",   CMF_PS1_SV1);
     1052            PM_SOURCES_READ_PSF("PS1_SV2",   CMF_PS1_SV2);
    10501053            PM_SOURCES_READ_PSF("PS1_DV1",   CMF_PS1_DV1);
    10511054            PM_SOURCES_READ_PSF("PS1_DV2",   CMF_PS1_DV2);
     1055            PM_SOURCES_READ_PSF("PS1_DV3",   CMF_PS1_DV3);
    10521056
    10531057            long *sourceIndex = NULL;
    … …  
    12541258        PM_SOURCES_READ_XSRC("PS1_V4",    CMF_PS1_V4);
    12551259        PM_SOURCES_READ_XSRC("PS1_SV1",   CMF_PS1_SV1);
     1260        PM_SOURCES_READ_XSRC("PS1_SV2",   CMF_PS1_SV2);
    12561261        PM_SOURCES_READ_XSRC("PS1_DV1",   CMF_PS1_DV1);
    12571262        PM_SOURCES_READ_XSRC("PS1_DV2",   CMF_PS1_DV2);
     1263        PM_SOURCES_READ_XSRC("PS1_DV3",   CMF_PS1_DV3);
    12581264    }
    12591265    psFree(tableHeader);
    … …  
    12901296        PM_SOURCES_READ_XFIT("PS1_V4",    CMF_PS1_V4);
    12911297        PM_SOURCES_READ_XFIT("PS1_SV1",   CMF_PS1_SV1);
     1298        PM_SOURCES_READ_XFIT("PS1_SV2",   CMF_PS1_SV2);
    12921299        PM_SOURCES_READ_XFIT("PS1_DV1",   CMF_PS1_DV1);
    12931300        PM_SOURCES_READ_XFIT("PS1_DV2",   CMF_PS1_DV2);
     1301        PM_SOURCES_READ_XFIT("PS1_DV3",   CMF_PS1_DV3);
    12941302    }
    12951303    psFree(tableHeader);
    … …  
    13251333        PM_SOURCES_READ_XRAD("PS1_V4",    CMF_PS1_V4);
    13261334        PM_SOURCES_READ_XRAD("PS1_SV1",   CMF_PS1_SV1);
     1335        PM_SOURCES_READ_XRAD("PS1_SV2",   CMF_PS1_SV2);
    13271336        PM_SOURCES_READ_XRAD("PS1_DV1",   CMF_PS1_DV1);
    13281337        PM_SOURCES_READ_XRAD("PS1_DV2",   CMF_PS1_DV2);
     1338        PM_SOURCES_READ_XRAD("PS1_DV3",   CMF_PS1_DV3);
    13291339    }
    13301340    psFree(tableHeader);

  • branches/eam_branches/ipp-20120805/psModules/src/objects/pmSourceIO.h

    r34367 r34383  
    3636MK_PROTO(CMF_PS1_V4);
    3737MK_PROTO(CMF_PS1_SV1);
     38MK_PROTO(CMF_PS1_SV2);
    3839MK_PROTO(CMF_PS1_DV1);
    3940MK_PROTO(CMF_PS1_DV2);

  • branches/eam_branches/ipp-20120805/psModules/src/objects/pmSourceIO_CMF.c.in

    r34381 r34383  
    142142        @>PS1_DV1,>PS1_V3,>PS1_SV1@ psMetadataAdd (row, PS_LIST_TAIL, "AP_FLUX_SIG",      PS_DATA_F32, "aperture flux error",                        source->apFluxErr);
    143143
    144         @<PS1_V3,PS1_SV?,PS1_DV?@ psMetadataAdd (row, PS_LIST_TAIL, "PEAK_FLUX_AS_MAG", PS_DATA_F32, "Peak flux expressed as magnitude",           outputs.peakMag);
     144        @<PS1_V3,PS1_SV1,PS1_DV?@ psMetadataAdd (row, PS_LIST_TAIL, "PEAK_FLUX_AS_MAG", PS_DATA_F32, "Peak flux expressed as magnitude",           outputs.peakMag);
    145145
    146146        @ALL@                     psMetadataAdd (row, PS_LIST_TAIL, "CAL_PSF_MAG",      PS_DATA_F32, "PSF Magnitude using supplied calibration",   outputs.calMag);
    … …  
    151151        @ALL,!PS1_V1@             psMetadataAdd (row, PS_LIST_TAIL, "DEC_PSF",          PS_DATA_F64, "PSF DEC coordinate (degrees)",               outputs.dec);
    152152
    153         @>=PS1_V3@                psMetadataAdd (row, PS_LIST_TAIL, "PEAK_FLUX_AS_MAG", PS_DATA_F32, "Peak flux expressed as magnitude",           outputs.peakMag);
     153        @>=PS1_V3,>PS1_SV1@       psMetadataAdd (row, PS_LIST_TAIL, "PEAK_FLUX_AS_MAG", PS_DATA_F32, "Peak flux expressed as magnitude",           outputs.peakMag);
    154154        @ALL@                     psMetadataAdd (row, PS_LIST_TAIL, "SKY",              PS_DATA_F32, "Sky level",                                  source->sky);
    155155        @ALL@                     psMetadataAdd (row, PS_LIST_TAIL, "SKY_SIGMA",        PS_DATA_F32, "Sigma of sky level",                         source->skyErr);
    … …  
    200200        @ALL@                      psMetadataAdd (row, PS_LIST_TAIL, "FLAGS",            PS_DATA_U32, "psphot analysis flags",                      source->mode);
    201201        @>PS1_V2,PS1_SV?,>PS1_DV1@ psMetadataAdd (row, PS_LIST_TAIL, "FLAGS2",           PS_DATA_U32, "psphot analysis flags",                      source->mode2);
    202         @>PS1_V2@                  psMetadataAdd (row, PS_LIST_TAIL, "PADDING2",         PS_DATA_S32, "more padding", 0);
     202        @PS1_V3,PS1_SV2@           psMetadataAdd (row, PS_LIST_TAIL, "PADDING2",         PS_DATA_S32, "more padding", 0);
    203203        @PS1_SV?@
    204204
Note: See TracChangeset for help on using the changeset viewer.