update prior files to include units

test/prior_files/binary_black_holes.prior

@@ -2,10 +2,10 @@
 # Note that you may wish to use more specific mass and distance parameters.
 # These commands are all known to tupak.gw.prior.
 # Lines beginning "#" are ignored.
-mass_1 = Uniform(name='mass_1', minimum=5, maximum=100)
-mass_2 = Uniform(name='mass_2', minimum=5, maximum=100)
-# chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100)
-# total_mass =  Uniform(name='total_mass', minimum=10, maximum=200)
+mass_1 = Uniform(name='mass_1', minimum=5, maximum=100, unit='$M_{\\odot}$')
+mass_2 = Uniform(name='mass_2', minimum=5, maximum=100, unit='$M_{\\odot}$')
+# chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100, unit='$M_{\\odot}$')
+# total_mass =  Uniform(name='total_mass', minimum=10, maximum=200, unit='$M_{\\odot}$')
 # mass_ratio =  Uniform(name='mass_ratio', minimum=0.125, maximum=1)
 # symmetric_mass_ratio =  Uniform(name='symmetric_mass_ratio', minimum=8 / 81, maximum=0.25)
 a_1 =  Uniform(name='a_1', minimum=0, maximum=0.8)
@@ -16,7 +16,7 @@ tilt_2 =  Sine(name='tilt_2')
 # cos_tilt_2 =  Uniform(name='cos_tilt_2', minimum=-1, maximum=1)
 phi_12 =  Uniform(name='phi_12', minimum=0, maximum=2 * np.pi)
 phi_jl =  Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi)
-luminosity_distance =  tupak.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=1e2, maximum=5e3)
+luminosity_distance =  tupak.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=1e2, maximum=5e3, unit='Mpc')
 dec =  Cosine(name='dec')
 ra =  Uniform(name='ra', minimum=0, maximum=2 * np.pi)
 iota =  Sine(name='iota')

tupak/gw/prior_files/GW150914.prior

 # These are the default priors for analysing GW150914.
-mass_1 = Uniform(name='mass_1', minimum=30, maximum=50)
-mass_2 = Uniform(name='mass_2', minimum=20, maximum=40)
+mass_1 = Uniform(name='mass_1', minimum=30, maximum=50, unit='$M_{\\odot}$')
+mass_2 = Uniform(name='mass_2', minimum=20, maximum=40, unit='$M_{\\odot}$')
 a_1 =  Uniform(name='a_1', minimum=0, maximum=0.8)
 a_2 =  Uniform(name='a_2', minimum=0, maximum=0.8)
 tilt_1 =  Sine(name='tilt_1')
 tilt_2 =  Sine(name='tilt_2')
 phi_12 =  Uniform(name='phi_12', minimum=0, maximum=2 * np.pi)
 phi_jl =  Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi)
-luminosity_distance =  tupak.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=1e2, maximum=1e3)
+luminosity_distance =  tupak.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=1e2, maximum=1e3, unit='Mpc')
 dec =  Cosine(name='dec')
 ra =  Uniform(name='ra', minimum=0, maximum=2 * np.pi)
 iota =  Sine(name='iota')
 psi =  Uniform(name='psi', minimum=0, maximum=np.pi)
 phase =  Uniform(name='phase', minimum=0, maximum=2 * np.pi)
-geocent_time = Uniform(1126259462.322, 1126259462.522, name='geocent_time')
+geocent_time = Uniform(1126259462.322, 1126259462.522, name='geocent_time', unit='$s$')
 # These are the calibration parameters as described in
 # https://journals.aps.org/prx/abstract/10.1103/PhysRevX.6.041015
 # recalib_H1_frequency_0 = 20

tupak/gw/prior_files/binary_black_holes.prior

@@ -2,10 +2,10 @@
 # Note that you may wish to use more specific mass and distance parameters.
 # These commands are all known to tupak.gw.prior.
 # Lines beginning "#" are ignored.
-mass_1 = Uniform(name='mass_1', minimum=5, maximum=100)
-mass_2 = Uniform(name='mass_2', minimum=5, maximum=100)
-# chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100)
-# total_mass =  Uniform(name='total_mass', minimum=10, maximum=200)
+mass_1 = Uniform(name='mass_1', minimum=5, maximum=100, unit='$M_{\\odot}$')
+mass_2 = Uniform(name='mass_2', minimum=5, maximum=100, unit='$M_{\\odot}$')
+# chirp_mass = Uniform(name='chirp_mass', minimum=25, maximum=100, unit='$M_{\\odot}$')
+# total_mass =  Uniform(name='total_mass', minimum=10, maximum=200, unit='$M_{\\odot}$')
 # mass_ratio =  Uniform(name='mass_ratio', minimum=0.125, maximum=1)
 # symmetric_mass_ratio =  Uniform(name='symmetric_mass_ratio', minimum=8 / 81, maximum=0.25)
 a_1 =  Uniform(name='a_1', minimum=0, maximum=0.8)
@@ -16,7 +16,7 @@ tilt_2 =  Sine(name='tilt_2')
 # cos_tilt_2 =  Uniform(name='cos_tilt_2', minimum=-1, maximum=1)
 phi_12 =  Uniform(name='phi_12', minimum=0, maximum=2 * np.pi)
 phi_jl =  Uniform(name='phi_jl', minimum=0, maximum=2 * np.pi)
-luminosity_distance =  tupak.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=1e2, maximum=5e3)
+luminosity_distance =  tupak.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=1e2, maximum=5e3, unit='Mpc')
 dec =  Cosine(name='dec')
 ra =  Uniform(name='ra', minimum=0, maximum=2 * np.pi)
 iota =  Sine(name='iota')

tupak/gw/prior_files/binary_neutron_stars.prior

@@ -2,15 +2,15 @@
 # Note that you may wish to use more specific mass and distance parameters.
 # These commands are all known to tupak.gw.prior.
 # Lines beginning "#" are ignored.
-mass_1 = Uniform(name='mass_1', minimum=1, maximum=2)
-mass_2 = Uniform(name='mass_2', minimum=1, maximum=2)
-# chirp_mass = Uniform(name='chirp_mass', minimum=0.87, maximum=1.74)
-# total_mass =  Uniform(name='total_mass', minimum=2, maximum=4)
+mass_1 = Uniform(name='mass_1', minimum=1, maximum=2, unit='$M_{\\odot}$')
+mass_2 = Uniform(name='mass_2', minimum=1, maximum=2, unit='$M_{\\odot}$')
+# chirp_mass = Uniform(name='chirp_mass', minimum=0.87, maximum=1.74, unit='$M_{\\odot}$')
+# total_mass =  Uniform(name='total_mass', minimum=2, maximum=4, unit='$M_{\\odot}$')
 # mass_ratio =  Uniform(name='mass_ratio', minimum=0.5, maximum=1)
 # symmetric_mass_ratio =  Uniform(name='symmetric_mass_ratio', minimum=0.22, maximum=0.25)
 a_1 =  Uniform(name='a_1', minimum= -0.05, maximum= 0.05)
 a_2 =  Uniform(name='a_2', minimum= -0.05, maximum= 0.05)
-luminosity_distance =  tupak.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=10, maximum=500)
+luminosity_distance =  tupak.gw.prior.UniformComovingVolume(name='luminosity_distance', minimum=10, maximum=500, unit='Mpc')
 dec =  Cosine(name='dec')
 ra =  Uniform(name='ra', minimum=0, maximum=2 * np.pi)
 iota =  Sine(name='iota')