Update home authored by Wei Changfeng's avatar Wei Changfeng
...@@ -81,8 +81,7 @@ In the plot, h_cross looks like a straight line, but it s correct. This is becau ...@@ -81,8 +81,7 @@ In the plot, h_cross looks like a straight line, but it s correct. This is becau
## (6)VItamin ## (6)VItamin
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
In the page above I can't tell the information of waveforms:
![test_waveform_and_pre-trained_model](uploads/1b016328f68833cec6ff50726bbe4b95/test_waveform_and_pre-trained_model.png)
Plot them : Plot them :
![data_0](uploads/e9a7fde77706c97e6937f44c4fe669d4/data_0.png) ![data_0](uploads/e9a7fde77706c97e6937f44c4fe669d4/data_0.png)
![data_1](uploads/d2383c7b84556e8d2ccbfd8d3bae3fb6/data_1.png) ![data_1](uploads/d2383c7b84556e8d2ccbfd8d3bae3fb6/data_1.png)
...@@ -91,39 +90,54 @@ Plot them : ...@@ -91,39 +90,54 @@ Plot them :
There are 4 corner plots due to 4 test data files. There are 4 corner plots due to 4 test data files.
![vitamin_corner_timeseries-0](uploads/c0b57a767538fae585290da8c8ba44d9/vitamin_corner_timeseries-0.png) ![vitamin_corner_timeseries-0](uploads/c0b57a767538fae585290da8c8ba44d9/vitamin_corner_timeseries-0.png)
mass_1=72.5678967 > mass_1=72.5678967
mass_2=39.7158247 > mass_2=39.7158247
luminosity_distance=1818.14522 > luminosity_distance=1818.14522
geocent_time=0.228755910 > geocent_time=0.228755910
theta_jn=4.67871413 > theta_jn=1.59060244
ra=1.59060244 > ra=2.25761658
dec=1.13254678 > dec=0.220253249
Basically, we can see that:
mass_1 is estimated much lower than real value.
distance is much higher.
ra has 2 peak, the left one is good, but the right one is higher.
dec can't be estimated clearly to the real value = 0.22
![vitamin_corner_timeseries-1](uploads/511671132510e8532e5b321690df9375/vitamin_corner_timeseries-1.png) ![vitamin_corner_timeseries-1](uploads/511671132510e8532e5b321690df9375/vitamin_corner_timeseries-1.png)
mass_1= > mass_1=79.5055181
mass_2= > mass_2=59.7295127
luminosity_distance= > luminosity_distance=1326.64890
geocent_time= > geocent_time=0.159704400
theta_jn= > theta_jn=0.694929922
ra= > ra=2.79268421
dec= > dec=-0.0544109042
Basically, we can see that:
mass_1 is estimated a bit lower than real value.
mass_2 is estimated a bit higher.
ra is estimated a bit lower.
![vitamin_corner_timeseries-2](uploads/d061df13078cb287369aeeb8f49b2d1e/vitamin_corner_timeseries-2.png) ![vitamin_corner_timeseries-2](uploads/d061df13078cb287369aeeb8f49b2d1e/vitamin_corner_timeseries-2.png)
mass_1= > mass_1=70.2724558
mass_2= > mass_2=63.5675168
luminosity_distance= > luminosity_distance=1881.73441
geocent_time= > geocent_time=0.158573090
theta_jn= > theta_jn=2.54178566
ra= > ra=1.96712040
dec= > dec=1.06435185
This one almost has a wonderful estimation. Just dec is estimated a bit higher.
![vitamin_corner_timeseries-3](uploads/a61d808b9ea479feef5a328ea2aae0f8/vitamin_corner_timeseries-3.png) ![vitamin_corner_timeseries-3](uploads/a61d808b9ea479feef5a328ea2aae0f8/vitamin_corner_timeseries-3.png)
mass_1= > mass_1=67.7930585
mass_2= > mass_2=50.8160513
luminosity_distance= > luminosity_distance=1493.29876
geocent_time= > geocent_time=0.291206202
theta_jn= > theta_jn=1.06098206
ra= > ra=4.05618476
dec= > dec=-0.171642340
Basically, we can see that:
mass_1 is estimated a bit lower than real value.
distance is much higher.
# 2.Data # 2.Data
## (1)Korean merger wavform ## (1)Korean merger wavform
...@@ -136,10 +150,12 @@ dec= ...@@ -136,10 +150,12 @@ dec=
## (2)Hunter’s BBH waveform generated from Bilby. ## (2)Hunter’s BBH waveform generated from Bilby.
"data_0.h5py" "data_0.h5py"
> Scaled waveforms (~10^-24-10^-23) > Scaled waveforms (~10^-24-10^-23)
> The hdf5 file has four directories: snrs, y_data_noisefree, y_data_noisy and rand_pars. > The hdf5 file has 4+ directories: snrs, y_data_noisefree, y_data_noisy and rand_pars.
> snrs == A list containing the optimal SNR of each waveform at each observatory (e.g. H1, L1 and V1 ... in that order). > snrs == A list containing the optimal SNR of each waveform at each observatory (e.g. H1, L1 and V1 ... in that order).
> y_data_noisy == The noisy waveform at each observatory with dimension (number of detectors, sampling frequency * duration). > y_data_noisy == The noisy waveform at each observatory with dimension (number of detectors, sampling frequency * duration).
> y_data_noisefree == The noisefree waveform at each observatory with dimension (number of detectors, sampling frequency * duration). > y_data_noisefree == The noisefree waveform at each observatory with dimension (number of detectors, sampling frequency * duration).
The parameter values are contained in the key 'x_data', where the values are in the following order:
['mass_1','mass_2','luminosity_distance','geocent_time','phase','theta_jn','psi','ra','dec']
# 3.Code # 3.Code
## (1)Adjust BBH waveform's amplitude the same to Korean waveform's ## (1)Adjust BBH waveform's amplitude the same to Korean waveform's
... ...
......