Added experimental capability to search through TICA FFIs

README.rst

@@ -175,6 +175,7 @@ Please include a self-contained example that fully demonstrates your problem or
 Changelog:
 ==========
 
+  - Added ability to use "TICA" FFIs. This is experimental and might be buggy.
   - Patch removes the un-needed `fitsio` dependency
   - Initial v1.0.0 release of `tesscube`.
 

pyproject.toml

@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "tesscube"
-version = "1.0.5"
+version = "1.1.0"
 description = ""
 authors = ["TESS Science Support Center <tesshelp@bigbang.gsfc.nasa.gov>", "Christina Hedges <christina.l.hedges@nasa.gov>"]
 license = "MIT"

src/tesscube/__init__.py

@@ -1,10 +1,21 @@
-__version__ = "1.0.5"
 # Standard library
 import os  # noqa
 import tempfile
 
 PACKAGEDIR = os.path.abspath(os.path.dirname(__file__))
 
+from importlib.metadata import version, PackageNotFoundError  # noqa
+
+
+def get_version():
+    try:
+        return version("tesscube")
+    except PackageNotFoundError:
+        return "unknown"
+
+
+__version__ = get_version()
+
 HDR_SIZE = 2880  # bytes
 BYTES_PER_PIX = 4  # float32
 MAX_CONCURRENT_DOWNLOADS = 10

src/tesscube/cube.py

@@ -7,7 +7,7 @@
 from astropy.time import Time
 from astropy.wcs.utils import fit_wcs_from_points
 
