def METHOD_NAME(self, toas):
    """Get the TOAs class for the TZRMJD.
    We are treating the TZRMJD as a special TOA.
    Note that any observatory clock corrections will be applied
    to this TOA, as with any other TOA. This does not affect the
    value of the TZRMJD parameter, however.
    """
    clkc_info = toas.clock_corr_info
    # If we have cached the TZR TOA and all the TZR* and clock info has not changed, then don't rebuild it
    if self.tz_cache is not None and (
        self.tz_clkc_info["include_bipm"] == clkc_info["include_bipm"]
        and self.tz_clkc_info["include_gps"] == clkc_info["include_gps"]
        and self.tz_planets == toas.planets
        and self.tz_ephem == toas.ephem
        and self.tz_hash
        == hash((self.TZRMJD.value, self.TZRSITE.value, self.TZRFRQ.value))
    ):
        return self.tz_cache
    # Otherwise we have to build the TOA and apply clock corrections
    # NOTE: Using TZRMJD.quantity.jd[1,2] so that the time scale can be properly
    # set to the TZRSITE default timescale (e.g. UTC for TopoObs and TDB for SSB)
    log.debug("Creating and dealing with the single TZR_toa for absolute phase")
    # TZR_toa = toa.TOA(
    #     (self.TZRMJD.quantity.jd1 - 2400000.5, self.TZRMJD.quantity.jd2),
    #     obs=self.TZRSITE.value,
    #     freq=self.TZRFRQ.quantity,
    # )
    # tz = toa.get_TOAs_list(
    #     [TZR_toa],
    #     include_bipm=clkc_info["include_bipm"],
    #     include_gps=clkc_info["include_gps"],
    #     ephem=toas.ephem,
    #     planets=toas.planets,
    # )
    tz = toa.get_TOAs_array(
        (self.TZRMJD.quantity.jd1 - 2400000.5, self.TZRMJD.quantity.jd2),
        obs=self.TZRSITE.value,
        freqs=self.TZRFRQ.quantity,
        include_bipm=clkc_info["include_bipm"],
        include_gps=clkc_info["include_gps"],
        ephem=toas.ephem,
        planets=toas.planets,
        tzr=True,
    )
    log.debug("Done with TZR_toa")
    self.tz_cache = tz
    self.tz_hash = hash((self.TZRMJD.value, self.TZRSITE.value, self.TZRFRQ.value))
    self.tz_clkc_info = clkc_info
    self.tz_planets = toas.planets
    self.tz_ephem = toas.ephem
    return tz