data1 = load('10hz.txt');

data1bin = hist(data1,0:25);

[data2 y] = poissonsim(11.4,100);



data = [data1 data2'];

hist(data,0:25);

hold on;

errorbar(0:25,data1bin,sqrt(data1bin),'.');

xlim([0 25]);

grid on;

xlabel('Frequency Bins');

ylabel('Number of Events');

title('Poisson data experimental (left); simulated (right)');