``````

In [ ]:

from pyoptools.all import *
from numpy import pi,sqrt

``````

## Crear un sistema sencillo con una lente

``````

In [ ]:

S=System(complist=[(L1,(0,0,100),(0,0,0))],n=1)
R=[Ray(pos=(0,0,0),dir=(0,.2,1),wavelength=.650),
Ray(pos=(0,0,0),dir=(0,-.2,1),wavelength=.650),
Ray(pos=(0,0,0),dir=(.2,0,1),wavelength=.650),
Ray(pos=(0,0,0),dir=(-.2,0,1),wavelength=.650),
Ray(pos=(0,0,0),dir=(0,0,1),wavelength=.650)]
S.propagate()
Plot3D(S,center=(0,0,100),size=(200,100),scale=4,rot=[(0,pi/2,0),(pi/20,-pi/10,0)])

``````

## Mirar algunos parametros de la lente

``````

In [ ]:

L1.paraxial_constants() #Distancia focal efectiva, punto focal anterior, punto focal posterior

``````

## Armar un sistema con 2 lentes

``````

In [ ]:

S=System(complist=[(L2,(0,0,100),(0,0,0)),(L3,(0,0,120),(0,0,0)) ],n=1)
R=[Ray(pos=(0,0,0),dir=(0,.2,1),wavelength=.650),
Ray(pos=(0,0,0),dir=(0,-.2,1),wavelength=.650),
Ray(pos=(0,0,0),dir=(.2,0,1),wavelength=.650),
Ray(pos=(0,0,0),dir=(-.2,0,1),wavelength=.650),
Ray(pos=(0,0,0),dir=(0,0,1),wavelength=.650)]
S.propagate()
Plot3D(S,center=(0,0,100),size=(200,100),scale=4,rot=[(0,pi/2,0),(pi/20,-pi/10,0)])

``````
``````

In [ ]:

S=System(complist=[(L2,(0,0,100),(0,0,0)),(L3,(0,0,120),(0,0,0)) ],n=1)
R=[]

for x in range(-30,40,10):
for y in range(-30,40,10):
R.append(Ray(pos=(0,0,0),dir=(x,y,100),wavelength=.650))
S.propagate()
Plot3D(S,center=(0,0,100),size=(200,100),scale=4,rot=[(0,pi/2,0),(pi/20,-pi/10,0)])

``````
``````

In [ ]:

C=CCD()

S=System(complist=[(L2,(0,0,100),(0,0,0)),(L3,(0,0,120),(0,0,0)),(C,(0,0,215),(0,0,0)) ],n=1)
R=[]

for x in range(-25,35,5):
for y in range(-25,35,5):
R.append(Ray(pos=(0,0,0),dir=(x,y,100),wavelength=.650))
S.propagate()
Plot3D(S,center=(0,0,100),size=(250,100),scale=4,rot=[(0,pi/2,0),(pi/20,-pi/10,0)])

``````
``````

In [ ]:

%pylab inline
spot_diagram(C)

``````

## Un espectroscopio con prisma

L=60 h=L/2.cos(pi/3.) h1=L/2sin(pi/3) S1=Plane(shape=Rectangular((L,L))) S2=Plane(shape=Rectangular((L,L))) S3=Plane(shape=Rectangular((L,L)))

T1=Plane(shape=Triangular(((2h,h1),(-2h,h1),(0,-h1)))) T2=Plane(shape=Triangular(((2h,h1),(-2h,h1),(0,-h1)))) P=Component(surflist=[(S1,(0,h1,0),(pi/2,0,0)), (S2,(0,0,h),(pi/6,0,0)), (S3,(0,0,-h),(-pi/6,0,0)), (T1,(L/2,0,0),(0,pi/2,0)), (T2,(-L/2,0,0),(0,pi/2,0)) ], material=material.schott["BK7"])

