Transverse momentum ($p_T$) and rapidity ($y$) spectra of $Z$ bosons and quarkonium states (some charmonium $c\overline{c}$ mesons such as $J/\psi$ and $\psi \left(2S\right)$ and some bottomonium $b\overline{b}$ mesons such as ${\rm Y}\left(1S\right)$, ${\rm Y}\left(2S\right)$, and ${\rm Y}\left(3S\right)$) produced in proton-proton (pp) and lead-lead (Pb-Pb) collisions at the large hadron collider (LHC) are uniformly described by a hybrid model of two-component Erlang distribution for $p_T$ spectrum and two-component Gaussian distribution for $y$ spectrum. The former distribution results from a multisource thermal model, and the latter one results from the revised Landau hydrodynamic model. The modelling results are in agreement with the experimental data measured in pp collisions at center-of-mass energies $\sqrt{s}=2.76$ and 7 TeV and in Pb-Pb collisions at center-of-mass energy per nucleon pair $\sqrt{s_{NN}}=2.76$ TeV. Based on the parameter values extracted from $p_T$ and $y$ spectra, the event patterns (particle scatter plots) in two-dimensional $p_T$-$y$ space and in three-dimensional velocity space are obtained.