root = "2017/05/04/20170504_182411_WireScan.h5|/"
wire = "x"
blms = ["",""]
cycles = 3


for i in range(len(blms)):
    path_back = root + "background/blm" + str(i+1)
    bg = get_attributes(path_back)["Mean"]    
    for cycle in range (cycles):    
        pos = load_data(root + wire+"_" + ("%04d" % (cycle+1)) + "/w_pos")
        sig = load_data(root + wire+"_" + ("%04d" % (cycle+1)) + "/blm" + str(j+1))         
        sig = blm_remove_spikes(sig)
        if bg is not None:
            sig = [v-bg for v in sig]        
        [com, rms] = profile_rms_stats(pos, sig,noise_std=0, n_sigma=3.5)
        #set_attribute(path, "RMS", rms)            
        #set_attribute(path, "RMS centroid", com)            
        print i, j, [com, rms]
        [off, amp, com, sigma] = profile_gauss_stats(pos, sig, off=None, amp=None, com=None, sigma=None)
        #set_attribute(path, "Gauss COM", com)            
        #set_attribute(path, "Gauss Sigma", sigma)                                            
        print i, j, [off, amp, com, sigma]