Commit 44da8d0a authored by Lucas Laplanche's avatar Lucas Laplanche
Browse files

nettoyage code

parent abe25178
...@@ -33,8 +33,6 @@ def refra(bypass_dbr=True, electric_field=0., wavelength=850e-9): ...@@ -33,8 +33,6 @@ def refra(bypass_dbr=True, electric_field=0., wavelength=850e-9):
plt.plot_refra(sl) plt.plot_refra(sl)
def al_doping(bypass_dbr=True): def al_doping(bypass_dbr=True):
sl = st.structure_eam_vcsel(eam_only=True, bypass_dbr=bypass_dbr) sl = st.structure_eam_vcsel(eam_only=True, bypass_dbr=bypass_dbr)
plt.plot_al_doping(sl) plt.plot_al_doping(sl)
...@@ -73,29 +71,6 @@ def reflectivity(bypass_dbr=True, ...@@ -73,29 +71,6 @@ def reflectivity(bypass_dbr=True,
return r return r
def dbr_reflectivity(start_wavelength=700e-9,
stop_wavelength=1000e-9,
electric_field=0.,
n_points=100,
plot=True):
wavelength = np.linspace(start_wavelength, stop_wavelength, num=n_points)
r = np.zeros(len(wavelength))
# wavelength in [m]
# wavelength must be a numpy array
for i in tqdm(range(len(wavelength))):
sl = st.structure_dbr(period=20, grading_type = 'linear slope')
sl = op.algaas_super_lattice_refractive_index(sl, electric_field, wavelength[i], lengyel=True)
n = sl['refractive_index'].to_numpy(dtype=np.complex128)
d = sl['thickness'].to_numpy(dtype=float)
r[i] = tmm.reflection(n, d, wavelength[i])
if plot:
plt.plot_reflectivity(wavelength, r)
else:
return
def reflectivity_heatmap(bypass_dbr=True, def reflectivity_heatmap(bypass_dbr=True,
...@@ -199,16 +174,6 @@ def clad_coupling_electromagnetic_amplitude(bypass_dbr=False, electric_field=0., ...@@ -199,16 +174,6 @@ def clad_coupling_electromagnetic_amplitude(bypass_dbr=False, electric_field=0.,
plt.plot_refra_clad_coupling_em(sl1, sl2, sl3) plt.plot_refra_clad_coupling_em(sl1, sl2, sl3)
def animation_clad_coupling(bypass_dbr=False, electric_field=0., wavelength=850e-9):
cladding = np.linspace(9.9e-9, 10.1e-9, 9)
for cl in cladding:
sl = st.structure_eam_vcsel(bypass_dbr=bypass_dbr, vcsel_only=True, grading_type='none', mqw_alloy_type='none', l_vcsel_clad=cl)
sl = op.algaas_super_lattice_refractive_index(sl, electric_field, wavelength)
sl = pt.cut_in_equal_layers_thickness(sl, 5e-9)
em = tmm.em_amplitude_scattering_matrix(sl, wavelength)
sl.insert(sl.shape[1], 'electromagnetic_amplitude', value=em)
plt.plot_refra_em(sl, str(cl))
def zandberg_electromagnetic_amplitude(electric_field=0., wavelength=1e-6): def zandberg_electromagnetic_amplitude(electric_field=0., wavelength=1e-6):
...@@ -273,7 +238,7 @@ def ftir_theory_comparison_reflectivity(filename='eam_2', electric_field=0., doc ...@@ -273,7 +238,7 @@ def ftir_theory_comparison_reflectivity(filename='eam_2', electric_field=0., doc
def mqw_psi(lz=0.01e-9): def mqw_psi(lz=0.01e-9):
sl = st.structure_eam_vcsel(eam_mqw_only=True, amount_eam_qw=5) sl = st.structure_eam_vcsel(amount_eam_qw=5)
sl = pt.cut_in_equal_layers_thickness(sl, lz) sl = pt.cut_in_equal_layers_thickness(sl, lz)
v_e, v_hh, eig_e, psi_e, eig_hh, psi_hh = qt.solve_schrodinger(sl) v_e, v_hh, eig_e, psi_e, eig_hh, psi_hh = qt.solve_schrodinger(sl)
np.savez_compressed('v_e', v_e) np.savez_compressed('v_e', v_e)
......
...@@ -35,6 +35,8 @@ def eam_heatmaps(): ...@@ -35,6 +35,8 @@ def eam_heatmaps():
r_file_name='heatmap_r_eam_2') r_file_name='heatmap_r_eam_2')
def reflectivity_eam_vcsel_clad_heatmap(): def reflectivity_eam_vcsel_clad_heatmap():
n_pts = 30 n_pts = 30
l_eam_clad =np.linspace(15, 20, num=n_pts) l_eam_clad =np.linspace(15, 20, num=n_pts)
...@@ -59,6 +61,8 @@ def reflectivity_eam_vcsel_clad_heatmap(): ...@@ -59,6 +61,8 @@ def reflectivity_eam_vcsel_clad_heatmap():
plt.plot_std_heatmap(l_eam_clad, l_vcsel_clad, r) plt.plot_std_heatmap(l_eam_clad, l_vcsel_clad, r)
def structure_for_matlab(vcsel_only=False, bypass_dbr=False, wavelength=850.5e-9): def structure_for_matlab(vcsel_only=False, bypass_dbr=False, wavelength=850.5e-9):
l_eam_clad = 5e-9 l_eam_clad = 5e-9
l_vcsel_clad = 15e-9 l_vcsel_clad = 15e-9
...@@ -70,5 +74,3 @@ def structure_for_matlab(vcsel_only=False, bypass_dbr=False, wavelength=850.5e-9 ...@@ -70,5 +74,3 @@ def structure_for_matlab(vcsel_only=False, bypass_dbr=False, wavelength=850.5e-9
n = sl['refractive_index'].apply(np.real).to_numpy(dtype=float) n = sl['refractive_index'].apply(np.real).to_numpy(dtype=float)
scipy.io.savemat('eam_vcsel_n_lz.mat', dict(n=n, lz=lz)) scipy.io.savemat('eam_vcsel_n_lz.mat', dict(n=n, lz=lz))
OBJECTIFS OBJECTIFS
Nettoyer le code de la structure python.
Traiter les courbes mesurées de la refléctivité temporelle. Traiter les courbes mesurées de la refléctivité temporelle.
Les comparer aux calculs théoriques. Les comparer aux calculs théoriques.
Faire l'épitaxie.
Faire le fit du FTIR. Faire le fit du FTIR.
Vérifier les recettes avec Alexandre Vérifier les recettes avec Alexandre
Faire l'épitaxie.
ACCOMPLI ACCOMPLI
Nettoyer le code de la structure python.
Faire les calculs VCSEL EAM avec oxydation eam. Faire les calculs VCSEL EAM avec oxydation eam.
Réinstaller crystal XE et faire les recettes EAM-VCSELs. Réinstaller crystal XE et faire les recettes EAM-VCSELs.
Corriger le tmm !!!!!!!!! Corriger le tmm !!!!!!!!!
...@@ -20,8 +17,4 @@ ACCOMPLI ...@@ -20,8 +17,4 @@ ACCOMPLI
ECHOUÉ ECHOUÉ
FTIR. FTIR.
Calcul carte de chaleur reflectivité Calcul carte de chaleur reflectivité
\ No newline at end of file
QUESTIONS
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment