# definition of optical system
s = OpticalSystem()
deg = math.pi/180.
lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="stop", decz=0.0), refname=s.rootcoordinatesystem.name)
lccomprism = s.addLocalCoordinateSystem(LocalCoordinates(name="prismcenter", decz=50.0), refname=lc0.name)
lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1", decz=-10.0, tiltx=30.*deg), refname=lccomprism.name) # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2", decz=10.0, tiltx=-30.*deg), refname=lccomprism.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=50.0), refname=lccomprism.name)
stopsurf = Surface(lc0)
frontsurf = Surface(lc1, shape=surfShape.Conic(lc1, curv=0), apert=CircularAperture(lc1, 20.0))
rearsurf = Surface(lc2, shape=surfShape.Conic(lc2, curv=0), apert=CircularAperture(lc2, 20.0))
image = Surface(lc3)
elem = OpticalElement(lc0, name="prism")
glass = ModelGlass(lc1)
elem.addMaterial("glass", glass)
elem.addSurface("stop", stopsurf, (None, None))
elem.addSurface("surf1", frontsurf, (None, "glass"))
elem.addSurface("surf2", rearsurf, ("glass", None))
elem.addSurface("image", image, (None, None))
lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="stop", decz=0.0),
refname=s.rootcoordinatesystem.name)
lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1", decz=-1.048),
refname=lc0.name) # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2", decz=4.0),
refname=lc1.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf3", decz=2.5),
refname=lc2.name)
lc4 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=97.2),
refname=lc3.name)
stopsurf = Surface(lc0, name="stopsurf")
frontsurf = Surface(lc1,
name="frontsurf",
shape=surfShape.Conic(lc1, curv=1. / 62.8, name='conic1'),
apert=CircularAperture(lc1, 12.7))
cementsurf = Surface(lc2,
name="cementsurf",
shape=surfShape.Conic(lc2, curv=-1. / 45.7,
name='conic2'),
apert=CircularAperture(lc2, 12.7))
rearsurf = Surface(lc3,
name="rearsurf",
shape=surfShape.Conic(lc3, curv=-1. / 128.2, name='conic3'),
apert=CircularAperture(lc3, 12.7))
image = Surface(lc4, name="imagesurf")
elem = OpticalElement(lc0, name="thorlabs_AC_254-100-A")
rnd_data1 = np.random.random((3, 3)) #np.eye(3)
lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="obj", decz=0.0),
refname=s.rootcoordinatesystem.name)
lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1",
decz=10.0,
tiltx=5. * math.pi / 180.0),
refname=lc0.name) # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2",
decz=20.0,
tiltx=10. * math.pi / 180.0),
refname=lc1.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=10.0),
refname=lc2.name)
stopsurf = Surface(lc0)
surf1 = Surface(lc1,
shape=surfShape.Conic(lc1, curv=1. / 24.0),
apert=CircularAperture(lc1, 5.0))
surf2 = Surface(lc2,
shape=surfShape.Conic(lc2, curv=-1. / 24.0),
apert=CircularAperture(lc2, 5.0))
image = Surface(lc3)
elem = OpticalElement(lc0, name="grinelement")
grin_strength = 0.5
def nfunc(x):
return grin_strength * np.exp(
-x[0]**2 - 4. * x[1]**2) + 1.0 #(2.5 - (x**2 + 100.0*y**4)/10.**2)
decy=-35),
refname=lc4.name)
lc5ap = s.addLocalCoordinateSystem(LocalCoordinates(name="oaparaap",
decz=0,
decy=35),
refname=lc5.name)
lc6 = s.addLocalCoordinateSystem(LocalCoordinates(name="image2",
decz=55,
tiltx=1 * math.pi / 32),
refname=lc5.name)
lc7 = s.addLocalCoordinateSystem(LocalCoordinates(name="image3", decz=5),
refname=lc6.name)
objectsurf = Surface(lc0)
m1surf = Surface(lc1,
shape=surfShape.Conic(lc1, curv=-0.01),
apert=CircularAperture(lc1, 20.))
m2surf = Surface(lc2,
shape=surfShape.Conic(lc2, curv=0.01),
apert=CircularAperture(lc2, 12.7))
m3surf = Surface(lc3,
shape=surfShape.Conic(lc3, curv=-0.006),
apert=CircularAperture(lc3, 12.7))
image1 = Surface(lc4)
oapara = Surface(lc3,
shape=surfShape.Conic(lc5, curv=0.01, cc=-1.),
apert=CircularAperture(lc5ap, 30.0))
image2 = Surface(lc6, apert=CircularAperture(lc6, 20.0))
image3 = Surface(lc7, apert=CircularAperture(lc7, 20.0))
elem = OpticalElement(lc0, name="TMA")
# definition of optical system
s = OpticalSystem()
lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="stop", decz=1.0),
refname=s.rootcoordinatesystem.name)
lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1", decz=10.0),
refname=lc0.name) # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2", decz=5.0),
refname=lc1.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=10.0),
refname=lc2.name)
stopsurf = Surface(lc0)
frontsurf = Surface(lc1,
shape=surfShape.Conic(lc1, curv=0),
apert=CircularAperture(lc1, 10.0))
rearsurf = Surface(lc2,
shape=surfShape.Conic(lc2, curv=0),
apert=CircularAperture(lc3, 10.0))
image = Surface(lc3)
elem = OpticalElement(lc0, name="crystalelem")
no = 1.5
neo = 1.8
myeps = np.array([[no, 0, 0], [0, no, 0], [0, 0, neo]])
crystal = AnisotropicMaterial(lc1, myeps)
logging.basicConfig(level=logging.DEBUG)
wavelength = 0.5876e-3
# definition of optical system
s = OpticalSystem()
lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="stop", decz=0.0), refname=s.rootcoordinatesystem.name)
lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1", decz=-1.048), refname=lc0.name) # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2", decz=4.0), refname=lc1.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf3", decz=2.5), refname=lc2.name)
lc4 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=97.2), refname=lc3.name)
stopsurf = Surface(lc0)
frontsurf = Surface(lc1, shape=surfShape.Conic(lc1, curv=1./62.8), apert=CircularAperture(lc1, 12.7))
cementsurf = Surface(lc2, shape=surfShape.Conic(lc2, curv=-1./45.7), apert=CircularAperture(lc2, 12.7))
rearsurf = Surface(lc3, shape=surfShape.Conic(lc3, curv=-1./128.2), apert=CircularAperture(lc3, 12.7))
image = Surface(lc4)
elem = OpticalElement(lc0, name="thorlabs_AC_254-100-A")
bk7 = ConstantIndexGlass(lc1, n=1.5168)
sf5 = ConstantIndexGlass(lc2, n=1.6727)
elem.addMaterial("BK7", bk7)
elem.addMaterial("SF5", sf5)
elem.addSurface("stop", stopsurf, (None, None))
elem.addSurface("front", frontsurf, (None, "BK7"))
decz=5.0,
tiltx=2.5 * degree),
refname=lc2.name)
lc4 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf4", decz=3.0),
refname=lc3.name)
lc5 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf5", decz=3.0),
refname=lc4.name)
lc6 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf6", decz=2.0),
refname=lc5.name)
lc7 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf7", decz=3.0),
refname=lc6.name)
lc8 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=19.0),
refname=lc7.name)
objectsurf = Surface(lc0)
surf1 = Surface(lc1, shape=surfShape.Conic(lc1, curv=1 / -5.922))
surf2 = Surface(lc2, shape=surfShape.Conic(lc2, curv=1 / -3.160))
surf3 = Surface(lc3, shape=surfShape.Conic(lc3, curv=1 / 15.884))
surf4 = Surface(lc4, shape=surfShape.Conic(lc4, curv=1 / -12.756))
stopsurf = Surface(lc5)
surf6 = Surface(lc6, shape=surfShape.Conic(lc6, curv=1 / 3.125))
surf7 = Surface(lc7, shape=surfShape.Conic(lc7, curv=0.1 * 1 / 1.479))
image = Surface(lc8)
elem = OpticalElement(lc0, name="lenssystem")
glass = ConstantIndexGlass(lc0, n=1.7)
glass2 = ConstantIndexGlass(lc0, n=1.5)
elem.addMaterial("glass", glass)
elem.addMaterial("glass2", glass2)
s = OpticalSystem()
lc0 = s.addLocalCoordinateSystem(LocalCoordinates(name="stop", decz=0.0),
refname=s.rootcoordinatesystem.name)
lc1 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf1", decz=-1.048),
refname=lc0.name) # objectDist
lc2 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf2", decz=4.0),
refname=lc1.name)
lc3 = s.addLocalCoordinateSystem(LocalCoordinates(name="surf3", decz=2.5),
refname=lc2.name)
lc4 = s.addLocalCoordinateSystem(LocalCoordinates(name="image", decz=97.2),
refname=lc3.name)
stopsurf = Surface(lc0)
frontsurf = Surface(lc1,
shape=surfShape.Conic(lc1, curv=1. / 62.8),
apert=CircularAperture(lc1, 12.7))
cementsurf = Surface(lc2,
shape=surfShape.Conic(lc2, curv=-1. / 45.7),
apert=CircularAperture(lc2, 12.7))
rearsurf = Surface(lc3,
shape=surfShape.Conic(lc3, curv=-1. / 128.2),
apert=CircularAperture(lc3, 12.7))
image = Surface(lc4)
elem = OpticalElement(lc0, name="thorlabs_AC_254-100-A")
bk7 = ConstantIndexGlass(lc1, n=1.5168)
sf5 = ConstantIndexGlass(lc2, n=1.6727)
elem.addMaterial("BK7", bk7)
decz=-100,
decy=-35),
refname=lc4.name)
lc5ap = s.addLocalCoordinateSystem(LocalCoordinates(name="oaparaap",
decz=0,
decy=35),
refname=lc5.name)
lc6 = s.addLocalCoordinateSystem(LocalCoordinates(name="image2",
decz=52.8,
tiltx=1 * math.pi / 32),
refname=lc5.name)
lc7 = s.addLocalCoordinateSystem(LocalCoordinates(name="image3", decz=5),
refname=lc6.name)
objectsurf = Surface(lc0)
m1surf = Surface(lc1, shape=surfShape.Conic(lc1, curv=-0.01))
m2surf = Surface(lc2, shape=surfShape.Conic(lc2, curv=0.01))
m3surf = Surface(lc3, shape=surfShape.Conic(lc3, curv=-0.006))
image1 = Surface(lc4)
oapara = Surface(lc3, shape=surfShape.Conic(lc5, curv=0.01, cc=-1.))
image2 = Surface(lc6, apert=CircularAperture(lc6, 20.0))
image3 = Surface(lc7, apert=CircularAperture(lc7, 20.0))
elem = OpticalElement(lc0, name="TMA")
elem.addMaterial("air", air)
elem.addSurface("object", objectsurf, (None, None))
elem.addSurface("m1", m1surf, (None, None))
elem.addSurface("m2", m2surf, (None, None))
elem.addSurface("m3", m3surf, (None, None))
