Changes
Page history
Update home
authored
Sep 18, 2020
by
Wei Changfeng
Show whitespace changes
Inline
Side-by-side
home.md
View page @
f45c1c25
...
@@ -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:

Plot them :
Plot them :




...
@@ -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.


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


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.


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.


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
...
...
...
...