Changeset 30596

Timestamp:

Feb 13, 2011, 10:48:38 AM (15 years ago)

Author:

eugene

Message:

adding option to enhance the noise (initially for BigBoss tests, from Daniel Farrow)

Location:

tags/ipp-20101215/ppImage/src

Files:

• Makefile.am (modified) (2 diffs)
• ppImage.h (modified) (3 diffs)
• ppImageAddNoise.c (added)
• ppImageLoop.c (modified) (2 diffs)
• ppImageOptions.c (modified) (2 diffs)
• ppImageRandomGaussian.c (added)

tags/ipp-20101215/ppImage/src/Makefile.am

@@ -2 +2 @@
 
 noinst_HEADERS = \
-        ppImage.h 
+        ppImage.h
 
 if HAVE_SVNVERSION
@@ -61 +61 @@
         ppImageFileCheck.c \
         ppImageVersion.c \
-        ppImageMemory.c
+        ppImageMemory.c \
+        ppImageAddNoise.c \
+        ppImageRandomGaussian.c
 
 CLEANFILES = *~

tags/ipp-20101215/ppImage/src/ppImage.h

@@ -53 +53 @@
     bool doCrosstalkMeasure;            // measure crosstalk signal
     bool doCrosstalkCorrect;            // apply crosstalk correction
+    bool addNoise;                      // Add noise to degrade MD image to 3pi
 
     // output files requested
@@ -155 +156 @@
 bool ppImageDetrendBias(pmReadout *inputReadout, pmReadout *bias, pmReadout *dark, ppImageOptions *options);
 
-//bool ppImageDetrendNonLinear(pmReadout *input, ppImageOptions *options);
 bool ppImageDetrendNonLinear(pmReadout *input, pmFPAview *linearity, pmConfig *config);
 bool ppImageDetrendNonLinearLookup(pmReadout *input, psMetadataItem *dataItem);
@@ -291 +291 @@
 void ppImageMemoryDump(const char *description);
 
+
+//Functions needed to degrade MD exposures to 3pi exposures
+
+bool ppImageAddNoise(pmConfig *config, ppImageOptions *options, pmFPAview *view, pmFPA *fpa) ;
+double ppImageRandomGaussian (const psRandom *rnd, double mean, double sigma);
+double ppImageRandomGaussianNorm (const psRandom *rnd);
+void ppImageRandomGaussianFree(void);
+
+
 #endif

tags/ipp-20101215/ppImage/src/ppImageLoop.c

@@ -110 +110 @@
                     ESCAPE("Unable to set bit masks");
                 }
+                // optionally degrade a MD image to 3pi exposure times
+                if (!ppImageAddNoise(config, options, view, input->fpa)){
+                    ESCAPE("Unable to degrade MD image to 3pi");
+                }
 
                 // perform the detrend analysis
@@ -152 +156 @@
             ESCAPE("Unable to free detrend images");
         }
-    
+
         // Apply the fringe correction
         if (options->doFringe) {

tags/ipp-20101215/ppImage/src/ppImageOptions.c

@@ -42 +42 @@
     options->applyParity     = false;   // Apply Cell parities
     options->doMaskStats     = false;   // Calculate mask fractions
-    
+    options->addNoise        = false;  //Degrade an MD image to a 3pi image
+
     // output files requested
     options->BaseFITS        = false;   // create output image
@@ -239 +240 @@
     options->doPatternCell = psMetadataLookupBool(NULL, recipe, "PATTERN.CELL");
     options->doMaskStats = psMetadataLookupBool(NULL, recipe, "MASK.STATS");
+    options->addNoise = psMetadataLookupBool(NULL, recipe, "ADDNOISE");
 
     options->doStats = false;