Skip to content
Snippets Groups Projects
Commit f75c2788 authored by Gregory Ashton's avatar Gregory Ashton
Browse files

Fix prior range when calculating the direct table

parent 1ef6d597
No related branches found
No related tags found
No related merge requests found
...@@ -7,7 +7,7 @@ label = 'directed_setup' ...@@ -7,7 +7,7 @@ label = 'directed_setup'
data_label = '{}_data'.format(label) data_label = '{}_data'.format(label)
# Properties of the GW data # Properties of the GW data
sqrtSX = 2e-23 sqrtSX = 1e-23
tstart = 1000000000 tstart = 1000000000
Tspan = 100*86400 Tspan = 100*86400
tend = tstart + Tspan tend = tstart + Tspan
...@@ -16,8 +16,8 @@ tend = tstart + Tspan ...@@ -16,8 +16,8 @@ tend = tstart + Tspan
F0_center = 30 F0_center = 30
F1_center = 1e-10 F1_center = 1e-10
F2 = 0 F2 = 0
Alpha = 5e-3 Alpha = np.radians(83.6292)
Delta = 6e-2 Delta = np.radians(22.0144)
tref = .5*(tstart+tend) tref = .5*(tstart+tend)
...@@ -52,11 +52,11 @@ data.make_data() ...@@ -52,11 +52,11 @@ data.make_data()
predicted_twoF = data.predict_fstat() predicted_twoF = data.predict_fstat()
theta_prior = {'F0': {'type': 'unif', theta_prior = {'F0': {'type': 'unif',
'lower': F0-DeltaF0/2., 'lower': F0-DeltaF0,
'upper': F0+DeltaF0/2.}, 'upper': F0+DeltaF0},
'F1': {'type': 'unif', 'F1': {'type': 'unif',
'lower': F1-DeltaF1/2., 'lower': F1-DeltaF1,
'upper': F1+DeltaF1/2.}, 'upper': F1+DeltaF1},
'F2': F2, 'F2': F2,
'Alpha': Alpha, 'Alpha': Alpha,
'Delta': Delta 'Delta': Delta
...@@ -73,4 +73,5 @@ mcmc = pyfstat.MCMCFollowUpSearch( ...@@ -73,4 +73,5 @@ mcmc = pyfstat.MCMCFollowUpSearch(
tref=tref, minStartTime=tstart, maxStartTime=tend, tref=tref, minStartTime=tstart, maxStartTime=tend,
nwalkers=nwalkers, ntemps=ntemps, nsteps=[nsteps, nsteps], nwalkers=nwalkers, ntemps=ntemps, nsteps=[nsteps, nsteps],
log10temperature_min=log10temperature_min) log10temperature_min=log10temperature_min)
#mcmc.run(Nsegs0=20, R=10)
mcmc.run(run_setup) mcmc.run(run_setup)
\begin{tabular}{c|cccc} \begin{tabular}{c|cccc}
Stage & $\Nseg$ & $\Tcoh^{\rm days}$ &$\Nsteps$ & $\Vpe$ \\ \hline Stage & $\Nseg$ & $\Tcoh^{\rm days}$ &$\Nsteps$ & $\Vpe$ \\ \hline
0 & 20 & 5.0 & 25 & 10.0 \\ 0 & 20 & 5.0 & 25 & 60.0 \\
1 & 7 & 14.3 & 25 & $1{\times}10^{2}$ \\ 1 & 7 & 14.3 & 25 & $5{\times}10^{2}$ \\
2 & 2 & 50.0 & 25 & $1{\times}10^{3}$ \\ 2 & 2 & 50.0 & 25 & $5{\times}10^{3}$ \\
3 & 1 & 100.0 & 25,25 & $2{\times}10^{3}$ \\ 3 & 1 & 100.0 & 25,25 & $1{\times}10^{4}$ \\
\end{tabular} \end{tabular}
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment