Multiplicity distributions as probes of quarkonia production mechanisms