from epsilons import epsilon_Au
wl = 800e-9
n_air = 1.0
n_glass = 1.5
high_na = 0.8
low_na = 0.2
widths = N.arange(100, 501, 10) * 1e-9
k0 = 2 * N.pi / wl
fig1 = P.figure(figsize=(plot_config.pagewidth, 2.5))
# fig2 = P.figure(figsize=(plot_config.pagewidth, 2.5))
calc = SlitSystem(slit_widths=widths,
wavelength=wl,
metal_epsilon=epsilon_Au(wl),
metal_thickness=200e-9,
incidence_index=n_glass,
slit_index=n_air,
outcoupling_index=n_air,
collection_index=n_air,
numerical_aperture=low_na,
angle_of_incidence=0.0,
caching=False)
ratio1 = calc.collection_eff_TM / calc.collection_eff_TE
# absolute1 = calc.collection_eff_TE
calc = SlitSystem(slit_widths=widths,
wavelength=wl,
metal_epsilon=epsilon_Au(wl),
wl = 800e-9
n_air = 1.0
n_glass = 1.5
high_na = 0.8
low_na = 0.2
widths = N.arange(100, 501, 10) * 1e-9
k0 = 2 * N.pi / wl
fig1 = P.figure(figsize=(plot_config.pagewidth, 2.5))
# fig2 = P.figure(figsize=(plot_config.pagewidth, 2.5))
calc = SlitSystem(
slit_widths=widths,
wavelength=wl,
metal_epsilon=epsilon_Au(wl),
metal_thickness=200e-9,
incidence_index=n_glass,
slit_index=n_air,
outcoupling_index=n_air,
collection_index=n_air,
numerical_aperture=low_na,
angle_of_incidence=0.0,
caching=False)
ratio1 = calc.collection_eff_TM / calc.collection_eff_TE
# absolute1 = calc.collection_eff_TE
calc = SlitSystem(
slit_widths=widths,
wavelength=wl,
import plot_config
import epsilons
from matplotlib import pyplot as P
from matplotlib import gridspec, patches
from diffraction import calculation_figure
# Free parameters
wl = 830e-9 # light wavelength in free space in meters
width = 26e-6 # width of view in meters
eps2 = 1.0 # dielectric medium is air, n = 1
w0 = 2.2e-6 # width parameter of incident beam in meters
L = 50e-6 # length of slits in meters (also size of simulation)
res = 1e-7 # size of one pixel in meters
# Calculated parameters
eps1 = epsilons.epsilon_Au(wl)
k0 = 2 * N.pi / wl
kSP = k0 * N.sqrt(eps1 * eps2 / (eps1 + eps2))
gs = gridspec.GridSpec(1, 6, width_ratios=[1, 1, 1, 1, 1, 0.05])
fig = P.figure(figsize=(plot_config.pagewidth, 1.21))
calc_args = {'k': kSP, 'w0': w0, 'L': L, 'res': res, 'exact': False,
'normalize': True, 'find_zeroes': True}
axis_args = {'aspect': 'equal', 'xticks': [], 'yticks': []}
circle_args = {'facecolor': 'none', 'linestyle': 'dashed'}
ax = fig.add_subplot(gs[0], **axis_args)
p = calculation_figure(ax, 3, 50e-6, **calc_args)
# Manually put the circle in since the goddamn contour plot doesn't work
ax.add_patch(patches.Circle((0, 0), 4.0, **circle_args))
import plot_config
import epsilons
from matplotlib import pyplot as P
from matplotlib import gridspec
from diffraction import calculation_figure, measurement_figure
# Free parameters
wl = 830e-9 # light wavelength in free space in meters
width = 20e-6 # width of view in meters
eps2 = 1.0 # dielectric medium is air, n = 1
w0 = 2.2e-6 # width parameter of incident beam in meters
L = 50e-6 # length of slits in meters (also size of simulation)
res = 1e-7 # size of one pixel in meters
# Calculated parameters
eps1 = epsilons.epsilon_Au(wl)
k0 = 2 * N.pi / wl
kSP = k0 * N.sqrt(eps1 * eps2 / (eps1 + eps2))
gs = gridspec.GridSpec(1, 6, width_ratios=[1, 1, 1, 1, 1, 0.05])
fig = P.figure(figsize=(plot_config.pagewidth, 1.21))
axis_args = {'aspect': 'equal', 'xticks': [], 'yticks': []}
ax = fig.add_subplot(gs[0], **axis_args)
reflection = misc.imread('../../data/20100429 l=3.5 plate/far field.png')[
150:350, 115:315]
ax.imshow(reflection, origin='upper')
ax.subfigure_label('(a)', pos='lower_right')
ax = fig.add_subplot(gs[1], **axis_args)
