cambios de mi notebook de casa

analisis/plots/20231123_CPTconmicromocion3/CPT_plotter_20231123.py

@@ -1163,7 +1163,6 @@ Temp_vec = [0.0017384071359855713, 0.0006380145223733723, 0.0007457923288975645,
 ErrorTemp_vec = [0.0004229476096548473, 0.00014439375508413987, 0.00013204015146487435, 9.307939678673377e-05, 0.000100129717662808, 0.0001841318633900307, 0.0003595040837509976, 0.0005950353892849986, 0.001866844309182069, 0.012656306714647434, 0.13143081065882864]
 
 """
-print(ErrorTemp_vec)
 
 def expo(x,tau,A,B):
     return A*np.exp(x/tau)+B
@@ -1177,7 +1176,7 @@ def InverseMicromotionSpectra(beta, A, det, x0, gamma, B):
     P = ((jv(0, beta)**2)/((((det-x0)**2)+(0.5*gamma)**2)**2))*(-2*(det-x0))
     i = 1
     #print(P)
-    while i <= 5:
+    while i <= 2:
         P = P + (-2*(det-x0))*((jv(i, beta))**2)/(((((det-x0)+i*ftrap)**2)+(0.5*gamma)**2)**2) + (-2*(det-x0))*(((jv(-i, beta))**2)/((((det-x0)-i*ftrap)**2)+(0.5*gamma)**2)**2)
         i = i + 1
         #print(P)
@@ -1191,7 +1190,7 @@ def MicromotionSpectra(beta,det, gamma):
     P = (jv(0, beta)**2)/(((det)**2)+(0.5*gamma)**2)
     i = 1
     #print(P)
-    while i <= 5:
+    while i <= 2:
         P = P + ((jv(i, beta))**2)/((((det)+i*ftrap)**2)+(0.5*gamma)**2) + ((jv(-i, beta))**2)/((((det)-i*ftrap)**2)+(0.5*gamma)**2) 
         i = i + 1
         #print(P)
@@ -1210,7 +1209,7 @@ def InverseDerivMicromotionSpectra(beta, det, gamma):
     P = ((jv(0, beta)**2)/((((det)**2)+(0.5*gamma)**2)**2))*(-2*(det))
     i = 1
     #print(P)
-    while i <= 5:
+    while i <= 2:
         P = P + (-2*(det))*((jv(i, beta))**2)/(((((det)+i*ftrap)**2)+(0.5*gamma)**2)**2) + (-2*(det))*(((jv(-i, beta))**2)/((((det)-i*ftrap)**2)+(0.5*gamma)**2)**2)
         i = i + 1
         #print(P)
@@ -1229,7 +1228,7 @@ def FinalTemp(beta,det, C,D):
     #det = -12e6
     #print(MicromotionSpectra(beta,det,gamma))
     #D=0
-    return ((hbar/(1*kb))*C*MicromotionSpectra(beta,2*np.pi*det,2*np.pi*gamma)+D*(beta**2))*InverseDerivMicromotionSpectra(beta, 2*np.pi*det, 2*np.pi*gamma)
+    return ((hbar/(kb))*C*MicromotionSpectra(beta,2*np.pi*det,2*np.pi*gamma)+D*(beta**2))*InverseDerivMicromotionSpectra(beta, 2*np.pi*det, 2*np.pi*gamma)
     #return (C*MicromotionSpectra(beta,det,gamma))*InverseDerivMicromotionSpectra(beta, det, gamma)
 
 
@@ -1258,8 +1257,8 @@ def FinalTemp_fixedall(beta, det, D):
 Temperatura vs beta con un ajuste exponencial
 """
 
-popt_rho22_balance, pcov_rho22_balance = curve_fit(FinalTemp,list(np.arange(0,0.3,0.001))+list(Betas_vec[:9]), list(0.5e-3*np.ones(len(np.arange(0,0.3,0.001))))+[t for t in Temp_vec[:9]],p0=(-10e6, 1,1e-27), bounds=((-30e6,0,0),(0e6,2,1e-25))) #esto ajusta muy bien
-#popt_rho22_balance, pcov_rho22_balance = curve_fit(FinalTemp,list(Betas_vec[:9]), [t for t in Temp_vec[:9]],p0=(-11e6, 1,1e-28), bounds=((-30e6,0,1e-30),(0e6,2,1e-25))) #esto ajusta muy bien
+#popt_rho22_balance, pcov_rho22_balance = curve_fit(FinalTemp,list(np.arange(0,0.3,0.001))+list(Betas_vec[:9]), list(0.5e-3*np.ones(len(np.arange(0,0.3,0.001))))+[t for t in Temp_vec[:9]],p0=(-10e6, 1,1e-27), bounds=((-30e6,0,0),(0e6,2,1e-25))) #esto ajusta muy bien
+popt_rho22_balance, pcov_rho22_balance = curve_fit(FinalTemp,list(Betas_vec[:9]), [t for t in Temp_vec[:9]],p0=(-11e6,0.8,4e-28)) #esto ajusta muy bien
 #popt_rho22_balance, pcov_rho22_balance = curve_fit(FinalTemp_fixedall,list(Betas_vec[:9]), [t*1e3 for t in Temp_vec[:9]],p0=(-10,1e-20), bounds=((-30,0),(0,1e6))) #esto ajusta muy bien
 
 
@@ -1291,6 +1290,28 @@ plt.ylabel('Temperature (mK)')
 plt.title('With RF heating')
 plt.grid()
 
+
+#%%
+
+hbar = 1.05e-34
+kb = 1.38e-23
+
+beta_vector = np.arange(0,1.5,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)
+    comparison.append(factor2/factor1)
+
+plt.figure()
+plt.plot(beta_vector,comparison)
+plt.xlabel('Beta')
+plt.ylabel('Ratio RF heat/Normal heat')
+plt.axhline(1)
+
+
 #%%
 """
 Calculo limite doppler para distintos detunings pero sin micromocion