Merge branch 'master' of ssh://code.df.uba.ar:5022/nnunez/artiq_experiments

analisis/plots/20231123_CPTconmicromocion3/CPT_plotter_20231123.py

@@ -60,6 +60,7 @@ for i, fname in enumerate(CPT_FILES.split()):
     UVCPTAmp.append(np.array(data['datasets']['UV_CPT_amp']))
     No_measures.append(np.array(data['datasets']['no_measures']))
     Voltages.append(np.array(data['datasets']['scanning_voltages']))
+    print(np.array(data['datasets']['no_measures']))
 
 def Split(array,n):
     length=len(array)/n
@@ -89,7 +90,7 @@ CountsSplit_2ions.append(Split(Counts[4],len(Freqs[4])))
 Ploteo curvas para ver que tal son
 """
 
-jvec = [2] # de la 1 a la 9 vale la pena, despues no
+jvec = [4] # de la 1 a la 9 vale la pena, despues no
 
 
 Frequencies = Freqs[0]
@@ -151,7 +152,9 @@ alpha = 0
 drivefreq = 2*np.pi*22.135*1e6
 
 
-SelectedCurveVec = [1,2,3,4,5,6,7,8,9,10,11]
+#SelectedCurveVec = [1,2,3,4,5,6,7,8,9,10,11]
+
+SelectedCurveVec = [4]
 
 popt_SA_vec = []
 pcov_SA_vec = []
@@ -218,7 +221,11 @@ for selectedcurve in SelectedCurveVec:
         pcov_SA_vec.append(pcov_3_SA)
 
         FittedEITpi_3_SA_short, Detunings_3_SA_short = FitEIT_MM_single(FreqsDR, *popt_3_SA, plot=True)
-        freqslong = np.arange(min(FreqsDR), max(FreqsDR)+FreqsDR[1]-FreqsDR[0], 0.1*(FreqsDR[1]-FreqsDR[0]))
+        freqslong = np.arange(min(FreqsDR)-1*5e6, max(FreqsDR)+FreqsDR[1]-FreqsDR[0]-1*5e6, 0.1*(FreqsDR[1]-FreqsDR[0]))
+        
+        popt_3_SA[4]=1
+        popt_3_SA[5]=0
+        popt_3_SA[6]=0
         FittedEITpi_3_SA_long, Detunings_3_SA_long = FitEIT_MM_single(freqslong, *popt_3_SA, plot=True)
 
     DetuningsUV_vec.append(popt_3_SA[1])
@@ -246,6 +253,23 @@ for selectedcurve in SelectedCurveVec:
 
     print(f'listo med {selectedcurve}')
     print(popt_3_SA)
+    
+    
+#%%
+
+plt.figure(figsize=(0.8*5,0.7*3))
+#plt.errorbar(Detunings_3_SA_short, CountsDR, yerr=2*np.sqrt(CountsDR), fmt='o', color='darkgreen', alpha=0.5, capsize=2, markersize=2)
+plt.plot(Detunings_3_SA_long, [100*f for f in FittedEITpi_3_SA_long], color='darkolivegreen', linewidth=3, label=f'med {selectedcurve}')
+#plt.title(f'Sdop: {round(popt[0], 2)}, Spr: {round(popt[1], 2)}, T: {round(popt[2]*1e3, 2)} mK, detDop: {DetDoppler} MHz')
+plt.xlabel(r'$\Delta_{\mathrm{IR}}$ (MHz)',fontsize=10)
+plt.ylabel('Fluorescence (arb. units)',fontsize=10)
+#plt.legend(loc='upper left', fontsize=20)
+plt.grid()
+plt.xlim(-60,18)
+plt.tight_layout()
+plt.savefig('/home/nico/Nextcloud/G_liaf/Publicaciones/Papers/2023 Micromotion CPT/Drafts/cpt_fig1_v2.pdf')
+
+
 
 
 #%%

artiq_master/last_rid.pyon

-16638
+16673