Merge branch 'inconsistent_prior_files' into 'master'

inconsistent_prior_files

Closes #490

See merge request !799

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))

bilby/gw/prior_files/aligned_spins_binary_black_holes.prior → 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))

bilby/gw/prior_files/aligned_spins_waveform_tides_on.prior → bilby/gw/prior_files/aligned_spins_bbh_tides_on.prior

-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')
+# 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')
-lambda_1 = Uniform(name='lambda_2', minimum=0, maximum=5000)
+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)

bilby/gw/prior_files/aligned_spins_bns.prior

+# 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))

bilby/gw/prior_files/aligned_spins_bns_tides_on.prior

+# 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)

bilby/gw/prior_files/aligned_spins_waveform_tides_off.prior

-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')

test/prior_files/precessing_spins_binary_black_holes.prior → 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')

bilby/gw/prior_files/precessing_spins_waveform_tides_on.prior → bilby/gw/prior_files/precessing_spins_bbh_tides_on.prior

-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')
+# 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')
-lambda_1 = Uniform(name='lambda_2', minimum=0, maximum=5000)
+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)

bilby/gw/prior_files/precessing_spins_bns.prior

+# 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')

bilby/gw/prior_files/precessing_spins_bns_tides_on.prior

+# 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)

bilby/gw/prior_files/precessing_spins_waveform_tides_off.prior

-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')

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():

test/prior_files/aligned_spins_bbh.prior

+# 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))

test/prior_files/aligned_spins_bbh_tides_on.prior

+# 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)

test/prior_files/aligned_spins_bns.prior

+# 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))

test/prior_files/aligned_spins_bns_tides_on.prior

+# 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)

test/prior_files/aligned_spins_waveform_tides_on.prior

-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)

bilby/gw/prior_files/precessing_spins_binary_black_holes.prior → 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')

test/prior_files/precessing_spins_bbh_tides_on.prior

+# 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)

test/prior_files/precessing_spins_bns.prior

+# 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')