Akvo now uses the full boostrap noise algorithm from GJI paper

akvo/gui/akvoGUI.py

             qv = np.array(qv)
             iQ = np.argsort(np.array(qq))
             qq = np.array(qq)[iQ]
-            print("qq", qq, iQ)
             qv = np.array(qv)[iQ]
             INFO["Pulses"][pulse] = {}
             INFO["Pulses"][pulse]["units"] = "A"
                 INFO["Gated"][pulse]["windows"] = VectorXr( self.RAWDataProc.GATEDWINDOW ) 
                 for ichan in self.RAWDataProc.DATADICT[pulse]["chan"]:
                     INFO["Gated"][pulse]["Chan. " + str(ichan)] = {} 
-                    INFO["Gated"][pulse]["Chan. " + str(ichan)]["STD"] =  VectorXr( np.std(self.RAWDataProc.GATED[ichan]["NR"], axis=0) )
+                    #INFO["Gated"][pulse]["Chan. " + str(ichan)]["STD"] =  VectorXr( np.std(self.RAWDataProc.GATED[ichan]["NR"], axis=0) )
+                    INFO["Gated"][pulse]["Chan. " + str(ichan)]["STD"] = VectorXr( np.average(self.GATED[chan]["BN"], axis=0) )
                     for ipm in range(self.RAWDataProc.DATADICT["nPulseMoments"]):     
                         INFO["Gated"][pulse]["Chan. " + str(ichan)]["Q-"+str(ipm) + " CA"] = VectorXr(self.RAWDataProc.GATED[ichan]["CA"][ipm])   
                         INFO["Gated"][pulse]["Chan. " + str(ichan)]["Q-"+str(ipm) + " RE"] = VectorXr(self.RAWDataProc.GATED[ichan]["RE"][ipm])