We study the distribution Ex[exp(-q∫0t1(a,b)(Xs)ds); Xt ∈ dy], where -∞ ≤ a < b < ∞, and where q, t > 0 and x ∈ R for a spectrally negative Lévy process X. More precisely, we identify the Laplace transform with respect to t of this measure in terms of the scale functions of the underlying process. Our results are then used to price step options and the particular case of an exponential spectrally negative Lévy jump-diffusion model is discussed.