data = lab3dataconv('Tuesday_F08_writetime2');
N = length(data(:,1));
M = max(data(:,2));
data(:,2) = 0.33/M*data(:,2);

subplot(211);
plot((1:(N-1))/10000,data(1:(N-1),2));
xlabel('Time (in minutes)');
ylabel('Power in Bragg-diffracted beam (in % rel. to main beam)');
title('Storage time');

% Round 2
% 
data = lab3dataconv('Tuesday_F08_erasetime2');
N = length(data(:,1));
M = max(data(:,2));
data(:,2) = 0.33/M*data(:,2);

subplot(212);
plot((1:(N-1))/10000,data(1:(N-1),2));
xlabel('Time (in minutes)');
ylabel('Power in Bragg-diffracted beam (in % rel. to main beam)');
title('Erase time');