elena/gaussianoptics
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Update ionBeam854nm.py

Browse Source
This commit is contained in:
Elena Arenskötter 2023-11-04 13:44:26 +01:00
parent b894e76e9f
commit 4abb5fbd44
1 changed files with 76 additions and 76 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

152
ionBeam854nm.py
View File

@ -1,76 +1,76 @@
# -*- coding: utf-8 -*- # -*- coding: utf-8 -*-
""" """
Created on Fri Aug 14 14:17:54 2020 Created on Fri Aug 14 14:17:54 2020
@author: Jan @author: Elena
""" """
import numpy as np import numpy as np
import matplotlib.pyplot as plt import matplotlib.pyplot as plt
import gaussianoptics as go import gaussianoptics as go
#Laser wavelength #Laser wavelength
lam = 854e-9 lam = 854e-9
#Number of points #Number of points
N =5001 N =5001
#Focal length and position definition #Focal length and position definition
lens_coupler = 12e-3 lens_coupler = 12e-3
f1_tele = 30e-3 f1_tele = 30e-3
pos_f1 = 100e-3 pos_f1 = 100e-3
f2_tele = 300e-3 f2_tele = 300e-3
pos_f2 = 330e-3 pos_f2 = 330e-3
haloPosDif393_854 = 1.4e-3 haloPosDif393_854 = 1.4e-3
f_halo1 = 25e-3 f_halo1 = 25e-3
pos_halo1 = 300e-3 pos_halo1 = 300e-3
f_halo2 = 25e-3 f_halo2 = 25e-3
pos_halo2 = 50e-3 - haloPosDif393_854 pos_halo2 = 50e-3 - haloPosDif393_854
f1_colTele = 300e-3; f1_colTele = 300e-3;
pos_colTele1= 400e-3; pos_colTele1= 400e-3;
f2_colTele = 25e-3; f2_colTele = 25e-3;
pos_colTele2= 497.5e-3; pos_colTele2= 497.5e-3;
#Definition of optical arangement #Definition of optical arangement
opticalbench = [['path',1,0,'f'], opticalbench = [['path',1,0,'f'],
['thinlens',12e-3,lens_coupler,'r'], ['thinlens',12e-3,lens_coupler,'r'],
['path',1,0,'f'], ['path',1,0,'f'],
['thinlens',f1_tele,pos_f1,'r'], ['thinlens',f1_tele,pos_f1,'r'],
['path',1,0,'f'], ['path',1,0,'f'],
['thinlens',f2_tele,pos_f2,'r'], ['thinlens',f2_tele,pos_f2,'r'],
['path',1,0,'f'], ['path',1,0,'f'],
['thinlens',f_halo1,pos_halo1,'r'], ['thinlens',f_halo1,pos_halo1,'r'],
['path',1,0,'f'], ['path',1,0,'f'],
['thinlens',f_halo2,pos_halo2,'r'], ['thinlens',f_halo2,pos_halo2,'r'],
['path',1,0,'f'], ['path',1,0,'f'],
['thinlens',f1_colTele,pos_colTele1,'r'], ['thinlens',f1_colTele,pos_colTele1,'r'],
['path',1,0,'f'], ['path',1,0,'f'],
['thinlens',f2_colTele,pos_colTele2,'r'], ['thinlens',f2_colTele,pos_colTele2,'r'],
['path',1,0.3,'r'] ['path',1,0.3,'r']
] ]
#Beamwaist inside fiber 780HP #Beamwaist inside fiber 780HP
w0 = 2.5e-6 w0 = 2.5e-6
z0 = np.pi*w0**2/lam z0 = np.pi*w0**2/lam
q0 = complex(0,z0) q0 = complex(0,z0)
[x,w,R,OB,M_tot] = go.gaussianoptics(opticalbench,lam,q0,0,N) [x,w,R,OB,M_tot] = go.gaussianoptics(opticalbench,lam,q0,0,N)
fig, ax = plt.subplots() fig, ax = plt.subplots()
ax.plot(x*1000, w*1000, color='black') ax.plot(x*1000, w*1000, color='black')
ax.plot(x*1000, -w*1000, color='black') ax.plot(x*1000, -w*1000, color='black')
plt.grid(True,which="both") plt.grid(True,which="both")
plt.show() plt.show()