
import os,glob
import numpy as np
import matplotlib.pyplot as plt
import scipy.constants as const

# enterprise
import enterprise
from enterprise.pulsar import Pulsar
from enterprise.signals import parameter, white_signals, gp_signals, signal_base, utils

# enterprise extensions
from enterprise_extensions import models, model_utils
from enterprise_extensions.frequentist import optimal_statistic as OS
from enterprise_extensions.deterministic import cw_block_circ
from enterprise_extensions.sampler import setup_sampler

# sampler
from PTMCMCSampler.PTMCMCSampler import PTSampler as ptmcmc

#plotting
import corner

#diagnostics
import la_forge.core as co
import la_forge.diagnostics as dg

/opt/homebrew/Caskroom/mambaforge/base/envs/cw_hunt/lib/python3.10/site-packages/enterprise/signals/utils.py:13: UserWarning: pkg_resources is deprecated as an API. See https://setuptools.pypa.io/en/latest/pkg_resources.html. The pkg_resources package is slated for removal as early as 2025-11-30. Refrain from using this package or pin to Setuptools<81.
  from pkg_resources import Requirement, resource_filename
libstempo not installed. PINT or libstempo are required to use par and tim files.
RuntimeWarning: This platform does not support extended precision floating-point, and PINT will run at reduced precision.

Optional mpi4py package is not installed.  MPI support is not available.

/opt/homebrew/Caskroom/mambaforge/base/envs/cw_hunt/lib/python3.10/site-packages/arviz/__init__.py:50: FutureWarning: 
ArviZ is undergoing a major refactor to improve flexibility and extensibility while maintaining a user-friendly interface.
Some upcoming changes may be backward incompatible.
For details and migration guidance, visit: https://python.arviz.org/en/latest/user_guide/migration_guide.html


DATA_DIR = './hunters/'
par_files = sorted(glob.glob(os.path.join(DATA_DIR, '*.par')))
tim_files = sorted(glob.glob(os.path.join(DATA_DIR, '*.tim')))


psrs = []
for par, tim in zip(par_files, tim_files):
    print(f"Loading {par} and {tim}...")
    psr = Pulsar(par, tim)
    psrs.append(psr)

2026-06-16 07:54:32.281 | DEBUG    | pint.toa:get_TOAs:195 - Using EPHEM = DE440 from the given model
2026-06-16 07:54:32.281 | DEBUG    | pint.toa:get_TOAs:211 - Using CLOCK = BIPM2021 from the given model
2026-06-16 07:54:32.307 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:32.311 | INFO     | pint.observatory:gps_correction:230 - Applying GPS to UTC clock correction (~few nanoseconds)
2026-06-16 07:54:32.311 | DEBUG    | pint.observatory:_load_gps_clock:108 - Loading global GPS clock file
2026-06-16 07:54:32.333 | INFO     | pint.observatory.global_clock_corrections:get_file:129 - File index.txt to be downloaded due to download policy if_expired: https://raw.githubusercontent.com/ipta/pulsar-clock-corrections/main/index.txt

Loading ./hunters/JPSR00.par and ./hunters/JPSR00.tim...

2026-06-16 07:54:37.146 | DEBUG    | pint.observatory.clock_file:__init__:812 - Global clock file gps2utc.clk saving kwargs={'bogus_last_correction': False, 'valid_beyond_ends': False}
2026-06-16 07:54:37.197 | DEBUG    | pint.observatory.clock_file:read_tempo2_clock_file:463 - Loading TEMPO2-format observatory clock correction file gps2utc.clk (/Users/ashokan/.astropy/cache/download/url/d3c81b5766f4bfb84e65504c8a453085/contents) with bogus_last_correction=False
2026-06-16 07:54:37.210 | INFO     | pint.observatory:find_clock_file:973 - Using global clock file for gps2utc.clk with bogus_last_correction=False
2026-06-16 07:54:37.211 | INFO     | pint.observatory:bipm_correction:245 - Applying TT(TAI) to TT(BIPM2021) clock correction (~27 us)
2026-06-16 07:54:37.211 | INFO     | pint.observatory:_load_bipm_clock:119 - Loading BIPM clock version bipm2021
2026-06-16 07:54:37.212 | DEBUG    | pint.observatory.clock_file:__init__:812 - Global clock file tai2tt_bipm2021.clk saving kwargs={'bogus_last_correction': False, 'valid_beyond_ends': False}
2026-06-16 07:54:37.220 | DEBUG    | pint.observatory.clock_file:read_tempo2_clock_file:463 - Loading TEMPO2-format observatory clock correction file tai2tt_bipm2021.clk (/Users/ashokan/.astropy/cache/download/url/e00edeef4edde217d65207a9abeb6a8c/contents) with bogus_last_correction=False
2026-06-16 07:54:37.224 | INFO     | pint.observatory:find_clock_file:973 - Using global clock file for tai2tt_bipm2021.clk with bogus_last_correction=False
2026-06-16 07:54:37.225 | DEBUG    | pint.observatory.clock_file:__init__:812 - Global clock file time_gbt.dat saving kwargs={'bogus_last_correction': False, 'valid_beyond_ends': False}
2026-06-16 07:54:37.231 | DEBUG    | pint.observatory.clock_file:read_tempo_clock_file:613 - Loading TEMPO-format observatory clock correction file time_gbt.dat (/Users/ashokan/.astropy/cache/download/url/599e3ebbfc317e090244ee1ef4c79374/contents) with bogus_last_correction=False
2026-06-16 07:54:37.232 | INFO     | pint.observatory.clock_file:read_tempo_clock_file:703 - Disregarding suspicious MJD -2271.5 in TEMPO clock file
2026-06-16 07:54:37.243 | INFO     | pint.observatory:find_clock_file:973 - Using global clock file for time_gbt.dat with bogus_last_correction=False
2026-06-16 07:54:37.244 | INFO     | pint.observatory.topo_obs:clock_corrections:340 - Applying observatory clock corrections for observatory='gbt'.
2026-06-16 07:54:37.273 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:37.273 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:37.282 | DEBUG    | pint.toa:get_TOAs:310 - Planet PosVels will be calculated.
2026-06-16 07:54:37.282 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:37.338 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.741 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [-1.16671432e+10  1.34133391e+11  5.81464894e+10] m
2026-06-16 07:54:37.741 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.743 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.743 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.744 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.745 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.745 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.746 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.747 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:37.751 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:101 - Creating and dealing with the single TZR_toa for absolute phase
2026-06-16 07:54:37.752 | DEBUG    | pint.toa:__init__:1377 - No pulse number flags found in the TOAs
2026-06-16 07:54:37.752 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:37.755 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:37.755 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:37.756 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:37.757 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [0. 0. 0.] m
2026-06-16 07:54:37.757 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.757 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.758 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.758 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.759 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.759 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.759 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.760 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:37.760 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:122 - Done with TZR_toa
WARNING: Could not find pulsar distance for PSR JPSR00. Setting value to 1 with 20% uncertainty.
2026-06-16 07:54:37.833 | DEBUG    | pint.toa:get_TOAs:195 - Using EPHEM = DE440 from the given model
2026-06-16 07:54:37.834 | DEBUG    | pint.toa:get_TOAs:211 - Using CLOCK = BIPM2021 from the given model
2026-06-16 07:54:37.857 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:37.862 | INFO     | pint.observatory:gps_correction:230 - Applying GPS to UTC clock correction (~few nanoseconds)
2026-06-16 07:54:37.862 | INFO     | pint.observatory:bipm_correction:245 - Applying TT(TAI) to TT(BIPM2021) clock correction (~27 us)
2026-06-16 07:54:37.863 | INFO     | pint.observatory.topo_obs:clock_corrections:340 - Applying observatory clock corrections for observatory='gbt'.
2026-06-16 07:54:37.883 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:37.883 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:37.892 | DEBUG    | pint.toa:get_TOAs:310 - Planet PosVels will be calculated.
2026-06-16 07:54:37.892 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:37.894 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.903 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [-1.16671432e+10  1.34133391e+11  5.81464894e+10] m
2026-06-16 07:54:37.904 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.904 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.905 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.905 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.905 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.906 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.906 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.907 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:37.911 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:101 - Creating and dealing with the single TZR_toa for absolute phase
2026-06-16 07:54:37.912 | DEBUG    | pint.toa:__init__:1377 - No pulse number flags found in the TOAs
2026-06-16 07:54:37.912 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:37.912 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:37.912 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:37.913 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:37.914 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [0. 0. 0.] m
2026-06-16 07:54:37.914 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.914 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.915 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.915 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.915 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.916 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.916 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:37.916 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:37.917 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:122 - Done with TZR_toa
WARNING: Could not find pulsar distance for PSR JPSR01. Setting value to 1 with 20% uncertainty.
2026-06-16 07:54:37.955 | DEBUG    | pint.toa:get_TOAs:195 - Using EPHEM = DE440 from the given model
2026-06-16 07:54:37.956 | DEBUG    | pint.toa:get_TOAs:211 - Using CLOCK = BIPM2021 from the given model
2026-06-16 07:54:37.975 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:37.980 | INFO     | pint.observatory:gps_correction:230 - Applying GPS to UTC clock correction (~few nanoseconds)
2026-06-16 07:54:37.981 | INFO     | pint.observatory:bipm_correction:245 - Applying TT(TAI) to TT(BIPM2021) clock correction (~27 us)
2026-06-16 07:54:37.981 | INFO     | pint.observatory.topo_obs:clock_corrections:340 - Applying observatory clock corrections for observatory='gbt'.

