VItamin-and-Bilby-with-Different-Noises.md

@@ -51,7 +51,7 @@ Fig.14  random seed = 10000
 ![BBH_03](uploads/aac24811cfec9b3141650aa726c66bca/BBH_03.png)   
 
 # 3. Bilby's tutorial    
-We did PE run provided by Bilby twice to see the difference between the results.    
+We did PE run provided by Bilby twice to see the consistency of the results.    
 ## (1) 4 parameters    
 ### (i) posterior     
 ![4](uploads/4149e24039801f12dfaa9d806691a5bf/4.png)    
@@ -66,7 +66,7 @@ JSD = 0.0462
 #### (c) distance     
 JSD = 0.0581     
 ![d](uploads/c3c3391f318b691c3214a5995f784ea8/d.png)     
-    
+It can be seen that the consistency of results is good.       
 ## (2) 15 parameters    
 ### (i) posterior    
 ![15cbc](uploads/a77d9f9030948fdd57cac91d6d52c0cf/15cbc.png)     
@@ -81,6 +81,48 @@ JSD = 0.0472
 #### (c) distance     
 JSD = 0.0359    
 ![d2](uploads/208d8cb343dbe5fd8d9a7f7630a1c3c8/d2.png)     
+It can be seen that the consistency of results is good.         
+Therefore, two tutorial scripts(4 and 15 parameters) can be accepted that they ran well on my laptop.      
+Then, in response to the suspicion that there may be errors in my previous scripts (I may add or delete some lines by mistake), I list the differences between the 15-parameter script and mine:          
+(i) Compared to the tutorial, my script misses the line:     
+'''     
+bilby.core.utils.setup_logger(outdir=outdir, label=label)     
+'''     
+However this doesn't matter to the results.     
+(ii) Priors are used differently. Priors used by tutorial are:(only different parts are shown)     
+'''     
+priors['geocent_time'] = bilby.core.prior.Uniform(
+    minimum=ref_geocent_time - 0.1,
+    maximum=ref_geocent_time + 0.1,
+    name='geocent_time', latex_label='$t_c$', unit='$s$') 
+'''     
+Priors used by me are:     
+'''       
+priors['luminosity_distance'] = bilby.core.prior.Uniform(name='luminosity_distance', minimum=1000, maximum=5000, unit='$Mpc_{\\odot}$')     
+priors['dec'] = bilby.core.prior.Cosine(name='dec')     
+priors['tilt_1'] = bilby.core.prior.Sine(name='tilt_1')      
+priors['tilt_2'] = bilby.core.prior.Sine(name='tilt_2')     
+priors['theta_jn'] = bilby.core.prior.Sine(name='theta_jn')     
+priors['geocent_time'] = bilby.core.prior.Uniform(
+    minimum=ref_geocent_time - 1,
+    maximum=ref_geocent_time + 1,
+    name='geocent_time', latex_label='$t_c$', unit='$s$')     
+'''        
+(iii) Result function.    
+Tutorial:    
+'''     
+result = bilby.run_sampler(
+    likelihood=likelihood, priors=priors, sampler='dynesty', npoints=2500,
+    outdir=outdir,
+    label=label, maxmcmc=5000,
+    conversion_function=bilby.gw.conversion.generate_all_bbh_parameters)     
+'''    
+Mine:    
+'''     
+result = bilby.run_sampler(
+    likelihood=likelihood, priors=priors, sampler='dynesty', npoints=1000,
+     outdir=outdir, label=label)     
+'''     
 ## (3) 15 parameters with (m1,m2) instead of (chirp mass, mass ratio)    
 We get results with good consistency in the sections above. This is to test if (m1,m2) will decrease the convergence.    
 ### (i) posterior    
@@ -96,7 +138,10 @@ JSD = 0.0459
 #### (a) distance     
 JSD = 0.0336       
 ![15d](uploads/0e97da3389404c5d4a673e4fb1423e26/15d.png)     
-## (4) BH enoucnter with marginalization(distance, time, phase)     
+It can be seen that the consistency of results is good.         
+    
+# BH encounter    
+## (1) BH enoucnter with marginalization(distance, time, phase)     
 ### (i) posterior     
 ![encounter](uploads/3681c0bbd9a90c7f5c9684e5ec5a891f/encounter.png)    
 ![encounter_cp](uploads/c6fdc02274cdb24f7556690946761aea/encounter_cp.png)     
@@ -110,7 +155,7 @@ JSD = 0.1759
 #### (a) distance     
 JSD = 0.1096    
 ![encounter_margil_d](uploads/81863c4571955f98fe5aa3f7d5112dba/encounter_margil_d.png)      
-## (5) BH enoucnter without marginalization(distance, time, phase)     
+## (2) BH enoucnter without marginalization(distance, time, phase)     
 ### (i) posterior     
 ![encounter_without_margil](uploads/155ebf8256d4bdb4f3aeeab20a974c9b/encounter_without_margil.png)    
 ![encounter_without_margil_cp](uploads/b5b138c6f0426eff581a95ed1201de1e/encounter_without_margil_cp.png)