Changeset 34531 for branches/czw_branch/20120906/ppStack/src/ppStackCombineFinal.c

Timestamp:

Oct 12, 2012, 4:12:45 PM (14 years ago)

Author:

watersc1

Message:

Finally was able to get a background model out of ppStack. Now to get a background restored model out as well.

File:

• branches/czw_branch/20120906/ppStack/src/ppStackCombineFinal.c (modified) (3 diffs)

branches/czw_branch/20120906/ppStack/src/ppStackCombineFinal.c

@@ -15 +15 @@
   
   psArray *inputs  = psArrayAlloc(options->num);
-  //  const psArray *cells = options->cells;
   for (int i = 0; i < options->num; i++) {
     ppStackFileActivationSingle(config, PPSTACK_FILES_BKG, true, i);
@@ -23 +22 @@
     inputs->data[i] = ro;
     pmFPAfileClose(file,view);
-    ppStackFileActivationSingle(config, PPSTACK_FILES_BKG, false, i);
   }
   if (!pmStackSimpleMedianCombine(bkgRO,inputs)) {
@@ -29 +27 @@
     return(false);
   }
-/*   for (int i = 0; i < options->num; i++) { */
-/*     //    ppStackFilesIterateUp(config); */
-/* /\*     pmFPAfile *file = pmFPAfileSelectSingle(config->files, "PPSTACK.INPUT.BKGMODEL", i); *\/ */
-/* /\*     pmFPAview *view = ppStackFilesIterateDown(config); *\/ */
-/*     pmFPAfile *file = pmFPAfileSelectSingle(config->files, "PPSTACK.INPUT.BKGMODEL", i); */
-    
-/*     ppStackFileActivationSingle(config, PPSTACK_FILES_BKG, false, i); */
-/*   }   */
+  if (!ppStackWriteImage("/tmp/test_forced.bkgmdl.fits",
+                         bkgRO->parent->parent->parent->hdu->header,
+                         bkgRO->image,
+                         config)) {
+    fprintf(stderr,"Failed to write image because fail.\n");
+  }
+  
+  for (int i = 0; i < options->num; i++) {
+    pmFPAfile *file = pmFPAfileSelectSingle(config->files, "PPSTACK.INPUT.BKGMODEL", i);
+    pmFPAview *view = ppStackFilesIterateDown(config);
+    bool success = pmFPAfileClose(file,view);
+    if (!success) {
+      psTrace("ppStack",5,"I failed at closing a file.\n");
+    }
+              
+    psFitsClose(file->fits);
+    file->fits = NULL;
+    file->header = NULL;
+    file->state = PM_FPA_STATE_CLOSED;
+    file->wrote_phu = false;
+    ppStackFileActivationSingle(config, PPSTACK_FILES_BKG, false, i);
+  }
   psFree(inputs);
-  //  psFree(stack->bkgFits);
   bkgRO->data_exists = true;
+  bkgRO->parent->data_exists = true;
+  bkgRO->parent->parent->data_exists = true;
   return(true);
 }