tuti

0a7cb61a · Nicolas Nunez Barreto · a92d95f5 · 0a7cb61a · 0a7cb61a · 0a7cb61a
Commit 0a7cb61a authored Apr 03, 2024 by Nicolas Nunez Barreto
7 changed files
--- a/analisis/plots/20230804_RotationalDopplerShift_v2/RDS_CPT2.py
+++ b/analisis/plots/20230804_RotationalDopplerShift_v2/RDS_CPT2.py
@@ -12,8 +12,8 @@ from scipy import interpolate

 #C:\Users\Usuario\Documents\artiq\artiq_experiments\analisis\plots\20220106_CPT_DosLaseres_v08_TISA_DR\Data

-os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20230804_RotationalDopplerShift_v2/Data')
-
+#os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20230804_RotationalDopplerShift_v2/Data')
+os.chdir('C://Users//nicon//Doctorado//artiq_experiments//analisis//plots//20230804_RotationalDopplerShift_v2//Data')


 """
@@ -201,6 +201,7 @@ plt.title('Variando potencia de IR2 para potencia de IR1 fija')
 #%%

 """
+FIG PARA PAPER
 Comparo resonancias DD con haces OAM en config colineal con desplazada
 """

@@ -270,7 +271,7 @@ plt.xticks([425,430,435,440,445,450,455],fontname='STIXgeneral',fontsize=15)
 plt.ylim(1.,3.2)
 plt.grid()
 plt.tight_layout()
-plt.savefig('/home/nico/Nextcloud/Nico/Doctorado/Charlas/2023 Europe/DDresonancesexperimental.pdf')
+#plt.savefig('/home/nico/Nextcloud/Nico/Doctorado/Charlas/2023 Europe/DDresonancesexperimental.pdf')
 #plt.legend()
 #plt.title('Comparacion de config colineal (insensible) con desplazada (sensible)')


--- a/analisis/plots/20230804_RotationalDopplerShift_v2/RDS_location.py
+++ b/analisis/plots/20230804_RotationalDopplerShift_v2/RDS_location.py
@@ -12,7 +12,7 @@ from scipy import interpolate

 #C:\Users\Usuario\Documents\artiq\artiq_experiments\analisis\plots\20220106_CPT_DosLaseres_v08_TISA_DR\Data

-os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20230804_RotationalDopplerShift_v2/Data')
+#os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20230804_RotationalDopplerShift_v2/Data')




--- a/analisis/plots/20230815_RotationalDopplerShift_v3/RDS_location2.py
+++ b/analisis/plots/20230815_RotationalDopplerShift_v3/RDS_location2.py
@@ -12,7 +12,7 @@ from scipy import interpolate

 #C:\Users\Usuario\Documents\artiq\artiq_experiments\analisis\plots\20220106_CPT_DosLaseres_v08_TISA_DR\Data

-os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20230815_RotationalDopplerShift_v3/Data')
+os.chdir('C://Users//nicon//Doctorado//artiq_experiments//analisis//plots//20230815_RotationalDopplerShift_v3//Data')




--- a/analisis/plots/20230817_RotationalDopplerShift_v5/RDS_piezo.py
+++ b/analisis/plots/20230817_RotationalDopplerShift_v5/RDS_piezo.py
@@ -12,9 +12,9 @@ from scipy import interpolate

 #C:\Users\Usuario\Documents\artiq\artiq_experiments\analisis\plots\20220106_CPT_DosLaseres_v08_TISA_DR\Data

-os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20230817_RotationalDopplerShift_v5/Data')
-
-
+#os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20230817_RotationalDopplerShift_v5/Data')
+os.chdir('C://Users//nicon//Doctorado//artiq_experiments//analisis//plots//20230817_RotationalDopplerShift_v5//Data')
+#os.chdir('C:\\Users\\nicon\\Nextcloud\\G_liaf\\Publicaciones\\Papers\\2024 Rotational Doppler Effect\\Figuras\\Fig1\\Data')
 """
 en este codigo ploteo espectros CPT de resonancias D-D para configuracion +2/+2 y +2/-2 (usando pentaprisma)
 """
@@ -453,7 +453,7 @@ plt.ylim(-0.1,1.1)
 plt.grid()
 #plt.axvline(3, color='salmon')
 plt.tight_layout()
-plt.savefig('/home/nico/Nextcloud/G_liaf/Publicaciones/Papers/2024 Rotational Doppler Effect/Figuras/Material/fig3_b.pdf')
+#plt.savefig('/home/nico/Nextcloud/G_liaf/Publicaciones/Papers/2024 Rotational Doppler Effect/Figuras/Material/fig3_b.pdf')


 #%%
@@ -472,7 +472,7 @@ plt.grid()
 #plt.axvline(3, color='salmon')
 plt.legend(prop={'family':'STIXgeneral', 'size':7},loc='upper center',fontsize=5,markerscale=1)
 plt.tight_layout()
-plt.savefig('/home/nico/Nextcloud/G_liaf/Publicaciones/Papers/2024 Rotational Doppler Effect/Figuras/Material/fig3_c.pdf')
+#plt.savefig('/home/nico/Nextcloud/G_liaf/Publicaciones/Papers/2024 Rotational Doppler Effect/Figuras/Material/fig3_c.pdf')


 #%%

--- a/analisis/plots/20230817_RotationalDopplerShift_v5/RDS_piezoazimut.py
+++ b/analisis/plots/20230817_RotationalDopplerShift_v5/RDS_piezoazimut.py
@@ -12,8 +12,8 @@ from scipy import interpolate

 #C:\Users\Usuario\Documents\artiq\artiq_experiments\analisis\plots\20220106_CPT_DosLaseres_v08_TISA_DR\Data

-os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20230817_RotationalDopplerShift_v5/Data')
-
+#os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20230817_RotationalDopplerShift_v5/Data')
+os.chdir('C://Users//nicon//Doctorado//artiq_experiments//analisis//plots//20230817_RotationalDopplerShift_v5//Data')

 """
 en este codigo ploteo espectros CPT de resonancias D-D para configuracion +2/+2 y +2/-2 (usando pentaprisma)
@@ -191,6 +191,6 @@ if FullCurve:
    #plt.axvline(3, color='salmon')
    plt.legend()

-    plt.savefig('/home/nico/Nextcloud/G_liaf/Publicaciones/Papers/2024 Rotational Doppler Effect/Figuras/Material/fig2_a.pdf')
+    #plt.savefig('/home/nico/Nextcloud/G_liaf/Publicaciones/Papers/2024 Rotational Doppler Effect/Figuras/Material/fig2_a.pdf')


--- a/analisis/plots/20230817_RotationalDopplerShift_v5/RDS_piezodirections.py
+++ b/analisis/plots/20230817_RotationalDopplerShift_v5/RDS_piezodirections.py
@@ -12,7 +12,8 @@ from scipy import interpolate

 #C:\Users\Usuario\Documents\artiq\artiq_experiments\analisis\plots\20220106_CPT_DosLaseres_v08_TISA_DR\Data

-os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20230817_RotationalDopplerShift_v5/Data')
+#os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20230817_RotationalDopplerShift_v5/Data')
+os.chdir('C://Users//nicon//Doctorado//artiq_experiments//analisis//plots//20230817_RotationalDopplerShift_v5//Data')


 """
@@ -94,10 +95,10 @@ def ErrorDRdepth(p, f, b):
    derivadap = 1/((f-b)**2)
    derivadaf = ((p-b)/((f-b)**2))**2
    derivadab = ((p-f)/((f-b)**2))**2
-    return 2*np.sqrt(derivadap*ep*ep + derivadaf*ef*ef + derivadab*eb*eb)
+    return 1*np.sqrt(derivadap*ep*ep + derivadaf*ef*ef + derivadab*eb*eb)

-def Lorentzian( x, A, B, x0, gam ):
-    return A * gam**2 / ( gam**2 + ( x - x0 )**2) + B
+def Lorentzian( x, A, B, x0, gam, C):
+    return A * gam**2 / ( gam**2 + ( x - x0 )**2) + B - C*(x-x0)

 #%%
 """
@@ -117,12 +118,11 @@ palette = sns.color_palette("tab10")
 pmlocmedvec = np.arange(0,len(PIEZOVER_FILES),1)


-#pmlocmedvec = [26]
-
+#pmlocmedvec = [1]

 plt.figure()

-bkg = np.min(PiezoVerCounts[5])
+bkg = np.mean(PiezoVerCounts[5][0:2])

 pmdepthsdrver=[]
 errorpmdepthsdrver=[]
@@ -136,7 +136,31 @@ 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)))
+
+    if med==1:
+        popt, pcov = curve_fit(Lorentzian, Freqs[0:-10], Counts[0:-10], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.5,0,-1000),(0,1e4, 436.5, 1,1000)))
+
+    elif med==2:
+        popt, pcov = curve_fit(Lorentzian, Freqs[0:-5], Counts[0:-5], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.75,0,-1000),(0,1e4, 435.85, 1,1000)))
+
+    elif med==3:
+        popt, pcov = curve_fit(Lorentzian, Freqs[0:-5], Counts[0:-5], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.5,0,-1000),(0,1e4, 436.1, 1,1000)))
+
+    elif med==5:
+        popt, pcov = curve_fit(Lorentzian, Freqs[5:], Counts[5:], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.5,0,-1000),(0,1e4, 436.1, 1,1000)))
+
+    elif med==9:
+        popt, pcov = curve_fit(Lorentzian, Freqs[:-8], Counts[:-8], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.5,0,-1000),(0,1e4, 436.1, 1,1000)))
+
+    elif med==10:
+        popt, pcov = curve_fit(Lorentzian, Freqs[20:], Counts[20:], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.5,0,-1000),(0,1e4, 436.1, 1,1000)))
+
+    elif med==17:
+        popt, pcov = curve_fit(Lorentzian, Freqs[:-9], Counts[:-9], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.5,0,-1000),(0,1e4, 436.1, 1,1000)))
+
+
+    else:
+        popt, pcov = curve_fit(Lorentzian, Freqs, Counts, p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.5,0,-1000),(0,1e4, 436.1, 1,1000)))

    pmdepthsdrver.append(1-(np.min(Lorentzian(Freqs,*popt))-bkg)/(popt[1]-bkg))
    errorpmdepthsdrver.append(ErrorDRdepth(np.min(Lorentzian(Freqs,*popt)),popt[1], bkg))
@@ -160,17 +184,18 @@ plt.grid()
 # #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()
+if len(pmlocmedvec)!=1:
+    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()

 #%%
 """
@@ -190,12 +215,12 @@ palette = sns.color_palette("tab10")
 pmlocmedvec = np.arange(0,len(PIEZODIAG_FILES),1)


-#pmlocmedvec = [0,1]
+#pmlocmedvec = [1]


 plt.figure()

-bkg = np.min(PiezoDiagCounts[5])
+bkg = np.mean(PiezoDiagCounts[5][130:140])

 pmdepthsdrdiag=[]
 errorpmdepthsdrdiag=[]
@@ -211,22 +236,31 @@ 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==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)))
+
+    if med==0:
+        popt, pcov = curve_fit(Lorentzian, Freqs, Counts, p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.7,0,-1000),(0,1e4, 435.9, 1,1000)))
+
    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)))
+        popt, pcov = curve_fit(Lorentzian, Freqs[:-5], Counts[:-5], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.5,0,-1000),(0,1e4, 435.9, 1,1000)))
+
+    elif med==2:
+        popt, pcov = curve_fit(Lorentzian, Freqs[:-25], Counts[:-25], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.5,0,-1000),(0,1e4, 436.1, 1,1000)))
+
+    elif med==3:
+        popt, pcov = curve_fit(Lorentzian, Freqs[:-4], Counts[:-4], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.7,0,-1000),(0,1e4, 436.1, 1,1000)))
 
    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)))
+        popt, pcov = curve_fit(Lorentzian, Freqs[10:-53]+Freqs[-30:-1], Counts[10:-53]+Counts[-30:-1], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.5,0,-1000),(0,1e4, 436.1, 1,1000)))
+
+    elif med==20:
+        popt, pcov = curve_fit(Lorentzian, Freqs[:-10], Counts[:-10], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.7,0,-1000),(0,1e4, 436.1, 1,1000)))
+
+    elif med==22:
+        popt, pcov = curve_fit(Lorentzian, Freqs[:-10], Counts[:-10], p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.7,0,-1000),(0,1e4, 436.1, 1,1000)))
+
        
    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,0.1), bounds=((-10000,0,435.5,0,-1000),(0,1e4, 436.1, 1,1000)))

    pmdepthsdrdiag.append(1-(np.min(Lorentzian(Freqs,*popt))-bkg)/(popt[1]-bkg))
    errorpmdepthsdrdiag.append(ErrorDRdepth(np.min(Lorentzian(Freqs,*popt)),popt[1], bkg))
@@ -251,18 +285,19 @@ plt.grid()
 # #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.plot(np.arange(0,len(Intensitydiag),1), [p for p in Anchos], '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()
+if len(pmlocmedvec)!=1:
+    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.plot(np.arange(0,len(Intensitydiag),1), [p for p in Anchos], '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()


 #%%
@@ -293,7 +328,7 @@ plt.errorbar([s*scale for s in np.arange(16,len(Intensity),1)-16], [p for p in p
 plt.xlabel('Ion position')
 plt.ylabel('Intensity / DR Relative depth')
 #plt.xticks([1,2,3,4,5])
-plt.xlim(-1,15)
+#plt.xlim(-1,15)
 plt.ylim(-0.1,1.1)
 plt.grid()
 #plt.axvline(3, color='salmon')
@@ -338,7 +373,7 @@ 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)))
+    popt, pcov = curve_fit(Lorentzian, Freqs, Counts, p0=(-200,2100,435.8,0.05,0.1), bounds=((-10000,0,435.5,0,-1000),(0,1e4, 436.1, 1,1000)))

    pmdepthsdrverd.append(1-(np.min(Lorentzian(Freqs,*popt))-bkg)/(popt[1]-bkg))
    errorpmdepthsdrverd.append(ErrorDRdepth(np.min(Lorentzian(Freqs,*popt)),popt[1], bkg))

--- a/analisis/plots/20240312_RotationalDopplerShift_news/RDS_piezobeamsizes_2024_invariance.py
+++ b/analisis/plots/20240312_RotationalDopplerShift_news/RDS_piezobeamsizes_2024_invariance.py
@@ -12,7 +12,7 @@ from scipy import interpolate

 #C:\Users\Usuario\Documents\artiq\artiq_experiments\analisis\plots\20220106_CPT_DosLaseres_v08_TISA_DR\Data

-os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20240312_RotationalDopplerShift_news/Data')
+#os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20240312_RotationalDopplerShift_news/Data')


 """
