make change sampled parameters a little more pathological

examples/injection_examples/change_sampled_parameters.py

@@ -21,8 +21,8 @@ outdir = 'outdir'
 np.random.seed(151226)
 
 injection_parameters = dict(
-    mass_1=36., mass_2=29., a_1=0.4, a_2=0.3, tilt_1=0.5, tilt_2=1.0,
-    phi_12=1.7, phi_jl=0.3, luminosity_distance=500, iota=0.4, psi=2.659,
+    total_mass=66., mass_ratio=0.9, a_1=0.4, a_2=0.3, tilt_1=0.5, tilt_2=1.0,
+    phi_12=1.7, phi_jl=0.3, luminosity_distance=2000, iota=0.4, psi=2.659,
     phase=1.3, geocent_time=1126259642.413, ra=1.375, dec=-1.2108)
 
 waveform_arguments = dict(waveform_approximant='IMRPhenomPv2',
@@ -34,8 +34,8 @@ waveform_generator = tupak.gw.waveform_generator.WaveformGenerator(
     frequency_domain_source_model=tupak.gw.source.lal_binary_black_hole,
     parameter_conversion=
         tupak.gw.conversion.convert_to_lal_binary_black_hole_parameters,
-    parameters=injection_parameters, waveform_arguments=waveform_arguments)
-hf_signal = waveform_generator.frequency_domain_strain()
+    waveform_arguments=waveform_arguments)
+# hf_signal = waveform_generator.frequency_domain_strain()
 
 # Set up interferometers.
 ifos = tupak.gw.detector.InterferometerList(['H1', 'L1', 'V1', 'K1'])
@@ -52,14 +52,16 @@ priors.pop('mass_2')
 priors.pop('luminosity_distance')
 priors['chirp_mass'] = tupak.prior.Uniform(
     name='chirp_mass', latex_label='$m_c$', minimum=13, maximum=45)
-priors['mass_ratio'] = tupak.prior.Uniform(
-    name='mass_ratio', latex_label='q', minimum=0.1, maximum=1)
+priors['symmetric_mass_ratio'] = tupak.prior.Uniform(
+    name='symmetric_mass_ratio', latex_label='q', minimum=0.1, maximum=0.25)
 priors['redshift'] = tupak.prior.Uniform(
     name='redshift', latex_label='$z$', minimum=0, maximum=0.5)
 # These parameters will not be sampled
 for key in ['a_1', 'a_2', 'tilt_1', 'tilt_2', 'phi_12', 'phi_jl', 'psi',
             'ra', 'dec', 'geocent_time', 'phase']:
     priors[key] = injection_parameters[key]
+priors.pop('iota')
+priors['cos_iota'] = np.cos(injection_parameters['iota'])
 print(priors)
 
 # Initialise GravitationalWaveTransient