JSD-in-VIatmin.md

+  
+
+# 1. Signals     
+## 1.1 BBH        
+> m1=78     
+> m2=72   
+> distance=1400         
+> ra=0.89         
+> dec=-0.94         
+> psi=1.54            
+> t0=0.22    
+> theta_jn=1.51      
+    
+SNR：    
+> H1=10.0   
+> L1=12.9   
+> V1=10.5   
+    
+Fig 1.1.1    
+![BBH](uploads/d2afb339e996771d9906dba04701b497/BBH.png)    
+Fig 1.1.2    
+![noisy_BBH](uploads/ce7dc90c7827266cf33037da0a209318/noisy_BBH.png)    
+    
+## 1.2 BH encounter       
+waveform: "h_m1_L0.9_l2m2_r300.dat"   
+> total_mass=150     
+> distance=5000         
+> ra=0.89         
+> dec=-0.94         
+> psi=1.54            
+> t0=0.22     
+    
+SNR：    
+> H1=10.6   
+> L1=10.6   
+> V1=12.6   
+    
+Fig 1.2.1    
+![encounter](uploads/46309280f992e4772a36bbde70d57e72/encounter.png)   
+Fig 1.2.2    
+![noisy_encounter](uploads/55aa1400011cced00b426fb90661497f/noisy_encounter.png)    
+# 2. Posteriors.    
+(1) BBH          
+![vitamin_corner_timeseries-0](uploads/25ddcbc085573513e718098e6678bf6d/vitamin_corner_timeseries-0.png)    
+(2) BH encounter     
+![vitamin_corner_timeseries-0](uploads/12d79365d5b244dea61370e95912f608/vitamin_corner_timeseries-0.png)   
+ 
+# 3. J-S Divergence     
+100 BBHs with different random noise were input to VItamin. We calculated JSD for each two samples and got 50 mean JSDs(each of them are from 8 JSDs for 8 parameters).     
+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       
+Fig 3.1.1  m1 JSD   
+![m1](uploads/ed4bfde3e2b073bea5fea9ab0799ee7e/m1.png)   
+Fig 3.1.2  m1 JSD distribution     
+![m1_jsd](uploads/d45c1854bc87bfdfece2bf0316f12e93/m1_jsd.png)    
+Fig 3.1.3  m2 JSD   
+![m2](uploads/62a0d80864c908bb09e7fbcd575fa07c/m2.png)    
+Fig 3.1.4  m2 JSD distribution     
+![m2_jsd](uploads/0dd831fe4083a85af20e96e739e7c62c/m2_jsd.png)   
+Fig 3.1.5  distance JSD   
+![d](uploads/65f64ac7fa55bf06f3c2edae5b00087e/d.png)    
+Fig 3.1.6  distance JSD distribution     
+![d_jsd](uploads/7da57b16762157f2cd2a447e76ce6cec/d_jsd.png)   
+
+## 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)     
\ No newline at end of file