@@ -380,18 +380,21 @@ plt.figure(figsize=(3.5,3))
 ms = 4
 cs = 2

-plt.errorbar(np.arange(x1,len(Intensityver1)*20+x1,20), [p for p in pmdepthsdrver1], yerr= errorpmdepthsdrver1, fmt='o',color='red', capsize=cs, markersize=ms)
-plt.plot([x*20 for x in xchicofinal],sf([i/np.max(IntensityChico) for i in IntensityChico],11,3),linestyle='dashed',color='red',alpha=0.5)
+scale = 140/6

-plt.errorbar(np.arange(x2,len(Intensityver2)*20+x2,20), [p for p in pmdepthsdrver2], yerr= errorpmdepthsdrver2, fmt='o',color='blue', capsize=cs, markersize=ms)
-plt.plot([x*20 for x in xgrandefinal],sf([i/np.max(IntensityGrande) for i in IntensityGrande],11,3),linestyle='dashed',color='blue',alpha=0.5)

-plt.xlim(-20,600+x2+10)
+plt.errorbar([f/scale for f in np.arange(x1,len(Intensityver1)*20+x1,20)], [p for p in pmdepthsdrver1], yerr= errorpmdepthsdrver1, fmt='o',color='red', capsize=cs, markersize=ms)
+plt.plot([f/scale for f in [x*20 for x in xchicofinal]],sf([i/np.max(IntensityChico) for i in IntensityChico],11,3),linestyle='dashed',color='red',alpha=0.5)
+
+plt.errorbar([f/scale for f in np.arange(x2,len(Intensityver2)*20+x2,20)], [p for p in pmdepthsdrver2], yerr= errorpmdepthsdrver2, fmt='o',color='blue', capsize=cs, markersize=ms)
+plt.plot([f/scale for f in [x*20 for x in xgrandefinal]],sf([i/np.max(IntensityGrande) for i in IntensityGrande],11,3),linestyle='dashed',color='blue',alpha=0.5)
+
+plt.xlim(-20/scale,(600+x2+10)/scale)
 plt.grid()

 plt.xlabel('Distance to center of beam (um)',fontname='STIXgeneral',fontsize=10)
 plt.ylabel('DR Relative depth',fontname='STIXgeneral',fontsize=10)
-plt.xticks([0,100,200,300,400,500,600,700],fontname='STIXgeneral',fontsize=10)
+plt.xticks([0,5,10,15,20,25,30,35],fontname='STIXgeneral',fontsize=10)
 plt.yticks([0,0.2,0.4,0.6,0.8,1],fontname='STIXgeneral',fontsize=10)
 plt.tight_layout()
 plt.savefig('/home/nico/Nextcloud/G_liaf/Publicaciones/Papers/2024 Rotational Doppler Effect/Figuras/Material/fig3_a.pdf')