agrego todo antes de irme

analisis/plots/20230804_RotationalDopplerShift_v2/RDS_CPT2.py

@@ -262,7 +262,7 @@ secular=False
 
 
 plt.plot([2*f*1e-6 for f in ExtraIR1_Freqs[1][1:]], [1e-3*c for c in ExtraCounts[1][1:]], '-o', markersize=2, label='Colineales')
-plt.xlabel('Frequency (MHz)', fontname='STIXgeneral',fontsize=15)
+plt.xlabel('IR1 Frequency (MHz)', fontname='STIXgeneral',fontsize=15)
 plt.ylabel('kCounts',fontname='STIXgeneral',fontsize=15)
 plt.xlim(425, 455)
 plt.yticks([1,1.5,2,2.5,3],fontname='STIXgeneral',fontsize=15)

analisis/plots/20230817_RotationalDopplerShift_v5/RDS_sidebands2.py

@@ -50,41 +50,14 @@ def SeeKeys(files):
 PiezoVerCounts = []
 PiezoVerFrequencies = []
 
-PIEZOVER_FILES = np.arange(922, 948,1)
+PIEZOVER_FILES = [167,168,169]
 
 for i in PIEZOVER_FILES:
     #print(str(i) + ' - ' + fname)
-    data = h5py.File(f'VaryingBeamlocation/Directions/Vertical/000014{i}-IR_Scan_withcal_optimized'+'.h5', 'r')
+    data = h5py.File(f'Sidebands/000015{i}-IR_Scan_withcal_optimized'+'.h5', 'r')
     PiezoVerCounts.append(np.array(data['datasets']['counts_spectrum']))
     PiezoVerFrequencies.append(np.array(data['datasets']['IR1_Frequencies']))
 
-PiezoVerDCounts = []
-PiezoVerDFrequencies = []
-
-PIEZOVERD_FILES = list(np.arange(70, 97,1))+list(np.arange(99,116))
-
-
-
-for i in PIEZOVERD_FILES:
-    #print(str(i) + ' - ' + fname)
-    if i <100:
-        data = h5py.File(f'VaryingBeamlocation/Directions/VerticalDense/0000150{i}-IR_Scan_withcal_optimized'+'.h5', 'r')
-    else:
-        data = h5py.File(f'VaryingBeamlocation/Directions/VerticalDense/000015{i}-IR_Scan_withcal_optimized'+'.h5', 'r')
-    PiezoVerDCounts.append(np.array(data['datasets']['counts_spectrum']))
-    PiezoVerDFrequencies.append(np.array(data['datasets']['IR1_Frequencies']))
-
-
-PiezoDiagCounts = []
-PiezoDiagFrequencies = []
-
-PIEZODIAG_FILES = list(np.arange(948, 976,1))
-
-for i in PIEZODIAG_FILES:
-    #print(str(i) + ' - ' + fname)
-    data = h5py.File(f'VaryingBeamlocation/Directions/Diagonal/000014{i}-IR_Scan_withcal_optimized'+'.h5', 'r')
-    PiezoDiagCounts.append(np.array(data['datasets']['counts_spectrum']))
-    PiezoDiagFrequencies.append(np.array(data['datasets']['IR1_Frequencies']))
 
 
 def ErrorDRdepth(p, f, b):
@@ -117,12 +90,12 @@ palette = sns.color_palette("tab10")
 pmlocmedvec = np.arange(0,len(PIEZOVER_FILES),1)
 
 
-#pmlocmedvec = [26]
+pmlocmedvec = [2]
 
 
 plt.figure()
 
-bkg = np.min(PiezoVerCounts[5])
+#bkg = np.min(PiezoVerCounts[5])
 
 pmdepthsdrver=[]
 errorpmdepthsdrver=[]
@@ -135,235 +108,37 @@ for med in pmlocmedvec:
 
     Freqs = [2*f*1e-6 for f in PiezoVerFrequencies[med][1:]]
     Counts = [c for c in PiezoVerCounts[med][1:]]
-
-    popt, pcov = curve_fit(Lorentzian, Freqs, Counts, p0=(-200,2100,435.8,0.05), bounds=((-10000,0,435.5,0),(0,1e4, 436.1, 1)))
-
-    pmdepthsdrver.append(1-(np.min(Lorentzian(Freqs,*popt))-bkg)/(popt[1]-bkg))
-    errorpmdepthsdrver.append(ErrorDRdepth(np.min(Lorentzian(Freqs,*popt)),popt[1], bkg))
-        
-    Intens = popt[1]
     
