Changeset 25069 for branches/simtest_nebulous_branches/ppSim/src/ppSimSequenceDark.c

Timestamp:

Aug 13, 2009, 4:08:06 PM (17 years ago)

Author:

giebink

Message:

Support for simtests using nebulous for SIMTEST and SIMMOSAIC cameras.

File:

• branches/simtest_nebulous_branches/ppSim/src/ppSimSequenceDark.c (modified) (3 diffs)

branches/simtest_nebulous_branches/ppSim/src/ppSimSequenceDark.c

@@ -1 +1 @@
 # include "ppSimSequence.h"
 
-bool ppSimSequenceDark (FILE *simfile, FILE *inject, psMetadata *sequence, int nSeq, psRandom *rng, const char *path, const char *basename, const char *ppSimCommand, const char *injectCommand) {
+bool ppSimSequenceDark (FILE *simfile, FILE *inject, psMetadata *sequence, int nSeq, psRandom *rng, const char *path, const char *basename, const char *ppSimCommand, const char *injectCommand, const char *camera) {
 
     bool status, setRate;
@@ -9 +9 @@
     min = psMetadataLookupF32 (&status, sequence, "DARK.MIN");
     if (status) {
-        max = psMetadataLookupF32 (&status, sequence, "DARK.MAX");
-        setRate = true;
+                max = psMetadataLookupF32 (&status, sequence, "DARK.MAX");
+                setRate = true;
     }
 
@@ -22 +22 @@
     for (int i = 0; i < nImages->n; i++) {
             
-        float exptime = exptimes->data.F32[i];
-        float n = nImages->data.S32[i];
+                float exptime = exptimes->data.F32[i];
+                float n = nImages->data.S32[i];
 
-        for (int j = 0; j < n; j++) {
+            for (int j = 0; j < n; j++) {
 
-            // XXX need to add output filename
-            psString filename = NULL;
-            if (path) {
-                psStringAppend (&filename, "%s/%s.%03d.%03d", path, basename, nSeq, nImage);
-            } else {
-                psStringAppend (&filename, "%s.%03d.%03d", basename, nSeq, nImage);
-            }
+                // XXX need to add output filename
+                psString filename = NULL;
+                if (path) {
+                        psStringAppend (&filename, "%s/%s.%03d.%03d", path, basename, nSeq, nImage);
+                } else {
+                        psStringAppend (&filename, "%s.%03d.%03d", basename, nSeq, nImage);
+                }
 
-            // define the ppSim command
-            psString command = NULL;
+                // define the ppSim command
+                psString command = NULL;
 
-            psStringAppend (&command, "%s -type DARK", ppSimCommand);
-      
-            if (setRate) {
-                double frnd = psRandomUniform(rng);
-                float rate = min + (max - min)*frnd;
-                psStringAppend (&command, " -darkrate %f", rate);
-            }
+                psStringAppend (&command, "%s -type DARK", ppSimCommand);
+          
+                if (setRate) {
+                                double frnd = psRandomUniform(rng);
+                                float rate = min + (max - min)*frnd;
+                                psStringAppend (&command, " -darkrate %f", rate);
+                }
 
-            psStringAppend (&command, " -exptime %f", exptime);
+                psStringAppend (&command, " -exptime %f", exptime);
 
-            psStringAppend (&command, " %s", filename);
+                psStringAppend (&command, " %s", filename);
 
-            fprintf (simfile, "%s\n", command);
-            psFree (command);
-                            
-            // define the inject command
-            // path should be dirname/filename
-            command = psStringCopy (injectCommand);
-            psStringAppend (&command, " %s*.fits",    filename);
-            fprintf (inject, "%s\n", command);
-            psFree (command);
-            psFree (filename);
+                fprintf (simfile, "%s\n", command);
+                psFree (command);
+                                
+                // define the inject command
+                    // path should be dirname/filename
+                    // for SIMTEST camera nebulous insertion append .fits
+                    // SIMMOSAIC nebulous insertions are appended with ChipXX.fits
+                    // non-nebulous insertions use *.fits
+                    const char *neb = "neb";
 
-            nImage ++;
-        }
+                    if (strncmp (path, neb, 3) == 0) {
+                                
+                if (strcmp (camera, "SIMTEST") == 0) {
+                    command = psStringCopy (injectCommand); 
+                    psStringAppend (&command, " %s.fits", filename);
+                    fprintf (inject, "%s\n", command);
+                                        
+                } else if (strcmp (camera, "SIMMOSAIC") == 0) {
+                    for (int i = 0; i < 2; i++) {
+                        for (int j = 0; j < 2; j++) {
+                            command = psStringCopy (injectCommand);
+                                        psStringAppend (&command, " %s.Chip%d%d.fits", filename, i, j);
+                                        fprintf (inject, "%s\n", command);
+                                            command = NULL;
+                                            }
+                    }
+                                            
+                } else {
+                    psLogMsg ("ppSimSequence", PS_LOG_WARN, "unsupported camera type for nebulous insertion");
+                    exit (1);
+                }
+                
+            } else {
+                command = psStringCopy (injectCommand);
+                psStringAppend (&command, " %s*.fits",    filename);
+                fprintf (inject, "%s\n", command);
+            }
+
+            psFree (command);
+            psFree (filename);
+            nImage ++;
+        }
     }
     return true;