Loading ./hunters/JPSR01.par and ./hunters/JPSR01.tim...
Loading ./hunters/JPSR02.par and ./hunters/JPSR02.tim...

2026-06-16 07:54:38.002 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:38.003 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:38.011 | DEBUG    | pint.toa:get_TOAs:310 - Planet PosVels will be calculated.
2026-06-16 07:54:38.011 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:38.013 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.023 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [-1.16671432e+10  1.34133391e+11  5.81464894e+10] m
2026-06-16 07:54:38.023 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.024 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.024 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.025 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.025 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.026 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.026 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.027 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:38.031 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:101 - Creating and dealing with the single TZR_toa for absolute phase
2026-06-16 07:54:38.032 | DEBUG    | pint.toa:__init__:1377 - No pulse number flags found in the TOAs
2026-06-16 07:54:38.032 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:38.032 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:38.033 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:38.033 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:38.034 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [0. 0. 0.] m
2026-06-16 07:54:38.034 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.035 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.035 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.036 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.036 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.037 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.037 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.037 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:38.038 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:122 - Done with TZR_toa
WARNING: Could not find pulsar distance for PSR JPSR02. Setting value to 1 with 20% uncertainty.
2026-06-16 07:54:38.078 | DEBUG    | pint.toa:get_TOAs:195 - Using EPHEM = DE440 from the given model
2026-06-16 07:54:38.078 | DEBUG    | pint.toa:get_TOAs:211 - Using CLOCK = BIPM2021 from the given model
2026-06-16 07:54:38.098 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:38.102 | INFO     | pint.observatory:gps_correction:230 - Applying GPS to UTC clock correction (~few nanoseconds)
2026-06-16 07:54:38.102 | INFO     | pint.observatory:bipm_correction:245 - Applying TT(TAI) to TT(BIPM2021) clock correction (~27 us)
2026-06-16 07:54:38.103 | INFO     | pint.observatory.topo_obs:clock_corrections:340 - Applying observatory clock corrections for observatory='gbt'.
2026-06-16 07:54:38.123 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:38.123 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:38.130 | DEBUG    | pint.toa:get_TOAs:310 - Planet PosVels will be calculated.
2026-06-16 07:54:38.131 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:38.133 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.143 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [-1.16671433e+10  1.34133391e+11  5.81464894e+10] m
2026-06-16 07:54:38.143 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.144 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.144 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.144 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.145 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.145 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.146 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.147 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:38.151 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:101 - Creating and dealing with the single TZR_toa for absolute phase
2026-06-16 07:54:38.152 | DEBUG    | pint.toa:__init__:1377 - No pulse number flags found in the TOAs
2026-06-16 07:54:38.152 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:38.153 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:38.153 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:38.154 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:38.155 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [0. 0. 0.] m
2026-06-16 07:54:38.155 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.156 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.156 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.157 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.157 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.157 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.158 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.158 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:38.159 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:122 - Done with TZR_toa
WARNING: Could not find pulsar distance for PSR JPSR03. Setting value to 1 with 20% uncertainty.
2026-06-16 07:54:38.200 | DEBUG    | pint.toa:get_TOAs:195 - Using EPHEM = DE440 from the given model
2026-06-16 07:54:38.200 | DEBUG    | pint.toa:get_TOAs:211 - Using CLOCK = BIPM2021 from the given model
2026-06-16 07:54:38.225 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:38.229 | INFO     | pint.observatory:gps_correction:230 - Applying GPS to UTC clock correction (~few nanoseconds)
2026-06-16 07:54:38.230 | INFO     | pint.observatory:bipm_correction:245 - Applying TT(TAI) to TT(BIPM2021) clock correction (~27 us)
2026-06-16 07:54:38.230 | INFO     | pint.observatory.topo_obs:clock_corrections:340 - Applying observatory clock corrections for observatory='gbt'.
2026-06-16 07:54:38.251 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:38.252 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:38.260 | DEBUG    | pint.toa:get_TOAs:310 - Planet PosVels will be calculated.

