Update VItamin and Bilby with Different Noises authored by Wei Changfeng's avatar Wei Changfeng
...@@ -54,6 +54,10 @@ Fig.14 random seed = 10000 ...@@ -54,6 +54,10 @@ Fig.14 random seed = 10000
We did PE run provided by Bilby twice to see the consistency of the results. We did PE run provided by Bilby twice to see the consistency of the results.
## (1) 4 parameters ## (1) 4 parameters
### (i) posterior ### (i) posterior
bbh_4parameters_marginalization.ipynb
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/blob/master/2021_05/bbh_4parameters_marginalization.ipynb
bbh_4parameters_marginalization_cp.ipynb
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/blob/master/2021_05/bbh_4parameters_marginalization_cp.ipynb
Fig.15 Fig.15
![4](uploads/4149e24039801f12dfaa9d806691a5bf/4.png) ![4](uploads/4149e24039801f12dfaa9d806691a5bf/4.png)
Fig.16 Fig.16
...@@ -73,6 +77,10 @@ Fig.19 ...@@ -73,6 +77,10 @@ Fig.19
![d](uploads/c3c3391f318b691c3214a5995f784ea8/d.png) ![d](uploads/c3c3391f318b691c3214a5995f784ea8/d.png)
It can be seen that the consistency of results is good. It can be seen that the consistency of results is good.
## (2) 15 parameters ## (2) 15 parameters
bbh_15parameters_marginalization.ipynb
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/blob/master/2021_05/bbh_15parameters_marginalization.ipynb
bbh_15parameters_marginalization_cp.ipynb
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/blob/master/2021_05/bbh_15parameters_marginalization_cp.ipynb
### (i) posterior ### (i) posterior
Fig.20 Fig.20
![15cbc](uploads/a77d9f9030948fdd57cac91d6d52c0cf/15cbc.png) ![15cbc](uploads/a77d9f9030948fdd57cac91d6d52c0cf/15cbc.png)
...@@ -161,6 +169,10 @@ result = bilby.run_sampler( ...@@ -161,6 +169,10 @@ result = bilby.run_sampler(
## (3) 15 parameters with (m1,m2) instead of (chirp mass, mass ratio) ## (3) 15 parameters with (m1,m2) instead of (chirp mass, mass ratio)
We get results with good consistency in the sections above. This is to test if (m1,m2) will decrease the convergence. We get results with good consistency in the sections above. This is to test if (m1,m2) will decrease the convergence.
bbh_15parameters_marginalization_m1m2.ipynb
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/blob/master/2021_05/bbh_15parameters_marginalization_m1m2.ipynb
bbh_15parameters_marginalization_m1m2_cp.ipynb
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/blob/master/2021_05/bbh_15parameters_marginalization_m1m2_cp.ipynb
### (i) posterior ### (i) posterior
Fig.25 Fig.25
![15cbc_n](uploads/3866e0c185dce78896dd1b18100c8a06/15cbc_n.png) ![15cbc_n](uploads/3866e0c185dce78896dd1b18100c8a06/15cbc_n.png)
...@@ -183,6 +195,10 @@ It can be seen that the consistency of results is still good. Hence we decide to ...@@ -183,6 +195,10 @@ It can be seen that the consistency of results is still good. Hence we decide to
# 4. BH encounter # 4. BH encounter
## (1) BH enoucnter with marginalization(distance, time, phase) ## (1) BH enoucnter with marginalization(distance, time, phase)
encounter_signal2_q1_marginalization.ipynb
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/blob/master/2021_05/encounter_signal2_q1_marginalization.ipynb
encounter_signal2_q1_marginalization_cp.ipynb
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/blob/master/2021_05/encounter_signal2_q1_marginalization_cp.ipynb
### (i) posterior ### (i) posterior
Fig.30 Fig.30
![encounter](uploads/3681c0bbd9a90c7f5c9684e5ec5a891f/encounter.png) ![encounter](uploads/3681c0bbd9a90c7f5c9684e5ec5a891f/encounter.png)
...@@ -202,7 +218,13 @@ JSD = 0.1096 ...@@ -202,7 +218,13 @@ JSD = 0.1096
Fig.34 Fig.34
![encounter_margil_d](uploads/81863c4571955f98fe5aa3f7d5112dba/encounter_margil_d.png) ![encounter_margil_d](uploads/81863c4571955f98fe5aa3f7d5112dba/encounter_margil_d.png)
Wrong template (using BBH template to model BH encounter) reduces the consistency of results. Wrong template (using BBH template to model BH encounter) reduces the consistency of results.
## (2) BH enoucnter without marginalization(distance, time, phase) ## (2) BH enoucnter without marginalization(distance, time, phase)
encounter_signal2_q1.ipynb
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/blob/master/2021_05/encounter_signal2_q1.ipynb
encounter_signal2_q1_cp.ipynb
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/blob/master/2021_05/encounter_signal2_q1_cp.ipynb
### (i) posterior ### (i) posterior
Fig.35 Fig.35
![encounter_without_margil](uploads/155ebf8256d4bdb4f3aeeab20a974c9b/encounter_without_margil.png) ![encounter_without_margil](uploads/155ebf8256d4bdb4f3aeeab20a974c9b/encounter_without_margil.png)
...@@ -230,11 +252,13 @@ https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/tree/ ...@@ -230,11 +252,13 @@ https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/tree/
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/tree/master/2021_05/test_2 https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/tree/master/2021_05/test_2
Though they produced the same number of samples this time, I checked the first ten lines of the sample files and they are not the same. Though they produced the same number of samples this time, I checked the first ten lines of the sample files and they are not the same.
### (i) test 1 ### (i) test 1
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/tree/master/2021_05/test_01
Fig.40 Fig.40
![test1](uploads/a0d28c07feda2275d7f2862d9ab32a8f/test1.png) ![test1](uploads/a0d28c07feda2275d7f2862d9ab32a8f/test1.png)
Fig.41 Fig.41
![corner1](uploads/5c4fb50dbd6d6a4ade2dfcd78f417ad9/corner1.png) ![corner1](uploads/5c4fb50dbd6d6a4ade2dfcd78f417ad9/corner1.png)
### (i) test 2 ### (i) test 2
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/tree/master/2021_05/test_02
Fig.42 Fig.42
![test2](uploads/db55b617fe1a2b00080c0c78f1fe3d27/test2.png) ![test2](uploads/db55b617fe1a2b00080c0c78f1fe3d27/test2.png)
Fig.43 Fig.43
...@@ -245,8 +269,10 @@ https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/tree/ ...@@ -245,8 +269,10 @@ https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/tree/
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/tree/master/2021_05/test_cpnest_2 https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/tree/master/2021_05/test_cpnest_2
They produced the different number of samples. I checked the first ten lines of the sample files and they are not the same. They produced the different number of samples. I checked the first ten lines of the sample files and they are not the same.
### (i) test 1 ### (i) test 1
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/tree/master/2021_05/test_cpnest_1
Fig.44 Fig.44
![corner_cpn_1](uploads/f1df89d85618e3a74432b06333aad372/corner_cpn_1.png) ![corner_cpn_1](uploads/f1df89d85618e3a74432b06333aad372/corner_cpn_1.png)
### (i) test 2 ### (i) test 2
https://gilsay.physics.gla.ac.uk/gitlab/wei.changfeng/weichangfeng_project/tree/master/2021_05/test_cpnest_2
Fig.45 Fig.45
![corner_cpn_2](uploads/81c603af077b79714bfd0dbbf8fcd384/corner_cpn_2.png) ![corner_cpn_2](uploads/81c603af077b79714bfd0dbbf8fcd384/corner_cpn_2.png)
\ No newline at end of file