actualizacion codigo modulacion rf

artiq_master/repository/Tests/Fluorescence_correlation.py

@@ -190,9 +190,13 @@ class FluorescenceCorrelation(EnvExperiment):
 
         ## Datos
         self.set_dataset("time_stamps", [], broadcast=True, archive=True)
-        self.set_dataset("tiem_stamps_corrected", [], broadcast=True, archive=True)
+        #self.set_dataset("time_stamps_corrected", [], broadcast=True, archive=True)
         self.set_dataset("counts", [], broadcast=True, archive=True)
-        self.set_dataset("counts_err", [], broadcast=True, archive=True)
+
+        # Datos postprocesado
+        self.set_dataset("fit_params", [], broadcast=True, archive=True)
+        self.set_dataset("x_data", [], broadcast=True, archive=True)
+        self.set_dataset("y_data", [], broadcast=True, archive=True)
 
         # Parámetros
         self.set_dataset("f_rf", self.f_rf, broadcast=False, archive=True)
@@ -543,12 +547,15 @@ class FluorescenceCorrelation(EnvExperiment):
 
         x_extended = np.concatenate([bin_edges[:-1] + i * self.tau_rf for i in range(3)])
         y_extended = np.tile(hist, 3)
-        x_extended_dense = np.arange(np.min(x_extended),np.max(x_extended),0.1*(x_extended[1]-x_extended[0]))
+        # x_extended_dense = np.arange(np.min(x_extended),np.max(x_extended),0.1*(x_extended[1]-x_extended[0]))
 
         params, _ = curve_fit(self.sinusoidal, x_extended, y_extended, p0=[max(y_extended), 2*np.pi/self.tau_rf, 0, np.mean(y_extended)])
-        # amps.append(np.abs(params[0]))
 
-        return x_extended, y_extended, params
+        self.append_to_dataset("fit_params", params )
+        self.append_to_dataset("x_data", x_extended )
+        self.append_to_dataset("y_data", y_extended )
+
+        return 
     
     @rpc
     def run(self):