maxTime = max(LP_Real_Event) + 5;
dt = min(diff(LP_Real_Event)) / 2
LP_Real_EventWave = 0;
current = 1;
numspikes = length(LP_Real_Event);
for i = 1 :  maxTime / dt,
  time = i * dt;
  LP_Real_EventWave(i, 1) = time;
  if (((LP_Real_Event(current) - time) >= -dt) & ((LP_Real_Event(current) - time) <= dt))
    LP_Real_EventWave(i, 2) = 1;
    if current < numspikes
      current = current + 1;
    end
  else
    LP_Real_EventWave(i, 2) = 0;
  end
end