Changeset 38315

Timestamp:

May 24, 2015, 10:18:49 AM (11 years ago)

Author:

eugene

Message:

more relphot testing work

File:

• branches/eam_branches/ipp-20150419/Ohana/src/addstar/test/relphot.flatcorr.dvo (modified) (7 diffs)

branches/eam_branches/ipp-20150419/Ohana/src/addstar/test/relphot.flatcorr.dvo

@@ -26 +26 @@
 $DO_SYNTH = 0
 $DO_TYCHO = 0
+$NC_CLOUDS = 0.08
+$UC_CLOUDS = 0.02
 
 # set various global variables
@@ -84 +86 @@
   $CHIP_DX = 1000
   $CHIP_DY = 1000
-  $DEFECT_FRAC = 0.03
+  $DEFECT_FRAC = 0.05
   $OFFSET_FRAC_UC = 0.0
   $OFFSET_FRAC_NC = 0.5
@@ -92 +94 @@
   # the number of exposures is defined by filt_uc:n
   list filt_uc -split g g g g r r r r
-  $filt_uc:n = 2
+  $filt_uc:n = 4
 
   # tmpseq is used to generate vectors in this function
@@ -109 +111 @@
   for i 0 $filt_uc:n
     klam_uc[$i] = $klam_nominal:$filt_uc:$i 
-    zpt_uc[$i]  = $zpt_nominal:$filt_uc:$i + 2.5*log(exptime_uc[$i]) + klam_uc[$i]*(airmass_uc[$i] - 1.0) + 0.04*rnd(klam_uc[$i]) - 0.02
-  # zpt_uc[$i]  = $zpt_nominal:$filt_uc:$i + 2.5*log(exptime_uc[$i]) + klam_uc[$i]*(airmass_uc[$i] - 1.0)
+    zpt_uc[$i]  = $zpt_nominal:$filt_uc:$i + 2.5*log(exptime_uc[$i]) + klam_uc[$i]*(airmass_uc[$i] - 1.0) + $UC_CLOUDS*(rnd(klam_uc[$i]) - 0.5)
   end    
   set mjd_uc = 55000.00 + 0.01*(tmpseq % 3) + 10.0*int(tmpseq / 3)
@@ -116 +117 @@
   # mjd and zpt values for not-ubercal'ed images
   list filt_nc -split g g g g r r r r
-  $filt_nc:n = 2
+  $filt_nc:n = 4
 
   # tmpseq is used to generate vectors in this function
@@ -130 +131 @@
   for i 0 $filt_nc:n
     klam_nc[$i] = $klam_nominal:$filt_nc:$i 
-    zpt_nc[$i]  = $zpt_nominal:$filt_nc:$i + 2.5*log(exptime_nc[$i]) + klam_nc[$i]*(airmass_nc[$i] - 1.0) + 0.08*rnd(klam_nc[$i]) - 0.04
-  # zpt_nc[$i]  = $zpt_nominal:$filt_nc:$i + 2.5*log(exptime_nc[$i]) + klam_nc[$i]*(airmass_nc[$i] - 1.0)
+    zpt_nc[$i]  = $zpt_nominal:$filt_nc:$i + 2.5*log(exptime_nc[$i]) + klam_nc[$i]*(airmass_nc[$i] - 1.0) + $NC_CLOUDS*(rnd(klam_nc[$i]) - 0.5)
   end    
   set mjd_nc = 55000.10 + 0.01*(tmpseq % 3) + 10.0*int(tmpseq / 3)
@@ -521 +521 @@
   # create deviant photometry (1%)
   # I'm setting the flag value 0x02 (EXTMODEL) if the object is deviated.  This is a simple way to pass the info
+  # note that outlier rejection only occurs for Nmeas > 5
   set stars_mi = stars_mi + 0.2*badPhot
   set stars_fl = zero(stars_mi) + 0x00000080*badFlag + 0x00000100*poorFlag + 0x00000002*badPhot