Loading ./hunters/JPSR03.par and ./hunters/JPSR03.tim...
Loading ./hunters/JPSR04.par and ./hunters/JPSR04.tim...

2026-06-16 07:54:38.260 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:38.263 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.273 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [-1.16671431e+10  1.34133391e+11  5.81464894e+10] m
2026-06-16 07:54:38.273 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.274 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.274 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.275 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.275 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.276 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.276 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.277 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:38.281 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:101 - Creating and dealing with the single TZR_toa for absolute phase
2026-06-16 07:54:38.283 | DEBUG    | pint.toa:__init__:1377 - No pulse number flags found in the TOAs
2026-06-16 07:54:38.283 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:38.283 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:38.284 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:38.284 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:38.285 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [0. 0. 0.] m
2026-06-16 07:54:38.286 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.286 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.287 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.287 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.287 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.288 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.288 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.289 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:38.289 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:122 - Done with TZR_toa
WARNING: Could not find pulsar distance for PSR JPSR04. Setting value to 1 with 20% uncertainty.
2026-06-16 07:54:38.334 | DEBUG    | pint.toa:get_TOAs:195 - Using EPHEM = DE440 from the given model
2026-06-16 07:54:38.334 | DEBUG    | pint.toa:get_TOAs:211 - Using CLOCK = BIPM2021 from the given model
2026-06-16 07:54:38.358 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:38.362 | INFO     | pint.observatory:gps_correction:230 - Applying GPS to UTC clock correction (~few nanoseconds)
2026-06-16 07:54:38.363 | INFO     | pint.observatory:bipm_correction:245 - Applying TT(TAI) to TT(BIPM2021) clock correction (~27 us)
2026-06-16 07:54:38.363 | INFO     | pint.observatory.topo_obs:clock_corrections:340 - Applying observatory clock corrections for observatory='gbt'.
2026-06-16 07:54:38.385 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:38.385 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:38.395 | DEBUG    | pint.toa:get_TOAs:310 - Planet PosVels will be calculated.
2026-06-16 07:54:38.395 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:38.397 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.407 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [-1.16671431e+10  1.34133391e+11  5.81464894e+10] m
2026-06-16 07:54:38.408 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.409 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.409 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.410 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.410 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.410 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.411 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.411 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:38.416 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:101 - Creating and dealing with the single TZR_toa for absolute phase
2026-06-16 07:54:38.417 | DEBUG    | pint.toa:__init__:1377 - No pulse number flags found in the TOAs
2026-06-16 07:54:38.417 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:38.417 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:38.418 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:38.418 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:38.419 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [0. 0. 0.] m
2026-06-16 07:54:38.419 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.420 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.420 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.421 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.421 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.421 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.422 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.422 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:38.423 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:122 - Done with TZR_toa
WARNING: Could not find pulsar distance for PSR JPSR05. Setting value to 1 with 20% uncertainty.
2026-06-16 07:54:38.465 | DEBUG    | pint.toa:get_TOAs:195 - Using EPHEM = DE440 from the given model
2026-06-16 07:54:38.465 | DEBUG    | pint.toa:get_TOAs:211 - Using CLOCK = BIPM2021 from the given model
2026-06-16 07:54:38.490 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:38.494 | INFO     | pint.observatory:gps_correction:230 - Applying GPS to UTC clock correction (~few nanoseconds)
2026-06-16 07:54:38.495 | INFO     | pint.observatory:bipm_correction:245 - Applying TT(TAI) to TT(BIPM2021) clock correction (~27 us)
2026-06-16 07:54:38.495 | INFO     | pint.observatory.topo_obs:clock_corrections:340 - Applying observatory clock corrections for observatory='gbt'.
2026-06-16 07:54:38.515 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.

Loading ./hunters/JPSR05.par and ./hunters/JPSR05.tim...
Loading ./hunters/JPSR06.par and ./hunters/JPSR06.tim...

