Merge branch 'master' of https://code.df.uba.ar/nnunez/artiq_experiments

analisis/plots/20250519_CPT_bladetrap/CPT_BLADETRAP_CARO.py

+# -*- coding: utf-8 -*-
+"""
+Created on Tue May 20 10:59:18 2025
+
+@author: Carolina
+"""
+
+# -*- coding: utf-8 -*-
+"""
+Created on Sun May 18 19:12:45 2025
+
+@author: Carolina
+"""
+
+import h5py
+import matplotlib.pyplot as plt
+import numpy as np
+import sys
+import re
+import ast
+from scipy.optimize import curve_fit
+import os
+from scipy import interpolate
+
+
+"""
+Primero tengo mediciones de espectros cpt de un ion variando la tension dc_A
+"""
+
+
+#os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20231123_CPTconmicromocion3/Data/')
+#os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20250516_CPT_bladetrap/Data')
+# os.chdir(r'C:\Users\Carolina\Documents\artiq_experiments\analisis\plots\20250519_CPT_bladetrap\Data')
+
+os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20250519_CPT_bladetrap/Data')
+CPT_FILES = os.listdir()
+
+
+def SeeKeys(files):
+    for i, fname in enumerate(files.split()):
+        data = h5py.File(fname+'.h5', 'r') # Leo el h5: Recordar que nuestros datos estan en 'datasets'
+        print(fname)
+        print(list(data['datasets'].keys()))
+
+Counts = []
+Freqs = []
+AmpTisa = []
+UVCPTAmp = []
+No_measures = []
+Voltages = []
+
+
+for i, fname in enumerate(CPT_FILES):
+    print(str(i) + ' - ' + fname)
+    #print(fname)
+    data = h5py.File(fname, 'r') # Leo el h5: Recordar que nuestros datos estan en 'datasets'
+
+    # Aca hago algo repugnante para poder levantar los strings que dejamos
+    # que además tenian un error de tipeo al final. Esto no deberá ser necesario
+    # cuando se solucione el error este del guardado.
+    Freqs.append(np.array(data['datasets']['IR1_Frequencies']))
+    Counts.append(np.array(data['datasets']['counts_spectrum']))
+    #AmpTisa.append(np.array(data['datasets']['TISA_CPT_amp']))
+    print(np.array(data['datasets']['no_measures']))
+    
+    
+#%%
+
+
+# index = [7, 1, 4, 6]
+index = [6]
+amp = [0.05, 0.04, 0.03, 0.02]
+
+for i in index:
+    plt.plot([2*f*1e-6 for f in Freqs[i]],Counts[i],'o-')
+    
\ No newline at end of file