function [intts] = timewindow(tsin,n)
%% generates intensities of photon count streams

length = uint32(size(tsin,1));
if n>length
    n=length;
end

if length<1000
% n ist Zeitfenster
% tsin ist die Photocounts Zeitreihe
% Erstellen der Multiplikationsmatrix für kurze Zeitreihen
sumMatrix = zeros(length-n+1,length);
   for i=1:((length)-n)
       sumVec = [zeros(1,(i-1)) ones(1,n) zeros(1,length-i-n+1)];
       sumMatrix(i,:) = sumVec;
   end

% neue Intensitäten durch multiplizieren   
intts(:,1) = sumMatrix*tsin(:,1);
intts(:,2) = sumMatrix*tsin(:,2);


% Gegen Speicherüberlauf auf kosten der Rechenzeit
else
    intts=zeros(length-n+1,2);
    for i=0:length-n
        intts(i+1,1) = sum(tsin(i+1:n+i,1));
        intts(i+1,2) = sum(tsin(i+1:n+i,2));
    end
end
 clear sumMatrix
clear sumVec

