an extension of the work by @Nug2002 [^2], based primarily on data from @Lev2005 for type Ibc, @DiC2002 for IIn and @Gil1999 for the remaining two types. We approximate the distribution of absolute magnitudes in the $B$ band by Gaussians with mean and scatter of $-17.50$ and 1.0 for type Ib/c, $-16.70$ and 1.0 for type IIP, $-17.9$ and 0.90 for type IIL, $-18.5$ and 1.0 (after discarding two outliers with $M_{B}=-15.1,-22.2$) for IIn [@Rich2014 for $h=0.73$]. The assumed luminosity functions for all classes of supernovae are shown in Fig. \[abs\_lum\]. We note that the borderline between different classes may not be entirely clear cut. For example, there may be a continuous transition between type IIP and type IIL supernovae [e.g., @2014ApJ...786...67A]. In this sense, these types can be seen as being representative for the fainter and more luminous ends of the type II supernova population.

Computation
-----------

First, we generate a random sample of lens galaxies at redshifts $z<1.2$. We find that it is sufficient to realize $10^{5}$ lens galaxies and then rescale the final counts of strongly lensed supernovae according to the actual number of lens galaxies contained in the assumed comoving volume. Then, we populate the observational 