preperation for the meeting of 20/12/2021

LISA_x_IceGiant.py

@@ -42,7 +42,7 @@ infile.close()
 def isclose(a, b, rel_tol=1e-04, abs_tol=0.0):
     return abs(a-b) <= max(rel_tol * max(abs(a), abs(b)), abs_tol)
 
-def comparison(mode='s', weight=1e-4):
+def comparison(mode='m', weight=1e-4):
     bins = binaries[mode]
     ig_bins = binaries_ig[mode]
     P = np.array([b['pars']['P'] for b in bins])
@@ -81,5 +81,5 @@ def comparison(mode='s', weight=1e-4):
     
     plt.savefig(fig+'rel_uncertainty_added.png')
     
-comparison(weight=1e-4)
+comparison(weight=0,mode='s')
 plt.show()
\ No newline at end of file

uncertainty_plot_Tamanini_ig.py

@@ -31,11 +31,11 @@ from ig_binary import ig_binary
 
 keys = enumerate(['K','P','phi_0', 'theta_S_Ec', 'phi_S_Ec', 'theta_L', 'phi_L', 'f_0', 'ln(A)'])
 
-n0 = 10
+n0 = 15
 mP = 10
-a = 1
-b = -2
-mode = 'm'
+a = -2
+b = 1.3
+mode = 's'
 
 Ps = np.logspace(a,b,n0)
 uncs = np.zeros((n0,3),np.float64)
@@ -77,7 +77,7 @@ plt.loglog(Ps,uncs2[:,1],'b--')
 plt.loglog(Ps,uncs2[:,2],'g--')
 '''
 
-plt.title(r'Positions as in [1], $M_{b1,2}=0.23M_\odot,r=1kpc,M_P=1M_J, T_{obs}=4 yr$')
+plt.title(r'Positions as in [1], $M_{b1,2}=0.23M_\odot,r=1kpc,M_P=10 M_J$')
 plt.legend()
 plt.grid()
 plt.ylabel(r'$\sqrt{(\Gamma^{-1})_{ii}}/\lambda_i\cdot$SNR$\cdot M_P/M_J$')