Update home authored by Wei Changfeng's avatar Wei Changfeng
...@@ -78,7 +78,8 @@ In the plot, h_cross looks like a straight line, but it s correct. This is becau ...@@ -78,7 +78,8 @@ In the plot, h_cross looks like a straight line, but it s correct. This is becau
## (5)Add white noise ## (5)Add white noise
## (6)VItamin ## (6)VItamin: Quickstart test
This step is to make sure that Vitamin works in the right way.
Run VItamin with test waveform provided by VItamin documentation: Run VItamin with test waveform provided by VItamin documentation:
https://hagabbar.github.io/vitamin_b/quickstart.html https://hagabbar.github.io/vitamin_b/quickstart.html
...@@ -185,6 +186,26 @@ This one almost has a wonderful estimation. Just dec is estimated a bit higher. ...@@ -185,6 +186,26 @@ This one almost has a wonderful estimation. Just dec is estimated a bit higher.
Basically, we can see that: Basically, we can see that:
mass_1 is estimated a bit lower than real value. mass_1 is estimated a bit lower than real value.
distance is much higher. distance is much higher.
## (7)VItamin: bilby test
This step is to test VItamin with BBH waveform generated by bilby.
The injection parameters are all set as median value of prior's range:
> mass_1=57.5
> mass_2=57.5
> luminosity_distance=2000
> geocent_time=0.25
> theta_jn=3.141592653589793/2
> ra=3.141592653589793
> dec=0.
Waveform is injected in 3 detectors(L1, H1 and V1), here we plot the one in L1:
![L1_BBH_bilby](uploads/a6499b451955b2afff2560646806c463/L1_BBH_bilby.png)
Then whiten it:
![whitened_L1_BBH_bilby](uploads/d89ee0ed55ed37a750f04907c81a8cb9/whitened_L1_BBH_bilby.png)
Input the waveform (in 3 detectors) to VItamin, and we get an estimation:
![vitamin_corner_timeseries-0](uploads/278f3e2351bae412a7a63644ccb09f10/vitamin_corner_timeseries-0.png)
The posterior doesn't look good, especially theta_jn.
What's the problem, is it an error from waveform's amplitude ? Because we can see that this nosiy data's amplitude (in L1) is 10^-22, the whitened one is 10^1.
# 2.Data # 2.Data
## (1)Korean merger wavform ## (1)Korean merger wavform
... ...
......