meds

artiq_master/last_rid.pyon

-19363
+19584

artiq_master/repository/Tests/loadingions.py

+
 """
-Created on Fri May  3 22:37:57 2024
+Created on 25/11/2024
 
 @author: juamp y rom
 """
@@ -219,8 +220,14 @@ class InOut(EnvExperiment):
     #El famoso código de reconocimiento
     
     @rpc
-    def func(self,i,s,mu,mu0,b,c,ruido)  -> TList(TInt64):#Necesaria porque por alguna razon a artiq no le gusta que "Recon" funcione directo.
-        recon_dev=self.Recon(i,s,self.t_med,self.tseg,mu,mu0,b,c,ruido,self.paso,self.pasom,self.med,self.medmax,self.N,self.fluo,self.atrapados,self.desviacion)
+    def func(self,i,s,mu,mu0,b,c,ruido,calibrate)  -> TList(TInt64):#Necesaria porque por alguna razon a artiq no le gusta que "Recon" funcione directo.
+        if calibrate==True:#Para que no se pierda valores cuando calibra
+            mu_calib=2*ruido
+            mu0_calib=2*ruido
+        else:
+            mu_calib=mu
+            mu0_calib=mu0
+        recon_dev=self.Recon(i,s,self.t_med,self.tseg,mu_calib,mu0_calib,b,c,ruido,self.paso,self.pasom,self.med,self.medmax,self.N,self.fluo,self.atrapados,self.desviacion)
         print(recon_dev[3])
         return recon_dev
 
@@ -245,10 +252,6 @@ class InOut(EnvExperiment):
                     mus=mu0
                 else:
                     mus=mu
-#                if s[-1]<(fluo[0]+ruido):
-#                    atrapados.append(z*N[-1])
-#                    desviacion.append(.1)
-#                else:
                 a=abs(fluo[-1]-fluo[-2])
                 print(a)
                 atrapados.append(z*(1+round(abs(mus-a)/mu)))
@@ -341,7 +344,7 @@ class InOut(EnvExperiment):
                     self.append_to_dataset("counts",rta)
                     self.append_to_dataset("time",kk*self.t_tir+self.t_exp*i/int(self.t_exp/self.t_med))
 
-                    recon_dev=self.func(i,s,100,100,b,c,ruido0)
+                    recon_dev=self.func(i,s,0,0,b,c,ruido0,True)#El true esta para que calibre. El 0,0 esta para poner algo, pero el valor es irrelevante
                     b,c,fl,n=recon_dev[0],recon_dev[1],recon_dev[2],recon_dev[3]
                     if fl!=0:
                         fluoresc=fluoresc+[fl]
@@ -429,7 +432,7 @@ class InOut(EnvExperiment):
             self.append_to_dataset("counts",rta)
             self.append_to_dataset("time",self.t_exp*i/int(self.t_exp/self.t_med))
 
-            recon_dev=self.func(i,sig,mu,mu0,b,c,ruido0)
+            recon_dev=self.func(i,sig,mu,mu0,b,c,ruido0,False)#Aca NO estamos calibrando.
 
             b,c,fl,n=recon_dev[0],recon_dev[1],recon_dev[2],recon_dev[3]
             #self.append_to_dataset("fluo",fl)