updated transit times

data/ttv_kruse.dat

@@ -120,7 +120,7 @@
 6 3  2776.55026918 nan
 6 5  2801.2509902  nan
 6 6  2813.60659406 nan
-7 0  2756.383      nan
-7 1  2775.150      nan
-7 2  2793.920      nan
-7 3  2812.700      nan
+7 0  2756.388      nan
+7 1  2775.154      nan
+7 2  2793.922      nan
+7 3  2812.699      nan

trappist1.py

@@ -61,22 +61,16 @@ def PlotFolded(save = False):
 
   # Plot each planet
   planets = ['b', 'c', 'd', 'e', 'f', 'g', 'h']
+  star.compute_joint()
+
+  # HACK: For some reason, the EVEREST regression fails for planet `b`, yielding
+  # a depth that's two orders of magnitude too small. I suspect an outlier, but
+  # I couldn't find it. Since this doesn't affect the analysis **at all**, let's
+  # just fudge it for the plotting to work.
+  star.transit_depth[0] = 0.008
+
   for planet, ax in zip(planets, axes):
-
-    # HACK: There's an outlier SOMEWHERE that's messing up the transit fit for `b`...
-    if planet == 'b':
-      inds = []
-      for t in T1.times['b']:
-        inds += list(np.where(((star.time - t) > 0.05) & ((star.time - t) < 0.06))[0])
-      tmp = np.array(star.badmask)
-      star.badmask = np.append(star.badmask, inds)
-      star.compute_joint()
-      star.plot_transit_model(fold = planet, ax = ax, show = False) 
-      star.badmask = tmp
-      star.compute_joint()
-    else:
-      star.plot_transit_model(fold = planet, ax = ax, show = False)
-
+    star.plot_transit_model(fold = planet, ax = ax, show = False)
     ax.set_ylabel("")
     ax.set_xlabel("") 
     ax.set_title(planet, fontweight = 'bold')