>>> import numpy as np
>>> from astropy import units as u
Create a quantity from a pre-existing sampled distribution. The MC direction is the first dimension.
>>> parr = np.random.poisson([1, 5, 30, 400],(1000, 4))
>>> distr = u.Distribution(parr, u.kpc)