L=SphericalLens(radius=25,curvature_s1=1./200.,curvature_s2=-1./200,thickness=10,material=material.schott["BK7"]) C=CCD() PCCD=243.5 S=System(complist=[(P,(0,0,0),(0,0,0)), (L,(0,.3560,60),(-0.357,0,0)), (C,(0,.35PCCD,PCCD),(-0.357,0,0))],n=1)

R=[]

for x in range (-10,15,5): for y in range (-10,15,5): R.append(Ray(pos=(x,y+35,-100),dir=(0,-.35,1),wavelength=.45)) R.append(Ray(pos=(x,y+35,-100),dir=(0,-.35,1),wavelength=.5)) R.append(Ray(pos=(x,y+35,-100),dir=(0,-.35,1),wavelength=.55)) R.append(Ray(pos=(x,y+35,-100),dir=(0,-.35,1),wavelength=.60)) R.append(Ray(pos=(x,y+35,-100),dir=(0,-.35,1),wavelength=.65)) S.ray_add(R) S.propagate()

Plot3D(S,center=(0,.35*60,60),size=(500,200),scale=2,rot=[(0,pi/2+.2,0),(-.1,0,0)])

spot_diagram_c(C)

Plot3D(S,center=(0,.35*PCCD,PCCD),size=(50,30),scale=10,rot=[(0,pi/2+.1,0)])

## Otro espectroscopio con prisma

L=60 h=L/2.cos(pi/3.) h1=L/2sin(pi/3) S1=Plane(shape=Rectangular((L,L))) S2=Plane(shape=Rectangular((L,L))) S3=Plane(shape=Rectangular((L,L)))

T1=Plane(shape=Triangular(((2h,h1),(-2h,h1),(0,-h1)))) T2=Plane(shape=Triangular(((2h,h1),(-2h,h1),(0,-h1)))) P=Component(surflist=[(S1,(0,h1,0),(pi/2,0,0)), (S2,(0,0,h),(pi/6,0,0)), (S3,(0,0,-h),(-pi/6,0,0)), (T1,(L/2,0,0),(0,pi/2,0)), (T2,(-L/2,0,0),(0,pi/2,0)) ], material=material.schott["BK7"])

L=Component(surflist=[(S5,(0,0,5),(0,0,pi/2)), (S4,(0,0,-5),(0,0,pi/2)) ], material=material.schott["BK7"])

C=CCD() PCCD=243.5 S=System(complist=[(P,(0,0,0),(0,0,0)), (L,(0,.3560,60),(-0.357,0,0)), (C,(0,.35PCCD,PCCD),(-0.357,0,0))],n=1)

R=[]

for x in range (-5,6,1): for y in range (-5,6,2): R.append(Ray(pos=(x,y+35,-100),dir=(0,-.35,1),wavelength=.45)) R.append(Ray(pos=(x,y+35,-100),dir=(0,-.35,1),wavelength=.5)) R.append(Ray(pos=(x,y+35,-100),dir=(0,-.35,1),wavelength=.55)) R.append(Ray(pos=(x,y+35,-100),dir=(0,-.35,1),wavelength=.60)) R.append(Ray(pos=(x,y+35,-100),dir=(0,-.35,1),wavelength=.65)) S.ray_add(R) S.propagate()

Plot3D(S,center=(0,.35*60,60),size=(500,200),scale=2,rot=[(0,pi/2+.2,0),(-.1,0,0)])

spot_diagram_c(C)

Plot3D(S,center=(0,.35*PCCD,PCCD),size=(50,30),scale=10,rot=[(0,pi/2+.1,0)])

C=CCD()

# AP=Aperture(shape=Rectangular(size=(50,50)),ap_shape=Triangular(coord=((0,15),(15,-15),(-15,-15))))

APC=Component(surflist=[(AP,(0,0,0),(0,0,0))])

S=System(complist=[(L2,(0,0,100),(0,0,0)), (L3,(0,0,120),(0,0,0)), (C,(0,0,215),(0,0,0)), (APC,(0,0,110),(0,0,0)), ],n=1) R=[]

for x in range(-25,35,2): for y in range(-25,35,2): R.append(Ray(pos=(0,0,0),dir=(x,y,100),wavelength=.650)) S.ray_add(R) S.propagate() Plot3D(S,center=(0,0,100),size=(250,100),scale=4,rot=[(0,pi/4,0)])

## Reflección total interna

P0=RightAnglePrism(width=50,height=50,material=material.schott["BK7"],reflectivity=0) P1=RightAnglePrism(width=50,height=50,material=material.schott["BK7"],reflectivity=0)

S=System(complist=[(P0,(0,0,20),(0,0,pi/2)),(P1,(0,80,20),(0,0,-pi/2))],n=1)

R=[] for x in range (-5,6,2): for y in range (-5,6,2): R.append(Ray(pos=(x,y,-100),dir=(0,0,1),wavelength=.45))

for x in range (-5,6,2): for z in range (35,45,2): R.append(Ray(pos=(x,100,z),dir=(0,-1,0),wavelength=.540))

## Multiples fuentes y distintos tipos de fuentes

r_b=point_source_r(origin=(0.,0.,0.),direction=(0.,0.,0),span=pi/256 ,num_rays=10,wavelength=0.470, label="blue")

r_b1=point_source_r(origin=(30.,0.,0.),direction=(0.,0.,0),span=pi/256 ,num_rays=10,wavelength=0.470, label="blue1")

# ,num_rays=(10,10),wavelength=0.470, label="blue1")

N_BK7=material.schott['BK7'] SF5=material.schott['SF5']

# Dobletes 32-327 Edmund Scientific

DB1=Doublet(radius=12.5, curvature_s1 =1./61.47, curvature_s2 =-1./44.64, curvature_s3 =-1./129.94, thickness_l1 = 6., thickness_l2 = 2.5, material_l1 = N_BK7, material_l2 = SF5)

DB2=Doublet(radius=12.5, curvature_s1 =1./61.47, curvature_s2 =-1./44.64, curvature_s3 =-1./129.94, thickness_l1 = 6., thickness_l2 = 2.5, material_l1 = N_BK7, material_l2 = SF5)

# Definition of a detector plane

ccd=CCD() ccd1=CCD()

# Place de tetectors at the focal planes of the lenses

os=System(complist=[(DB1,(30,0,200),(0,0,0)), (DB2,(0,0,200),(pi,0,0)), (ccd,(30,0,float(400)),(0,0,0)), (ccd1,(0,0,400),(0,0,0)) ],n=1)

os.propagate()

display(Plot3D(os,center=(0,0,200), size=(500,100),scale=2,rot=[(0,pi/2+.1,0),(-pi/4,0,0)]))

display(spot_diagram(ccd)) figure() display(spot_diagram(ccd1))

## Par de dobletes

L1=library.Edmund.get("45408") ##f20,D5 L2=library.Edmund.get("45407")

ccd=CCD(size=(15,15))

spx=10.0pi/180. spy=10.0pi/180.

R1=point_source_c(origin=(0, 0., -20), direction=(0., 0, 0), span=(spy, spx),num_rays=(10,10),wavelength=0.600 )

S=System(complist=[(L1, (0, 0, 0), (0, 0, 0)), (L2, (0, 0, 7), (0, 0, 0)), (ccd, (0, 0, 25), (0, 0, 0)),

``````               ], n=1)

``````

## Lente Asférica

r_b= parallel_beam_c(size=(2,2),num_rays=(5,5), wavelength=.650)

R=2.75 k=-0.6139160 A2=0 A4=5.8891900E-04 A6=-1.7660200E-05 A8=1.0102500E-05 A10=-3.9148700E-06

r2=poly2d((0,0,0,1.,0,1.)) r4=r2r2 r6=r4r2 r8=r4r4 r10=r8r2

poly=A2r2+A4r4+ A6r6 +A8r8 +A10*r10

R=-3.1885400 k=-12.6638600

A2=0 A4=1.2458340e-02 A6=-3.7119450e-03 A8=5.1223910e-04 A10=-3.1085780e-05 poly=A2r2+A4r4+ A6r6 +A8r8 +A10*r10

asf1=Aspherical(Kx=k, Ky=k, Ax=1./R,Ay=1./R, shape=Circular(radius=2.5), poly=poly, reflectivity=.5)

oc=Component(surflist=[(asf2, (0, 0, 0), (0, 0, 0)), (asf1, (0, 0, 2.8+.35), (0,0, 0)) ], material=1.58913)

ccd=CCD(size=(3,3))

os=System(complist=[(oc,(0,0,20),(0,0,0)), (ccd,(0,0,20+2.8+2.14),(0,0,0)), ],n=1)

display(Plot3D(os,center=(0,0,10), size=(50,20),scale=16,rot=[(0,-3*pi/8,0)]))

spot_diagram(ccd)

N_BK7=material.schott["BK7"] N_BAK4=material.schott["BAK4"] N_SF10=material.schott["SF10"]

# r_b= parallel_beam_c(size=(10,10),num_rays=(10,10), wavelength=.470)

r_b=point_source_r(origin=(0.,0.,0.),direction=(0.,0.,0),span=pi/512 ,num_rays=100,wavelength=0.470)

# r_g= parallel_beam_c(size=(10,10),num_rays=(10,10), wavelength=.540)

r_g=point_source_r(origin=(0.,0.,0.),direction=(0.,0.,0),span=pi/512 ,num_rays=100,wavelength=0.540)

# r_r= parallel_beam_c(size=(10,10),num_rays=(10,10), wavelength=.670)

r_r=point_source_r(origin=(0.,0.,0.),direction=(0.,0.,0),span=pi/512 ,num_rays=100,wavelength=0.670)

# Definition

oc=Doublet(radius=25, curvature_s1 =1./162.59, curvature_s2 =-1./123.82, curvature_s3 =-1./402.58, thickness_l1 = 9.75, thickness_l2 = 3.50, material_l1 = N_BAK4, material_l2 = N_SF10)

bs=BeamSplitingCube(size=50,material=N_BK7,reflectivity=0.5)

ccd=CCD()

# Place de tetectors at the focal planes of the lenses

os=System(complist=[(oc,(0,0,500),(.50,0,0)), (ccd,(0,0,float(990)),(0,0,0)), (bs,(0,0,float(750)),(0,0,0)), ],n=1)

display(Plot3D(os,center=(0,0,500), size=(1200,200),scale=1,rot=[(0,pi/2,0),(.3,0,0)]))

spot_diagram_c(ccd)

## Sistema con divisores de haz

N_BK7=material.schott['BK7']

bs=BeamSplitingCube(size=8,material=N_BK7,reflectivity=0.5)

L1=library.Edmund.get("45408") ##f20,D5 L2=library.Edmund.get("31861") ccd=CCD(size=(10,10), transparent=False)

spx=10.0pi/180. spy=10.0pi/180.

R1=point_source_c(origin=(0, 0., -15), direction=(0., 0, 0), span=(spy, spx),num_rays=(10,10),wavelength=0.450 )

S=System(complist=[(L1, (0, 0, 0), (0, 0, 0)), (bs, (0, 0, 15), (0, 0, 0)), (L2, (-10, 0, 15), (0, 90pi/180., 0)), (bs, (-20, 0, 15), (90pi/180., 0, 0))

``````               ], n=1)

``````

S.propagate() Plot3D(S,center=(-15,0,7), size=(50,50),scale=16,rot=[(0,-pi/2,0),(-pi/2+.1,-pi/4,0)])

# Ejercicio:

Simular un expansor de haz fabricado con 2 lentes esfericas de las siguientes caracteristicas:

1. Lente convergente fabricada en BK7 de 50 mm de diametro con radios de curvatura 100 y -100 mm
2. Lente convergente fabricada en BK7 de 50 mm de diametro con radios de curvatura 100 y -100 mm
``````

In [ ]:

``````