The infinitely divisible distributions of positive random variables C_t, S_t and T_t with Laplace transforms in x ( 1 / \cosh \sqrt{ 2 x } )^t, (\sqrt{2 x / \sinh \sqrt{ 2 x } )^t, and ( (\tanh \sqrt{ 2 x} ) / \sqrt{ 2 x } )^t respectively are characterized for various t >0 in a number of different ways: by simple relations between their moments and cumulants, by corresponding relations between the distributions and their Levy measures, by recursions for their Mellin transforms, and by differential equations satisfied by their Laplace transforms. Some of these results are interpreted probabilistically via known appearances of these distributions for t =1 or 2 in the description of the laws of various functionals of Brownian motion and Bessel processes, such as the heights and lengths of excursions of a one-dimensional Brownian motion. The distributions of C_1 and S_2 are also known to appear in the Mellin representations of two important functions in analytic number theory, the Riemann zeta function and the Dirichlet L-function associated with the quadratic character modulo 4. Related families of infinitely divisible laws, including the gamma, logistic and generalized hyperbolic secant distributions, are derived from S_t and C_t by operations such as Brownian subordination, exponential tilting, and weak limits, and characterized in various ways.