run("Diagnostics/sig_process_wrapper")

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

stats = []


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

        samples[0].append(ć);samples[1].append(rms);samples[2].append(com);samples[3].append(sigma)
        
        from mathutils import Gaussian
        g = Gaussian(amp, com, sigma)
        gauss = [g.value(v)+off for v in pos]
        plot([data, sp, sig, gauss], ["data", "sp", "signal", "gauss", ], xdata = pos, title="Fit blm" + str(i+1) + " - " + str(cycle+1))   
    stats.append([])
    for sample in samples:
        sample =  [v for v in sample if v is not None]
        stats[i].append( (mean(sample), stdev(sample)) if len(sample)>0 else (float("nan"), float("nan")) )
print stats

msg = ""
for i in range(len(blms)):  
    msg += "BLM " + str(i)+ ":"
    msg += "    RMS:" + str(stats[i][0])
    msg += "    RMS COM:" + str(stats[i][1])
    msg += "    Gauss COM:" + str(stats[i][2])
    msg += "    Gauss Sigma:" + str(stats[i][3])

print msg