inconsistent_prior_files

Meant to close #490 (closed) . Basically just copy pasting the 2 existing binary_black_holes.prior files and add tides priors at the end. I started "fixing" @nicola.delillo 's file but then out of consistency it just seemed the more sensible thing to do.

Do we want a bunch of similar files for bns, nsbh etc or is that sufficient?

bilby/gw/prior_files/aligned_spins_binary_black_holes_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)
+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')
-lambda_1 = Uniform(name='lambda_2', minimum=0, maximum=5000)
+lambda_1 = Uniform(name='lambda_1', minimum=0, maximum=5000)
 lambda_2 = Uniform(name='lambda_2', minimum=0, maximum=5000)