Changeset 28589
- Timestamp:
- Jul 1, 2010, 4:05:29 PM (16 years ago)
- File:
-
- 1 edited
-
trunk/ippToPsps/src/ippToPspsBatchDetection.c (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ippToPsps/src/ippToPspsBatchDetection.c
r28579 r28589 83 83 84 84 // stuff to keep from psf.hdr header 85 int sourceId = -1; 86 int imageId = -1; 87 float fwhmMaj; 88 float fwhmMin; 89 float momentMaj; 90 float momentMin; 91 float psfFwhm; 92 float momentFwhm; 85 int sourceId = -1, imageId = -1; 86 float fwhmMaj, fwhmMin, momentMaj, momentMin, psfFwhm, momentFwhm; 93 87 long pspsImageId = -1; 94 88 … … 96 90 SkyList *skylist = NULL; 97 91 dvoDetection *dvoDetections = NULL; 98 int maxDvoDetId = -1; 99 int numDvoDetections = -1; 92 int maxDvoDetId = -1, numDvoDetections = -1; 100 93 Image *image = NULL; 101 94 … … 108 101 strftime(timeStr, sizeof(timeStr), "%Y-%m-%d", ts); 109 102 110 uint32_t totalDetections = 0; 111 uint32_t s,d; 112 uint32_t invalidDvoRows; 113 uint32_t smfJumps; 114 uint32_t unmatched; 103 uint32_t s,d, invalidDvoRows, smfJumps, unmatched, totalDetections = 0; 115 104 116 105 long longnull = -999; … … 118 107 int anynull = 0; 119 108 120 char ccdNumber[3] ;121 char extensionName[15]; 122 109 char ccdNumber[3], extensionName[15]; 110 111 // for storing FITS column data 123 112 long ippIDet[MAXDETECT]; 124 float instMag[MAXDETECT]; 125 float instMagErr[MAXDETECT]; 126 float peakMag[MAXDETECT]; 113 float instMag[MAXDETECT], instMagErr[MAXDETECT], peakMag[MAXDETECT]; 127 114 uint64_t flags[MAXDETECT]; 128 long objID[MAXDETECT]; 129 long detectID[MAXDETECT]; 130 long ippObjID[MAXDETECT]; 131 long ippDetectID[MAXDETECT]; 132 long imageID[MAXDETECT]; 115 long objID[MAXDETECT], detectID[MAXDETECT], ippObjID[MAXDETECT], ippDetectID[MAXDETECT], imageID[MAXDETECT]; 133 116 double obsTimes[MAXDETECT]; 134 float instFlux[MAXDETECT]; 135 float instFluxErr[MAXDETECT]; 136 float peakFlux[MAXDETECT]; 137 int8_t filterIDs[MAXDETECT]; 138 int8_t surveyIDs[MAXDETECT]; 117 float instFlux[MAXDETECT], instFluxErr[MAXDETECT], peakFlux[MAXDETECT]; 118 int8_t filterIDs[MAXDETECT], surveyIDs[MAXDETECT]; 139 119 140 120 char** assocDate = (char**)calloc(MAXDETECT, sizeof(char**)); … … 149 129 } 150 130 151 long maxObjID = LONG_MIN; 152 long minObjID = LONG_MAX; 131 long maxObjID = LONG_MIN, minObjID = LONG_MAX; 153 132 uint64_t imageFlags; 154 133 short nOta = 0; 155 134 long i; 156 135 bool isDuplicate; 157 uint32_t numOfDuplicates; 158 uint32_t numInvalidFlux; 159 long numDetectionsOut; 160 long totalDetectionsOut = 0; 161 long numPhotoRef, totalNumPhotoRef = 0; 136 uint32_t numOfDuplicates, numInvalidFlux; 137 long numDetectionsOut, totalDetectionsOut = 0, numPhotoRef, totalNumPhotoRef = 0; 162 138 long removeList[MAXDETECT]; 163 139 long thisObjId; … … 201 177 status=0; fits_read_key(fitsIn, TLONG, "SOURCEID", &sourceId, NULL, &status); 202 178 status=0; fits_read_key(fitsIn, TLONG, "NASTRO", &numPhotoRef, NULL, &status); 203 totalNumPhotoRef += numPhotoRef; 204 psLogMsg("ippToPsps", PS_LOG_INFO, "FFF %ld %ld\n", numPhotoRef, totalNumPhotoRef); 179 totalNumPhotoRef += numPhotoRef; // total up for storing in FrameMeta 180 205 181 // access DVO database 206 182 skylist = dvoSkyListByExternID(this->dvoConfig, sourceId, imageId, &image); 207 183 if (skylist == NULL) { 208 184 psError(PS_ERR_IO, false, 209 "DVO: can't find SkyList for sourceId='%d' imageId='%d' (CCD = XY%s): skipping \n",185 "DVO: can't find SkyList for sourceId='%d' imageId='%d' (CCD = XY%s): skipping this chip\n", 210 186 sourceId, imageId, ccdNumber); 211 187 continue;
Note:
See TracChangeset
for help on using the changeset viewer.
