meds de rde invariante de escala

analisis/plots/20230912_RotationalDopplerShift_v6/RDS_piezobeamsizes.py

@@ -551,16 +551,16 @@ if depthscurve:
     plt.grid()
     #plt.axvline(3, color='salmon')
     plt.legend()
-
+#%%
     plt.figure()
-    #plt.plot(np.arange(0,len(Intensityver1),1), [i/np.max(Intensityver1) for i in Intensityver1], '-o',markersize=8)
-    #plt.errorbar(np.arange(0,len(Intensityver1),1), [p for p in pmdepthsdrver1], yerr= errorpmdepthsdrver1, fmt='o', capsize=3, markersize=8)
-    plt.plot(np.arange(0,len(Intensityver2),1), [i/np.max(Intensityver2) for i in Intensityver2], '-o',markersize=8)
-    plt.errorbar(np.arange(0,len(Intensityver2),1), [p for p in pmdepthsdrver2], yerr= errorpmdepthsdrver2, fmt='o', capsize=3, markersize=8)
+    # plt.plot(np.arange(0,len(Intensityver1),1), [i/np.max(Intensityver1) for i in Intensityver1], '-o',markersize=8)
+    # plt.errorbar(np.arange(0,len(Intensityver1),1), [p for p in pmdepthsdrver1], yerr= errorpmdepthsdrver1, fmt='o', capsize=3, markersize=8)
+    # plt.plot(np.arange(0,len(Intensityver2),1), [i/np.max(Intensityver2) for i in Intensityver2], '-o',markersize=8)
+    # plt.errorbar(np.arange(0,len(Intensityver2),1), [p for p in pmdepthsdrver2], yerr= errorpmdepthsdrver2, fmt='o', capsize=3, markersize=8)
     plt.plot(np.arange(0,len(Intensityver3),1), [i/np.max(Intensityver3) for i in Intensityver3], '-o',color='darkgreen',markersize=8)
     plt.errorbar(np.arange(0,len(Intensityver3),1), [p for p in pmdepthsdrver3], yerr= errorpmdepthsdrver3, color='darkgreen',fmt='o', alpha=0.7, capsize=3, markersize=8)
-    #plt.plot([m for m in np.arange(0,len(Intensityver4),1)], [i/np.max(Intensityver4) for i in Intensityver4], '-o',color='navy',markersize=8)
-    #plt.errorbar([m for m in np.arange(0,len(Intensityver4),1)], [p for p in pmdepthsdrver4], yerr=list(errorpmdepthsdrver4), fmt='o', color='navy', alpha=0.7,capsize=3, markersize=8)
+    plt.plot([m for m in np.arange(0,len(Intensityver4),1)], [i/np.max(Intensityver3) for i in Intensityver4], '-o',color='navy',markersize=8)
+    plt.errorbar([m for m in np.arange(0,len(Intensityver4),1)], [p for p in pmdepthsdrver4], yerr=list(errorpmdepthsdrver4), fmt='o', color='navy', alpha=0.7,capsize=3, markersize=8)
     #plt.xlim(-0.5, 12.7)
     plt.ylim(-0.1,1.1)
     plt.grid()
@@ -657,7 +657,7 @@ 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.savefig('/home/nico/Nextcloud/Nico/Doctorado/Charlas/2023 Europe/DDresonancesexperimental_fine.pdf')
 # plt.legend()
 print(f'Ancho: {round(1e3*popt[3],2)} kHz')
 

analisis/plots/20230912_RotationalDopplerShift_v6/RDS_piezobeamsizes_vert.py

@@ -447,13 +447,18 @@ if plotcurvita:
 
 #%%
 
+a=0.3
 
 plt.figure()