2026-06-16 07:54:38.516 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:38.524 | DEBUG    | pint.toa:get_TOAs:310 - Planet PosVels will be calculated.
2026-06-16 07:54:38.524 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:38.526 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.535 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [-1.16671432e+10  1.34133391e+11  5.81464894e+10] m
2026-06-16 07:54:38.536 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.537 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.537 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.538 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.538 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.539 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.539 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.540 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:38.544 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:101 - Creating and dealing with the single TZR_toa for absolute phase
2026-06-16 07:54:38.545 | DEBUG    | pint.toa:__init__:1377 - No pulse number flags found in the TOAs
2026-06-16 07:54:38.545 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:38.546 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:38.546 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:38.547 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:38.548 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [0. 0. 0.] m
2026-06-16 07:54:38.548 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.549 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.549 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.550 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.550 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.550 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.551 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.551 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:38.552 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:122 - Done with TZR_toa
WARNING: Could not find pulsar distance for PSR JPSR06. Setting value to 1 with 20% uncertainty.
2026-06-16 07:54:38.596 | DEBUG    | pint.toa:get_TOAs:195 - Using EPHEM = DE440 from the given model
2026-06-16 07:54:38.597 | DEBUG    | pint.toa:get_TOAs:211 - Using CLOCK = BIPM2021 from the given model
2026-06-16 07:54:38.630 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:38.635 | INFO     | pint.observatory:gps_correction:230 - Applying GPS to UTC clock correction (~few nanoseconds)
2026-06-16 07:54:38.636 | INFO     | pint.observatory:bipm_correction:245 - Applying TT(TAI) to TT(BIPM2021) clock correction (~27 us)
2026-06-16 07:54:38.636 | INFO     | pint.observatory.topo_obs:clock_corrections:340 - Applying observatory clock corrections for observatory='gbt'.
2026-06-16 07:54:38.657 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:38.658 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:38.666 | DEBUG    | pint.toa:get_TOAs:310 - Planet PosVels will be calculated.
2026-06-16 07:54:38.666 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:38.668 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.678 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [-1.16671433e+10  1.34133391e+11  5.81464894e+10] m
2026-06-16 07:54:38.678 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.679 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.679 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.680 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.680 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.681 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.681 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.682 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:38.686 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:101 - Creating and dealing with the single TZR_toa for absolute phase
2026-06-16 07:54:38.687 | DEBUG    | pint.toa:__init__:1377 - No pulse number flags found in the TOAs
2026-06-16 07:54:38.688 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:38.689 | DEBUG    | pint.toa:compute_TDBs:2278 - Computing TDB columns.
2026-06-16 07:54:38.689 | DEBUG    | pint.toa:compute_TDBs:2299 - Using EPHEM = DE440 for TDB calculation.
2026-06-16 07:54:38.690 | DEBUG    | pint.toa:compute_posvels:2374 - Computing PosVels of observatories, Earth and planets, using DE440
2026-06-16 07:54:38.691 | DEBUG    | pint.toa:compute_posvels:2432 - SSB obs pos [0. 0. 0.] m
2026-06-16 07:54:38.691 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.692 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.692 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.692 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.693 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.693 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.694 | INFO     | pint.solar_system_ephemerides:load_kernel:194 - Set solar system ephemeris to de440 through astropy
2026-06-16 07:54:38.694 | DEBUG    | pint.toa:compute_posvels:2446 - Adding columns ssb_obs_pos ssb_obs_vel obs_sun_pos obs_jupiter_pos obs_saturn_pos obs_venus_pos obs_uranus_pos obs_neptune_pos obs_earth_pos
2026-06-16 07:54:38.695 | DEBUG    | pint.models.absolute_phase:get_TZR_toa:122 - Done with TZR_toa
WARNING: Could not find pulsar distance for PSR JPSR07. Setting value to 1 with 20% uncertainty.
2026-06-16 07:54:38.741 | DEBUG    | pint.toa:get_TOAs:195 - Using EPHEM = DE440 from the given model
2026-06-16 07:54:38.742 | DEBUG    | pint.toa:get_TOAs:211 - Using CLOCK = BIPM2021 from the given model
2026-06-16 07:54:38.767 | DEBUG    | pint.toa:apply_clock_corrections:2232 - Applying clock corrections (include_bipm = True)
2026-06-16 07:54:38.772 | INFO     | pint.observatory:gps_correction:230 - Applying GPS to UTC clock correction (~few nanoseconds)

Loading ./hunters/JPSR07.par and ./hunters/JPSR07.tim...
Loading ./hunters/JPSR08.par and ./hunters/JPSR08.tim...


len(psrs)

20


print("Pulsar names:")
for psr in psrs:
    print(psr.name)

Pulsar names:
JPSR00
JPSR01
JPSR02
JPSR03
JPSR04
JPSR05
JPSR06
JPSR07
JPSR08
JPSR09
JPSR10
JPSR11
JPSR12
JPSR13
JPSR14
JPSR15
JPSR16
JPSR17
JPSR18
JPSR19


for psr in psrs:
    print(len(psr.toas))

147
147
147
147
147
147
147
147
147
147
147
147
147
147
147
147
147
147
147
147


for psr in psrs:
    end_toa = psr.toas.max()
    start_toa = psr.toas.min()
    timespan_years = (end_toa - start_toa)/(86400*365.25)
    print("Pulsar ",psr.name, " has a timespan of ", timespan_years, " years of data.")

