Подгонка функции Гаусса

У меня есть гистограмма (, см. ниже ), и я пытаюсь найти среднее значение и стандартное отклонение вместе с кодом, который соответствует кривой моей гистограмме. Я думаю, что в SciPy или matplotlib есть что-то, что может помочь, но каждый пример, который я пробовал, не работает.

import matplotlib.pyplot as plt
import numpy as np

with open('gau_b_g_s.csv') as f:
    v = np.loadtxt(f, delimiter= ',', dtype="float", skiprows=1, usecols=None)

fig, ax = plt.subplots()

plt.hist(v, bins=500, color='#7F38EC', histtype='step')

plt.title("Gaussian")
plt.axis([-1, 2, 0, 20000])

plt.show()