-    Intensityver.append(Intens)
-    # errorIntensity.append(2*np.sqrt(np.mean(Piezo1Counts[med][1:][0:20]))+np.sqrt(bkg))
-        
-    
-    if med not in [800]:
-        plt.plot([2*f*1e-6 for f in PiezoVerFrequencies[med][1:]], [c for c in PiezoVerCounts[med][1:]], '-o', markersize=2, alpha=0.7)
-        plt.plot(Freqs,Lorentzian(Freqs,*popt))
-    
-    jj=jj+1
-# plt.xlabel('Frecuencia (MHz)')
-# plt.ylabel('Counts')
-#plt.xlim(435.2, 436.5)
-plt.grid()
-# plt.legend()
-# #plt.title('Espectros para distintas geometrías')
-
-
-plt.figure()
-plt.plot(np.arange(0,len(Intensityver),1), [i/np.max(Intensityver) for i in Intensityver], '-o',markersize=8)
-plt.plot(np.arange(0,len(Intensityver),1), [p for p in pmdepthsdrver], 'o',markersize=8)
-plt.xlabel('Ion position')
-plt.ylabel('Intensity / DR Relative depth')
-#plt.xticks([1,2,3,4,5])
-#plt.xlim(200,3200)
-plt.ylim(-0.1,1.1)
-plt.grid()
-#plt.axvline(3, color='salmon')
-plt.legend()
-
-#%%
-"""
-Ahora voy a intentar ajustarlas con una lorentziana que es mejor
-"""
-import seaborn as sns
-"""
-Resonancias DD configuracion +2/-2 colineal variando la ubicacion del ion en los haces
-
-Moviendo diagonalmente el haz
-"""
-
-
-
-palette = sns.color_palette("tab10")
-
-pmlocmedvec = np.arange(0,len(PIEZODIAG_FILES),1)
-
-
-#pmlocmedvec = [0,1]
-
-
-plt.figure()
-
-bkg = np.min(PiezoDiagCounts[5])
-
-pmdepthsdrdiag=[]
-errorpmdepthsdrdiag=[]
-
-Intensitydiag = []
-errorIntensitydiag = []
-
-jj=0
-for med in pmlocmedvec:
-
-    Freqs = [2*f*1e-6 for f in PiezoDiagFrequencies[med][1:]]
-    Counts = [c for c in PiezoDiagCounts[med][1:]]
+    if med==1:
+        Counts[119]=2100
+        Counts[240]=2050
 
     if med==2:
-        Freqs = Freqs[1:-30]
-        Counts = Counts[1:-30]
-        popt, pcov = curve_fit(Lorentzian, Freqs, Counts, p0=(-200,2100,435.8,0.05), bounds=((-10000,0,435.5,0),(0,1e4, 436.1, 1)))
-    elif med==1:
-        Freqs = Freqs[10:-30]
-        Counts = Counts[10:-30]
-        popt, pcov = curve_fit(Lorentzian, Freqs, Counts, p0=(-200,2100,435.8,0.05), bounds=((-10000,0,435.7,0),(0,1e4, 436.1, 1)))
- 
-    elif med==5:
-        Freqs = Freqs[10:-55]+Freqs[-30:-1]
-        Counts = Counts[10:-55]+Counts[-30:-1]
-        popt, pcov = curve_fit(Lorentzian, Freqs, Counts, p0=(-200,2100,435.8,0.05), bounds=((-10000,0,435.5,0),(0,1e4, 436.1, 1)))
+        Counts[425]=1400
+        Counts[469]=1300
+        Counts[470]=1300
+        Counts[471]=1300
         
-    else:
-        popt, pcov = curve_fit(Lorentzian, Freqs, Counts, p0=(-200,2100,435.8,0.05), bounds=((-10000,0,435.5,0),(0,1e4, 436.1, 1)))
+    #popt, pcov = curve_fit(Lorentzian, Freqs, Counts, p0=(-200,2100,435.8,0.05), bounds=((-10000,0,435.5,0),(0,1e4, 436.1, 1)))
 
-    pmdepthsdrdiag.append(1-(np.min(Lorentzian(Freqs,*popt))-bkg)/(popt[1]-bkg))
-    errorpmdepthsdrdiag.append(ErrorDRdepth(np.min(Lorentzian(Freqs,*popt)),popt[1], bkg))
+    #pmdepthsdrver.append(1-(np.min(Lorentzian(Freqs,*popt))-bkg)/(popt[1]-bkg))
+    #errorpmdepthsdrver.append(ErrorDRdepth(np.min(Lorentzian(Freqs,*popt)),popt[1], bkg))
         
-    Intens = popt[1]
+    #Intens = popt[1]
     
-    Intensitydiag.append(Intens)
+    #Intensityver.append(Intens)
     # errorIntensity.append(2*np.sqrt(np.mean(Piezo1Counts[med][1:][0:20]))+np.sqrt(bkg))
         
     
     if med not in [800]:
