function htt = makechirplet(simId, ampl, startSamp, sampFreq)

[id, f0, tau, durationSec, alpha] = simparse(simId);
if isempty(alpha)
    alpha = 0;
end

% Deterimine length of time vector in samples (floor trucated)
durationSamp = floor(durationSec*sampFreq);
% Check that durationSamp is not < 1
if (durationSamp < 1)
    msgId = 'makeh:durationSamp';
    error(msgId,'%s: %d',msgId,durationSamp);
end
% Time vector = times at which IFO samples waveform
ti = ((1:durationSamp)./sampFreq)+((1-mod(startSamp,1))/sampFreq);
% N.B. This accounts for intersample start times

% Generate the Gaussian waveform that is also the base for SG and CG
% Offset times for Gaussians
tg = ti-(ti(end)/2);
% Generate Gaussian
gauss = exp(-(tg/tau).^2);

htt(1,:) = (ampl*(2*pi*tau^2)^(-1/4)).*exp(-(tg).^2/(4*tau^2)).*cos(2*pi*f0.*(tg)+(alpha.*(tg).^2)/(4*tau.*tg.^2));
htt(2,:) = (ampl*(2*pi*tau^2)^(-1/4)).*exp(-(tg).^2/(4*tau^2)).*sin(2*pi*f0.*(tg)+(alpha.*(tg).^2)/(4*tau.*tg.^2));

return