IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 38635


Ignore:
Timestamp:
Jul 24, 2015, 5:32:13 AM (11 years ago)
Author:
eugene
Message:

adding dvo methods for galphot

Location:
trunk/Ohana/src/libdvo
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/Ohana/src/libdvo/include/dvo.h

    r38553 r38635  
    12161216# undef LENSFIELD
    12171217
     1218# define GALPHOT_FIELD(NAME, VALUE) float GalphotValue_##NAME (PhotCode *code, GalPhot *galphot, int Ngalphot);
     1219
     1220GALPHOT_FIELD(GAL_MAG,         mag);           
     1221GALPHOT_FIELD(GAL_MAG_ERR,     magErr);   
     1222GALPHOT_FIELD(GAL_MAJ,         majorAxis);             
     1223GALPHOT_FIELD(GAL_MAJ_ERR,     majorAxisErr);
     1224GALPHOT_FIELD(GAL_MIN,         minorAxis);             
     1225GALPHOT_FIELD(GAL_MIN_ERR,     minorAxisErr);   
     1226GALPHOT_FIELD(GAL_THETA,       theta);     
     1227GALPHOT_FIELD(GAL_THETA_ERR,   thetaErr); 
     1228GALPHOT_FIELD(GAL_INDEX,       index);     
     1229GALPHOT_FIELD(GAL_CHISQ,       chisq);     
     1230GALPHOT_FIELD(GAL_NPIX,        Npix);       
     1231GALPHOT_FIELD(GAL_TYPE,        modelType);       
     1232
     1233# undef GALPHOT_FIELD
     1234
    12181235# endif // DVO_H

  • trunk/Ohana/src/libdvo/include/dvodb.h

    r38601 r38635  
    9797  MAG_OPTION_E1,
    9898  MAG_OPTION_E2,
     99
     100  MAG_OPTION_GAL_MAG,
     101  MAG_OPTION_GAL_MAG_ERR,
     102  MAG_OPTION_GAL_MAJ,
     103  MAG_OPTION_GAL_MAJ_ERR,
     104  MAG_OPTION_GAL_MIN,
     105  MAG_OPTION_GAL_MIN_ERR,
     106  MAG_OPTION_GAL_THETA,
     107  MAG_OPTION_GAL_THETA_ERR,
     108  MAG_OPTION_GAL_INDEX,
     109  MAG_OPTION_GAL_CHISQ,
     110  MAG_OPTION_GAL_NPIX,
     111  MAG_OPTION_GAL_TYPE,       
    99112} dvoMagOptionType;
    100113
    … …  
    304317      AVE_REF_COLOR_RED,
    305318      AVE_PHOT, // photometry class of values
    306       AVE_GAL_MAG,
    307       AVE_GAL_MAG_ERR,
    308       AVE_GAL_MAJ,
    309       AVE_GAL_MAJ_ERR,
    310       AVE_GAL_MIN,
    311       AVE_GAL_MIN_ERR,
    312       AVE_GAL_THETA,
    313       AVE_GAL_THETA_ERR,
    314       AVE_GAL_INDEX,
    315       AVE_GAL_CHISQ,
    316       AVE_GAL_NPIX,
    317       AVE_GAL_TYPE,       
    318319      AVE_E_BV,               // extinction (mags)
    319320      AVE_E_BV_ERR,

  • trunk/Ohana/src/libdvo/src/dbExtractAverages.c

    r38601 r38635  
    407407        case MAG_OPTION_NONE:
    408408          break;
     409
     410        case MAG_OPTION_GAL_MAG:       { value.Flt = GalphotValue_GAL_MAG      (field->photcode, galphot, average->Ngalphot); break; }
     411        case MAG_OPTION_GAL_MAG_ERR:   { value.Flt = GalphotValue_GAL_MAG_ERR  (field->photcode, galphot, average->Ngalphot); break; }
     412        case MAG_OPTION_GAL_MAJ:       { value.Flt = GalphotValue_GAL_MAJ      (field->photcode, galphot, average->Ngalphot); break; }
     413        case MAG_OPTION_GAL_MAJ_ERR:   { value.Flt = GalphotValue_GAL_MAJ_ERR  (field->photcode, galphot, average->Ngalphot); break; }
     414        case MAG_OPTION_GAL_MIN:       { value.Flt = GalphotValue_GAL_MIN      (field->photcode, galphot, average->Ngalphot); break; }
     415        case MAG_OPTION_GAL_MIN_ERR:   { value.Flt = GalphotValue_GAL_MIN_ERR  (field->photcode, galphot, average->Ngalphot); break; }
     416        case MAG_OPTION_GAL_THETA:     { value.Flt = GalphotValue_GAL_THETA    (field->photcode, galphot, average->Ngalphot); break; }
     417        case MAG_OPTION_GAL_THETA_ERR: { value.Flt = GalphotValue_GAL_THETA_ERR(field->photcode, galphot, average->Ngalphot); break; }
     418        case MAG_OPTION_GAL_INDEX:     { value.Flt = GalphotValue_GAL_INDEX    (field->photcode, galphot, average->Ngalphot); break; }
     419        case MAG_OPTION_GAL_CHISQ:     { value.Flt = GalphotValue_GAL_CHISQ    (field->photcode, galphot, average->Ngalphot); break; }
     420        case MAG_OPTION_GAL_NPIX:      { value.Flt = GalphotValue_GAL_NPIX     (field->photcode, galphot, average->Ngalphot); break; }
     421        case MAG_OPTION_GAL_TYPE:      { value.Flt = GalphotValue_GAL_TYPE     (field->photcode, galphot, average->Ngalphot); break; }
    409422
    410423        default:
    … …  
    427440    case AVE_DEC_GALMODEL:     { value.Flt = starpar->D;        break; }
    428441
    429     case AVE_GAL_MAG:          { value.Flt = NAN;               break; }
    430     case AVE_GAL_MAG_ERR:      { value.Flt = NAN;               break; }
    431     case AVE_GAL_MAJ:          { value.Flt = NAN;               break; }
    432     case AVE_GAL_MAJ_ERR:      { value.Flt = NAN;               break; }
    433     case AVE_GAL_MIN:          { value.Flt = NAN;               break; }
    434     case AVE_GAL_MIN_ERR:      { value.Flt = NAN;               break; }
    435     case AVE_GAL_THETA:        { value.Flt = NAN;               break; }
    436     case AVE_GAL_THETA_ERR:    { value.Flt = NAN;               break; }
    437     case AVE_GAL_INDEX:        { value.Flt = NAN;               break; }
    438     case AVE_GAL_CHISQ:        { value.Flt = NAN;               break; }
    439     case AVE_GAL_NPIX:         { value.Flt = NAN;               break; }
    440     case AVE_GAL_TYPE:         { value.Flt = NAN;               break; }
    441 
    442442    case AVE_TYPE:
    443443      break;

  • trunk/Ohana/src/libdvo/src/dbExtractMeasures.c

    r38153 r38635  
    284284        case MAG_OPTION_E2:
    285285
     286        case MAG_OPTION_GAL_MAG:       
     287        case MAG_OPTION_GAL_MAG_ERR:   
     288        case MAG_OPTION_GAL_MAJ:       
     289        case MAG_OPTION_GAL_MAJ_ERR:   
     290        case MAG_OPTION_GAL_MIN:       
     291        case MAG_OPTION_GAL_MIN_ERR:   
     292        case MAG_OPTION_GAL_THETA:     
     293        case MAG_OPTION_GAL_THETA_ERR:
     294        case MAG_OPTION_GAL_INDEX:     
     295        case MAG_OPTION_GAL_CHISQ:     
     296        case MAG_OPTION_GAL_NPIX:     
     297        case MAG_OPTION_GAL_TYPE:     
     298
    286299        case MAG_OPTION_NONE:
    287300          break;

  • trunk/Ohana/src/libdvo/src/dbFields.c

    r38601 r38635  
    164164  int i;
    165165  for (i = 0; i < Nfields; i++) {
    166     if (fields[i].ID == AVE_GAL_MAG         ) return TRUE;
    167     if (fields[i].ID == AVE_GAL_MAG_ERR     ) return TRUE;
    168     if (fields[i].ID == AVE_GAL_MAJ         ) return TRUE;
    169     if (fields[i].ID == AVE_GAL_MAJ_ERR     ) return TRUE;
    170     if (fields[i].ID == AVE_GAL_MIN         ) return TRUE;
    171     if (fields[i].ID == AVE_GAL_MIN_ERR     ) return TRUE;
    172     if (fields[i].ID == AVE_GAL_THETA       ) return TRUE;
    173     if (fields[i].ID == AVE_GAL_THETA_ERR   ) return TRUE;
    174     if (fields[i].ID == AVE_GAL_INDEX       ) return TRUE;
    175     if (fields[i].ID == AVE_GAL_CHISQ       ) return TRUE;
    176     if (fields[i].ID == AVE_GAL_NPIX        ) return TRUE;
    177     if (fields[i].ID == AVE_GAL_TYPE        ) return TRUE;
     166    if (fields[i].magOption == MAG_OPTION_NONE) continue; // non-measure fields do not have a photcode
     167
     168    if (fields[i].magOption == MAG_OPTION_GAL_MAG         ) return TRUE;
     169    if (fields[i].magOption == MAG_OPTION_GAL_MAG_ERR     ) return TRUE;
     170    if (fields[i].magOption == MAG_OPTION_GAL_MAJ         ) return TRUE;
     171    if (fields[i].magOption == MAG_OPTION_GAL_MAJ_ERR     ) return TRUE;
     172    if (fields[i].magOption == MAG_OPTION_GAL_MIN         ) return TRUE;
     173    if (fields[i].magOption == MAG_OPTION_GAL_MIN_ERR     ) return TRUE;
     174    if (fields[i].magOption == MAG_OPTION_GAL_THETA       ) return TRUE;
     175    if (fields[i].magOption == MAG_OPTION_GAL_THETA_ERR   ) return TRUE;
     176    if (fields[i].magOption == MAG_OPTION_GAL_INDEX       ) return TRUE;
     177    if (fields[i].magOption == MAG_OPTION_GAL_CHISQ       ) return TRUE;
     178    if (fields[i].magOption == MAG_OPTION_GAL_NPIX        ) return TRUE;
     179    if (fields[i].magOption == MAG_OPTION_GAL_TYPE        ) return TRUE;
    178180  }
    179181  return FALSE;
    … …  
    265267  if (!strcasecmp (string, "E1"))            return MAG_OPTION_E1;
    266268  if (!strcasecmp (string, "E2"))            return MAG_OPTION_E2;
     269
     270  if (!strcasecmp (string, "GAL_MAG"    ))   return MAG_OPTION_GAL_MAG       ;
     271  if (!strcasecmp (string, "GAL_MAG_ERR"))   return MAG_OPTION_GAL_MAG_ERR   ;
     272  if (!strcasecmp (string, "GAL_MAJ"    ))   return MAG_OPTION_GAL_MAJ       ;
     273  if (!strcasecmp (string, "GAL_MAJ_ERR"))   return MAG_OPTION_GAL_MAJ_ERR   ;
     274  if (!strcasecmp (string, "GAL_MIN"    ))   return MAG_OPTION_GAL_MIN       ;
     275  if (!strcasecmp (string, "GAL_MIN_ERR"))   return MAG_OPTION_GAL_MIN_ERR   ;
     276  if (!strcasecmp (string, "GAL_THETA"  ))   return MAG_OPTION_GAL_THETA     ;
     277  if (!strcasecmp (string, "GAL_THETA_ERR")) return MAG_OPTION_GAL_THETA_ERR ;
     278  if (!strcasecmp (string, "GAL_INDEX"  ))   return MAG_OPTION_GAL_INDEX     ;
     279  if (!strcasecmp (string, "GAL_CHISQ"  ))   return MAG_OPTION_GAL_CHISQ     ;
     280  if (!strcasecmp (string, "GAL_NPIX"   ))   return MAG_OPTION_GAL_NPIX      ;
     281  if (!strcasecmp (string, "GAL_TYPE"   ))   return MAG_OPTION_GAL_TYPE      ;
    267282
    268283  return MAG_OPTION_NONE;
    … …  
    709724  if (!strcasecmp (fieldName, "REF_COLOR_RED"))  ESCAPE (AVE_REF_COLOR_RED,  OPIHI_FLT);
    710725
    711   if (!strcasecmp (fieldName, "GAL_MAG"    ))      ESCAPE (AVE_GAL_MAG       ,   OPIHI_FLT);
    712   if (!strcasecmp (fieldName, "GAL_MAG_ERR"))      ESCAPE (AVE_GAL_MAG_ERR   ,   OPIHI_FLT);
    713   if (!strcasecmp (fieldName, "GAL_MAJ"    ))      ESCAPE (AVE_GAL_MAJ       ,   OPIHI_FLT);
    714   if (!strcasecmp (fieldName, "GAL_MAJ_ERR"))      ESCAPE (AVE_GAL_MAJ_ERR   ,   OPIHI_FLT);
    715   if (!strcasecmp (fieldName, "GAL_MIN"    ))      ESCAPE (AVE_GAL_MIN       ,   OPIHI_FLT);
    716   if (!strcasecmp (fieldName, "GAL_MIN_ERR"))      ESCAPE (AVE_GAL_MIN_ERR   ,   OPIHI_FLT);
    717   if (!strcasecmp (fieldName, "GAL_THETA"  ))      ESCAPE (AVE_GAL_THETA     ,   OPIHI_FLT);
    718   if (!strcasecmp (fieldName, "GAL_THETA_ERR"))    ESCAPE (AVE_GAL_THETA_ERR ,   OPIHI_FLT);
    719   if (!strcasecmp (fieldName, "GAL_INDEX"  ))      ESCAPE (AVE_GAL_INDEX     ,   OPIHI_FLT);
    720   if (!strcasecmp (fieldName, "GAL_CHISQ"  ))      ESCAPE (AVE_GAL_CHISQ     ,   OPIHI_FLT);
    721   if (!strcasecmp (fieldName, "GAL_NPIX"   ))      ESCAPE (AVE_GAL_NPIX      ,   OPIHI_FLT);
    722   if (!strcasecmp (fieldName, "GAL_TYPE"   ))      ESCAPE (AVE_GAL_TYPE      ,   OPIHI_FLT);
    723 
    724726  if (!strcasecmp (fieldName, "E_BV"))             ESCAPE (AVE_E_BV          ,   OPIHI_FLT);
    725727  if (!strcasecmp (fieldName, "E_BV_ERR"))         ESCAPE (AVE_E_BV_ERR      ,   OPIHI_FLT);

  • trunk/Ohana/src/libdvo/src/dvo_photcode_ops.c

    r38553 r38635  
    18541854LENSFIELD(E1);
    18551855LENSFIELD(E2);
     1856
     1857# define GALPHOT_FIELD(NAME, VALUE)                                             \
     1858  float GalphotValue_##NAME (PhotCode *code, GalPhot *galphot, int Ngalphot) { \
     1859    int n;                                                              \
     1860    if (code == NULL) return NAN;                                       \
     1861    for (n = 0; n < Ngalphot; n++) {                                    \
     1862      short equivCode = GetPhotcodeEquivCodebyCode (galphot[n].photcode); \
     1863      if (!equivCode) continue;                                         \
     1864      if (equivCode != code->code) continue;                            \
     1865      float value = galphot[n].VALUE;                                   \
     1866      return (value);                                                   \
     1867    }                                                                   \
     1868    return NAN;                                                         \
     1869  }
     1870
     1871GALPHOT_FIELD(GAL_MAG,         mag);           
     1872GALPHOT_FIELD(GAL_MAG_ERR,     magErr);   
     1873GALPHOT_FIELD(GAL_MAJ,         majorAxis);             
     1874GALPHOT_FIELD(GAL_MAJ_ERR,     majorAxisErr);
     1875GALPHOT_FIELD(GAL_MIN,         minorAxis);             
     1876GALPHOT_FIELD(GAL_MIN_ERR,     minorAxisErr);   
     1877GALPHOT_FIELD(GAL_THETA,       theta);     
     1878GALPHOT_FIELD(GAL_THETA_ERR,   thetaErr); 
     1879GALPHOT_FIELD(GAL_INDEX,       index);     
     1880GALPHOT_FIELD(GAL_CHISQ,       chisq);     
     1881GALPHOT_FIELD(GAL_NPIX,        Npix);       
     1882GALPHOT_FIELD(GAL_TYPE,        modelType);       
Note: See TracChangeset for help on using the changeset viewer.