More post-process tweaks

git-svn-id: https://svn.ligo.caltech.edu/svn/bayeswave/trunk@49 c56465c9-8126-4a4f-9d7d-ac845eff4865

postprocess/plot_waveform.py

@@ -78,19 +78,19 @@ for ifo in ifoList:
         print "   "+str(count)+" waveforms found for IFO "+str(int(ifo)+1)
         median_waveform = np.median(strain_master, axis=0)
 
-        # down_vec = np.zeros(median_waveform.size)
-        # up_vec = np.zeros(median_waveform.size)
-        # for sample in range(median_waveform.size):
-        #    vector = np.transpose(strain_master)[sample]
-        #    sort_vec = np.sort(vector)
-        #    up_quart = sort_vec[ int(0.84*sort_vec.size) ]
-        #    down_quart = sort_vec[ int(0.16*sort_vec.size) ]
-        #    down_vec[sample] = down_quart
-        #    up_vec[sample] = up_quart
+        down_vec = np.zeros(median_waveform.size)
+        up_vec = np.zeros(median_waveform.size)
+        for sample in range(median_waveform.size):
+            vector = np.transpose(strain_master)[sample]
+            sort_vec = np.sort(vector)
+            up_quart = sort_vec[ int(0.84*sort_vec.size) ]
+            down_quart = sort_vec[ int(0.16*sort_vec.size) ]
+            down_vec[sample] = down_quart
+            up_vec[sample] = up_quart
 
         # -- Get axis info
-        # wave = up_vec
-        wave = median_waveform
+        wave = up_vec
+        # wave = median_waveform
         wave_max = np.maximum(wave.max(), d_s.max())
         sig_times = time[wave > 0.05*wave.max()]
         axisList.append([sig_times.min(), sig_times.max(), -wave_max*1.1, wave_max*1.1])