-plt.plot(np.arange(0,len(Intensityver1),1), [i/np.max(Intensityver1) for i in Intensityver1], '-o', color='navy',alpha=0.5,markersize=8)
-plt.errorbar(np.arange(0,len(Intensityver1),1), [p for p in pmdepthsdrver1], yerr= errorpmdepthsdrver1, fmt='o', color='navy',alpha=1, capsize=3, markersize=8)
+plt.plot([x*a for x in np.arange(0,len(Intensityver1),1)], [i/np.max(Intensityver2) for i in Intensityver1], '-o', color='navy',alpha=0.5,markersize=8)
+#plt.plot(np.arange(0,len(Intensityver1),1), [i/np.max(Intensityver1) for i in Intensityver1], '-o', color='navy',alpha=0.5,markersize=8)
+plt.errorbar([x*a for x in np.arange(0,len(Intensityver1),1)], [p for p in pmdepthsdrver1], yerr= errorpmdepthsdrver1, fmt='o', color='navy',alpha=1, capsize=3, markersize=8)
+
+
 plt.plot(np.arange(0,len(Intensityver2),1), [i/np.max(Intensityver2) for i in Intensityver2], '-o',color='firebrick',markersize=8,alpha=0.5)
+#plt.plot(np.arange(0,len(Intensityver2),1), [i/np.max(Intensityver2) for i in Intensityver2], '-o',color='firebrick',markersize=8,alpha=0.5)
 plt.errorbar(np.arange(0,len(Intensityver2),1), [p for p in pmdepthsdrver2], yerr= errorpmdepthsdrver2, fmt='o', color='firebrick',alpha=1,capsize=3, markersize=8)
-plt.xlim(-0.5, 36)
+plt.xlim(-0.5, 10)
 plt.ylim(-0.1,1.1)
 plt.grid()
 
@@ -960,3 +965,26 @@ plt.ylim(-0.1,1.1)
 plt.grid()
 
 
+
+#%%
+
+"""
+Analuisis posterior
+"""
+
+k = 0.4
+#ahora la anterior con haz grande y alguna de las 4 potencias con haz chico
+plt.figure()
+plt.plot([k*f for f in np.arange(0,len(Intensityver1),1)], [i/np.max(Intensityver31) for i in Intensityver1], '-o', color='navy',alpha=0.5,markersize=8)
+plt.errorbar([k*f for f in np.arange(0,len(Intensityver1),1)], [p for p in pmdepthsdrver1], yerr= errorpmdepthsdrver1, fmt='o', color='navy',alpha=1, capsize=3, markersize=8)
+# plt.plot(np.arange(0,len(Intensityver2),1), [i/np.max(Intensityver2) for i in Intensityver2], '-o', color='navy',alpha=0.5,markersize=8)
+# plt.errorbar(np.arange(0,len(Intensityver2),1), [p for p in pmdepthsdrver2], yerr= errorpmdepthsdrver2, fmt='o', color='navy',alpha=1, capsize=3, markersize=8)
+plt.plot(np.arange(0,len(Intensityver31),1), [i/np.max(Intensityver31) for i in Intensityver31], '-o',color='green',markersize=8,alpha=0.5)
+plt.errorbar(np.arange(0,len(Intensityver31),1), [0.8*p for p in pmdepthsdrver31], yerr= errorpmdepthsdrver31, fmt='o', color='green',alpha=1,capsize=3, markersize=8)
+plt.xlim(-0.5, 10)
+plt.ylim(-0.1,1.1)
+plt.grid()
+
+
+
+

analisis/plots/20230912_RotationalDopplerShift_v7/teoryanalysis.py

@@ -81,8 +81,8 @@ fact = (np.max(LG1)/np.max(LG2))
 
 
 plt.figure()
-plt.plot(r/w1,LG1)
-plt.plot(r/w2,LG2*fact)
+plt.plot(r,LG1)
+#plt.plot(r/w2,LG2*fact)
 plt.title('Scaled ratio between LG and RDE, x axis scaled with LG waist')
 
 

analisis/plots/20231214_CPTconmicromocioncristals/Data/EITfit/MM_eightLevel_2repumps_python_scripts.py

@@ -129,7 +129,7 @@ def LtempCalculus(beta, drivefreq, forma=1):
     if forma==2:
         deltaKro = np.diag([1, 1, 1, 1, 1, 1, 1, 1])
         Ltemp = (-1j)*(np.kron(Hint, deltaKro) - np.kron(deltaKro, Hint))
-
+    print(np.matrix(Ltemp))
     Omega = np.zeros((64, 64), dtype=np.complex_)
     
     for i in range(64):
