X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=files%2Fmatplotlib_hist.py;h=f8a1105f0b2ef716594bb901f5f8caffb40e4a67;hb=633e6ce55fd8ea49c4c07f93be9dc30b9e11e229;hp=c37a5f6b7bc379f7bce3717b629960f609564d39;hpb=4edf019d31293b50d5003ea7ea02dba7f4959596;p=django_unres.git

diff --git a/files/matplotlib_hist.py b/files/matplotlib_hist.py
index c37a5f6..f8a1105 100755
--- a/files/matplotlib_hist.py
+++ b/files/matplotlib_hist.py
@@ -19,7 +19,7 @@ else:
  x,y,s,r= np.loadtxt('remd_all.stat',usecols=(7,3,0,9),unpack=True)
  x0,s0,r0= np.loadtxt('remd_all0.stat',usecols=(7,0,9),unpack=True) 
 
-hall,binall=np.histogram(y,bins='auto',density=False)
+hall,binall=np.histogram(y,bins=40,density=False)
 
 plt.xlim(min(binall), max(binall[hall>4]))
 #plt.ylim(0,max(hall)/4)
@@ -34,7 +34,7 @@ Tremd=map(float,sys.argv[1].split())
 colors = cm.rainbow(np.linspace(0, 1, len(Tremd)))
 for T,c in zip(Tremd,colors):
  yt=y[x==T]
- h,bin=np.histogram(yt,bins='auto',density=False)
+ h,bin=np.histogram(yt,bins=40,range=(min(binall),max(binall[hall>4])),density=False)
  center = (bin[:-1] + bin[1:]) / 2
  plt.plot(center,h,'-',color=c)
 # plt.bar(bin[:-1], h, width = bin[2]-bin[1],color=c)