Fixing bug in mideclipse time of WASP-189

bmorris3 · Apr 29, 2020 · faf03c8 · faf03c8
1 parent dc98021
commit faf03c8
Show file tree

Hide file tree

Showing 7 changed files with 2,494 additions and 2,798 deletions.
diff --git a/example/WASP-189_build.ipynb b/example/WASP-189_build.ipynb
@@ -49,7 +49,6 @@
     "n_visits = 4\n",
     "\n",
     "for i in range(n_visits):\n",
-    "\n",
     "    cheops_orbit_min = 99.5  # minutes\n",
     "\n",
     "    all_times = np.linspace(p.t0 + p.per * (0.5 + i) - 0.1 * p.per, \n",
@@ -78,8 +77,8 @@
     "    centroid_x = location_x[0] + 0.2 * np.random.randn(n_points)\n",
     "    centroid_y = location_y[0] + 0.2 * np.random.randn(n_points)\n",
     "\n",
-    "    p.t_secondary = 0.5\n",
     "    p.fp = 80e-6\n",
+    "\n",
     "    model = TransitModel(p, bjd_time, \n",
     "                         supersample_factor=3,\n",
     "                         transittype='secondary',\n",
@@ -101,6 +100,7 @@
     "\n",
     "    fluxerr = np.std(flux) * np.ones(len(bjd_time))\n",
     "\n",
+    "\n",
     "    n_outliers = 50\n",
     "\n",
     "    flux[np.random.randint(0, flux.shape[0], size=n_outliers)] += 3.5 * flux.std() * np.random.randn(n_outliers)\n",
@@ -134,9 +134,16 @@
     "\n",
     "    for name, fmt, arr in formatter: \n",
     "        ra[name] = arr\n",
-    "\n",
+    "        \n",
     "    fits.writeto(f'../linea/data/example_wasp189_visit{i}.fits', ra, overwrite=True)"
    ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": []
   }
  ],
  "metadata": {

diff --git a/linea/data/example_wasp189_visit0.fits b/linea/data/example_wasp189_visit0.fits
diff --git a/linea/data/example_wasp189_visit1.fits b/linea/data/example_wasp189_visit1.fits
diff --git a/linea/data/example_wasp189_visit2.fits b/linea/data/example_wasp189_visit2.fits
diff --git a/linea/data/example_wasp189_visit3.fits b/linea/data/example_wasp189_visit3.fits
diff --git a/linea/data/planets.json b/linea/data/planets.json
@@ -24,7 +24,7 @@
         "per": 2.7240338,
         "rp": 0.06099180272790763,
         "t0": 2456706.4558,
-        "t_secondary": 0.5,
+        "t_secondary": 2456707.8178169,
         "u": [
             0,
             0

diff --git a/linea/tests/test_core_eclipse.py b/linea/tests/test_core_eclipse.py
@@ -52,7 +52,6 @@ def generate_recarrays_WASP189(depth_ppm=80, seed=42, n_outliers=50,
         centroid_x = location_x[0] + 0.2 * np.random.randn(n_points)
         centroid_y = location_y[0] + 0.2 * np.random.randn(n_points)
 
-        p.t_secondary = 0.5
         p.fp = depth_ppm * 1e-6
         model = TransitModel(p, bjd_time,
                              supersample_factor=3,