  We introduce the numbers of hot and cold spots, $n_h$ and $n_c$, of excursion
sets of the CMB temperature anisotropy maps as statistical observables that can
discriminate different non-Gaussian models. We numerically compute them from
simulations of non-Gaussian CMB temperature fluctuation maps. The first kind of
non-Gaussian model we study is the local type primordial non-Gaussianity. The
second kind of models have some specific form of the probability distribution
function from which the temperature fluctuation value at each pixel is drawn,
obtained using HEALPIX. We find the characteristic non-Gaussian deviation
shapes of $n_h$ and $n_c$, which is distinct for each of the models under
consideration. We further demonstrate that $n_h$ and $n_c$ carry additional
information compared to the genus, which is just their linear combination,
making them valuable additions to the Minkowski Functionals in constraining
non-Gaussianity.
