Long transient MCMC search

MCMC search for a long transient CW signal.

 import pyfstat
 import os
 import numpy as np

 outdir = os.path.join("PyFstat_example_data", "PyFstat_example_long_transient_search")
 if not os.path.isdir(outdir) or not np.any(
     [f.endswith(".sft") for f in os.listdir(outdir)]
 ):
     raise RuntimeError(
         "Please first run PyFstat_example_make_data_for_long_transient_search.py !"
     )

 tstart = 1000000000
 duration = 200 * 86400

 inj = {
     "tref": tstart,
     "F0": 30.0,
     "F1": -1e-10,
     "F2": 0,
     "Alpha": 0.5,
     "Delta": 1,
     "transient_tstart": tstart + 0.25 * duration,
     "transient_duration": 0.5 * duration,
 }

 DeltaF0 = 6e-7
 DeltaF1 = 1e-13

 # to make the search cheaper, we exactly target the transientStartTime
 # to the injected value and only search over TransientTau
 theta_prior = {
     "F0": {
         "type": "unif",
         "lower": inj["F0"] - DeltaF0 / 2.0,
         "upper": inj["F0"] + DeltaF0 / 2.0,
     },
     "F1": {
         "type": "unif",
         "lower": inj["F1"] - DeltaF1 / 2.0,
         "upper": inj["F1"] + DeltaF1 / 2.0,
     },
     "F2": inj["F2"],
     "Alpha": inj["Alpha"],
     "Delta": inj["Delta"],
     "transient_tstart": tstart + 0.25 * duration,
     "transient_duration": {
         "type": "halfnorm",
         "loc": 0.001 * duration,
         "scale": 0.5 * duration,
     },
 }

 ntemps = 2
 log10beta_min = -1
 nwalkers = 100
 nsteps = [100, 100]

 mcmc = pyfstat.MCMCTransientSearch(
     label="transient_search",
     outdir=outdir,
     sftfilepattern=os.path.join(outdir, "*simulated_transient_signal*sft"),
     theta_prior=theta_prior,
     tref=inj["tref"],
     nsteps=nsteps,
     nwalkers=nwalkers,
     ntemps=ntemps,
     log10beta_min=log10beta_min,
     transientWindowType="rect",
 )
 mcmc.run(walker_plot_args={"plot_det_stat": True, "injection_parameters": inj})
 mcmc.print_summary()
 mcmc.plot_corner(add_prior=True, truths=inj)
 mcmc.plot_prior_posterior(injection_parameters=inj)