Merge branch 'nanograv:master' into implicitpar

CHANGELOG-unreleased.md

@@ -38,6 +38,7 @@ the released changes.
 - `pint.models.chromatic_model.Chromatic` as the base class for variable-index chromatic delays.
 - `pint.models.chromatic_model.ChromaticCM` for a Taylor series representation of the variable-index chromatic delay.
 - Whitened residuals (`white-res`) as a plotting axis in `pintk`
+- `TOAs.get_Tspan()` method
 ### Fixed
 - `pint.utils.split_swx()` to use updated `SolarWindDispersionX()` parameter naming convention 
 - Fix #1759 by changing order of comparison

src/pint/simulation.py

@@ -131,10 +131,10 @@ def make_fake_toas(
     name: str = "fake",
     subtract_mean: bool = True,
 ) -> pint.toa.TOAs:
-    """Make toas from an array of times
+    """Make toas from an array of TOAs
 
-    Can include alternating frequencies if fed an array of frequencies,
-    only works with one observatory at a time
+    Uses the input TOAs as a starting place, but then adjusts to fit the given model
+    and then optionally adds noise.
 
     Parameters
     ----------

src/pint/toa.py

@@ -429,7 +429,7 @@ def get_TOAs_list(
         )
     else:
         log.debug(
-            f"Not applying clock corrections since at least some TOAs already have clkcorr in their flags."
+            "Not applying clock corrections since at least some TOAs already have clkcorr in their flags."
         )
     if "tdb" not in t.table.colnames:
         t.compute_TDBs(method=tdb_method, ephem=ephem)
@@ -1690,6 +1690,10 @@ def get_mjds(self, high_precision: bool = False) -> u.Quantity:
         else:
             return self.table["mjd_float"].quantity
 
+    def get_Tspan(self) -> u.Quantity:
+        mjds = self.get_mjds()
+        return mjds.max() - mjds.min()
+
     def get_errors(self) -> u.Quantity:
         """Return a numpy array of the TOA errors in us."""
         return self.table["error"].quantity

tests/test_toa.py

@@ -113,6 +113,7 @@ def test_make_toas(self):
         assert toas.table["error"].unit == u.us
         assert toas.table["mjd"][0].precision == 9
         assert toas.table["mjd"][0].location is not None
+        assert toas.get_Tspan() == toas.table["mjd"].max() - toas.table["mjd"].min()
 
     def test_multiple_observatories_stay_attached(self):
         obs1 = "gbt"