Update JSD in VIatmin authored by Wei Changfeng's avatar Wei Changfeng
# 1. Signals
The parameters of BBH signal were adjusted until it has a good SNR on each detector.
Then we injected a BH encounter with mass ratio of one, fixing its sky location and total mass, adjusting distance to make its amplitude similar with the BBH. Also it SNRs should look good.
## 1.1 BBH
> m1=78
> m2=72
......@@ -51,55 +54,15 @@ The operation below were also adopted for 100 BH encounters.
> For VItamin, it makes more sense to calculate the JSD between BBH and BH encounter. It is possible to create reference posteriors with BBH waveforms since we know they exist and are a good approximation of what is actually observed, but we do not have any clear observations of BH encounters. So if we have a set of BBH posteriors to use as reference, then we can compute the JSD with any new BBH observation to test whether it is a BBH or BH encoutner.
## 3.1 VItamin JSD calculation
### 3.1.1 test
Fig 3.1.1.1 m1 JSD
![m1](uploads/ed4bfde3e2b073bea5fea9ab0799ee7e/m1.png)
Fig 3.1.1.2 m1 JSD distribution
![m1_jsd](uploads/d45c1854bc87bfdfece2bf0316f12e93/m1_jsd.png)
Fig 3.1.1.3 m2 JSD
![m2](uploads/62a0d80864c908bb09e7fbcd575fa07c/m2.png)
Fig 3.1.1.4 m2 JSD distribution
![m2_jsd](uploads/0dd831fe4083a85af20e96e739e7c62c/m2_jsd.png)
Fig 3.1.1.5 distance JSD
![d](uploads/65f64ac7fa55bf06f3c2edae5b00087e/d.png)
Fig 3.1.1.6 distance JSD distribution
![d_jsd](uploads/7da57b16762157f2cd2a447e76ce6cec/d_jsd.png)
Fig 3.1.1.7 BBH & BH encounter t0 JSD
![jsd_t0](uploads/61236ed49b935e8019e57173ec36a942/jsd_t0.png)
Fig 3.1.1.8 BBH & BH encounter t0 JSD distribution
![pdf_t0](uploads/344c1bb00f3d53aef4526b53d72833fd/pdf_t0.png)
### 3.1.2 JSD of t0
Fig 3.1.2.1 t0 JSD distribution
## 3.1 VItamin JSD of t0 calculation
Fig 3.1.1 t0 JSD distribution
![jsd_ddd](uploads/62d5d54bfde18c592678073be115eb46/jsd_ddd.png)
We also injected BBH using BH encounter peak value given by VItamin. From four BH encounters(mass ratio = 1, 4, 8, 16), we got four BBHs. Then these BBHs' VItamin posteriors were calculated JSDs with corresponding BH encounters.
> JSD of t0 are 0.1732, 0.2894, 0.0110, and 0.0983, for mass ratio of 1, 4, 8 and 16 respectively.
They're blue lines on the Fig 3.1.1.
I think these four JSD values don't have an obvious distance from the BBH JSD distribution. Remember in the previous test, a big difference between BBH and BH encounter posterior is BH encounter has double peaks on t0. However, new injected BBH (using BH encounter peak values) also have double peaks on t0 expect m4. This makes me confusing.
Therefore, I think it may work if we perform more this kind of tests. We have 100 random noise, but only used one for each BH encounter waveform already.
## 3.2 Bilby and VItamin JSD calculation
Fig 3.2.1 BBH m1 JSD
![BBH_m1](uploads/edde8c260119e8ed794b1a07b7934008/BBH_m1.png)
Fig 3.2.2 BH encounter m1 JSD
![encounter_m1](uploads/649d4c834d0e7a67a711b69322a18d07/encounter_m1.png)
Fig 3.2.3 BBH m2 JSD
![BBH_m2](uploads/18e7d025e8867eca8ac785417892860e/BBH_m2.png)
Fig 3.2.4 BH encounter m2 JSD
![encounter_m2](uploads/d7b6364e26c14761562d5bba0db7867b/encounter_m2.png)
Fig 3.2.5 BBH distance JSD
![BBH_distance](uploads/69dadcfe7168e2bc3a9702c045f92634/BBH_distance.png)
Fig 3.2.6 BH encounter distance JSD
![encounter_d](uploads/83b397b55e53199afdc5f8ccb6329eaa/encounter_d.png)
Fig 3.2.7 BBH dec JSD
![BBH_dec](uploads/b8fa2bfb8435dabac8af2d376638a964/BBH_dec.png)
Fig 3.2.8 BH encounter dec JSD
![encounter_dec](uploads/a563525b073d88d7491b4d11657e43de/encounter_dec.png)
Fig 3.2.9 BBH ra JSD
![BBH_ra](uploads/05b33dda2b281108dff69e2dfe3ace4b/BBH_ra.png)
Fig 3.2.10 BH encounter ra JSD
![encounter_ra](uploads/e281f62d24444d05a9a045d4061e34ba/encounter_ra.png)
Fig 3.2.11 BBH t0 JSD
![BBH_t0](uploads/1742bdb1f9ce411ddde561abfc47c1d1/BBH_t0.png)
Fig 3.2.12 BH encounter t0 JSD
![encounter_t0](uploads/3f7aa9a4e948554f9226d15c2af921fc/encounter_t0.png)
Fig 3.2.13 Comparison
![comparison](uploads/65f956ab4b9c933b6835eaf8a208e91f/comparison.png)