agrego meds xdxd

analisis/plots/20230720_EspectrosCristal2iones/Espectroscpt_cristal.py

@@ -20,6 +20,7 @@ MOTIONAL_FILES = """000013489-IR_Scan_withcal_optimized_andor
 000013491-IR_Scan_withcal_optimized_andor
 000013493-IR_Scan_withcal_optimized_andor
 000013494-IR_Scan_withcal_optimized_andor
+000013479-IR_Scan_withcal_optimized_andor
 """
 
 
@@ -62,18 +63,35 @@ for i, fname in enumerate(MOTIONAL_FILES.split()):
 En cristal de 2 iones veo espectros cpt. 
 """
 
-medvec = [0,1,2,3,4]
+medvec = [5]
 
 plt.figure()
 
+ftrap = 22.1
 
 for med in medvec:
-    #CountsRois = [CountsRoi1[med], CountsRoi2[med]]
-    CountsRois = [CountsRoi2[med]]
+    CountsRois = [CountsRoi1[med], CountsRoi2[med]]
+    #CountsRois = [CountsRoi1[med]]    
+    #CountsRois = [CountsRoi2[med]]
     i=0
     for counts in CountsRois:
-        plt.plot(IR1_Freqs[0][1:], [c for c in counts[1:]], '-o', markersize=2)
+        plt.plot([2*f*1e-6 for f in IR1_Freqs[0][1:]], [c for c in counts[1:]], '-o', markersize=2)
         i=i+1
+    plt.axvline(415.5)
+    plt.axvline(423)
+    plt.axvline(429)
+    plt.axvline(438.5)
+    
+    plt.axvline(415.5+ftrap, color='red')
+    plt.axvline(423+ftrap, color='red')
+    plt.axvline(429+ftrap, color='red')
+    plt.axvline(438.5+ftrap, color='red')
+
+    plt.axvline(415.5-ftrap, color='red')
+    plt.axvline(423-ftrap, color='red')
+    plt.axvline(429-ftrap, color='red')
+    plt.axvline(438.5-ftrap, color='red')
+    
 plt.xlabel('Frecuencia')
 plt.ylabel('Cuentas ROI')
 #plt.xlim(0.05,0.23)