tempCirc=circ_total1[:,30]
statsG.append(lister+" Tubule Diameter="+str(np.mean(tempCirc))+" stdev="+str(np.std(tempCirc)))
diameter_Glob.append(np.mean(tempCirc))
plt.hist(tempCirc,bins=100,normed=True)
plt.xlim((min(tempCirc), max(tempCirc)))
mean = np.mean(tempCirc)
variance = np.var(tempCirc)
sigma = np.sqrt(variance)
xtemp = np.linspace(min(tempCirc), max(tempCirc),100)
plt.plot(xtemp,mlab.normpdf(xtemp,mean,sigma))