astro-otter 0.1.0__py3-none-any.whl → 0.3.0__py3-none-any.whl

{astro_otter-0.1.0.dist-info → astro_otter-0.3.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
-Metadata-Version: 2.1
+Metadata-Version: 2.4
 Name: astro-otter
-Version: 0.1.0
+Version: 0.3.0
 Author-email: Noah Franz <nfranz@arizona.edu>
 License: MIT License
         
@@ -36,14 +36,15 @@ Classifier: Programming Language :: Python :: 3
 Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
 Classifier: Development Status :: 2 - Pre-Alpha
-Requires-Python: <3.12,>=3.10
+Requires-Python: >=3.9
 Description-Content-Type: text/markdown
 License-File: LICENSE
-Requires-Dist: numpy <2,>=1.20
-Requires-Dist: astropy >=6
+Requires-Dist: numpy<2,>=1.20
+Requires-Dist: astropy>=6
 Requires-Dist: pandas
 Requires-Dist: synphot
 Requires-Dist: typing-extensions
+Requires-Dist: pyarango
 Requires-Dist: matplotlib
 Requires-Dist: plotly
 Requires-Dist: astroquery
@@ -53,18 +54,21 @@ Requires-Dist: fundamentals
 Requires-Dist: astro-datalab
 Requires-Dist: sparclclient
 Requires-Dist: astro-ghost
-Requires-Dist: pyreadline3 ; platform_system == "Windows"
-Provides-Extra: dev
-Requires-Dist: ruff ; extra == 'dev'
-Requires-Dist: pre-commit ; extra == 'dev'
+Requires-Dist: pydantic
+Requires-Dist: pyreadline3; platform_system == "Windows"
 Provides-Extra: docs
-Requires-Dist: Sphinx >=3.0.0 ; extra == 'docs'
-Requires-Dist: myst-parser >=0.13 ; extra == 'docs'
-Requires-Dist: nbsphinx >=0.9.1 ; extra == 'docs'
-Requires-Dist: sphinx-book-theme >=0.0.33 ; extra == 'docs'
-Requires-Dist: sphinx-copybutton ; extra == 'docs'
-Requires-Dist: autodoc ; extra == 'docs'
-Requires-Dist: ipykernel ; extra == 'docs'
+Requires-Dist: Sphinx>=3.0.0; extra == "docs"
+Requires-Dist: myst-parser>=0.13; extra == "docs"
+Requires-Dist: nbsphinx>=0.9.1; extra == "docs"
+Requires-Dist: sphinx-book-theme>=0.0.33; extra == "docs"
+Requires-Dist: sphinx_copybutton; extra == "docs"
+Requires-Dist: autodoc; extra == "docs"
+Requires-Dist: ipykernel; extra == "docs"
+Requires-Dist: autodoc_pydantic; extra == "docs"
+Provides-Extra: dev
+Requires-Dist: ruff; extra == "dev"
+Requires-Dist: pre-commit; extra == "dev"
+Dynamic: license-file
 
 # OTTER API
 ### **O**pen mul**T**iwavelength **T**ransient **E**vent **R**epository
@@ -104,6 +108,39 @@ To install the OTTER API use
 python3 -m pip install astro-otter
 ```
 
+## Developer Instructions
+1. Set the `OTTER_ROOT` environment variable
+   ```
+   export OTTER_ROOT=/path/to/where/to/clone
+   ```
+2. Clone the relevant repos:
+   ```
+   git clone https://github.com/astro-otter/otter.git $OTTER_ROOT/otter
+   git clone https://github.com/astro-otter/otterdb.git $OTTER_ROOT/otterdb
+   ```
+3. Install the NASA ADS Python API by following the instructions at https://ads.readthedocs.io/en/latest/#getting-started
+4. Install otter, the API for this database. From
+   the root directory where you installed these repos:
+   ```
+   cd $OTTER_ROOT/otter
+   python -m pip install -e .
+   ```
+5. Process the data to build the local "database" (although it is really just a directory).
+   Then, you can build the "database" by running the
+   following commands:
+   ```
+   cd $OTTER_ROOT/otter/scripts/
+   python3 gen_summary_table.py --otterroot $OTTER_ROOT
+   ```
+6. Easily access the data using the Otter code! In python:
+  ```
+  import os
+  from otter import Otter
+  otter = Otter(os.path.join(os.environ['OTTER_ROOT'], 'otterdb', '.otter'))
+  res = otter.query(names='AT2018hyz')
+  print(res)
+  ```
+
 ## Installation from Source
 To install the OTTER API from the source code use
 ```

astro_otter-0.3.0.dist-info/RECORD

@@ -0,0 +1,18 @@
+astro_otter-0.3.0.dist-info/licenses/LICENSE,sha256=s9IPE8A3CAMEaZpDhj4eaorpmfLYGB0mIGphq301PUY,1067
+otter/__init__.py,sha256=pvX-TN7nLVmvKpkDi89Zxe-jMfHNiVMD3zsd_bPEK9Y,535
+otter/_version.py,sha256=XDCAY_avI84JEWPKHKBF9kXkCGBGi2CjNGEItozNsk4,76
+otter/exceptions.py,sha256=3lQF4AXVTfs9VRsVePQoIrXnramsPZbUL5crvf1s9Ng,1702
+otter/schema.py,sha256=eOxlrtp9TTbiENy38ueE2HcYlI-M56g0Ohg7zQeTxjk,10631
+otter/util.py,sha256=xKsNkkxGajML1rZZnR9d5rJV1_z1KeAlIdAs8t0814M,22440
+otter/io/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+otter/io/data_finder.py,sha256=v3jZCOhvysHUQG1FOwHdeJ-psTT-MsdO_GZalBpMBGc,38218
+otter/io/host.py,sha256=xv_SznZuvMoMVsZLqlcmlOyaqKCMZqlTQ_gkN4VBSTw,7139
+otter/io/otter.py,sha256=HiYQ5hgAvPYSHmU6taF6NumpElYzFgjjA0oIgLEo9_4,47471
+otter/io/transient.py,sha256=b_gOQxjk8z3VylZ9ZrAv-Dhqp_d2iXPICL65Hq3kYGg,41846
+otter/plotter/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+otter/plotter/otter_plotter.py,sha256=OQhuLgnMSzgtAjJF8SYBuQOyYcu7Pr0uia5P0G_7z5Q,2144
+otter/plotter/plotter.py,sha256=z09NwQVJS2tuwH3sv95DZv8xogjvf-7Gvj6iWCEx-gQ,9635
+astro_otter-0.3.0.dist-info/METADATA,sha256=7qokfMceKHN1Ai0FNKqY95LYiil9jf3BzAX_7_A5KVI,7127
+astro_otter-0.3.0.dist-info/WHEEL,sha256=ooBFpIzZCPdw3uqIQsOo4qqbA4ZRPxHnOH7peeONza0,91
+astro_otter-0.3.0.dist-info/top_level.txt,sha256=Wth72sCwBRUk3KZGknSKvLQDMFuJk6qiaAavMDOdG5k,6
+astro_otter-0.3.0.dist-info/RECORD,,

{astro_otter-0.1.0.dist-info → astro_otter-0.3.0.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (70.3.0)
+Generator: setuptools (80.0.1)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

otter/__init__.py

@@ -13,4 +13,7 @@ from .io.transient import Transient
 from .io.host import Host
 from .io.data_finder import DataFinder
 from .plotter.otter_plotter import OtterPlotter
-from .plotter.plotter import plot_light_curve, plot_sed
+from .plotter.plotter import plot_light_curve, plot_sed, quick_view, query_quick_view
+from . import util
+from . import schema
+from . import exceptions

otter/_version.py

@@ -2,4 +2,4 @@
 Just define the package version in one place
 """
 
-__version__ = "0.1.0"
+__version__ = "0.3.0"

otter/io/data_finder.py

@@ -11,15 +11,17 @@ import re
 import time
 import math
 from urllib.request import urlopen
+import requests
 
 from astropy import units as u
 from astropy.coordinates import SkyCoord
 from astropy.time import Time
 from astropy.table import Table
 from astropy.io.votable import parse_single_table
+from astropy.io import ascii
 
+import numpy as np
 import pandas as pd
-import requests
 import logging
 
 from fundamentals.stats import rolling_window_sigma_clip
@@ -382,23 +384,84 @@ class DataFinder(object):
         )
         return light_curve.data
 
-    def query_wise(self, radius: float = 5, **kwargs) -> Table:
+    def query_wise(
+        self,
+        radius: float = 5,
+        datadir: str = "ipac/",
+        overwrite: bool = False,
+        verbose=False,
+        **kwargs,
+    ) -> pd.DataFrame:
         """
         Query NEOWISE for their multiepoch photometry
 
+        The method used to query wise here was taken from this github repo:
+        https://github.com/HC-Hwang/wise_light_curves/tree/master
+        and you should cite this other paper that the authors of this code developed
+        it for: https://ui.adsabs.harvard.edu/abs/2020MNRAS.493.2271H/abstract
+
+        This will download the ipac data files to the "datadir" argument. by default,
+        these will go into os.getcwd()/ipac
+
         Args:
             radius (float) : The cone search radius in arcseconds
+            overwrite (bool) : Overwrite the existing datasets downloaded from wise
             **kwargs : Other optional arguments for the astroquery query_region
         Returns:
             An astropy Table of the multiepoch wise data for this host
         """
-        from astroquery.ipac.irsa import Irsa
+        # from https://www.cambridge.org/core/journals/
+        # publications-of-the-astronomical-society-of-australia/article/
+        # recalibrating-the-widefield-infrared-survey-explorer-wise-w4-filter/
+        # B238BFFE19A533A2D2638FE88CCC2E89
+        band_vals = {"w1": 3.4, "w2": 4.6, "w3": 12, "w4": 22}  # in um
 
-        wise_catalogs = "neowiser_p1bs_psd"
-        res = DataFinder._wrap_astroquery(
-            Irsa, self.coord, radius="5 arcsec", catalog=wise_catalogs, **kwargs
+        ra, dec = self.coord.ra.value, self.coord.dec.value
+
+        fbasename = f"wise_{self.name}"
+        allwise_name = f"{fbasename}_allwise.ipac"
+        neowise_name = f"{fbasename}_neowise.ipac"
+
+        if not os.path.exists(datadir):
+            os.makedirs(datadir)
+
+        self._download_single_data(
+            name=fbasename,
+            ra=ra,
+            dec=dec,
+            root_path=datadir,
+            radius=radius,
+            overwrite=overwrite,
         )
-        return res
+
+        allwise = ascii.read(f"ipac/{allwise_name}", format="ipac")
+        neowise = ascii.read(f"ipac/{neowise_name}", format="ipac")
+
+        allwise, neowise = self._only_good_data(allwise, neowise, verbose=verbose)
+        if verbose and (allwise is None or neowise is None):
+            print(f"Limited good infrared data for {self.name}, skipping!")
+
+        mjd, mag, mag_err, filts = self._make_full_lightcurve_multibands(
+            allwise, neowise, bands=["w1", "w2", "w3", "w4"]
+        )
+
+        df = pd.DataFrame(
+            dict(
+                name=[self.name] * len(mjd),
+                date_mjd=mjd,
+                filter=filts,
+                filter_eff=[band_vals[f] for f in filts],
+                filter_eff_unit=["um"] * len(mjd),
+                flux=mag,
+                flux_err=mag_err,
+                flux_unit=["mag(AB)"] * len(mjd),
+                upperlimit=[False] * len(mjd),
+            )
+        )
+
+        # clean up the wise data by filtering out negative flux
+        wise = df[df.flux > 0].reset_index(drop=True)
+        return wise
 
     def query_alma(self, radius: float = 5, **kwargs) -> Table:
         """
@@ -415,6 +478,11 @@ class DataFinder(object):
             An astropy Table of the multiepoch wise data for this host
         """
 
+        logger.warn(
+            "This method may not work if you are using a conda environment!\
+            This is a known issue in setuptools that is not resolved!"
+        )
+
         from astroquery.alma import Alma
 
         res = DataFinder._wrap_astroquery(
@@ -473,9 +541,7 @@ class DataFinder(object):
         )
         return res
 
-    def query_heasarc(
-        self, radius: u.Quantity = 5 * u.arcsec, heasarc_table: str = "xray", **kwargs
-    ) -> Table:
+    def query_heasarc(self, radius: u.Quantity = 5 * u.arcsec, **kwargs) -> Table:
         """
         Query Heasarc by the argument "heasarc_key" for the ra/dec associated with this
         DataLoader object.
@@ -494,9 +560,7 @@ class DataFinder(object):
         """
         from astroquery.heasarc import Heasarc
 
-        res = DataFinder._wrap_astroquery(
-            Heasarc, self.coord, mission=heasarc_table, radius=radius, **kwargs
-        )
+        res = DataFinder._wrap_astroquery(Heasarc, self.coord, radius=radius, **kwargs)
 
         return res
 
@@ -742,3 +806,232 @@ class DataFinder(object):
         print("completed the ``stack_photometry`` method")
 
         return alldata
+
+    """
+    The following code was taken and modified for the purposes of this package from
+    https://github.com/HC-Hwang/wise_light_curves/blob/master/wise_light_curves.py
+
+    Original Authors:
+    - Matthew Hill
+    - Hsiang-Chih Hwang
+
+    Update Author:
+    - Noah Franz
+    """
+
+    @staticmethod
+    def _get_by_position(ra, dec, radius=2.5):
+        allwise_cat = "allwise_p3as_mep"
+        neowise_cat = "neowiser_p1bs_psd"
+        query_url = "http://irsa.ipac.caltech.edu/cgi-bin/Gator/nph-query"
+        payload = {
+            "catalog": allwise_cat,
+            "spatial": "cone",
+            "objstr": " ".join([str(ra), str(dec)]),
+            "radius": str(radius),
+            "radunits": "arcsec",
+            "outfmt": "1",
+        }
+        r = requests.get(query_url, params=payload)
+        allwise = ascii.read(r.text)
+        payload = {
+            "catalog": neowise_cat,
+            "spatial": "cone",
+            "objstr": " ".join([str(ra), str(dec)]),
+            "radius": str(radius),
+            "radunits": "arcsec",
+            "outfmt": "1",
+            "selcols": "ra,dec,sigra,sigdec,sigradec,glon,glat,elon,elat,w1mpro,w1sigmpro,w1snr,w1rchi2,w2mpro,w2sigmpro,w2snr,w2rchi2,rchi2,nb,na,w1sat,w2sat,satnum,cc_flags,det_bit,ph_qual,sso_flg,qual_frame,qi_fact,saa_sep,moon_masked,w1frtr,w2frtr,mjd,allwise_cntr,r_allwise,pa_allwise,n_allwise,w1mpro_allwise,w1sigmpro_allwise,w2mpro_allwise,w2sigmpro_allwise,w3mpro_allwise,w3sigmpro_allwise,w4mpro_allwise,w4sigmpro_allwise",  # noqa: E501
+        }
+        r = requests.get(query_url, params=payload)
+
+        neowise = ascii.read(r.text, guess=False, format="ipac")
+
+        return allwise, neowise
+
+    @staticmethod
+    def _download_single_data(
+        name, ra, dec, root_path="ipac/", radius=2.5, overwrite=False
+    ):
+        # ra, dec: in degree
+        # name, ra, dec = row['Name'], row['RAJ2000'], row['DEJ2000']
+        # name = 'J' + ra + dec
+        if root_path[-1] != "/":
+            root_path += "/"
+        if (
+            not overwrite
+            and os.path.isfile(root_path + name + "_allwise.ipac")
+            and os.path.isfile(root_path + name + "_neowise.ipac")
+        ):
+            pass
+        else:
+            allwise, neowise = DataFinder._get_by_position(ra, dec, radius=radius)
+            allwise.write(
+                root_path + name + "_allwise.ipac", format="ascii.ipac", overwrite=True
+            )
+            neowise.write(
+                root_path + name + "_neowise.ipac", format="ascii.ipac", overwrite=True
+            )
+
+    @staticmethod
+    def _get_data_arrays(table, t, mag, magerr):
+        """Get the time series from a potentially masked astropy table"""
+        if table.masked:
+            full_mask = table[t].mask | table[mag].mask | table[magerr].mask
+            t = table[t].data
+            mag = table[mag].data
+            magerr = table[magerr].data
+
+            t.mask = full_mask
+            mag.mask = full_mask
+            magerr.mask = full_mask
+
+            return t.compressed(), mag.compressed(), magerr.compressed()
+
+        else:
+            return table[t].data, table[mag].data, table[magerr].data
+
+    @staticmethod
+    def _make_full_lightcurve(allwise, neowise, band):
+        """band = 'w1', 'w2', 'w3', or 'w4'"""
+        """Get a combined AllWISE and NEOWISE lightcurve from their Astropy tables"""
+
+        if band not in ["w1", "w2", "w3", "w4"]:
+            raise ValueError("band can only be w1, w2, w3, or w4")
+
+        use_neowise = band in {"w1", "w2"}
+        use_allwise = allwise is not None
+
+        if use_neowise and use_allwise:
+            t, m, e = DataFinder._get_data_arrays(
+                allwise, "mjd", band + "mpro_ep", band + "sigmpro_ep"
+            )
+            t_n, m_n, e_n = DataFinder._get_data_arrays(
+                neowise, "mjd", band + "mpro", band + "sigmpro"
+            )
+            t, m, e = (
+                np.concatenate((t, t_n)),
+                np.concatenate((m, m_n)),
+                np.concatenate((e, e_n)),
+            )
+
+        elif use_neowise and not use_allwise:
+            t, m, e = DataFinder._get_data_arrays(
+                neowise, "mjd", band + "mpro", band + "sigmpro"
+            )
+
+        elif not use_neowise and use_allwise:
+            t, m, e = DataFinder._get_data_arrays(
+                allwise, "mjd", band + "mpro_ep", band + "sigmpro_ep"
+            )
+
+        else:
+            raise Exception("No good allwise or neowise data!")
+
+        t_index = t.argsort()
+        t, m, e = map(lambda e: e[t_index], [t, m, e])
+
+        return t, m, e
+
+    @staticmethod
+    def _make_full_lightcurve_multibands(allwise, neowise, bands=["w1", "w2"]):
+        t, m, e = DataFinder._make_full_lightcurve(allwise, neowise, bands[0])
+        filts = [bands[0] for i in range(len(t))]
+        for band in bands[1:]:
+            try:
+                t_tmp, m_tmp, e_tmp = DataFinder._make_full_lightcurve(
+                    allwise, neowise, band
+                )
+            except Exception:
+                continue
+            t = np.concatenate((t, t_tmp))
+            m = np.concatenate((m, m_tmp))
+            e = np.concatenate((e, e_tmp))
+            filts += [band for i in range(len(t_tmp))]
+        return t, m, e, np.array(filts)
+
+    @staticmethod
+    def _cntr_to_source_id(cntr):
+        cntr = str(cntr)
+
+        # fill leanding 0s
+        if len(cntr) < 19:
+            num_leading_zeros = 19 - len(cntr)
+            cntr = "0" * num_leading_zeros + cntr
+
+        pm = "p"
+        if cntr[4] == "0":
+            pm = "m"
+
+        t = chr(96 + int(cntr[8:10]))
+
+        return "%s%s%s_%cc%s-%s" % (
+            cntr[0:4],
+            pm,
+            cntr[5:8],
+            t,
+            cntr[11:13],
+            cntr[13:19],
+        )
+
+    @staticmethod
+    def _only_good_data(allwise, neowise, verbose=False):
+        """
+        Select good-quality data. The criteria include:
+        - matching the all-wise ID
+
+        To be done:
+        - deal with multiple cntr
+
+        This filtering is described here:
+        https://wise2.ipac.caltech.edu/docs/release/neowise/expsup/sec2_3.html
+        """
+
+        neowise_prefilter_n = len(neowise)
+        neowise = neowise[
+            (neowise["qual_frame"] > 0.0)
+            * (neowise["qi_fact"] > 0.9)
+            * (neowise["saa_sep"] > 0)
+            * (neowise["moon_masked"] == "00")
+        ]
+        neowise_postfilter_n = len(neowise)
+        if verbose:
+            print(
+                f"Filtered out {neowise_prefilter_n-neowise_postfilter_n} neowise \
+                points, leaving {neowise_postfilter_n}"
+            )
+
+        cntr_list = []
+        for data in neowise:
+            if data["allwise_cntr"] not in cntr_list and data["allwise_cntr"] > 10.0:
+                cntr_list.append(data["allwise_cntr"])
+
+        if len(cntr_list) >= 2:
+            print("multiple cntr:")
+            print(cntr_list)
+            return None, neowise
+
+        if len(cntr_list) == 0:
+            # import pdb; pdb.set_trace()
+            # raise Exception('No center!')
+            return None, neowise
+
+        cntr = cntr_list[0]
+
+        source_id = DataFinder._cntr_to_source_id(cntr)
+
+        allwise_prefilter_n = len(allwise)
+        allwise = allwise[
+            (allwise["source_id_mf"] == source_id)
+            * (allwise["saa_sep"] > 0.0)
+            * (allwise["moon_masked"] == "0000")
+            * (allwise["qi_fact"] > 0.9)
+        ]
+        allwise_postfilter_n = len(neowise)
+        if verbose:
+            print(
+                f"Filtered out {allwise_prefilter_n-allwise_postfilter_n} allwise \
+                points, leaving {allwise_postfilter_n}"
+            )
+
+        return allwise, neowise

otter/io/host.py

@@ -4,6 +4,10 @@ for pulling in data corresponding to that host
 """
 
 from __future__ import annotations
+
+from urllib.request import urlopen
+import json
+
 import numpy as np
 from astropy.coordinates import SkyCoord
 from astropy import units as u
@@ -11,6 +15,10 @@ from astropy import units as u
 from .data_finder import DataFinder
 from ..exceptions import OtterLimitationError
 
+import logging
+
+logger = logging.getLogger(__name__)
+
 
 class Host(DataFinder):
     def __init__(
@@ -21,6 +29,7 @@ class Host(DataFinder):
         host_dec_units: str | u.Unit,
         host_name: str = None,
         host_redshift: float = None,
+        redshift_type: str = None,
         reference: list[str] = None,
         transient_name: str = None,
         **kwargs,
@@ -41,6 +50,8 @@ class Host(DataFinder):
                                                       SkyCoord
             host_name (str) : The name of the host galaxy
             host_redshift (float) : The redshift of the host galaxy
+            redshift_type (str) : Either "phot" or "spec", tells you if the redshift
+                                  is a phot-z or spec-z.
             reference (list[str]) : a list of bibcodes that found this to be the host
             transient_name (str) : the name of the transient associated with this host
             kwargs : Just here so we can pass **Transient['host'] into this constructor
@@ -52,6 +63,7 @@ class Host(DataFinder):
         self.redshift = host_redshift  # just here for ease of use
         self.bibcodes = reference
         self.transient_name = transient_name
+        self.redshift_type = redshift_type
 
     def pcc(self, transient_coord: SkyCoord, mag: float = None):
         """
@@ -104,3 +116,71 @@ class Host(DataFinder):
         prob = 1 - np.exp(-eta)
 
         return prob
+
+    ###################################################################################
+    ######### METHODS FOR FINDING HOSTS  ###########################
+    ###################################################################################
+    @staticmethod
+    def query_blast(tns_name: str) -> dict:
+        """
+        Query the BLAST host galaxy service
+
+        Args:
+            tns_name (str) : The TNS target name to grab the data from BLAST for
+
+        Returns:
+            best match BLAST ID, host ra, host dec, host redshift. Redshift will be
+            spectroscopic if available, otherwise photometric.
+        """
+
+        # clean up the input name a little
+        if tns_name[:2] == "AT":
+            tns_name = tns_name.replace("AT", "")
+
+        failed_query_res = None
+
+        # do the query
+        blast_base_url = "https://blast.scimma.org"
+        blast_query_url = f"{blast_base_url}/api/transient/?name={tns_name}&format=json"
+        with urlopen(blast_query_url) as response:
+            if response.status != 200:
+                logger.warn(f"BLAST query failed with response code {response.status}!")
+                return failed_query_res
+            else:
+                blast_data = json.loads(response.read())
+        if len(blast_data) == 0:
+            logger.warn("BLAST query returned no results!")
+            return failed_query_res
+
+        blast_data = blast_data[0]
+        if "host" not in blast_data:
+            logger.warn("BLAST query found the object but it has no host associated!")
+            return failed_query_res
+
+        blast_host = blast_data["host"]
+
+        if blast_host["redshift"] is not None:
+            # prefer spec-z over phot-z
+            z = blast_host["redshift"]
+            z_type = "spec"
+        else:
+            # well I guess we need to use phot-z
+            z = blast_host["photometric_redshift"]
+            z_type = "phot"
+
+        refs = [
+            "2021ApJ...908..170G",  # GHOST citation
+            "2024arXiv241017322J",  # BLAST citation
+        ]
+
+        return Host(
+            host_ra=blast_host["ra_deg"],
+            host_dec=blast_host["dec_deg"],
+            host_ra_units="deg",
+            host_dec_units="deg",
+            host_name=blast_host["id"],
+            host_redshift=z,
+            redshift_type=z_type,
+            reference=refs,
+            transient_name=tns_name,
+        )