Pulsar  JPSR00  has a timespan of  11.991786387164314  years of data.
Pulsar  JPSR01  has a timespan of  11.99178729202828  years of data.
Pulsar  JPSR02  has a timespan of  11.99178664988047  years of data.
Pulsar  JPSR03  has a timespan of  11.99178733920596  years of data.
Pulsar  JPSR04  has a timespan of  11.991786742696146  years of data.
Pulsar  JPSR05  has a timespan of  11.991787142014326  years of data.
Pulsar  JPSR06  has a timespan of  11.991785764288629  years of data.
Pulsar  JPSR07  has a timespan of  11.99178573576619  years of data.
Pulsar  JPSR08  has a timespan of  11.991785940050242  years of data.
Pulsar  JPSR09  has a timespan of  11.991785780323271  years of data.
Pulsar  JPSR10  has a timespan of  11.991786491849181  years of data.
Pulsar  JPSR11  has a timespan of  11.991786948521243  years of data.
Pulsar  JPSR12  has a timespan of  11.991786741615718  years of data.
Pulsar  JPSR13  has a timespan of  11.991786649727405  years of data.
Pulsar  JPSR14  has a timespan of  11.99178584261107  years of data.
Pulsar  JPSR15  has a timespan of  11.991786726815233  years of data.
Pulsar  JPSR16  has a timespan of  11.991785793120162  years of data.
Pulsar  JPSR17  has a timespan of  11.991786812230094  years of data.
Pulsar  JPSR18  has a timespan of  11.991785929103578  years of data.
Pulsar  JPSR19  has a timespan of  11.991787134767453  years of data.


# timing model
tm = gp_signals.TimingModel()

#Noise can be pre-estimated from the single-pulsar noise analysis. Here, lets try to estimate white-noise along with the gravitational wave signal.
efac = parameter.Uniform(0.1, 5.0)

# white noise
ef = white_signals.MeasurementNoise(efac=efac)


cw = cw_block_circ(amp_prior='log-uniform',#
                   skyloc=None,        # None = search full sky
                   log10_fgw=None,     # None = search over frequency
                   psrTerm=False,       # include pulsar term?
                   tref=53000*86400)   # reference time in seconds


# full model = white noise + timing model + CW
model = ef + tm + cw

# initialize PTA
pulsar_models = []
for psr in psrs:
    print(f"Setting up model for {psr.name}...")
    model(psr)  # This initializes the model for the pulsar
    pulsar_models.append(model(psr))

pta = signal_base.PTA(pulsar_models)

Setting up model for JPSR00...
Setting up model for JPSR01...
Setting up model for JPSR02...
Setting up model for JPSR03...
Setting up model for JPSR04...
Setting up model for JPSR05...
Setting up model for JPSR06...
Setting up model for JPSR07...
Setting up model for JPSR08...
Setting up model for JPSR09...
Setting up model for JPSR10...
Setting up model for JPSR11...
Setting up model for JPSR12...
Setting up model for JPSR13...
Setting up model for JPSR14...
Setting up model for JPSR15...
Setting up model for JPSR16...
Setting up model for JPSR17...
Setting up model for JPSR18...
Setting up model for JPSR19...


print(pta.summary())

enterprise v3.4.4,  Python v3.10.20
==========================================================================================

Signal Name                              Signal Class                   no. Parameters      
==========================================================================================
JPSR00_measurement_noise                 MeasurementNoise               1                   

params:
JPSR00_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR00_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR00_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR01_measurement_noise                 MeasurementNoise               1                   

params:
JPSR01_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR01_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR01_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR02_measurement_noise                 MeasurementNoise               1                   

params:
JPSR02_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR02_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR02_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR03_measurement_noise                 MeasurementNoise               1                   

params:
JPSR03_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR03_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR03_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR04_measurement_noise                 MeasurementNoise               1                   

params:
JPSR04_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR04_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR04_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR05_measurement_noise                 MeasurementNoise               1                   

params:
JPSR05_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR05_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR05_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR06_measurement_noise                 MeasurementNoise               1                   

params:
JPSR06_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR06_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR06_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR07_measurement_noise                 MeasurementNoise               1                   

params:
JPSR07_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR07_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR07_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR08_measurement_noise                 MeasurementNoise               1                   

params:
JPSR08_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR08_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR08_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR09_measurement_noise                 MeasurementNoise               1                   

params:
JPSR09_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR09_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR09_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR10_measurement_noise                 MeasurementNoise               1                   

