Update JSD in VIatmin authored by Wei Changfeng's avatar Wei Changfeng
......@@ -114,21 +114,21 @@ Therefore, I think it may work if we perform more this kind of tests. We have 10
## 3.3 More JSD calculation
Then I used 100 detector noises for JSD calculation. I calculated JSD for three parameters: m1, m2, and t0.
For each parameter, there are three kinds of JSD.
> (1) BBH reference, which is JSD between two BBHs with different detector noise.
> (2) BH encounter & BBH reference, which is JSD between BH encounter and BBH with the same detector noise.
> (3) BH encounter & injected BBH, which is JSD between BH encounter, and BBH injected by peak values of BH encounter's posterior, with the same detector noise.
> JSD(1): BBH reference, which is JSD between two BBHs with different detector noise.
> JSD(2): BH encounter & BBH reference, which is JSD between BH encounter and BBH with the same detector noise.
> JSD(3): BH encounter & injected BBH, which is JSD between BH encounter, and BBH injected by peak values of BH encounter's posterior, with the same detector noise.
Fig 3.3.1 JSD of m1: BH encounter & BBH reference
Fig 3.3.1 JSD(1) adn JSD(2) of m1
![m1_be_ref](uploads/818b862afbcf767a9b1ee8ec0ccebb66/m1_be_ref.png)
Fig 3.3.2 JSD of m1: BH encounter & injected BBH
Fig 3.3.2 JSD(1) adn JSD(3) of m1
![m1_be_inj](uploads/f9a51278beda36fd751aa653c4de95b4/m1_be_inj.png)
Fig 3.3.3 JSD of m2: BH encounter & BBH reference
Fig 3.3.3 JSD(1) adn JSD(2) of m2
![m2_be_ref](uploads/e51de0c0a22e377322cf05ed75fd5b94/m2_be_ref.png)
Fig 3.3.4 JSD of m2: BH encounter & injected BBH
Fig 3.3.4 JSD(1) adn JSD(3) of m2
![m2_be_inj](uploads/4ae8006cd7094dd7fcbfb8952d3e8ebc/m2_be_inj.png)
Fig 3.3.5 JSD of t0: BH encounter & BBH reference
Fig 3.3.5 JSD(1) adn JSD(2) of t0
![t0_be_ref](uploads/27faea38277d25115bde2dfff25f36cf/t0_be_ref.png)
Fig 3.3.6 JSD of t0: BH encounter & injected BBH
Fig 3.3.6 JSD(1) adn JSD(3) of t0
![t0_be_inj](uploads/360556eec165bde83e76acc75851a42e/t0_be_inj.png)
Here is the script which shows all JSD calculation. More Plots are included, such as comparison between JSD(1),(2),(3) for a single waveform, which can be seen from line 63.
......
......