function [Startime, EndTime,lb,ub] = Boundary(file,InstCurveNum,Start);

raw = load(file);
IRF = raw(:,InstCurveNum);
Decay = raw(:,Start);

fig = figure;
semilogy(IRF); hold on; semilogy(Decay,'g');
pause(0.001);
rectangle = getrect(fig);

Startime = round(rectangle(1));
EndTime = round(rectangle(1) + rectangle(3));



%%%%%%%%%%%%%% boundary conditions %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
AmpL = [0;0;0];               AmpU = [1;1;1];     % amplitude
LifeL = [1e-4;1e-4;1e-4];     LifeU = [5;5;10];  % lifetime
TShiftL = [-0.3];                   TShiftU = [0.3];          % time shift (ns)
BackL = [-5000;-5000];                BackU = [5000;5000];        % bkg shift
%%%%%%%%%%%%%%%%%%%   End %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

lb = [AmpL;LifeL;BackL;TShiftL]; % lower boundary
ub = [AmpU;LifeU;BackU;TShiftU]; % upper boundary