2

analisis/plots/20230912_RotationalDopplerShift_v7/teoryanalysis.py

@@ -26,6 +26,8 @@ amp2=2*np.pi/w2
 plt.figure()
 plt.plot(r,1*OAM(r,w1,amp1))
 plt.plot(r,1*OAM(r,w2,amp2))
+plt.plot(r,1*OAM(r,w2*1.5,amp2*1.5))
+plt.plot(r,1*OAM(r,w2*1.7,amp2*1.7))
 #%%
 plt.figure()
 plt.plot(r,OAM(r,w1,amp1)/RDE(r,2,2))

analisis/plots/20231123_CPTconmicromocion3/CPT_plotter_20231123.py

@@ -1355,8 +1355,8 @@ beta_vector = np.arange(0,2,0.01)
 comparison = []
 
 for beta_target in beta_vector:
-    factor1 = popt_rho22_balance[1]*(hbar/kb)*MicromotionSpectra(beta_target,2*np.pi*popt_rho22_balance[0],2*np.pi*23e6)
-    factor2 = popt_rho22_balance[2]*(beta_target**2)
+    factor1 = 1*(hbar/kb)*MicromotionSpectra(beta_target,2*np.pi*popt_rho22_balance[0],2*np.pi*23e6)
+    factor2 = popt_rho22_balance[1]*(beta_target**2)
     comparison.append(factor2/factor1)
 
 plt.figure()