-        plt.plot([2*f*1e-6 for f in PiezoDiagFrequencies[med][1:]], [c for c in PiezoDiagCounts[med][1:]], '-o', markersize=2, alpha=0.7)
-        plt.plot(Freqs,Lorentzian(Freqs,*popt))
-    
-    jj=jj+1
-# plt.xlabel('Frecuencia (MHz)')
-# plt.ylabel('Counts')
-#plt.xlim(435.2, 436.5)
-plt.grid()
-# plt.legend()
-# #plt.title('Espectros para distintas geometrías')
-
-
-plt.figure()
-plt.plot(np.arange(0,len(Intensitydiag),1), [i/np.max(Intensitydiag) for i in Intensitydiag], '-o',markersize=8)
-plt.plot(np.arange(0,len(Intensitydiag),1), [p for p in pmdepthsdrdiag], 'o',markersize=8)
-plt.xlabel('Ion position')
-plt.ylabel('Intensity / DR Relative depth')
-#plt.xticks([1,2,3,4,5])
-#plt.xlim(200,3200)
-plt.ylim(-0.1,1.1)
-plt.grid()
-#plt.axvline(3, color='salmon')
-plt.legend()
-
-
-#%%
-
-"""
-Ploteo en conjunto. La horizontal sale de RDS_piezo.py
-"""
-cap=3
-
-plt.figure()
-plt.plot(np.arange(0,len(Intensitydiag),1), [i/np.max(Intensitydiag) for i in Intensitydiag], '-o',markersize=8)
-#plt.plot(np.arange(0,len(Intensitydiag),1), [p for p in pmdepthsdrdiag], 'o',markersize=8, label='Diagonal')
-plt.errorbar(np.arange(0,len(Intensitydiag),1), [p for p in pmdepthsdrdiag], yerr=errorpmdepthsdrdiag, fmt='o',capsize=cap,markersize=8, label='Diagonal')
-
-
-#plt.plot(np.arange(0,len(Intensityver),1), [i/np.max(Intensityver) for i in Intensityver], '-o',markersize=8)
-#plt.plot(np.arange(0,len(Intensityver),1), [p for p in pmdepthsdrver], 'o',markersize=8, label='Vertical')
-plt.errorbar(np.arange(0,len(Intensityver),1), [p for p in pmdepthsdrver], yerr=errorpmdepthsdrver, fmt='o', capsize=cap,markersize=8, label='Vertical')
-
-
-scale = 1.6
-
-#plt.plot([s*scale for s in np.arange(16,len(Intensity),1)-16], [i/np.max(Intensity) for i in Intensity[16:]], '-o',markersize=8)
-#plt.plot([s*scale for s in np.arange(16,len(Intensity),1)-16], [p for p in pmdepthsdr[16:]], 'o',markersize=8, label='Horizontal')
-plt.errorbar([s*scale for s in np.arange(16,len(Intensity),1)-16], [p for p in pmdepthsdr[16:]], yerr=errorpmdepthsdr[16:], fmt='o', capsize=cap,markersize=8, label='Horizontal')
-
-
-plt.xlabel('Ion position')
-plt.ylabel('Intensity / DR Relative depth')
-#plt.xticks([1,2,3,4,5])
-plt.xlim(-1,15)
-plt.ylim(-0.1,1.1)
-plt.grid()
-#plt.axvline(3, color='salmon')
-plt.legend()
-
-
-#%%
-
-import seaborn as sns
-"""
-Resonancias DD configuracion +2/-2 colineal variando la ubicacion del ion en los haces
-
-Moviendo verticalmente el haz pero con mayor densidad de puuntos
-"""
-
-
-
-palette = sns.color_palette("tab10")
-
-pmlocmedvec = np.arange(0,len(PIEZOVERD_FILES),1)
-
-
-#pmlocmedvec = [6]
-
-
-plt.figure()
-
-#bkg = np.min(PiezoVerDCounts[5])
-bkg=150
-
-pmdepthsdrverd=[]
-errorpmdepthsdrverd=[]
-
-Intensityverd = []
-errorIntensityverd = []
-
-jj=0
-for med in pmlocmedvec:
-
-    Freqs = [2*f*1e-6 for f in PiezoVerDFrequencies[med][1:]]
-    Counts = [c for c in PiezoVerDCounts[med][1:]]
-
-    popt, pcov = curve_fit(Lorentzian, Freqs, Counts, p0=(-200,2100,435.8,0.05), bounds=((-10000,0,435.5,0),(0,1e4, 436.1, 1)))
-
-    pmdepthsdrverd.append(1-(np.min(Lorentzian(Freqs,*popt))-bkg)/(popt[1]-bkg))
-    errorpmdepthsdrverd.append(ErrorDRdepth(np.min(Lorentzian(Freqs,*popt)),popt[1], bkg))
-        
-    Intens = popt[1]
-    
-    Intensityverd.append(Intens)
-    errorIntensityverd.append(2*np.sqrt(Intens)+np.sqrt(bkg))
-        
-    
-    if med not in [800]:
-        plt.plot([2*f*1e-6 for f in PiezoVerDFrequencies[med][1:]], [c for c in PiezoVerDCounts[med][1:]], '-o', markersize=2, alpha=0.7)
-        plt.plot(Freqs,Lorentzian(Freqs,*popt))
+        plt.plot(Freqs, Counts, '-o', markersize=2, alpha=0.7,color=paleta[med])
+        #plt.plot(Freqs,Lorentzian(Freqs,*popt))
     
     jj=jj+1
