home.md

@@ -250,32 +250,71 @@ This is due to different normalization items.
 ![signal_freq_Weichangfeng_method](uploads/8e27887501d8aba9a1631c05342d6a3f/signal_freq_Weichangfeng_method.png)    
 ![signal_freq_Hunter_method](uploads/c5b668283727bdc9f544f81ea23465ee/signal_freq_Hunter_method.png)    
 
-## (12)Hunter's test data   
-![Hunter_s_test_data](uploads/4b5ab36ac71da7c423f9806519d925bf/Hunter_s_test_data.png)   
-![vitamin_corner_timeseries-1](uploads/613fbc3a8058783bf9c65c672449fbd1/vitamin_corner_timeseries-1.png)   
+
+## (13)VItamin: remake Hunter's test data  
+### (i) How to remake:   
+In this test we get frequency domain noise-free signals at detector directly.   
+What I did: To make my data exactly the same to Hunter's test data, I delete 'adding white Gaussian noise' step, which doesn't show up in Hunter's code.(The method of making test data is written in ‘gen_benchmark_pe.py’) And if I add this step, my data will look so different from Hunter's.   
+    
+Fig.1  Hunter's test data (whitened):  
+![01](uploads/9e822df0dd23fc5a8e40fe2bec2cc5d4/01.png)   
+Fig.2  my data (whitened):   
+![01_remake](uploads/e4f9fe82dfa6446fa8a28a7df775bb01/01_remake.png)   
+Fig.3  my data if added white Gaussian noise:  
+![01_nosiy](uploads/951819818736bdd1e33c3796bdbdd443/01_nosiy.png)   
+### (ii) What my original data looks like (before whitening and adding Gaussian white nosie):  
+Fig.4  The signal's plot in time domain:    
+![bilby_bbh_signal](uploads/34ca57b1a23a37b876deff13324a6dc8/bilby_bbh_signal.png)  
+I fail to generate a time-domain waveform whose amplitude comes to zero quickly after merging. But, I made waveform the same to Hunter's way and the whitened waveforms do look similar.   
+In Hunter's test data, the amplitude also doesn't come to zero after merging(as you can see in Fig. 1), and there is a vibration, although it is not obvious. I think whitening step reduces it.   
+
+### (iii) Comparison of their performances in VItamin:   
+Except the first one, all the corner plots from my data are very similar to Hunter's.    
+(a)   
+Fig.4 & Fig.5    
+![0](uploads/cdac31b163fa6022ac91f1d5e74f405b/0.png)   
+![0_remake_2](uploads/6aa4486b90e3f667e117e856aa531cca/0_remake_2.png)   
+> mass_1=72.6     
+> mass_2=39.7     
+> luminosity_distance=1818    
+> geocent_time=0.23+1126259642.5   
+> theta_jn=1.59     
+> ra=2.26      
+> dec=0.22    
+(b)   
+Fig.6 & Fig.7    
+![1](uploads/3e94af076a7a04f8fe633cc081c998b8/1.png)    
+![1_remake_2](uploads/a1de775dafd6b1864b1d1752158c4754/1_remake_2.png)   
 > mass_1=79.5     
 > mass_2=59.7     
 > luminosity_distance=1327      
-> geocent_time=0.16     
+> geocent_time=0.16+1126259642.5     
 > theta_jn=0.69     
 > ra=2.79    
 > dec=-0.05    
+(c)   
+Fig.8 & Fig.9    
+![2](uploads/d889d5e82185914b6d02aed30a323fc5/2.png)   
+![2_remake_2](uploads/7d9982185bdf8101d293c4d35c6aec36/2_remake_2.png)   
+> mass_1=70.3     
+> mass_2=63.6     
+> luminosity_distance=1882    
+> geocent_time=0.16+1126259642.5     
+> theta_jn=2.54     
+> ra=1.97     
+> dec=1.06       
+(d)  
+Fig.10 & Fig.11    
+![3](uploads/949727cfc1fc44a842a443f4186c5b40/3.png)   
+![3_remake_2](uploads/74e88b9fc572edd10fec40294344f124/3_remake_2.png)   
+> mass_1=67.8     
+> mass_2=50.8     
+> luminosity_distance=1493    
+> geocent_time=0.29+1126259642.5     
+> theta_jn=1.06     
+> ra=4.06    
+> dec=-0.17    
 
-## (13)VItamin: Bilby BBH waveform new test    
-In this test we get frequency domain noise-free signals at detector directly.     
-The signal's plot in time domain:    
-![bilby_bbh_signal](uploads/276a14a48f4c6800411974f55fcf3c69/bilby_bbh_signal.png)    
-The whitened noisy signal's plot in time domain:   
-![bilby_whitened_noisy_signal](uploads/187555468c782ac04c6380bbe313912f/bilby_whitened_noisy_signal.png)   
-I believe Hunter already wrote the option for truth line function. But the way to activate it just isn't showed yet.   
-![vitamin_corner_timeseries-0](uploads/edf32cc566eff4ff900fe269badc97ba/vitamin_corner_timeseries-0.png)    
-> mass_1=57.5     
-> mass_2=57.5     
-> luminosity_distance=1800    
-> geocent_time=0.25     
-> theta_jn=2.62      
-> ra=3.14    
-> dec=0.0    
 
 
 # 2.Data