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
## (6)VItamin
Run VItamin with test waveform provided by VItamin documentation:
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 :
![data_0](uploads/e9a7fde77706c97e6937f44c4fe669d4/data_0.png)
![data_1](uploads/d2383c7b84556e8d2ccbfd8d3bae3fb6/data_1.png)
......@@ -91,39 +90,54 @@ Plot them :
There are 4 corner plots due to 4 test data files.
![vitamin_corner_timeseries-0](uploads/c0b57a767538fae585290da8c8ba44d9/vitamin_corner_timeseries-0.png)
mass_1=72.5678967
mass_2=39.7158247
luminosity_distance=1818.14522
geocent_time=0.228755910
theta_jn=4.67871413
ra=1.59060244
dec=1.13254678
> mass_1=72.5678967
> mass_2=39.7158247
> luminosity_distance=1818.14522
> geocent_time=0.228755910
> theta_jn=1.59060244
> ra=2.25761658
> 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)
mass_1=
mass_2=
luminosity_distance=
geocent_time=
theta_jn=
ra=
dec=
> mass_1=79.5055181
> mass_2=59.7295127
> luminosity_distance=1326.64890
> geocent_time=0.159704400
> theta_jn=0.694929922
> ra=2.79268421
> 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)
mass_1=
mass_2=
luminosity_distance=
geocent_time=
theta_jn=
ra=
dec=
> mass_1=70.2724558
> mass_2=63.5675168
> luminosity_distance=1881.73441
> geocent_time=0.158573090
> theta_jn=2.54178566
> ra=1.96712040
> 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)
mass_1=
mass_2=
luminosity_distance=
geocent_time=
theta_jn=
ra=
dec=
> mass_1=67.7930585
> mass_2=50.8160513
> luminosity_distance=1493.29876
> geocent_time=0.291206202
> theta_jn=1.06098206
> ra=4.05618476
> dec=-0.171642340
Basically, we can see that:
mass_1 is estimated a bit lower than real value.
distance is much higher.
# 2.Data
## (1)Korean merger wavform
......@@ -136,10 +150,12 @@ dec=
## (2)Hunter’s BBH waveform generated from Bilby.
"data_0.h5py"
> 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).
> 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).
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
## (1)Adjust BBH waveform's amplitude the same to Korean waveform's
......
......