IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 39251 for trunk


Ignore:
Timestamp:
Dec 11, 2015, 4:57:24 PM (11 years ago)
Author:
watersc1
Message:

Change to calculate magnitudes for average flux values in the ForcedMeanObject table

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/Ohana/src/dvopsps/src/insert_FWobjects_dvopsps_catalog.c

    r39092 r39251  
    148148                 "FLUX_PSF_WRP_ERR    FLOAT,  "         
    149149                 "FLUX_PSF_WRP_STDEV  FLOAT,  "         
    150 
     150                 "MAG_PSF_WRP         FLOAT,  "
     151                 "MAG_PSF_WRP_ERR     FLOAT,  "
     152                 
    151153                 "FLUX_KRON_WRP       FLOAT,  "       
    152154                 "FLUX_KRON_WRP_ERR   FLOAT,  "   
    153155                 "FLUX_KRON_WRP_STDEV FLOAT,  "   
    154 
     156                 "MAG_KRON_WRP        FLOAT,  "
     157                 "MAG_KRON_WRP_ERR    FLOAT,  "
     158                 
    155159                 "FLUX_AP_WRP         FLOAT,  "             
    156160                 "FLUX_AP_WRP_ERR     FLOAT,  "             
    157                  "FLUX_AP_WRP_STDEV   FLOAT,  "             
     161                 "FLUX_AP_WRP_STDEV   FLOAT,  "
     162                 "MAG_AP_WRP          FLOAT,  "
     163                 "MAG_AP_WRP_ERR      FLOAT,  "
    158164
    159165                 "FLAGS               INT "     
    … …  
    188194                 "FLUX_STD_AP_R5     FLOAT, "
    189195                 "FLUX_FIL_AP_R5     FLOAT, "
     196                 "MAG_AP_R5          FLOAT, "
     197                 "MAG_ERR_AP_R5      FLOAT, "
    190198
    191199                 "FLUX_AP_R6         FLOAT, "
    … …  
    193201                 "FLUX_STD_AP_R6     FLOAT, "
    194202                 "FLUX_FIL_AP_R6     FLOAT, "
     203                 "MAG_AP_R6          FLOAT, "
     204                 "MAG_ERR_AP_R6      FLOAT, "
    195205
    196206                 "FLUX_AP_R7         FLOAT, "
    … …  
    198208                 "FLUX_STD_AP_R7     FLOAT, "
    199209                 "FLUX_FIL_AP_R7     FLOAT, "
     210                 "MAG_AP_R7          FLOAT, "
     211                 "MAG_ERR_AP_R7      FLOAT, "
    200212
    201213                 "X11_SM_OBJ         FLOAT, "
    … …  
    262274                 "FLUX_PSF_WRP, "             
    263275                 "FLUX_PSF_WRP_ERR, "         
    264                  "FLUX_PSF_WRP_STDEV, "         
     276                 "FLUX_PSF_WRP_STDEV, "
     277                 "MAG_PSF_WRP,  "
     278                 "MAG_PSF_WRP_ERR, "
    265279
    266280                 "FLUX_KRON_WRP, "       
    267281                 "FLUX_KRON_WRP_ERR, "   
    268                  "FLUX_KRON_WRP_STDEV, "   
     282                 "FLUX_KRON_WRP_STDEV, "
     283                 "MAG_KRON_WRP, "
     284                 "MAG_KRON_WRP_ERR, "
    269285
    270286                 "FLUX_AP_WRP, "             
    271287                 "FLUX_AP_WRP_ERR, "             
    272                  "FLUX_AP_WRP_STDEV, "             
     288                 "FLUX_AP_WRP_STDEV, "
     289                 "MAG_AP_WRP, "
     290                 "MAG_AP_WRP_ERR, "
    273291
    274292                 "FLAGS"     
    … …  
    287305                 "FLUX_STD_AP_R5, "
    288306                 "FLUX_FIL_AP_R5, "
     307                 "MAG_AP_R5, "
     308                 "MAG_ERR_AP_R5, "
    289309
    290310                 "FLUX_AP_R6, "
    … …  
    292312                 "FLUX_STD_AP_R6, "
    293313                 "FLUX_FIL_AP_R6, "
     314                 "MAG_AP_R6, "
     315                 "MAG_ERR_AP_R6, "
    294316
    295317                 "FLUX_AP_R7, "
    … …  
    297319                 "FLUX_STD_AP_R7, "
    298320                 "FLUX_FIL_AP_R7, "
     321                 "MAG_AP_R7, "
     322                 "MAG_ERR_AP_R7, "
    299323
    300324                 "X11_SM_OBJ, "
    … …  
    348372    float meanPSFFluxErr = NAN;
    349373    float meanPSFFluxStd = NAN;
     374    float meanPSFMag     = NAN;
     375    float meanPSFMagErr  = NAN;
    350376    if (isfinite(secfilt->dFpsfWrp) && isfinite(secfilt->FpsfWrp)) {
    351377      meanPSFFlux    = secfilt-> FpsfWrp;
    352378      meanPSFFluxErr = secfilt->dFpsfWrp;
    353379      meanPSFFluxStd = secfilt->sFpsfWrp;
     380      if (meanPSFFlux > 0.0) {
     381        meanPSFMag    = -2.5 * log10(meanPSFFlux / 3631.0);
     382        meanPSFMagErr = (2.5 * meanPSFFluxErr) / (meanPSFFlux * log(10));
     383      }
    354384    }
    355385
    … …  
    357387    float meanKronFluxErr = NAN;
    358388    float meanKronFluxStd = NAN;
     389    float meanKronMag     = NAN;
     390    float meanKronMagErr  = NAN;
    359391    if (isfinite(secfilt->dFkronWrp) && isfinite(secfilt->FkronWrp)) {
    360392      meanKronFlux    = secfilt-> FkronWrp;
    361393      meanKronFluxErr = secfilt->dFkronWrp;
    362394      meanKronFluxStd = secfilt->sFkronWrp;
     395      if (meanKronFlux > 0.0) {
     396        meanKronMag    = -2.5 * log10(meanKronFlux / 3631.0);
     397        meanKronMagErr = (2.5 * meanKronFluxErr) / (meanKronFlux * log(10));
     398      }     
    363399    }
    364400
    … …  
    366402    float meanApFluxErr = NAN;
    367403    float meanApFluxStd = NAN;
     404    float meanApMag     = NAN;
     405    float meanApMagErr  = NAN;
    368406    if (isfinite(secfilt->dFapWrp) && isfinite(secfilt->FapWrp)) {
    369407      meanApFlux    = secfilt-> FapWrp;
    370408      meanApFluxErr = secfilt->dFapWrp;
    371409      meanApFluxStd = secfilt->sFapWrp;
     410      if (meanApFlux > 0.0) {
     411        meanApMag    = -2.5 * log10(meanApFlux / 3631.0);
     412        meanApMagErr = (2.5 * meanApFluxErr) / (meanApFlux * log(10));
     413      }
    372414    }
    373415
    … …  
    382424    PRINT_FLOAT(sec_buffer, meanPSFFluxErr,   "%.6e, "); // uflux precision
    383425    PRINT_FLOAT(sec_buffer, meanPSFFluxStd,   "%.6e, "); // uflux precision
    384 
     426    PRINT_FLOAT(sec_buffer, meanPSFMag,       "%.6f, ");
     427    PRINT_FLOAT(sec_buffer, meanPSFMagErr,    "%.6f, ");   
     428   
    385429    PRINT_FLOAT(sec_buffer, meanKronFlux,     "%.6e, "); // uflux precision
    386430    PRINT_FLOAT(sec_buffer, meanKronFluxErr,  "%.6e, "); // uflux precision
    387431    PRINT_FLOAT(sec_buffer, meanKronFluxStd,  "%.6e, "); // uflux precision
     432    PRINT_FLOAT(sec_buffer, meanKronMag,      "%.6f, ");
     433    PRINT_FLOAT(sec_buffer, meanKronMagErr,   "%.6f, ");
    388434
    389435    PRINT_FLOAT(sec_buffer, meanApFlux,     "%.6e, "); // uflux precision
    390436    PRINT_FLOAT(sec_buffer, meanApFluxErr,  "%.6e, "); // uflux precision
    391437    PRINT_FLOAT(sec_buffer, meanApFluxStd,  "%.6e, "); // uflux precision
     438    PRINT_FLOAT(sec_buffer, meanApMag,      "%.6f, ");
     439    PRINT_FLOAT(sec_buffer, meanApMagErr,   "%.6f, ");
    392440
    393441    PrintIOBuffer (sec_buffer, "%u ", secfilt->flags);       
    … …  
    402450    PrintIOBuffer (cpy_buffer, "%hd, ", lensobj->photcode);         
    403451
     452    float magApR5    = NAN;
     453    float magErrApR5 = NAN;
     454    float magApR6    = NAN;
     455    float magErrApR6 = NAN;
     456    float magApR7    = NAN;
     457    float magErrApR7 = NAN;
     458
     459    if (lensobj->F_ApR5 > 0.0) {
     460      magApR5    = -2.5 * log10(lensobj->F_ApR5 / 3631.0);
     461      magErrApR5 = (2.5 * lensobj->dF_ApR5) / (lensobj->F_ApR5 * log(10));
     462    }
     463    if (lensobj->F_ApR6 > 0.0) {
     464      magApR6    = -2.5 * log10(lensobj->F_ApR6 / 3631.0);
     465      magErrApR6 = (2.5 * lensobj->dF_ApR6) / (lensobj->F_ApR6 * log(10));
     466    }
     467    if (lensobj->F_ApR7 > 0.0) {
     468      magApR7    = -2.5 * log10(lensobj->F_ApR7 / 3631.0);
     469      magErrApR7 = (2.5 * lensobj->dF_ApR7) / (lensobj->F_ApR7 * log(10));
     470    }
     471   
    404472    PRINT_FLOAT(cpy_buffer, lensobj-> F_ApR5,     "%.6e, ");
    405473    PRINT_FLOAT(cpy_buffer, lensobj->dF_ApR5,     "%.6e, ");
    406474    PRINT_FLOAT(cpy_buffer, lensobj->sF_ApR5,     "%.6e, ");
    407475    PRINT_FLOAT(cpy_buffer, lensobj->fF_ApR5,     "%.6e, ");
     476    PRINT_FLOAT(cpy_buffer, magApR5,              "%.6f, ");
     477    PRINT_FLOAT(cpy_buffer, magErrApR5,           "%.6f, ");
    408478    PRINT_FLOAT(cpy_buffer, lensobj-> F_ApR6,     "%.6e, ");
    409479    PRINT_FLOAT(cpy_buffer, lensobj->dF_ApR6,     "%.6e, ");
    410480    PRINT_FLOAT(cpy_buffer, lensobj->sF_ApR6,     "%.6e, ");
    411481    PRINT_FLOAT(cpy_buffer, lensobj->fF_ApR6,     "%.6e, ");
     482    PRINT_FLOAT(cpy_buffer, magApR6,              "%.6f, ");
     483    PRINT_FLOAT(cpy_buffer, magErrApR6,           "%.6f, ");
    412484    PRINT_FLOAT(cpy_buffer, lensobj-> F_ApR7,     "%.6e, ");
    413485    PRINT_FLOAT(cpy_buffer, lensobj->dF_ApR7,     "%.6e, ");
    414486    PRINT_FLOAT(cpy_buffer, lensobj->sF_ApR7,     "%.6e, ");
    415487    PRINT_FLOAT(cpy_buffer, lensobj->fF_ApR7,     "%.6e, ");
     488    PRINT_FLOAT(cpy_buffer, magApR7,              "%.6f, ");
     489    PRINT_FLOAT(cpy_buffer, magErrApR7,           "%.6f, ");
    416490    PRINT_FLOAT(cpy_buffer, lensobj->X11_sm_obj,  "%.6e, ");
    417491    PRINT_FLOAT(cpy_buffer, lensobj->X12_sm_obj,  "%.6e, ");
Note: See TracChangeset for help on using the changeset viewer.