diff --git a/Paper/DirectedMC/generate_table.py b/Paper/DirectedMC/generate_table.py index 63cf7015739d0962f4017c3942204d567ea2c032..c98489f8562c1205c8a6fd60b817c685258f1433 100644 --- a/Paper/DirectedMC/generate_table.py +++ b/Paper/DirectedMC/generate_table.py @@ -7,7 +7,7 @@ label = 'directed_setup' data_label = '{}_data'.format(label) # Properties of the GW data -sqrtSX = 2e-23 +sqrtSX = 1e-23 tstart = 1000000000 Tspan = 100*86400 tend = tstart + Tspan @@ -16,8 +16,8 @@ tend = tstart + Tspan F0_center = 30 F1_center = 1e-10 F2 = 0 -Alpha = 5e-3 -Delta = 6e-2 +Alpha = np.radians(83.6292) +Delta = np.radians(22.0144) tref = .5*(tstart+tend) @@ -52,11 +52,11 @@ data.make_data() predicted_twoF = data.predict_fstat() theta_prior = {'F0': {'type': 'unif', - 'lower': F0-DeltaF0/2., - 'upper': F0+DeltaF0/2.}, + 'lower': F0-DeltaF0, + 'upper': F0+DeltaF0}, 'F1': {'type': 'unif', - 'lower': F1-DeltaF1/2., - 'upper': F1+DeltaF1/2.}, + 'lower': F1-DeltaF1, + 'upper': F1+DeltaF1}, 'F2': F2, 'Alpha': Alpha, 'Delta': Delta @@ -73,4 +73,5 @@ mcmc = pyfstat.MCMCFollowUpSearch( tref=tref, minStartTime=tstart, maxStartTime=tend, nwalkers=nwalkers, ntemps=ntemps, nsteps=[nsteps, nsteps], log10temperature_min=log10temperature_min) +#mcmc.run(Nsegs0=20, R=10) mcmc.run(run_setup) diff --git a/Paper/directed_setup_run_setup.tex b/Paper/directed_setup_run_setup.tex index dbb27e8eb2d1a9310228ff4d720e32734115ee87..d639e825dd7eca4330ecf6ac763f887320223e8f 100644 --- a/Paper/directed_setup_run_setup.tex +++ b/Paper/directed_setup_run_setup.tex @@ -1,7 +1,7 @@ \begin{tabular}{c|cccc} Stage & $\Nseg$ & $\Tcoh^{\rm days}$ &$\Nsteps$ & $\Vpe$ \\ \hline -0 & 20 & 5.0 & 25 & 10.0 \\ -1 & 7 & 14.3 & 25 & $1{\times}10^{2}$ \\ -2 & 2 & 50.0 & 25 & $1{\times}10^{3}$ \\ -3 & 1 & 100.0 & 25,25 & $2{\times}10^{3}$ \\ +0 & 20 & 5.0 & 25 & 60.0 \\ +1 & 7 & 14.3 & 25 & $5{\times}10^{2}$ \\ +2 & 2 & 50.0 & 25 & $5{\times}10^{3}$ \\ +3 & 1 & 100.0 & 25,25 & $1{\times}10^{4}$ \\ \end{tabular}