dist_names = ['gamma', 'beta', 'rayleigh', 'norm', 'pareto']
for dist_name in dist_names:
    dist = getattr(scipy.stats, dist_name)
How to find probability distribution and parameters for real data? (Python 3)