A Bayesian analysis pipeline for continuous GW sources in the PTA band

The direct detection of gravitational waves (GWs) by pulsar timing arrays (PTAs) is very likely within the next decade. While the stochastic GW background is a promising candidate for detection it is also possible that single resolvable sources may be detectable as well. In this work we will focus on the detection and characterization of single GW sources from supermassive black hole binaries. We introduce a fully Bayesian data analysis pipeline that is meant to carry out a search, characterization, and evaluation phase. This will allow us to rapidly locate the global maxima in parameter space, map out the posterior, and finally weigh the evidence of a GW detection through a Bayes factor. Here we will make use of an adaptive metropolis algorithm and parallel tempering. We test this algorithm on realistic simulated data that are representative of modern PTAs. (paper)