@@ -399,7 +399,10 @@ if __name__ == "__main__":
     freqMax = 50
     freqStep = 5e-2
     
-    Frequencyvector, Fluovector = CPTspectrum8levels_MM(rabG, rabR, rabP, gPS, gPD, Detg, Detr, u, lwg, lwr, lwp, Temp, alpha, phidoppler, titadoppler, phiprobe, titaprobe, phirepump, titarepump, freqMin=freqMin, freqMax=freqMax, freqStep=freqStep, plot=plotCPT, solvemode=1)
+        
+    Freq, Fluo = CPTspectrum8levels_MM(sg, sp, gPS, gPD, Detg, u, lwg, lwp, Temp, alpha, phidoppler, titadoppler, phiprobe, titaprobe, Circularityprobe, beta, drivefreq, freqMin=-100, freqMax=100, freqStep=1e-1, plot=False, solvemode=1)
+    
+    #Frequencyvector, Fluovector = CPTspectrum8levels_MM(rabG, rabR, rabP, gPS, gPD, Detg, Detr, u, lwg, lwr, lwp, Temp, alpha, phidoppler, titadoppler, phiprobe, titaprobe, phirepump, titarepump, freqMin=freqMin, freqMax=freqMax, freqStep=freqStep, plot=plotCPT, solvemode=1)
     
     plt.plot(Frequencyvector, [100*f for f in Fluovector], label=str(titaprobe) + 'º, T: ' + str(Temp*1e3) + ' mK')
     plt.xlabel('Probe detuning (MHz)')

analisis/plots/20231214_CPTconmicromocioncristals/Data/EITfit/lolo_modelo_full_8niveles.py

@@ -267,6 +267,7 @@ def LtempCalculus(beta:float, drivefreq:float, forma=1):
     for i in range(64):
         Omega[i, i] = (1j)*drivefreq
 
+
     return Ltemp, Omega
 
 # LtempCalculus(0,1)

analisis/plots/20231226_CPTconmicromocion4/Data/EITfit/MM_eightLevel_2repumps_CPTPlotter.py

@@ -67,7 +67,7 @@ sp = 8 #nice range: 0.1 to 20 #p is for probe but is the repump
 
 drivefreq=2*np.pi*22.135*1e6 #ignore it
 #betavec = np.arange(0,1.1,0.1) #ignore it
-betavec=[0] #ignore it
+betavec=[10] #ignore it
 alphavec = [0] #ignore it
 
 
@@ -81,7 +81,7 @@ for sg in sgvec:
     for T in TempVec:
         for alpha in alphavec:
             for beta in betavec:
-                Frequencies, Fluorescence = PerformExperiment_8levels(sg, sp, gPS, gPD, DetDoppler, u, DopplerLaserLinewidth, ProbeLaserLinewidth, T, alpha, phidoppler, titadoppler, phiprobe, titaprobe,  beta, drivefreq, freqMin, freqMax, freqStep, circularityprobe=CircPr, plot=False, solvemode=1, detpvec=None)
+                Frequencies, Fluorescence = PerformExperiment_8levels_MM(sg, sp, gPS, gPD, DetDoppler, u, DopplerLaserLinewidth, ProbeLaserLinewidth, T, alpha, phidoppler, titadoppler, phiprobe, titaprobe,  beta, drivefreq, freqMin, freqMax, freqStep, circularityprobe=CircPr, plot=False, solvemode=1, detpvec=None)
                 
                 FrequenciesVec.append(Frequencies)
                 FluorescencesVec.append(Fluorescence)

analisis/plots/20231226_CPTconmicromocion4/Data/EITfit/MM_eightLevel_2repumps_python_scripts.py

@@ -136,6 +136,8 @@ def LtempCalculus(beta, drivefreq, forma=1):
         Omega[i, i] = (1j)*drivefreq
     
         
+    #print(np.allclose(np.matrix(Ltemp), np.matrix(Ltemp).T, rtol=1e-5, atol=1e-8))
+    print(np.max(np.matrix(Ltemp)))
     return np.matrix(Ltemp), np.matrix(Omega)