-# plt.xlabel('Frecuencia (MHz)')
-# plt.ylabel('Counts')
-#plt.xlim(435.2, 436.5)
+plt.xlabel('Frequency (MHz)')
+plt.ylabel('Counts')
+#plt.ylim(300,2600)
 plt.grid()
 # plt.legend()
 # #plt.title('Espectros para distintas geometrías')
 
-
-plt.figure()
-plt.errorbar(np.arange(0,len(Intensityverd),1), [i/np.max(Intensityverd) for i in Intensityverd], yerr=[e/np.max(Intensityverd) for e in errorIntensityverd], fmt='o', capsize=3, markersize=8)
-plt.errorbar(np.arange(0,len(Intensityverd),1), [p for p in pmdepthsdrverd], yerr=errorpmdepthsdrverd, fmt='o',color='limegreen',capsize=3,markersize=8)
-plt.xlabel('Screw step')
-plt.ylabel('Intensity / DR Relative depth')
-#plt.xticks([1,2,3,4,5])
-plt.xlim(-2,35)
-plt.ylim(-0.1,1.1)
-plt.grid()
-#plt.axvline(3, color='salmon')
-plt.legend()
\ No newline at end of file

analisis/plots/20230912_RotationalDopplerShift_v6/RDS_piezobeamsizes.py

@@ -46,8 +46,8 @@ def SeeKeys(files):
         print(fname)
         print(list(data['datasets'].keys()))
 
-def Lorentzian( x, A, B, x0, gam,C):
-    #C=0
+def Lorentzian( x, A, B, x0, gam):
+    C=0
     return A * gam**2 / ( gam**2 + ( x - x0 )**2) + B - C*(x - x0)
 
 
@@ -598,9 +598,12 @@ Resonancias DD configuracion +2/-2 colineal variando la ubicacion del ion en los
 Alta estadistica
 """
 
+def Lorentzian( x, A, B, x0, gam,C):
 
+    return A * gam**2 / ( gam**2 + ( x - x0 )**2) + B - C*(x - x0)
 
-palette = sns.color_palette("tab10")
+
+palette = sns.color_palette("rocket")
 
 pmlocmedvec = np.arange(0,len(HS_FILES),1)
 
@@ -644,17 +647,17 @@ for med in pmlocmedvec:
         
     fi=775
     if med not in [800]:
-        plt.plot([1e3*(2*f*1e-6-435)-fi for f in PiezoHSFrequencies[med][1:]], [1e-3*c for c in PiezoHSCounts[med][1:]], '-o', markersize=2, alpha=0.7)
-        plt.plot([1e3*(f-435)-fi for f in Freqs],[l*1e-3 for l in Lorentzian(Freqs,*popt)],linewidth=5)
+        plt.plot([1e3*(2*f*1e-6-435)-fi for f in PiezoHSFrequencies[med][1:]], [1e-3*c for c in PiezoHSCounts[med][1:]], '-o', markersize=2, alpha=0.7,color=palette[0])
+        plt.plot([1e3*(f-435)-fi for f in Freqs],[l*1e-3 for l in Lorentzian(Freqs,*popt)],linewidth=5,color=palette[0])
     
     jj=jj+1
-plt.xlabel('Frequency (kHz)',fontname='STIXgeneral',fontsize=50)
-plt.ylabel('kCounts',fontname='STIXgeneral',fontsize=50)
+plt.xlabel('IR1 Frequency (kHz)',fontname='STIXgeneral',fontsize=40)
+plt.ylabel('kCounts',fontname='STIXgeneral',fontsize=40)
 plt.xticks([-50,-25,0,25,50],fontname='STIXgeneral',fontsize=35)
 plt.yticks([3.5,4,4.5],fontname='STIXgeneral',fontsize=35)
 plt.tight_layout()
 plt.grid()
 plt.savefig('/home/nico/Nextcloud/Nico/Doctorado/Charlas/2023 Europe/DDresonancesexperimental_fine.pdf')
 # plt.legend()
-#plt.title(f'Ancho: {round(1e3*popt[3],2)} kHz')
+print(f'Ancho: {round(1e3*popt[3],2)} kHz')