d_T = np.linspace(-4, 4, n_T)
t_m = np.zeros(n_T)
t_a = np.zeros(n_T)
sdr = np.zeros(n_T)
fsrs = np.zeros(n_T)
opds = np.zeros(n_T)
theta_t_0 = theta_t
for i in range(n_T):
n1 = h1.generate_n_glass(h1.t + d_T[i])
d1 = h1.d_glass_thermal_expansion(h1.t + d_T[i])
print('glass', 2 * (h1.d_glass * h1.n_glass - n1 * d1) / lam)
n2 = h1.generate_n_air(h1.t + d_T[i])
d2 = h1.d_air_thermal_expansion(h1.t + d_T[i])
print('air', 2 * (h1.n_air * h1.d_air - n2 * d2) / lam)
opd = h1.opd_exact_pure(theta_t, n1, d1, n2, d2)
opds[i] = opd
fsrs[i] = h1.fsr(opd)
t_m[i] = h1.overall_transmittance(theta_d, f, h1.gamma_m, h1.fsr(opd))
t_a[i] = h1.overall_transmittance(theta_d, f, h1.gamma_a, h1.fsr(opd))
sdr[i] = t_m[i] / t_a[i]
fig, ax = plt.subplots()
# ax.plot(d_T, t_a, color='black')
ax.plot(d_T, sdr, color='black')
# ax.plot(d_T, (opds - fopd) / lam, color='black')
# ax.plot(d_d1 * 1000, fsrs, color='black')
# ax.set_ylim([0, 400])
ax.grid(True)
ax.set_xlabel(r"$\Delta\theta_t$ (degrees)")
ax.set_ylabel(r"SDR")
switch_lam = 2
lam = 532e-9 * switch_lam
theta_t = 1.5 * np.pi / 180
gamma_m = 1.40e9 / switch_lam # molecular signal spectral width
gamma_a = 50e6 / switch_lam # aerosol signal spectral width
fopd = 0.15 * switch_lam
t_ref = 20
t = 20
p = 1
f = 0.1
h1 = Copper(fopd, theta_t, gamma_m, gamma_a, lam, t, t_ref, p, d_opd_d_t=lam / 5, glass=LITHOSIL_Q)
# print(h1.d1, h1.d2, h1.d3)
# dt = np.linspace(-0.5, 0.5, 25)
dt = np.linspace(-5, 5, 25)
opd = h1.opd_exact_pure(theta_t,
h1.generate_n_glass(t + dt),
h1.d_glass_thermal_expansion(t + dt),
h1.generate_n_air(t + dt),
h1.d_air_thermal_expansion(t + dt))
opd_variation = (opd - fopd) / lam
fig, ax = plt.subplots()
ax.plot(dt, opd_variation, color='black')
ax.grid(True)
ax.set_xlabel(r"$\Delta$T ($^\circ$C)")
ax.set_ylabel(r"$\Delta$OPD ($\lambda$)")
ax.set_title("OPD Variation v. Temperature")
ax.text(0.21, 0.00001, r"T($\Delta$T=0) = 20 $^\circ$C")
plt.show()
p,
d_opd_d_t=lam / 5,
glass=LITHOSIL_Q)
n_theta = 400
theta_d = 0.002
d_theta_t = np.linspace(-6, 6, n_theta) * np.pi / 180
t_m = np.zeros(n_theta)
t_a = np.zeros(n_theta)
sdr = np.zeros(n_theta)
fsrs = np.zeros(n_theta)
theta_t_0 = theta_t
for i in range(n_theta):
h1.theta_t = theta_t_0 + d_theta_t[i]
opd = h1.opd_exact_pure(theta_t, h1.n_glass, h1.d_glass, h1.n_air,
h1.d_air)
fsrs[i] = h1.fsr(opd)
t_m[i] = h1.overall_transmittance(theta_d, f, h1.gamma_m, h1.fsr(opd))
t_a[i] = h1.overall_transmittance(theta_d, f, h1.gamma_a, h1.fsr(opd))
sdr[i] = t_m[i] / t_a[i]
fig, ax = plt.subplots()
ax.plot(d_theta_t * 1000, sdr, color='black')
# ax.plot(d_d1 * 1000, fsrs, color='black')
# ax.set_ylim([0, 400])
ax.grid(True)
ax.set_xlabel(r"$\Delta\theta_t$ (mrad)")
ax.set_ylabel(r"SDR")
ax.set_title("SDR v. Tilt Angle Variation")
ax.text(
30, 110, r"$\theta_t$($\Delta\theta_t$=0) = {0} mrad".format(