structures corrigées, interaction matlab ajoutée

calculation.py

@@ -66,11 +66,11 @@ def al_doping(bypass_dbr=True):
 
 
 def reflectivity(bypass_dbr=False,
-                 start_wavelength=849e-9,
-                 stop_wavelength=852e-9,
+                 start_wavelength=850e-9,
+                 stop_wavelength=851e-9,
                  electric_field=0.,
                  n_points=100,
-                 l_eam_clad=15e-9,
+                 l_eam_clad=5e-9,
                  l_vcsel_clad=15e-9,
                  plot=True):
     wavelength = np.linspace(start_wavelength, stop_wavelength, num=n_points)
@@ -79,7 +79,7 @@ def reflectivity(bypass_dbr=False,
     # wavelength in [m]
     # wavelength must be a numpy array
     for i in tqdm(range(len(wavelength))):
-        sl = st.structure_eam_vcsel(bypass_dbr=bypass_dbr, vcsel_only=True, eam_only=False, grading_type='linear digital', mqw_alloy_type='digital',
+        sl = st.structure_eam_vcsel(bypass_dbr=bypass_dbr, vcsel_only=False, eam_only=False, grading_type='linear digital', mqw_alloy_type='digital',
                                     l_eam_clad=l_eam_clad, l_vcsel_clad=l_vcsel_clad)
 
         sl = op.algaas_super_lattice_refractive_index(sl, electric_field, wavelength[i], lengyel=True)

macro.py

@@ -46,10 +46,10 @@ def reflectivity_eam_vcsel_clad_heatmap():
     plt.plot_std_heatmap(l_eam_clad, l_vcsel_clad, r)
 
 
-def structure_for_matlab(vcsel_only=False, bypass_dbr=False, wavelength=850e-9):
-    l_eam_clad = 8e-9
+def structure_for_matlab(vcsel_only=False, bypass_dbr=False, wavelength=850.5e-9):
+    l_eam_clad = 5e-9
     l_vcsel_clad = 15e-9
-    sl = st.structure_eam_vcsel(bypass_dbr=bypass_dbr, vcsel_only=vcsel_only, grading_type='mean', mqw_alloy_type='none',
+    sl = st.structure_eam_vcsel(bypass_dbr=bypass_dbr, vcsel_only=vcsel_only, grading_type='linear digital', mqw_alloy_type='digital',
                                 l_eam_clad=l_eam_clad, l_vcsel_clad=l_vcsel_clad)
     sl = op.algaas_super_lattice_refractive_index(sl, 0., wavelength)
 

xlsx_export.py

@@ -76,7 +76,7 @@ def eam_structures_description(document_folder=True):
 
 def vcsel_eam_structures_description(document_folder=False):
     # create the structures
-    vcsel_eam = st.structure_eam_vcsel(bypass_dbr=False, vcsel_only=False, grading_type='mean', mqw_alloy_type='none',
+    vcsel_eam = st.structure_eam_vcsel(bypass_dbr=False, vcsel_only=False, eam_only=True, grading_type='linear digital', mqw_alloy_type='digital',
                                 l_eam_clad=10e-9, l_vcsel_clad=15e-9)
     vcsel_eam_dbr_bypass = st.structure_eam_vcsel(bypass_dbr=True, grading_type='none', mqw_alloy_type='none')