params:
JPSR10_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR10_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR10_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR11_measurement_noise                 MeasurementNoise               1                   

params:
JPSR11_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR11_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR11_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR12_measurement_noise                 MeasurementNoise               1                   

params:
JPSR12_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR12_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR12_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR13_measurement_noise                 MeasurementNoise               1                   

params:
JPSR13_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR13_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR13_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR14_measurement_noise                 MeasurementNoise               1                   

params:
JPSR14_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR14_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR14_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR15_measurement_noise                 MeasurementNoise               1                   

params:
JPSR15_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR15_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR15_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR16_measurement_noise                 MeasurementNoise               1                   

params:
JPSR16_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR16_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR16_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR17_measurement_noise                 MeasurementNoise               1                   

params:
JPSR17_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR17_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR17_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR18_measurement_noise                 MeasurementNoise               1                   

params:
JPSR18_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR18_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR18_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
JPSR19_measurement_noise                 MeasurementNoise               1                   

params:
JPSR19_efac:Uniform(pmin=0.1, pmax=5.0)                                                   
__________________________________________________________________________________________
JPSR19_linear_timing_model               TimingModel                    0                   

params:
__________________________________________________________________________________________
JPSR19_cw                                CWSignal                       8                   

params:
cw_costheta:Uniform(pmin=-1, pmax=1)                                                      
cw_phi:Uniform(pmin=0, pmax=6.283185307179586)                                            
cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)                                                    
cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)                                                  
cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)                                                
cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)                                                
cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)                                       
cw_psi:Uniform(pmin=0, pmax=3.141592653589793)                                            
__________________________________________________________________________________________
==========================================================================================
Total params: 28
Varying params: 28
Common params: 160
Fixed params: 0
Number of pulsars: 20


print(f"Number of parameters: {len(pta.params)}")
for p in pta.params:
    print(f"  {p}")

Number of parameters: 28
  JPSR00_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR01_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR02_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR03_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR04_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR05_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR06_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR07_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR08_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR09_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR10_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR11_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR12_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR13_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR14_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR15_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR16_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR17_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR18_efac:Uniform(pmin=0.1, pmax=5.0)
  JPSR19_efac:Uniform(pmin=0.1, pmax=5.0)
  cw_cosinc:Uniform(pmin=-1.0, pmax=1.0)
  cw_costheta:Uniform(pmin=-1, pmax=1)
  cw_log10_Mc:Uniform(pmin=6.0, pmax=10.0)
  cw_log10_fgw:Uniform(pmin=-9.0, pmax=-7.0)
  cw_log10_h:Uniform(pmin=-18.0, pmax=-11.0)
  cw_phase0:Uniform(pmin=0.0, pmax=6.283185307179586)
  cw_phi:Uniform(pmin=0, pmax=6.283185307179586)
  cw_psi:Uniform(pmin=0, pmax=3.141592653589793)


#set up an initial starting point for the sampler
#make a dictionary of parameter names and random samples from their priors
xs = {}
for par in pta.params:
    xs[par.name] = par.sample()

#flatten the dictionary into a single array of parameter values
x0 = np.hstack([xs[par.name] for par in pta.params])
ndim = len(x0)


cov = np.diag(np.ones(ndim) * 0.1**2)

# Group parameters for more efficient sampling
# Group 1: all parameters together
groups = [list(range(0, ndim))]

# Group 2: global CW parameters only (last 8)
groups.append(list(range(ndim-8, ndim)))


sampler = ptmcmc(ndim, pta.get_lnlikelihood, pta.get_lnprior, cov,
                 groups=groups,
                 outDir='./results/',
                 resume=True)

N = 30000  #short run
#x0 = np.hstack([p.sample() for p in pta.params])
sampler.sample(x0, N, SCAMweight=30, AMweight=15, DEweight=50)

Finished 0.00 percent (0.00 percent of new work) in 0.001435 s Acceptance rate = 0

/opt/homebrew/Caskroom/mambaforge/base/envs/cw_hunt/lib/python3.10/site-packages/enterprise/signals/parameter.py:70: RuntimeWarning: divide by zero encountered in log

Finished 33.33 percent (33.33 percent of new work) in 25.210956 s Acceptance rate = 0.445711Adding DE jump with weight 50
Finished 100.00 percent (100.00 percent of new work) in 62.125574 s Acceptance rate = 0.5265
Run Complete


chain = np.loadtxt('./results/chain_1.txt')
pars = sorted(xs.keys())

# Find indices of global CW parameters (not per-pulsar)
cw_par_names = [p for p in pars if p.startswith('cw_')]
cw_par_idx = [pars.index(p) for p in cw_par_names]

print("CW parameters and their indices:")
for name, idx in zip(cw_par_names, cw_par_idx):
    print(f"  [{idx}] {name}")

# Plot corner
burn = int(0.25 * chain.shape[0])
cw_samples = chain[burn:, cw_par_idx]

fig = corner.corner(cw_samples,
                    labels=cw_par_names,
                    smooth=True,
                    show_titles=True,
                    title_kwargs={'fontsize': 9})
plt.show()

CW parameters and their indices:
  [20] cw_cosinc
  [21] cw_costheta
  [22] cw_log10_Mc
  [23] cw_log10_fgw
  [24] cw_log10_h
  [25] cw_phase0
  [26] cw_phi
  [27] cw_psi


#import la_forge.core as co
#import la_forge.diagnostics as dg
np.savetxt('./results/pars.txt', pars, fmt='%s')
chaindir = './results/'
c0 = co.Core(chaindir, burn=0.25, label='cw_hunt')
print('Number of parameters: \t',len(c0.params))
print('Chain shape: \t\t',c0.chain.shape)
print('Burn: \t\t\t' , c0.burn)

Number of parameters: 	 32
Chain shape: 		 (3001, 32)
Burn: 			 750


#check trace pltos
dg.plot_chains(c0,
               hist=False,
               pars=c0.params[:],
               ncols=4)


# Look at the Gelman Rubin diagnostic
dg.plot_grubin(c0)


# sampler = ptmcmc(ndim, pta.get_lnlikelihood, pta.get_lnprior, cov,
#                  groups=groups,
#                  outDir='./results/',
#                  resume=True)

# N = 5000000  # 3M samples
# x0 = np.hstack([p.sample() for p in pta.params])
# sampler.sample(x0, N, SCAMweight=30, AMweight=15, DEweight=50)

Resuming run from chain file ./results/cw_hunt_nopsrterm_9p5//chain_1.txt
Resuming with 386201 samples from file representing 3862001 original samples
Adding DE jump with weight 50

/opt/homebrew/Caskroom/mambaforge/base/envs/cw_hunt/lib/python3.10/site-packages/enterprise/signals/parameter.py:70: RuntimeWarning: divide by zero encountered in log

Finished 100.00 percent (100.00 percent of new work) in 3368.140394 s Acceptance rate = 0.154316
Run Complete


chain = np.loadtxt('./cw_search_predone/chain_1.txt')
pars = sorted(xs.keys())

cw_par_names = [p for p in pars if p.startswith('cw_')]# and 'p_dist' not in p and 'p_phase' not in p]
cw_par_idx = [pars.index(p) for p in cw_par_names]

# Plot corner
burn = int(0.5 * chain.shape[0])
cw_samples = chain[burn:, cw_par_idx]

fig = corner.corner(cw_samples,
                    #bins=40,
                    labels=cw_par_names,
                    smooth=True,
                    show_titles=True,
                    title_kwargs={'fontsize': 9})
plt.show()

THE IPTA 2026 CW HUNT STARTS HERE!¶

Assemble hunters¶

Make Enterprise Pulsar objects¶

Understand the data¶

Homework : Where are the pulsars located on the sky?¶

Set up PTA!¶

We will search for circular orbit CW signals. So we will use the `cw_block_circ` function from enterprise_extensions to create a signal model.¶

Always make sure you are happy with how your PTA looks!¶

What are the 28 parameters we are varying?!¶

What are the 160 common params?¶

How might this look like if we decided to search over pulsar-terms?¶

Prepare to sampler!¶

Plot posteriors of CW parameters!¶

Time to Diagnose the search.¶

What do you think?¶

Run longer!¶

In the interest of time, load this pre-done search.¶

THE IPTA 2026 CW HUNT STARTS HERE!¶

Assemble hunters¶

Make Enterprise Pulsar objects¶

Understand the data¶

Homework : Where are the pulsars located on the sky?¶

Set up PTA!¶

We will search for circular orbit CW signals. So we will use the cw_block_circ function from enterprise_extensions to create a signal model.¶

Always make sure you are happy with how your PTA looks!¶

What are the 28 parameters we are varying?!¶

What are the 160 common params?¶

How might this look like if we decided to search over pulsar-terms?¶

Prepare to sampler!¶

Plot posteriors of CW parameters!¶

Time to Diagnose the search.¶

What do you think?¶

Run longer!¶

In the interest of time, load this pre-done search.¶

We will search for circular orbit CW signals. So we will use the `cw_block_circ` function from enterprise_extensions to create a signal model.¶