-from . import BYTES_PER_PIX, DATA_OFFSET, log
+from . import BYTES_PER_PIX, DATA_OFFSET, log, HDR_SIZE
 from .fits import (
     get_header_dict,
     get_output_first_extention_header,
@@ -38,9 +38,17 @@ class TESSCube(QueryMixin, WCSMixin):
         The CCD number (1-4).
     """
 
-    def __init__(self, sector: int, camera: int, ccd: int):
+    def __init__(self, sector: int, camera: int, ccd: int, tica=False):
         self.sector, self.camera, self.ccd = sector, camera, ccd
-        self.object_key = f"tess/public/mast/tess-s{sector:04}-{camera}-{ccd}-cube.fits"
+        self.tica = tica
+        if self.tica:
+            self.object_key = (
+                f"tess/public/mast/tica/tica-s{sector:04}-{camera}-{ccd}-cube.fits"
+            )
+        else:
+            self.object_key = (
+                f"tess/public/mast/tess-s{sector:04}-{camera}-{ccd}-cube.fits"
+            )
         (
             self.nsets,
             self.nframes,
@@ -115,6 +123,8 @@ def last_hdu(self):
             DATA_OFFSET
             + (self.ncolumns * self.nframes * self.nsets * self.nrows) * BYTES_PER_PIX
         )
+        if self.tica:
+            end += HDR_SIZE * 4
         return get_last_hdu(object_key=self.object_key, end=end)
 
     @cached_property
@@ -124,15 +134,18 @@ def ffi_names(self):
 
     @cached_property
     def tstart(self):
-        return self.last_hdu.data["TSTART"]
+        return self.last_hdu.data["TSTART" if not self.tica else "STARTTJD"]
 
     @cached_property
     def tstop(self):
-        return self.last_hdu.data["TSTOP"]
+        return self.last_hdu.data["TSTOP" if not self.tica else "ENDTJD"]
 
     @cached_property
     def telapse(self):
-        return self.last_hdu.data["TELAPSE"]
+        if not self.tica:
+            return self.last_hdu.data["TELAPSE"]
+        else:
+            return self.tstop - self.tstart
 
     @lru_cache(maxsize=4)
     def get_ffi(
@@ -474,12 +487,18 @@ def time(self):
     @cached_property
     def timecorr(self):
         """Barycentric time correction for the center of the FFI."""
-        return self.last_hdu.data["BARYCORR"]
+        if not self.tica:
+            return self.last_hdu.data["BARYCORR"]
+        else:
+            return self.time * 0
 
     @cached_property
     def quality(self):
         """SPOC provided quality flags for each cadence."""
-        return self.last_hdu.data["DQUALITY"]
+        if not self.tica:
+            return self.last_hdu.data["DQUALITY"]
+        else:
+            return self.last_hdu.data["QUAL_BIT"]
 
     @cached_property
     def cadence_number(self):
@@ -492,7 +511,10 @@ def cadence_number(self):
     @cached_property
     def exposure_time(self):
         """Exposure time in days"""
-        return self.last_hdu.data["EXPOSURE"][0]
+        if not self.tica:
+            return self.last_hdu.data["EXPOSURE"][0]
+        else:
+            return self.last_hdu.data["INT_TIME"] / 86400
 
     @property
     def shape(self):

src/tesscube/fits.py

@@ -59,23 +59,29 @@ def get_header_dict(cube):
     )
     header_dict["PROCVER"] = fits.Card("PROCVER", __version__, "software version")
     header_dict["TSTART"] = fits.Card(
-        "TSTART",
-        cube.last_hdu.data["TSTART"][0],
+        "TSTART" if not cube.tica else "STARTTJD",
+        cube.last_hdu.data["TSTART" if not cube.tica else "STARTTJD"][0],
         "observation start time in TJD of first FFI",
     )
     header_dict["TSTOP"] = fits.Card(
-        "TSTOP",
-        cube.last_hdu.data["TSTOP"][-1],
+        "TSTOP" if not cube.tica else "ENDTJD",
+        cube.last_hdu.data["TSTOP" if not cube.tica else "ENDTJD"][-1],
         "observation stop time in TJD of last FFI",
     )
     header_dict["DATE-OBS"] = fits.Card(
         "DATE-OBS",
-        Time(cube.last_hdu.data["TSTART"][0] + 2457000, format="jd").isot,
+        Time(
+            cube.last_hdu.data["TSTART" if not cube.tica else "STARTTJD"][0] + 2457000,
+            format="jd",
+        ).isot,
         "TSTART as UTC calendar date",
     )
     header_dict["DATE-END"] = fits.Card(
         "DATE-END",
-        Time(cube.last_hdu.data["TSTOP"][-1] + 2457000, format="jd").isot,
+        Time(
+            cube.last_hdu.data["TSTOP" if not cube.tica else "ENDTJD"][-1] + 2457000,
+            format="jd",
+        ).isot,
         "TSTOP as UTC calendar date",
     )
     return header_dict
@@ -118,8 +124,10 @@ def get_output_primary_hdu(cube):
         "TICVER",
         "CRMITEN",
         "CRBLKSZ",
-        "CRSPOC",
     ]
+    if not cube.tica:
+        for a in ["CRSPOC"]:
+            keys.append(a)
     header_dict = cube.header_dict
     return fits.PrimaryHDU(header=fits.Header(cards=[header_dict[key] for key in keys]))
 
@@ -135,45 +143,50 @@ def get_output_first_extention_header(cube):
         "RA_OBJ",
         "DEC_OBJ",
         "EQUINOX",
-        "EXPOSURE",
         "TIMEREF",
         "TASSIGN",
         "TIMESYS",
         "BJDREFI",
         "BJDREFF",
         "TIMEUNIT",
         "TELAPSE",
-        "LIVETIME",
         "TSTART",
         "TSTOP",
         "DATE-OBS",
         "DATE-END",
-        "DEADC",
-        "TIMEPIXR",
-        "TIERRELA",
-        "INT_TIME",
-        "READTIME",
-        "FRAMETIM",
-        "NUM_FRM",
-        "TIMEDEL",
-        "BACKAPP",
-        "DEADAPP",
-        "VIGNAPP",
-        "GAINA",
-        "GAINB",
-        "GAINC",
-        "GAIND",
-        "READNOIA",
-        "READNOIB",
-        "READNOIC",
-        "READNOID",
-        "NREADOUT",
-        "FXDOFF",
-        "MEANBLCA",
-        "MEANBLCB",
-        "MEANBLCC",
-        "MEANBLCD",
     ]
+    if not cube.tica:
+        for a in [
+            "EXPOSURE",
+            "LIVETIME",
+            "DEADC",
+            "TIMEPIXR",
+            "TIERRELA",
+            "INT_TIME",
+            "READTIME",
+            "FRAMETIM",
+            "NUM_FRM",
+            "TIMEDEL",
+            "BACKAPP",
+            "DEADAPP",
+            "VIGNAPP",
+            "GAINA",
+            "GAINB",
+            "GAINC",
+            "GAIND",
+            "READNOIA",
+            "READNOIB",
+            "READNOIC",
+            "READNOID",
+            "NREADOUT",
+            "FXDOFF",
+            "MEANBLCA",
+            "MEANBLCB",
+            "MEANBLCC",
+            "MEANBLCD",
+        ]:
+            keys.append(a)
+
     header_dict = cube.header_dict
     return fits.Header(cards=[header_dict[key] for key in keys])
 
@@ -371,7 +384,7 @@ def _fix_primary_hdu(hdu):
     hdr["TICID"] = (None, "unique tess target identifier")
 
     delete_kwds_wildcards = [
-        "NAXIS*" "SC_*",
+        "NAXIS*SC_*",
         "RMS*",
         "A_*",
         "AP_*",