IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Jun 12, 2015, 10:28:36 PM (11 years ago)
Author:
eugene
Message:

merge changes from EAM branch ipp-20150419

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/Ohana/src/relphot/include/relphot.h

    r38441 r38466  
    9696
    9797typedef struct {
     98  double *flxlist;           // list of measure.mag values for a given star
     99  double *errlist;            // mag errors for a star
     100  double *wgtlist;            // weights to use for mean mags
     101  int    *ranking;            // weights to use for mean mags
     102  int    *measSeq;            // weights to use for mean mags
     103  int     Nlist;
     104} StatDataSet;
     105
     106typedef struct {
    98107  int Nfew;
    99108  int Ncode;
    … …  
    105114
    106115  // NOTE: the following arrays are (possibly) pre-allocated and carried down to each
    107   // thread.  The first 3 (list, dlist, wlist) are used in all relphot analyses; the
    108   // others are only used on the final output steps.
    109 
    110   double *Mpsflist;           // list of measure.mag values for a given star
    111   double *dpsflist;           // mag errors for a star
    112   double *wpsflist;           // weights to use for mean mags
    113 
    114   double *Maplist;            // ap mags for a star
    115   double *daplist;            // ap mags for a star
    116   double *waplist;            // ap mags for a star
    117 
    118   double *Mkronlist;          // kron mags for a star
    119   double *dkronlist;          // kron mag errors
    120   double *wkronlist;          // kron mag errors
     116  // thread.  The psfData are used in all relphot analyses; the others are only used on
     117  // the final output steps.
     118
     119  int Nsecfilt;
     120
     121  StatDataSet  *psfData; // one is allocated for each primary (average) photcode
     122  StatDataSet *aperData;
     123  StatDataSet *kronData;
    121124
    122125  double *psfqf_list;         // psfqf for all filters
    123126  double *psfqfperf_list;     // psfqfperf for all filters
    124127  double *stargal_list;       // stargal for all filters
     128
     129  int   *havePS1;          // this secfilt has synthetic mags
     130  int   *haveSYN;          // this secfilt has synthetic mags
     131  int   *measSYN;          // this measurement is the synthetic mag for this secfilt
     132  int   *needSYN;          // this secfilt mag should use synthetic mags
     133  float *minSYN;           // minimum synthetic mag below which synthetic should be forced
     134 
     135  float *psfQfMax;            // max psfQf value for this secfilt
     136  float *psfQfPerfMax;        // max psfQfperf value for this secfilt
     137
     138  int *Nmeas;                 // count of PS1 exposure (chip) measurements for this secfilt
     139  int *NmeasGood;             // count of PS1 exposure (chip) measurements for this secfilt
     140  int *Next;                  // count of PS1 exposure (chip) measurements for this secfilt
     141  int *NexpPS1;               // count of PS1 exposure (chip) measurements for this secfilt
     142  int *haveUbercal;           // does this secfilt have any ubercal data?
     143
     144  int *tessID;                // tess,proj,skycell to use for warp and diff analysis
     145  int *projID;
     146  int *skycellID;
     147
     148  float *minUbercalDist;
    125149
    126150  StatType psfstats;
    … …  
    375399void          liststats_setmode   PROTO((StatType *stats, char *strmode));
    376400int           liststats           PROTO((double *value, double *dvalue, double *wvalue, int N, StatType *stats));
     401int           liststats_init      PROTO((StatType *stats));
    377402Catalog      *load_catalogs       PROTO((SkyList *skylist, int *Ncatalog, int hostID, char *hostpath, char *syncfile));
    378403Catalog      *load_catalogs_parallel PROTO((SkyList *sky, int *Ncatalog, char *syncfile));
    … …  
    412437int           setMrel             PROTO((Catalog *catalog, int Ncatalog, FlatCorrectionTable *flatcorr));
    413438void          setMrelFinal        PROTO((Catalog *catalog, FlatCorrectionTable *flatcorr, int simpleAverage));
    414 int           setMrelOutput       PROTO((Catalog *catalog, int Ncatalog, int pass, FlatCorrectionTable *flatcorr));
     439int           setMrelOutput       PROTO((Catalog *catalog, int Ncatalog, FlatCorrectionTable *flatcorr));
    415440int           setMave             PROTO((Catalog *catalog, int Ncatalog));
    416441void          set_ZP              PROTO((double ZERO));
    … …  
    474499
    475500int print_measure_set_alt (Average *average, SecFilt *secfilt, Measure *measure);
    476 int setMrel_catalog_alt (Catalog *catalog, int Nc, int pass, FlatCorrectionTable *flatcorr, SetMrelInfo *results, int Nsecfilt);
    477 int setMrelAverageExposure (off_t meas, int cat, int pass, FlatCorrectionTable *flatcorr, SetMrelInfo *results, Average *average, AverageTiny *averageT, SecFilt *secfilt, Measure *measure, MeasureTiny *measureT, off_t *found);
    478 int setMrelAverageStack (off_t measureOffset, int cat, FlatCorrectionTable *flatcorr, SetMrelInfo *results, Average *average, SecFilt *secfilt, Measure *measure, off_t Nmeasure);
    479 int setMrelAverageForcedWarp (off_t measureOffset, int cat, int pass, SetMrelInfo *results, Average *average, SecFilt *secfilt, Measure *measure, off_t *found, off_t Nmeasure);
     501int setMrel_catalog_alt (Catalog *catalog, int Nc, int isSetMrelFinal, FlatCorrectionTable *flatcorr, SetMrelInfo *results, int Nsecfilt);
     502int setMrelAverageExposure (Catalog *catalog, int cat, off_t ave, int Nsecfilt, int isSetMrelFinal, FlatCorrectionTable *flatcorr, SetMrelInfo *results);
     503int setMrelAverageStack (Catalog *catalog, int cat, off_t ave, int Nsecfilt, FlatCorrectionTable *flatcorr);
     504int setMrelAverageForcedWarp (Catalog *catalog, int cat, off_t ave, int Nsecfilt, FlatCorrectionTable *flatcorr, SetMrelInfo *results);
    480505
    481506int setGlobalObjStats (Average *average, Measure *measure);
    482507
    483508void SetMrelInfoInit (SetMrelInfo *results, int allocLists);
     509void SetMrelInfoFree (SetMrelInfo *results);
     510void SetMrelInfoReset (SetMrelInfo *results);
     511void SetMrelInfoResetObject (SetMrelInfo *results);
     512
     513void StatDataSetFree (StatDataSet *dataset, int Nsecfilt);
     514StatDataSet *StatDataSetAlloc (int Nsecfilt, int Nmax);
    484515
    485516int init_synthetic_mags ();
    … …  
    533564int isGPC1synth (int photcode);
    534565int whichGPC1filter (int photcode);
    535 
    536 
    537 SynthZeroPoints *SynthZeroPointsLoad (char *filename);
     566int is2MASS (int photcode);
     567int isTYCHO (int photcode);
     568
     569int magStatsByRanking (StatDataSet *dataset, StatType *stats);
     570
     571int SynthZeroPointsLoad (char *filename);
     572SynthZeroPoints *SynthZeroPointsGet ();
    538573
    539574int relphot_synthphot (int hostID, char *hostpath);
    … …  
    542577int relphot_synthphot_catalog (Catalog *catalog, SynthZeroPoints *zpts);
    543578int relphot_synthphot_average (Average *average, SecFilt *secfilt, Measure *measure, SynthZeroPoints *zpts);
     579
     580void setMeasureRank (Catalog *catalog);
     581int getImageFlags (off_t meas, int cat);
     582int getMosaicFlags (off_t meas, int cat);
     583
Note: See TracChangeset for help on using the changeset viewer.