data = [6176, 11046, 670, 6146, 7945, 6864, 767, 7623, 7212, 9040, 3213, 
    6302, 10044, 10195, 9386, 7230, 4602, 6282, 8619, 7903, 6318, 13294, 
    6990, 5515, 9157]
mean_of_distribution = numpy.mean(data)
Fitting a gamma distribution with (python) Scipy