In [1]:

    

%pylab inline


    

    




/Users/hantke/Library/Python/2.7/lib/python/site-packages/matplotlib/font_manager.py:273: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.
  warnings.warn('Matplotlib is building the font cache using fc-list. This may take a moment.')






    




Populating the interactive namespace from numpy and matplotlib

In [2]:

    

import h5py,os,sys
from flash_tools import *

root_dir = os.path.abspath(".") + "/../"
sys.path.append(root_dir)
import params


    

In [19]:

    

#runs = range(110, 113+1) + [115,117,118,119]
#runs = [121,122,123]
#runs = [128,129]
#runs = range(146, 154+1)
#run = [161,164,165,166,168,169]
#runs = range(13, 16+1) + [18,19] + [21,22,23] + [34]
#runs = [43,44]
#runs = range(183,189)
#runs = [195] + range(197,201+1) + [203,205,206,207,208]
#runs = [203,205,206,207,208]
#runs = range(210, 215+1)
#runs = [60,179,180]
#runs = range(252, 267+1)
runs = [270,271]
fig, axs = subplots(len(runs), 1, figsize=(15, 20))
for ax, r in zip(axs, runs):
    #print r
    hitscores = read_hitscores(run_nr=r, sub_folder="hummingbird", ol=1)
    T = params.read_params(params_fn, run_nr=r)["hitscoreThreshold"]
    #print mean(hitscores[:5000])
    H = hitscores[len(hitscores)/2:len(hitscores)/2+1000]
    ax.plot(H, ".")
    #ax.axhline(2.*median(H), color="red")
    ax.axhline(T, color="red")
    ax.set_ylim(0, 50000)
    ax.set_title(r)
    #print median(H)


    

In [ ]: