From 9de7f3a82006cb642b709862b03e46e284c3e7c5 Mon Sep 17 00:00:00 2001
From: Marc Arene <arene.marc@gmail.com>
Date: Thu, 11 Jun 2020 19:33:19 -0500
Subject: [PATCH] inconsistent_prior_files
---
bilby/gw/prior.py | 8 +++----
...ck_holes.prior => aligned_spins_bbh.prior} | 4 ++--
.../aligned_spins_bbh_tides_on.prior | 18 +++++++++++++++
bilby/gw/prior_files/aligned_spins_bns.prior | 16 ++++++++++++++
.../aligned_spins_bns_tides_on.prior | 18 +++++++++++++++
.../aligned_spins_waveform_tides_off.prior | 12 ----------
.../aligned_spins_waveform_tides_on.prior | 14 ------------
.../gw/prior_files/precessing_spins_bbh.prior | 12 +++++-----
.../precessing_spins_bbh_tides_on.prior | 22 +++++++++++++++++++
.../gw/prior_files/precessing_spins_bns.prior | 20 +++++++++++++++++
.../precessing_spins_bns_tides_on.prior | 22 +++++++++++++++++++
.../precessing_spins_waveform_tides_off.prior | 16 --------------
.../precessing_spins_waveform_tides_on.prior | 18 ---------------
test/gw_prior_test.py | 8 +++----
test/prior_files/aligned_spins_bbh.prior | 16 ++++++++++++++
.../aligned_spins_bbh_tides_on.prior | 18 +++++++++++++++
test/prior_files/aligned_spins_bns.prior | 16 ++++++++++++++
.../aligned_spins_bns_tides_on.prior | 18 +++++++++++++++
.../aligned_spins_waveform_tides_on.prior | 14 ------------
.../prior_files/precessing_spins_bbh.prior | 12 +++++-----
.../precessing_spins_bbh_tides_on.prior | 22 +++++++++++++++++++
test/prior_files/precessing_spins_bns.prior | 20 +++++++++++++++++
.../precessing_spins_bns_tides_on.prior | 22 +++++++++++++++++++
test/prior_test.py | 8 +++----
24 files changed, 274 insertions(+), 100 deletions(-)
rename bilby/gw/prior_files/{aligned_spins_binary_black_holes.prior => aligned_spins_bbh.prior} (100%)
create mode 100644 bilby/gw/prior_files/aligned_spins_bbh_tides_on.prior
create mode 100644 bilby/gw/prior_files/aligned_spins_bns.prior
create mode 100644 bilby/gw/prior_files/aligned_spins_bns_tides_on.prior
delete mode 100644 bilby/gw/prior_files/aligned_spins_waveform_tides_off.prior
delete mode 100644 bilby/gw/prior_files/aligned_spins_waveform_tides_on.prior
rename test/prior_files/precessing_spins_binary_black_holes.prior => bilby/gw/prior_files/precessing_spins_bbh.prior (100%)
create mode 100644 bilby/gw/prior_files/precessing_spins_bbh_tides_on.prior
create mode 100644 bilby/gw/prior_files/precessing_spins_bns.prior
create mode 100644 bilby/gw/prior_files/precessing_spins_bns_tides_on.prior
delete mode 100644 bilby/gw/prior_files/precessing_spins_waveform_tides_off.prior
delete mode 100644 bilby/gw/prior_files/precessing_spins_waveform_tides_on.prior
create mode 100644 test/prior_files/aligned_spins_bbh.prior
create mode 100644 test/prior_files/aligned_spins_bbh_tides_on.prior
create mode 100644 test/prior_files/aligned_spins_bns.prior
create mode 100644 test/prior_files/aligned_spins_bns_tides_on.prior
delete mode 100644 test/prior_files/aligned_spins_waveform_tides_on.prior
rename bilby/gw/prior_files/precessing_spins_binary_black_holes.prior => test/prior_files/precessing_spins_bbh.prior (100%)
create mode 100644 test/prior_files/precessing_spins_bbh_tides_on.prior
create mode 100644 test/prior_files/precessing_spins_bns.prior
create mode 100644 test/prior_files/precessing_spins_bns_tides_on.prior
diff --git a/bilby/gw/prior.py b/bilby/gw/prior.py
index 25c97b847..e3523a36f 100644
--- a/bilby/gw/prior.py
+++ b/bilby/gw/prior.py
@@ -411,10 +411,10 @@ class BBHPriorDict(CBCPriorDict):
"""
if dictionary is None and filename is None:
if aligned_spin:
- fname = 'aligned_spins_binary_black_holes.prior'
+ fname = 'aligned_spins_bbh.prior'
logger.info('Using aligned spin prior')
else:
- fname = 'precessing_spins_binary_black_holes.prior'
+ fname = 'precessing_spins_bbh.prior'
filename = os.path.join(DEFAULT_PRIOR_DIR, fname)
logger.info('No prior given, using default BBH priors in {}.'.format(filename))
elif filename is not None:
@@ -517,9 +517,9 @@ class BNSPriorDict(CBCPriorDict):
BNSPriorDict.default_conversion_function
"""
if aligned_spin:
- default_file = 'aligned_spins_waveform_tides_on.prior'
+ default_file = 'aligned_spins_bns_tides_on.prior'
else:
- default_file = 'precessing_spins_waveform_tides_on.prior'
+ default_file = 'precessing_spins_bns_tides_on.prior'
if dictionary is None and filename is None:
filename = os.path.join(DEFAULT_PRIOR_DIR, default_file)
logger.info('No prior given, using default BNS priors in {}.'.format(filename))
diff --git a/bilby/gw/prior_files/aligned_spins_binary_black_holes.prior b/bilby/gw/prior_files/aligned_spins_bbh.prior
similarity index 100%
rename from bilby/gw/prior_files/aligned_spins_binary_black_holes.prior
rename to bilby/gw/prior_files/aligned_spins_bbh.prior
index 9dda5d5bb..4a15f7a3a 100644
--- a/bilby/gw/prior_files/aligned_spins_binary_black_holes.prior
+++ b/bilby/gw/prior_files/aligned_spins_bbh.prior
@@ -6,11 +6,11 @@ mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100)
-chi_1 = bilby.gw.prior.AlignedSpin(name='chi_1', a_prior=Uniform(minimum=0, maximum=0.99))
-chi_2 = bilby.gw.prior.AlignedSpin(name='chi_2', a_prior=Uniform(minimum=0, maximum=0.99))
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
theta_jn = Sine(name='theta_jn')
psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+chi_1 = bilby.gw.prior.AlignedSpin(name='chi_1', a_prior=Uniform(minimum=0, maximum=0.99))
+chi_2 = bilby.gw.prior.AlignedSpin(name='chi_2', a_prior=Uniform(minimum=0, maximum=0.99))
diff --git a/bilby/gw/prior_files/aligned_spins_bbh_tides_on.prior b/bilby/gw/prior_files/aligned_spins_bbh_tides_on.prior
new file mode 100644
index 000000000..4fc71953e
--- /dev/null
+++ b/bilby/gw/prior_files/aligned_spins_bbh_tides_on.prior
@@ -0,0 +1,18 @@
+# These are the default priors we use for BBH systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
+mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+chi_1 = bilby.gw.prior.AlignedSpin(name='chi_1', a_prior=Uniform(minimum=0, maximum=0.99))
+chi_2 = bilby.gw.prior.AlignedSpin(name='chi_2', a_prior=Uniform(minimum=0, maximum=0.99))
+lambda_1 = Uniform(name='lambda_1', minimum=0, maximum=5000)
+lambda_2 = Uniform(name='lambda_2', minimum=0, maximum=5000)
diff --git a/bilby/gw/prior_files/aligned_spins_bns.prior b/bilby/gw/prior_files/aligned_spins_bns.prior
new file mode 100644
index 000000000..a524de777
--- /dev/null
+++ b/bilby/gw/prior_files/aligned_spins_bns.prior
@@ -0,0 +1,16 @@
+# These are the default priors we use for BNS systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5)
+mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=0.4, maximum=4.4)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+chi_1 = bilby.gw.prior.AlignedSpin(name='chi_1', a_prior=Uniform(minimum=0, maximum=0.99))
+chi_2 = bilby.gw.prior.AlignedSpin(name='chi_2', a_prior=Uniform(minimum=0, maximum=0.99))
diff --git a/bilby/gw/prior_files/aligned_spins_bns_tides_on.prior b/bilby/gw/prior_files/aligned_spins_bns_tides_on.prior
new file mode 100644
index 000000000..dcb30189b
--- /dev/null
+++ b/bilby/gw/prior_files/aligned_spins_bns_tides_on.prior
@@ -0,0 +1,18 @@
+# These are the default priors we use for BNS systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5)
+mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=0.4, maximum=4.4)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+chi_1 = bilby.gw.prior.AlignedSpin(name='chi_1', a_prior=Uniform(minimum=0, maximum=0.99))
+chi_2 = bilby.gw.prior.AlignedSpin(name='chi_2', a_prior=Uniform(minimum=0, maximum=0.99))
+lambda_1 = Uniform(name='lambda_1', minimum=0, maximum=5000)
+lambda_2 = Uniform(name='lambda_2', minimum=0, maximum=5000)
diff --git a/bilby/gw/prior_files/aligned_spins_waveform_tides_off.prior b/bilby/gw/prior_files/aligned_spins_waveform_tides_off.prior
deleted file mode 100644
index b15bfa6c5..000000000
--- a/bilby/gw/prior_files/aligned_spins_waveform_tides_off.prior
+++ /dev/null
@@ -1,12 +0,0 @@
-chirp_mass = Uniform(name='chirp_mass', minimum=3.56, maximum=3.68)
-mass_ratio = Uniform(name='mass_ratio', minimum=0.05, maximum=1)
-mass_1 = Constraint(name='mass_1', minimum=1.00, maximum=22.0)
-mass_2 = Constraint(name='mass_2', minimum=1.00, maximum=2.95)
-chi_1 = bilby.gw.prior.AlignedSpin(name='chi_1', a_prior=Uniform(minimum=0, maximum=0.5))
-chi_2 = bilby.gw.prior.AlignedSpin(name='chi_2', a_prior=Uniform(minimum=0, maximum=0.05))
-luminosity_distance = bilby.core.prior.PowerLaw(alpha=2, name='luminosity_distance', minimum=1, maximum=750, unit='Mpc')
-ra = Uniform(minimum=0, maximum=2 * np.pi, name='ra', boundary='periodic')
-dec = Cosine(name='dec')
-theta_jn = Sine(name='theta_jn')
-psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
-phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
diff --git a/bilby/gw/prior_files/aligned_spins_waveform_tides_on.prior b/bilby/gw/prior_files/aligned_spins_waveform_tides_on.prior
deleted file mode 100644
index 6c318f1b3..000000000
--- a/bilby/gw/prior_files/aligned_spins_waveform_tides_on.prior
+++ /dev/null
@@ -1,14 +0,0 @@
-chirp_mass = Uniform(name='chirp_mass', minimum=3.56, maximum=3.68)
-mass_ratio = Uniform(name='mass_ratio', minimum=0.05, maximum=1)
-mass_1 = Constraint(name='mass_1', minimum=1.00, maximum=22.0)
-mass_2 = Constraint(name='mass_2', minimum=1.00, maximum=2.95)
-chi_1 = bilby.gw.prior.AlignedSpin(name='chi_1', a_prior=Uniform(minimum=0, maximum=0.5))
-chi_2 = bilby.gw.prior.AlignedSpin(name='chi_2', a_prior=Uniform(minimum=0, maximum=0.05))
-luminosity_distance = bilby.core.prior.PowerLaw(alpha=2, name='luminosity_distance', minimum=1, maximum=750, unit='Mpc')
-ra = Uniform(minimum=0, maximum=2 * np.pi, name='ra', boundary='periodic')
-dec = Cosine(name='dec')
-theta_jn = Sine(name='theta_jn')
-psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
-phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
-lambda_1 = Uniform(name='lambda_2', minimum=0, maximum=5000)
-lambda_2 = Uniform(name='lambda_2', minimum=0, maximum=5000)
diff --git a/test/prior_files/precessing_spins_binary_black_holes.prior b/bilby/gw/prior_files/precessing_spins_bbh.prior
similarity index 100%
rename from test/prior_files/precessing_spins_binary_black_holes.prior
rename to bilby/gw/prior_files/precessing_spins_bbh.prior
index 06549e63e..27094d689 100644
--- a/test/prior_files/precessing_spins_binary_black_holes.prior
+++ b/bilby/gw/prior_files/precessing_spins_bbh.prior
@@ -6,15 +6,15 @@ mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100)
-a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
-a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
-tilt_1 = Sine(name='tilt_1')
-tilt_2 = Sine(name='tilt_2')
-phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
-phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
theta_jn = Sine(name='theta_jn')
psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
+a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
+tilt_1 = Sine(name='tilt_1')
+tilt_2 = Sine(name='tilt_2')
+phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
+phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
diff --git a/bilby/gw/prior_files/precessing_spins_bbh_tides_on.prior b/bilby/gw/prior_files/precessing_spins_bbh_tides_on.prior
new file mode 100644
index 000000000..81da0482f
--- /dev/null
+++ b/bilby/gw/prior_files/precessing_spins_bbh_tides_on.prior
@@ -0,0 +1,22 @@
+# These are the default priors we use for BBH systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
+mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
+a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
+tilt_1 = Sine(name='tilt_1')
+tilt_2 = Sine(name='tilt_2')
+phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
+phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
+lambda_1 = Uniform(name='lambda_1', minimum=0, maximum=5000)
+lambda_2 = Uniform(name='lambda_2', minimum=0, maximum=5000)
diff --git a/bilby/gw/prior_files/precessing_spins_bns.prior b/bilby/gw/prior_files/precessing_spins_bns.prior
new file mode 100644
index 000000000..cb663dbbe
--- /dev/null
+++ b/bilby/gw/prior_files/precessing_spins_bns.prior
@@ -0,0 +1,20 @@
+# These are the default priors we use for BNS systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5)
+mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=0.4, maximum=4.4)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
+a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
+tilt_1 = Sine(name='tilt_1')
+tilt_2 = Sine(name='tilt_2')
+phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
+phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
diff --git a/bilby/gw/prior_files/precessing_spins_bns_tides_on.prior b/bilby/gw/prior_files/precessing_spins_bns_tides_on.prior
new file mode 100644
index 000000000..d0c7e875c
--- /dev/null
+++ b/bilby/gw/prior_files/precessing_spins_bns_tides_on.prior
@@ -0,0 +1,22 @@
+# These are the default priors we use for BNS systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5)
+mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=0.4, maximum=4.4)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
+a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
+tilt_1 = Sine(name='tilt_1')
+tilt_2 = Sine(name='tilt_2')
+phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
+phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
+lambda_1 = Uniform(name='lambda_1', minimum=0, maximum=5000)
+lambda_2 = Uniform(name='lambda_2', minimum=0, maximum=5000)
diff --git a/bilby/gw/prior_files/precessing_spins_waveform_tides_off.prior b/bilby/gw/prior_files/precessing_spins_waveform_tides_off.prior
deleted file mode 100644
index 3bc8ebab2..000000000
--- a/bilby/gw/prior_files/precessing_spins_waveform_tides_off.prior
+++ /dev/null
@@ -1,16 +0,0 @@
-chirp_mass = Uniform(name='chirp_mass', minimum=3.56, maximum=3.68)
-mass_ratio = Uniform(name='mass_ratio', minimum=0.05, maximum=1)
-mass_1 = Constraint(name='mass_1', minimum=1.00, maximum=22.0)
-mass_2 = Constraint(name='mass_2', minimum=1.00, maximum=22.0)
-a_1 = Uniform(name='a_1', minimum=0, maximum=0.89)
-a_2 = Uniform(name='a_2', minimum=0, maximum=0.05)
-tilt_1 = Sine(name='tilt_1')
-tilt_2 = Sine(name='tilt_2')
-phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
-phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
-luminosity_distance = bilby.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=1, maximum=750, unit='Mpc')
-ra = Uniform(minimum=0, maximum=2 * np.pi, name='ra', boundary='periodic')
-dec = Cosine(name='dec')
-theta_jn = Sine(name='theta_jn')
-psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
-phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
diff --git a/bilby/gw/prior_files/precessing_spins_waveform_tides_on.prior b/bilby/gw/prior_files/precessing_spins_waveform_tides_on.prior
deleted file mode 100644
index d82a4f173..000000000
--- a/bilby/gw/prior_files/precessing_spins_waveform_tides_on.prior
+++ /dev/null
@@ -1,18 +0,0 @@
-chirp_mass = Uniform(name='chirp_mass', minimum=3.56, maximum=3.68)
-mass_ratio = Uniform(name='mass_ratio', minimum=0.05, maximum=1)
-mass_1 = Constraint(name='mass_1', minimum=1.00, maximum=22.0)
-mass_2 = Constraint(name='mass_2', minimum=1.00, maximum=22.0)
-a_1 = Uniform(name='a_1', minimum=0, maximum=0.89)
-a_2 = Uniform(name='a_2', minimum=0, maximum=0.05)
-tilt_1 = Sine(name='tilt_1')
-tilt_2 = Sine(name='tilt_2')
-phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
-phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
-luminosity_distance = bilby.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=1, maximum=750, unit='Mpc')
-ra = Uniform(minimum=0, maximum=2 * np.pi, name='ra', boundary='periodic')
-dec = Cosine(name='dec')
-theta_jn = Sine(name='theta_jn')
-psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
-phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
-lambda_1 = Uniform(name='lambda_2', minimum=0, maximum=5000)
-lambda_2 = Uniform(name='lambda_2', minimum=0, maximum=5000)
diff --git a/test/gw_prior_test.py b/test/gw_prior_test.py
index 171977652..7e864bd65 100644
--- a/test/gw_prior_test.py
+++ b/test/gw_prior_test.py
@@ -26,7 +26,7 @@ class TestBBHPriorDict(unittest.TestCase):
)
self.filename = os.path.join(
os.path.dirname(os.path.realpath(__file__)),
- "prior_files/precessing_spins_binary_black_holes.prior",
+ "prior_files/precessing_spins_bbh.prior",
)
self.bbh_prior_dict = bilby.gw.prior.BBHPriorDict(filename=self.filename)
for key, value in self.bbh_prior_dict.items():
@@ -267,13 +267,13 @@ class TestPackagedPriors(unittest.TestCase):
""" Test that the prepackaged priors load """
def test_aligned(self):
- filename = "aligned_spins_binary_black_holes.prior"
+ filename = "aligned_spins_bbh.prior"
prior_dict = bilby.gw.prior.BBHPriorDict(filename=filename)
self.assertTrue("chi_1" in prior_dict)
self.assertTrue("chi_2" in prior_dict)
def test_binary_black_holes(self):
- filename = "precessing_spins_binary_black_holes.prior"
+ filename = "precessing_spins_bbh.prior"
prior_dict = bilby.gw.prior.BBHPriorDict(filename=filename)
self.assertTrue("a_1" in prior_dict)
@@ -298,7 +298,7 @@ class TestBNSPriorDict(unittest.TestCase):
)
self.filename = os.path.join(
os.path.dirname(os.path.realpath(__file__)),
- "prior_files/aligned_spins_waveform_tides_on.prior",
+ "prior_files/aligned_spins_bns_tides_on.prior",
)
self.bns_prior_dict = bilby.gw.prior.BNSPriorDict(filename=self.filename)
for key, value in self.bns_prior_dict.items():
diff --git a/test/prior_files/aligned_spins_bbh.prior b/test/prior_files/aligned_spins_bbh.prior
new file mode 100644
index 000000000..4a15f7a3a
--- /dev/null
+++ b/test/prior_files/aligned_spins_bbh.prior
@@ -0,0 +1,16 @@
+# These are the default priors we use for BBH systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
+mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+chi_1 = bilby.gw.prior.AlignedSpin(name='chi_1', a_prior=Uniform(minimum=0, maximum=0.99))
+chi_2 = bilby.gw.prior.AlignedSpin(name='chi_2', a_prior=Uniform(minimum=0, maximum=0.99))
diff --git a/test/prior_files/aligned_spins_bbh_tides_on.prior b/test/prior_files/aligned_spins_bbh_tides_on.prior
new file mode 100644
index 000000000..4fc71953e
--- /dev/null
+++ b/test/prior_files/aligned_spins_bbh_tides_on.prior
@@ -0,0 +1,18 @@
+# These are the default priors we use for BBH systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
+mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+chi_1 = bilby.gw.prior.AlignedSpin(name='chi_1', a_prior=Uniform(minimum=0, maximum=0.99))
+chi_2 = bilby.gw.prior.AlignedSpin(name='chi_2', a_prior=Uniform(minimum=0, maximum=0.99))
+lambda_1 = Uniform(name='lambda_1', minimum=0, maximum=5000)
+lambda_2 = Uniform(name='lambda_2', minimum=0, maximum=5000)
diff --git a/test/prior_files/aligned_spins_bns.prior b/test/prior_files/aligned_spins_bns.prior
new file mode 100644
index 000000000..a524de777
--- /dev/null
+++ b/test/prior_files/aligned_spins_bns.prior
@@ -0,0 +1,16 @@
+# These are the default priors we use for BNS systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5)
+mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=0.4, maximum=4.4)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+chi_1 = bilby.gw.prior.AlignedSpin(name='chi_1', a_prior=Uniform(minimum=0, maximum=0.99))
+chi_2 = bilby.gw.prior.AlignedSpin(name='chi_2', a_prior=Uniform(minimum=0, maximum=0.99))
diff --git a/test/prior_files/aligned_spins_bns_tides_on.prior b/test/prior_files/aligned_spins_bns_tides_on.prior
new file mode 100644
index 000000000..dcb30189b
--- /dev/null
+++ b/test/prior_files/aligned_spins_bns_tides_on.prior
@@ -0,0 +1,18 @@
+# These are the default priors we use for BNS systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5)
+mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=0.4, maximum=4.4)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+chi_1 = bilby.gw.prior.AlignedSpin(name='chi_1', a_prior=Uniform(minimum=0, maximum=0.99))
+chi_2 = bilby.gw.prior.AlignedSpin(name='chi_2', a_prior=Uniform(minimum=0, maximum=0.99))
+lambda_1 = Uniform(name='lambda_1', minimum=0, maximum=5000)
+lambda_2 = Uniform(name='lambda_2', minimum=0, maximum=5000)
diff --git a/test/prior_files/aligned_spins_waveform_tides_on.prior b/test/prior_files/aligned_spins_waveform_tides_on.prior
deleted file mode 100644
index 6c318f1b3..000000000
--- a/test/prior_files/aligned_spins_waveform_tides_on.prior
+++ /dev/null
@@ -1,14 +0,0 @@
-chirp_mass = Uniform(name='chirp_mass', minimum=3.56, maximum=3.68)
-mass_ratio = Uniform(name='mass_ratio', minimum=0.05, maximum=1)
-mass_1 = Constraint(name='mass_1', minimum=1.00, maximum=22.0)
-mass_2 = Constraint(name='mass_2', minimum=1.00, maximum=2.95)
-chi_1 = bilby.gw.prior.AlignedSpin(name='chi_1', a_prior=Uniform(minimum=0, maximum=0.5))
-chi_2 = bilby.gw.prior.AlignedSpin(name='chi_2', a_prior=Uniform(minimum=0, maximum=0.05))
-luminosity_distance = bilby.core.prior.PowerLaw(alpha=2, name='luminosity_distance', minimum=1, maximum=750, unit='Mpc')
-ra = Uniform(minimum=0, maximum=2 * np.pi, name='ra', boundary='periodic')
-dec = Cosine(name='dec')
-theta_jn = Sine(name='theta_jn')
-psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
-phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
-lambda_1 = Uniform(name='lambda_2', minimum=0, maximum=5000)
-lambda_2 = Uniform(name='lambda_2', minimum=0, maximum=5000)
diff --git a/bilby/gw/prior_files/precessing_spins_binary_black_holes.prior b/test/prior_files/precessing_spins_bbh.prior
similarity index 100%
rename from bilby/gw/prior_files/precessing_spins_binary_black_holes.prior
rename to test/prior_files/precessing_spins_bbh.prior
index 06549e63e..27094d689 100644
--- a/bilby/gw/prior_files/precessing_spins_binary_black_holes.prior
+++ b/test/prior_files/precessing_spins_bbh.prior
@@ -6,15 +6,15 @@ mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100)
-a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
-a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
-tilt_1 = Sine(name='tilt_1')
-tilt_2 = Sine(name='tilt_2')
-phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
-phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
dec = Cosine(name='dec')
ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
theta_jn = Sine(name='theta_jn')
psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
+a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
+tilt_1 = Sine(name='tilt_1')
+tilt_2 = Sine(name='tilt_2')
+phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
+phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
diff --git a/test/prior_files/precessing_spins_bbh_tides_on.prior b/test/prior_files/precessing_spins_bbh_tides_on.prior
new file mode 100644
index 000000000..81da0482f
--- /dev/null
+++ b/test/prior_files/precessing_spins_bbh_tides_on.prior
@@ -0,0 +1,22 @@
+# These are the default priors we use for BBH systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=5, maximum=100)
+mass_2 = Constraint(name='mass_2', minimum=5, maximum=100)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
+a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
+tilt_1 = Sine(name='tilt_1')
+tilt_2 = Sine(name='tilt_2')
+phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
+phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
+lambda_1 = Uniform(name='lambda_1', minimum=0, maximum=5000)
+lambda_2 = Uniform(name='lambda_2', minimum=0, maximum=5000)
diff --git a/test/prior_files/precessing_spins_bns.prior b/test/prior_files/precessing_spins_bns.prior
new file mode 100644
index 000000000..cb663dbbe
--- /dev/null
+++ b/test/prior_files/precessing_spins_bns.prior
@@ -0,0 +1,20 @@
+# These are the default priors we use for BNS systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5)
+mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=0.4, maximum=4.4)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
+a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
+tilt_1 = Sine(name='tilt_1')
+tilt_2 = Sine(name='tilt_2')
+phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
+phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
diff --git a/test/prior_files/precessing_spins_bns_tides_on.prior b/test/prior_files/precessing_spins_bns_tides_on.prior
new file mode 100644
index 000000000..d0c7e875c
--- /dev/null
+++ b/test/prior_files/precessing_spins_bns_tides_on.prior
@@ -0,0 +1,22 @@
+# These are the default priors we use for BNS systems.
+# Note that you may wish to use more specific mass and distance parameters.
+# These commands are all known to bilby.gw.prior.
+# Lines beginning "#" are ignored.
+mass_1 = Constraint(name='mass_1', minimum=0.5, maximum=5)
+mass_2 = Constraint(name='mass_2', minimum=0.5, maximum=5)
+mass_ratio = Uniform(name='mass_ratio', minimum=0.125, maximum=1)
+chirp_mass = Uniform(name='chirp_mass', minimum=0.4, maximum=4.4)
+luminosity_distance = bilby.gw.prior.UniformSourceFrame(name='luminosity_distance', minimum=1e2, maximum=5e3)
+dec = Cosine(name='dec')
+ra = Uniform(name='ra', minimum=0, maximum=2 * np.pi, boundary='periodic')
+theta_jn = Sine(name='theta_jn')
+psi = Uniform(name='psi', minimum=0, maximum=np.pi, boundary='periodic')
+phase = Uniform(name='phase', minimum=0, maximum=2 * np.pi, boundary='periodic')
+a_1 = Uniform(name='a_1', minimum=0, maximum=0.99)
+a_2 = Uniform(name='a_2', minimum=0, maximum=0.99)
+tilt_1 = Sine(name='tilt_1')
+tilt_2 = Sine(name='tilt_2')
+phi_12 = Uniform(name='phi_12', minimum=0, maximum=2 * np.pi, boundary='periodic')
+phi_jl = Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi, boundary='periodic')
+lambda_1 = Uniform(name='lambda_1', minimum=0, maximum=5000)
+lambda_2 = Uniform(name='lambda_2', minimum=0, maximum=5000)
diff --git a/test/prior_test.py b/test/prior_test.py
index a0c57c493..6222752d4 100644
--- a/test/prior_test.py
+++ b/test/prior_test.py
@@ -911,7 +911,7 @@ class TestPriorDict(unittest.TestCase):
self.prior_set_from_dict = bilby.core.prior.PriorDict(dictionary=self.priors)
self.default_prior_file = os.path.join(
os.path.dirname(os.path.realpath(__file__)),
- "prior_files/precessing_spins_binary_black_holes.prior",
+ "prior_files/precessing_spins_bbh.prior",
)
self.prior_set_from_file = bilby.core.prior.PriorDict(
filename=self.default_prior_file
@@ -1252,7 +1252,7 @@ class TestFillPrior(unittest.TestCase):
self.priors = bilby.core.prior.PriorDict(dictionary=self.priors)
self.default_prior_file = os.path.join(
os.path.dirname(os.path.realpath(__file__)),
- "prior_files/precessing_spins_binary_black_holes.prior",
+ "prior_files/precessing_spins_bbh.prior",
)
self.priors.fill_priors(self.likelihood, self.default_prior_file)
@@ -1292,7 +1292,7 @@ class TestCreateDefaultPrior(unittest.TestCase):
def test_bbh_params(self):
prior_file = os.path.join(
os.path.dirname(os.path.realpath(__file__)),
- "prior_files/precessing_spins_binary_black_holes.prior",
+ "prior_files/precessing_spins_bbh.prior",
)
prior_set = bilby.core.prior.PriorDict(filename=prior_file)
for prior in prior_set:
@@ -1306,7 +1306,7 @@ class TestCreateDefaultPrior(unittest.TestCase):
def test_unknown_prior(self):
prior_file = os.path.join(
os.path.dirname(os.path.realpath(__file__)),
- "prior_files/precessing_spins_binary_black_holes.prior",
+ "prior_files/precessing_spins_bbh.prior",
)
self.assertIsNone(
bilby.core.prior.create_default_prior(
--
GitLab