From f75c2788deb7d255d734e5ea3a0abe5e2724fddc Mon Sep 17 00:00:00 2001
From: "gregory.ashton" <gregory.ashton@ligo.org>
Date: Thu, 22 Dec 2016 14:26:26 +0100
Subject: [PATCH] Fix prior range when calculating the direct table
---
Paper/DirectedMC/generate_table.py | 15 ++++++++-------
Paper/directed_setup_run_setup.tex | 8 ++++----
2 files changed, 12 insertions(+), 11 deletions(-)
diff --git a/Paper/DirectedMC/generate_table.py b/Paper/DirectedMC/generate_table.py
index 63cf701..c98489f 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 dbb27e8..d639e82 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}
--
GitLab