From 250343d76cf9c04978de58abed5ea682d240eb6e Mon Sep 17 00:00:00 2001 From: "Lucia F. de la Bella" Date: Thu, 29 Oct 2020 18:49:23 +0000 Subject: [PATCH 1/2] Add plot for three mergers' --- .../examples/merger_rates.yml | 37 +++++++++++++++++++ docs/gravitational_wave/index.rst | 31 ++++++++++++++++ 2 files changed, 68 insertions(+) create mode 100644 docs/gravitational_wave/examples/merger_rates.yml diff --git a/docs/gravitational_wave/examples/merger_rates.yml b/docs/gravitational_wave/examples/merger_rates.yml new file mode 100644 index 000000000..02ff2e5da --- /dev/null +++ b/docs/gravitational_wave/examples/merger_rates.yml @@ -0,0 +1,37 @@ +cosmology: !astropy.cosmology.default_cosmology.get +z_range: !numpy.arange [0, 2.01, 0.1] +mag_lim: 30 +sky_area: 1 deg2 +M_star_red: !astropy.modeling.models.Linear1D [-0.70798041, -20.37196157] +phi_star_red: !astropy.modeling.models.Exponential1D [0.0035097, -1.41649] +alpha_red: -0.5 +M_star_blue: !astropy.modeling.models.Linear1D [-0.9408582, -20.40492365] +phi_star_blue: !astropy.modeling.models.Exponential1D [0.00370253, -9.73858] +alpha_blue: -1.3 +redshift: !skypy.galaxy.redshift.schechter_lf_redshift + redshift: $z_range + M_star: $M_star_red + phi_star: $phi_star_red + alpha: $alpha_red + m_lim: $mag_lim + sky_area: $sky_area +magnitude: !skypy.galaxy.luminosity.schechter_lf_magnitude + redshift: $redshift + M_star: $M_star_red + alpha: $alpha_red + m_lim: $mag_lim +luminosity: !skypy.galaxy.luminosity.luminosity_from_absolute_magnitude + absolute_magnitude: $magnitude + zeropoint: Lsun_B +ns_ns_rate: !skypy.gravitational_wave.b_band_merger_rate + population: NS-NS + optimism: low + luminosity: $luminosity +ns_bh_rate: !skypy.gravitational_wave.b_band_merger_rate + population: NS-BH + optimism: low + luminosity: $luminosity +bh_bh_rate: !skypy.gravitational_wave.b_band_merger_rate + population: BH-BH + optimism: low + luminosity: $luminosity diff --git a/docs/gravitational_wave/index.rst b/docs/gravitational_wave/index.rst index 0bd0f0db1..11a6e85c5 100644 --- a/docs/gravitational_wave/index.rst +++ b/docs/gravitational_wave/index.rst @@ -2,4 +2,35 @@ Gravitational Waves (`skypy.gravitational_wave`) ************************************************ +SkyPy provides methods to model physical processes involving gravitational waves. +Here we demonstrate calculating the merger rates for three +different types of compact binary mergers. +For this we use `~skypy.gravitational_wave.b_band_merger_rate` and +the configuration file provided in the `~skypy.examples` directory: + + +.. plot:: + :include-source: true + + from skypy.pipeline import Pipeline + pipeline = Pipeline.read('examples/merger_rates.yml') + pipeline.execute() + + import matplotlib.pyplot as plt + + plt.figure() + plt.hist(pipeline['ns_ns_rate'], histtype='step', bins=np.logspace(-4,2,50), label='NS-NS') + plt.hist(pipeline['ns_bh_rate'], histtype='step', bins=np.logspace(-4,2,50), label='NS-BH') + plt.hist(pipeline['bh_bh_rate'], histtype='step', bins=np.logspace(-4,2,50), label='BH-BH') + plt.legend() + plt.xscale('log') + plt.yscale('log') + plt.xlabel('$Merger rate, R\,$[yr$^{-1}$]') + + plt.show() + + +‘NS-NS’ means neutron star - neutron star, ‘NS-BH’ is neutron star +- black hole and ‘BH-BH’ is black hole - black hole. + .. automodule:: skypy.gravitational_wave From b174b20f4231ede5c270d3c4b51fdb43040a18cc Mon Sep 17 00:00:00 2001 From: "Lucia F. de la Bella" Date: Tue, 3 Nov 2020 15:35:03 +0000 Subject: [PATCH 2/2] Adding table to config file and fix text --- .../examples/merger_rates.yml | 56 ++++++++++--------- docs/gravitational_wave/index.rst | 13 +++-- 2 files changed, 36 insertions(+), 33 deletions(-) diff --git a/docs/gravitational_wave/examples/merger_rates.yml b/docs/gravitational_wave/examples/merger_rates.yml index 02ff2e5da..d8d96cc75 100644 --- a/docs/gravitational_wave/examples/merger_rates.yml +++ b/docs/gravitational_wave/examples/merger_rates.yml @@ -8,30 +8,32 @@ alpha_red: -0.5 M_star_blue: !astropy.modeling.models.Linear1D [-0.9408582, -20.40492365] phi_star_blue: !astropy.modeling.models.Exponential1D [0.00370253, -9.73858] alpha_blue: -1.3 -redshift: !skypy.galaxy.redshift.schechter_lf_redshift - redshift: $z_range - M_star: $M_star_red - phi_star: $phi_star_red - alpha: $alpha_red - m_lim: $mag_lim - sky_area: $sky_area -magnitude: !skypy.galaxy.luminosity.schechter_lf_magnitude - redshift: $redshift - M_star: $M_star_red - alpha: $alpha_red - m_lim: $mag_lim -luminosity: !skypy.galaxy.luminosity.luminosity_from_absolute_magnitude - absolute_magnitude: $magnitude - zeropoint: Lsun_B -ns_ns_rate: !skypy.gravitational_wave.b_band_merger_rate - population: NS-NS - optimism: low - luminosity: $luminosity -ns_bh_rate: !skypy.gravitational_wave.b_band_merger_rate - population: NS-BH - optimism: low - luminosity: $luminosity -bh_bh_rate: !skypy.gravitational_wave.b_band_merger_rate - population: BH-BH - optimism: low - luminosity: $luminosity +tables: + merger_rates: + redshift: !skypy.galaxy.redshift.schechter_lf_redshift + redshift: $z_range + M_star: $M_star_red + phi_star: $phi_star_red + alpha: $alpha_red + m_lim: $mag_lim + sky_area: $sky_area + magnitude: !skypy.galaxy.luminosity.schechter_lf_magnitude + redshift: $merger_rates.redshift + M_star: $M_star_red + alpha: $alpha_red + m_lim: $mag_lim + luminosity: !skypy.galaxy.luminosity.luminosity_from_absolute_magnitude + absolute_magnitude: $merger_rates.magnitude + zeropoint: Lsun_B + ns_ns_rate: !skypy.gravitational_wave.b_band_merger_rate + population: NS-NS + optimism: low + luminosity: $merger_rates.luminosity + ns_bh_rate: !skypy.gravitational_wave.b_band_merger_rate + population: NS-BH + optimism: low + luminosity: $merger_rates.luminosity + bh_bh_rate: !skypy.gravitational_wave.b_band_merger_rate + population: BH-BH + optimism: low + luminosity: $merger_rates.luminosity diff --git a/docs/gravitational_wave/index.rst b/docs/gravitational_wave/index.rst index 11a6e85c5..6e0811a76 100644 --- a/docs/gravitational_wave/index.rst +++ b/docs/gravitational_wave/index.rst @@ -10,27 +10,28 @@ the configuration file provided in the `~skypy.examples` directory: .. plot:: - :include-source: true + :include-source: false from skypy.pipeline import Pipeline pipeline = Pipeline.read('examples/merger_rates.yml') pipeline.execute() + merger = pipeline['merger_rates'] import matplotlib.pyplot as plt plt.figure() - plt.hist(pipeline['ns_ns_rate'], histtype='step', bins=np.logspace(-4,2,50), label='NS-NS') - plt.hist(pipeline['ns_bh_rate'], histtype='step', bins=np.logspace(-4,2,50), label='NS-BH') - plt.hist(pipeline['bh_bh_rate'], histtype='step', bins=np.logspace(-4,2,50), label='BH-BH') + plt.hist(merger['ns_ns_rate'], histtype='step', bins=np.logspace(-4,2,50), label='NS-NS') + plt.hist(merger['ns_bh_rate'], histtype='step', bins=np.logspace(-4,2,50), label='NS-BH') + plt.hist(merger['bh_bh_rate'], histtype='step', bins=np.logspace(-4,2,50), label='BH-BH') plt.legend() plt.xscale('log') plt.yscale('log') - plt.xlabel('$Merger rate, R\,$[yr$^{-1}$]') + plt.xlabel('Merger rate, $R\, [year^{-1}$]') plt.show() -‘NS-NS’ means neutron star - neutron star, ‘NS-BH’ is neutron star +Note: ‘NS-NS’ means neutron star - neutron star, ‘NS-BH’ is neutron star - black hole and ‘BH-BH’ is black hole - black hole. .. automodule:: skypy.gravitational_wave