We study the distribution Ex[exp(-q∫0t 1(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.
"Joint distribution of a spectrally negative Lévy process and its occupation time, with step option pricing in view." Adv. in Appl. Probab. 48 (1) 274 - 297, March 2016.