diff --git a/Paper/Examples/fully_coherent_search_using_MCMC.py b/Paper/Examples/fully_coherent_search_using_MCMC.py index bb53c79937bf7920bf185e590fca0db1bdf192a4..ec34696d281e3f31bb69fa7b3d296511dbfd8eca 100644 --- a/Paper/Examples/fully_coherent_search_using_MCMC.py +++ b/Paper/Examples/fully_coherent_search_using_MCMC.py @@ -53,7 +53,7 @@ nsteps = [50, 50] mcmc = pyfstat.MCMCSearch( label='fully_coherent_search_using_MCMC', outdir='data', - sftfilepath='data/*'+data_label+'*sft', theta_prior=theta_prior, tref=tref, + sftfilepath='data/*'+data_label+'-*sft', theta_prior=theta_prior, tref=tref, minStartTime=tstart, maxStartTime=tend, nsteps=nsteps, nwalkers=nwalkers, ntemps=ntemps, log10temperature_min=log10temperature_min) mcmc.run(context='paper', subtractions=[30, -1e-10]) diff --git a/Paper/Examples/transient_search_using_MCMC.py b/Paper/Examples/transient_search_using_MCMC.py index 7eefea0b53e01dda410a7998fc8bb40f5526228d..b45fcf545819dc63d21792e55358caefe929f944 100644 --- a/Paper/Examples/transient_search_using_MCMC.py +++ b/Paper/Examples/transient_search_using_MCMC.py @@ -55,7 +55,10 @@ mcmc = pyfstat.MCMCSearch( nwalkers=nwalkers, ntemps=ntemps, log10temperature_min=log10temperature_min) mcmc.run() -mcmc.plot_cumulative_max() +fig, ax = plt.subplots() +mcmc.plot_cumulative_max(ax=ax) +ax.set_xlabel('Days from $t_\mathrm{start}$') +fig.savefig('data/transient_search_initial_stage_twoFcumulative') mcmc.print_summary() theta_prior = {'F0': {'type': 'unif', @@ -69,9 +72,9 @@ theta_prior = {'F0': {'type': 'unif', 'Delta': Delta, 'transient_tstart': {'type': 'unif', 'lower': data_tstart, - 'upper': data_tend}, + 'upper': data_tend-0.2*duration}, 'transient_duration': {'type': 'halfnorm', - 'loc': 0, + 'loc': 0.01*duration, 'scale': 0.5*duration} } diff --git a/Paper/follow_up_run_setup.tex b/Paper/follow_up_run_setup.tex index 711cb402182e5f06fc668f4b400c092986be34f3..b7b46fff1e8bd0096a67689a0216ea6cfd8f6299 100644 --- a/Paper/follow_up_run_setup.tex +++ b/Paper/follow_up_run_setup.tex @@ -1,11 +1,11 @@ \begin{tabular}{c|cccccc} Stage & $\Nseg$ & $\Tcoh^{\rm days}$ &$\Nsteps$ & $\V$ & $\Vsky$ & $\Vpe$ \\ \hline -0 & 100 & 1.0 & 200 & $1{\times}10^{2}$ & 7.0 & 10.0 \\ -1 & 56 & 1.8 & 200 & $1{\times}10^{3}$ & 20.0 & 40.0 \\ -2 & 31 & 3.2 & 200 & $1{\times}10^{4}$ & 70.0 & $1{\times}10^{2}$ \\ -3 & 17 & 5.9 & 200 & $1{\times}10^{5}$ & $2{\times}10^{2}$ & $5{\times}10^{2}$ \\ -4 & 9 & 11.1 & 200 & $1{\times}10^{6}$ & $8{\times}10^{2}$ & $2{\times}10^{3}$ \\ -5 & 5 & 20.0 & 200 & $1{\times}10^{7}$ & $2{\times}10^{3}$ & $5{\times}10^{3}$ \\ -6 & 2 & 50.0 & 200 & $1{\times}10^{8}$ & $3{\times}10^{3}$ & $3{\times}10^{4}$ \\ -7 & 1 & 100.0 & 200,200 & $3{\times}10^{8}$ & $4{\times}10^{3}$ & $6{\times}10^{4}$ \\ +0 & 100 & 1.0 & 200 & 95.0 & 6.8 & 14.0 \\ +1 & 56 & 1.8 & 200 & $9.6{\times}10^{2}$ & 22.0 & 45.0 \\ +2 & 31 & 3.2 & 200 & $1.0{\times}10^{4}$ & 69.0 & $1.5{\times}10^{2}$ \\ +3 & 17 & 5.9 & 200 & $1.1{\times}10^{5}$ & $2.3{\times}10^{2}$ & $4.8{\times}10^{2}$ \\ +4 & 9 & 11.1 & 200 & $1.3{\times}10^{6}$ & $7.5{\times}10^{2}$ & $1.7{\times}10^{3}$ \\ +5 & 5 & 20.0 & 200 & $1.1{\times}10^{7}$ & $2.0{\times}10^{3}$ & $5.5{\times}10^{3}$ \\ +6 & 2 & 50.0 & 200 & $1.0{\times}10^{8}$ & $3.3{\times}10^{3}$ & $3.1{\times}10^{4}$ \\ +7 & 1 & 100.0 & 200,200 & $2.7{\times}10^{8}$ & $4.4{\times}10^{3}$ & $6.2{\times}10^{4}$ \\ \end{tabular} diff --git a/Paper/follow_up_walkers.png b/Paper/follow_up_walkers.png index 9f8511b98ebfde60cae5ffaa875c3d909c54582d..5a6c9fdc23cc6fdfbd6b40751041fe5171f8d2be 100644 Binary files a/Paper/follow_up_walkers.png and b/Paper/follow_up_walkers.png differ diff --git a/Paper/fully_coherent_search_using_MCMC_walkers.png b/Paper/fully_coherent_search_using_MCMC_walkers.png index d7f82cc39211b04ecaf3a421e705103e1ef1b97d..e50c8205bf3a63fd851b3fab2f1f1fc5461db0a1 100644 Binary files a/Paper/fully_coherent_search_using_MCMC_walkers.png and b/Paper/fully_coherent_search_using_MCMC_walkers.png differ diff --git a/Paper/paper_cw_mcmc.tex b/Paper/paper_cw_mcmc.tex index 354041d1b4fa16bbfd58e05147de22f9f40b7200..eb5c4a694ce22125213ecd529d8926365b5eafa2 100644 --- a/Paper/paper_cw_mcmc.tex +++ b/Paper/paper_cw_mcmc.tex @@ -53,15 +53,15 @@ \title{MCMC follow-up of continuous gravitational wave candidates} - \author{G. Ashton} - \email[E-mail: ]{gregory.ashton@ligo.org} - \affiliation{Max Planck Institut f{\"u}r Gravitationsphysik - (Albert Einstein Institut) and Leibniz Universit\"at Hannover, - 30161 Hannover, Germany} - \author{R. Prix} - \affiliation{Max Planck Institut f{\"u}r Gravitationsphysik - (Albert Einstein Institut) and Leibniz Universit\"at Hannover, - 30161 Hannover, Germany} +% \author{G. Ashton} +% \email[E-mail: ]{gregory.ashton@ligo.org} +% \affiliation{Max Planck Institut f{\"u}r Gravitationsphysik +% (Albert Einstein Institut) and Leibniz Universit\"at Hannover, +% 30161 Hannover, Germany} +% \author{R. Prix} +% \affiliation{Max Planck Institut f{\"u}r Gravitationsphysik +% (Albert Einstein Institut) and Leibniz Universit\"at Hannover, +% 30161 Hannover, Germany} \date{\today} @@ -1215,13 +1215,13 @@ the improvement in fit due to the inclustion of the transient parameters was sufficient to compensate for the greater prior volume and produce an improvement in the evidence for the model. -\begin{figure}[htb] +\begin{figure*}[htb] \centering -\includegraphics[width=0.45\textwidth]{transient_search_corner} +\includegraphics[width=0.8\textwidth]{transient_search_corner} \caption{Posterior distributions for a targeted search of data containing a simulated transient signal and Gaussian noise.} \label{fig_transient_posterior} -\end{figure} +\end{figure*} diff --git a/Paper/single_glitch_glitchSearch_corner.png b/Paper/single_glitch_glitchSearch_corner.png index de502c11e3c407a0136e6a32009f9a3c5922bea0..95bfceb917d3b8f061e121ed78f57c442addc411 100644 Binary files a/Paper/single_glitch_glitchSearch_corner.png and b/Paper/single_glitch_glitchSearch_corner.png differ diff --git a/Paper/transient_search_corner.png b/Paper/transient_search_corner.png index cc829501625520d77e5917b6f6e6be7460f962e8..bff7f78e921bf2684c07f91e1201084b15587867 100644 Binary files a/Paper/transient_search_corner.png and b/Paper/transient_search_corner.png differ diff --git a/Paper/transient_search_initial_stage_twoFcumulative.png b/Paper/transient_search_initial_stage_twoFcumulative.png index 987c54ae52727f6cf61c2136a1545ed95fb00fc9..85fb01c8e8d865d7262a9b0102ec2798cf221b4e 100644 Binary files a/Paper/transient_search_initial_stage_twoFcumulative.png and b/Paper/transient_search_initial_stage_twoFcumulative.png differ