modelSpectrum (peakT, peakI, lengthT, widthT)

function mspec = modelSpectrum (peakT, peakI, lengthT, widthT)

%----------------------------------------------------
% Model spectrum for the peak list in resampled time
% (assuming constant peak width and Lorentzian shape)
%
% INPUT:
% peakT - peak locations;
% peakI - peak intensities (maxima);
% lengthT - length of TOF record for model spectrum;
% widthT - half-width of a peak in time points;
%
% OUTPUT:
% mspec - model spectrum;
%
% USAGE:
% mspec = modelSpectrum (peakT, peakI, lengthT, widthT);
% Dependency: none
%

t1=1:lengthT; t1=t1'; tau=1/widthT;
mspec=t1*0;

for i=1:length(peakT)
    pki=lorf(t1, peakT(i), tau);
    pki=peakI(i)*pki/max(pki);
    mspec=mspec+pki;
end;


return;