Commit b377f9a3 authored by Marcus Haberland's avatar Marcus Haberland
Browse files

included mode=='m' support for comparison(mode,weight) to play around with

parent 69560b36
......@@ -51,7 +51,11 @@ def comparison(mode='m', weight=1e-4):
F_fin = []
Tam_fin = []
for i, Pi in enumerate(P):
vals = np.array([bins[i]['pars']['K'],Pi*yr,1.])
if mode == 's':
vals = np.array([bins[i]['pars']['K'],Pi*yr,1.])
if mode == 'm':
b = bins[i]['binary']
vals = np.array([b.K,b.P,1.,1.,1.,1.,1.,np.log(b.a0),b.f_1])
for j, Pj in enumerate(P_ig):
if isclose(Pi,Pj):
P_fin.append(Pi)
......@@ -66,7 +70,7 @@ def comparison(mode='m', weight=1e-4):
plt.figure(dpi=300)
plt.tight_layout()
plt.title(r'LISA x IceGiant for $S_n^{LISA}(10 mHZ)/S_n^{IG}(10 mHz) = 10^{-4}$')
plt.title(r'LISA x IceGiant for $S_n^{LISA}(10 mHZ)/S_n^{IG}(10 mHz) = $'+'{:.1e}'.format(weight))
plt.loglog(P_fin,Tam_fin[:,0],'rx',label=r'$\sigma_K/K$')
plt.loglog(P_fin,Tam_fin[:,1],'bx',label=r'$\sigma_P/P$')
......@@ -81,5 +85,5 @@ def comparison(mode='m', weight=1e-4):
plt.savefig(fig+'rel_uncertainty_added.png')
comparison(weight=0,mode='s')
comparison(weight=1e-4,mode='m')
plt.show()
\ No newline at end of file
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment