Statistical properties of neuron firing are studied in the framework of a nonlinear leaky integrate-and-fire model that is driven by a slow periodic subthreshold signal. The firing events are characterized by first passage time densities. The experimentally better accessible interspike interval density generally depends on the sojourn times in a refractory state of the neuron. This aspect is not part of the integrate-and-fire model and must be modelled additionally. For a large class of refractory dynamics, a general expression for the interspike interval density is given and further evaluated for the two cases with an instantaneous resetting (i.e. no refractory state) and a refractory state possessing a deterministic lifetime. First passage time densities and interspike interval densities following from the proposed theory compare favorably with precise numerical simulations.