F.lens(fr2, ) # relay:2 lens:3
F.forvard(fr1 + fr2)
F.lens(fr1, ) # relay:2 lens:1
F.forvard(fr1) # MOT (2nd pass)
#F16= propagate( fslm, 0.1*cm, F16, 0, 0, 4)
imshow((F.value * F.value.conj()).real)
if False:
##### With the knife edge between the SLM (at SLM focus) relay:1,lens:1.
F = Field(N, side_length, wavelength) # beginning field
F.gaussian_aperture(gaussian_size)
F.value[:] *= lg # SLM LG phase
F.lens(fslm) # 'SLM' lens
F.forvard(fslm)
F.rectangular_screen(50, 100, 25 + 100 * um, 0) # knife edge
F.forvard(fslmR - fslm)
F.lens(fslmR) # relay:1 lens:1
F.forvard(2 * fslmR)
F.lens(fslmR) # relay:1 lens:2
F.forvard(fslmR + fslm) # MOT (1st pass)
ishow((F.value * F.value.conj()).real)
fprintf(stderr, 'press enter to continue')
raw_input()
F.forvard(fr1)
F.lens(fr1) # relay:2 lens:1
F.forvard(fr1 + fr2)
F.lens(fr2) # relay:2 lens:2
F.forvard(2 * fr2)
#!/usr/bin/env python
"""
LightPipes for Python Optical Toolbox
One Hole Diffraction.
"""
import common
from lightpipes import Field
from pylab import *
m=1
nm=1e-9*m
mm=1e-3*m
cm=1e-2*m
wavelength = 550*nm;
size = 5*mm
N = 100
R = 1*mm
z = 25*cm
F = Field(N,size,wavelength)
F.circular_aperture(R)
F.rectangular_screen(R*2, R/8, 0, 0, -45)
F.circular_screen(R/4)
F.fresnel(z)
imshow( abs(F.value)**2 ) # plot intensity
title('Intensity Distribution in the somewhat far-field');
show()
