  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 (SMBHBs). 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 (AM)
algorithm and parallel tempering. We test this algorithm on realistic simulated
data that are representative of modern PTAs.
