|
|
@@ -72,13 +72,23 @@ def analyze(data, dist, z_vec, wl_idx):
|
|
|
|
|
|
H1 = data_in_air[:,6,wl_idx]
|
|
|
H2 = data_in_gold[:,6,wl_idx]
|
|
|
+ E1 = data_in_air[:,4,wl_idx]
|
|
|
+ E2 = data_in_gold[:,4,wl_idx]
|
|
|
for i in range(len(z_vec)):
|
|
|
z = z_vec[i]*1e-9
|
|
|
print("z =",z)
|
|
|
H1_0 = H1[i]/np.exp(-kappa1[wl_idx]*z)
|
|
|
H2_0 = H2[i]/np.exp(-kappa2[wl_idx]*z)
|
|
|
- print("H0 air ",H1_0," from H1",H1[i])
|
|
|
- print("H0 gold",H2_0," from H2",H2[i])
|
|
|
+ E1_0 = E1[i]/np.exp(-kappa1[wl_idx]*z)
|
|
|
+ E2_0 = E2[i]/np.exp(-kappa2[wl_idx]*z)*eps2[wl_idx]
|
|
|
+ print("H0 air (%5.4g %+5.4gj)"%(np.real(H1_0), np.imag(H1_0)),
|
|
|
+ " from H1 (%5.4g %+5.4gj)"%(np.real(H1[i]), np.imag(H1[i])))
|
|
|
+ print("H0 gold (%5.4g %+5.4gj)"%(np.real(H2_0), np.imag(H2_0)),
|
|
|
+ " from H2 (%5.4g %+5.4gj)"%(np.real(H2[i]), np.imag(H2[i])))
|
|
|
+ print("E0 air (%5.4g %+5.4gj)"%(np.real(E1_0), np.imag(E1_0)),
|
|
|
+ " from E1 (%5.4g %+5.4gj)"%(np.real(E1[i]), np.imag(E1[i])))
|
|
|
+ print("E0*eps2 (%5.4g %+5.4gj)"%(np.real(E2_0), np.imag(E2_0)),
|
|
|
+ " from E2 (%5.4g %+5.4gj)"%(np.real(E2[i]), np.imag(E2[i])))
|
|
|
# H1_0 = H1/np.exp(-kappa1*
|
|
|
# print(H1[0], H2[0],H1[0]- H2[0])
|
|
|
# pl_data = (np.absolute(data_gold[:,2,wl_idx]*np.sqrt(dist)))
|
|
|
@@ -105,11 +115,11 @@ def main ():
|
|
|
|
|
|
#WLs=[300,350,400,450,600,700,800]
|
|
|
#WLs=[600,700, 800, 450]
|
|
|
- WLs=[800]#, 450]
|
|
|
+ WLs=[600]#, 450]
|
|
|
WLs_idx = get_WLs_idx(WLs, data)
|
|
|
|
|
|
|
|
|
- dist = 8 #mkm
|
|
|
+ dist = 10 #mkm
|
|
|
wl_idx = WLs_idx[0]
|
|
|
|
|
|
analyze(data, dist, z_vec, wl_idx)
|
