agrego cambios a ion statistics

analisis/plots/20221017_IonStatistics/IonStatistics.py

@@ -6,7 +6,7 @@ import re
 import ast
 from scipy.optimize import curve_fit
 import os
-
+import scipy.stats as sts
 # Solo levanto algunos experimentos
 Stat_files = [8731]
 
@@ -53,8 +53,17 @@ for i, fname in enumerate(Stat_files):
 #plt.figure()
 #plt.plot(Stat_Bins[0][:-1], Stat_Heigths[0])
 
+bins = np.arange(100,350, 1)
+
 plt.figure()
-plt.hist(Stat_Heigths[0], bins=np.arange(100,350, 1), histtype='step')
+plt.hist(Stat_Heigths[0], bins=bins, histtype='step',density = True)
+
+poisson = sts.poisson.pmf(bins,np.mean(Stat_Heigths))
+
+plt.plot(bins+0.5,poisson)
+
+poissoneidad =  np.var(Stat_Heigths)/np.mean(Stat_Heigths)
+plt.title('Varianza/media = {:.4f}'.format(poissoneidad))