X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=files%2Fmatplotlib_hist.py;h=f8a1105f0b2ef716594bb901f5f8caffb40e4a67;hb=3cd10cb73439c774e87c98f6873b12eff8109ffa;hp=2c54ebedd2dc538a7897e555830e3d7b3e3c0d26;hpb=a6e756bc874c05ea14f7087df64561ec1d8a98bc;p=django_unres.git diff --git a/files/matplotlib_hist.py b/files/matplotlib_hist.py index 2c54ebe..f8a1105 100755 --- a/files/matplotlib_hist.py +++ b/files/matplotlib_hist.py @@ -22,7 +22,7 @@ else: hall,binall=np.histogram(y,bins=40,density=False) plt.xlim(min(binall), max(binall[hall>4])) -plt.ylim(0,max(hall)/4) +#plt.ylim(0,max(hall)/4) plt.ylabel('number of samples') plt.xlabel('potential energy [kcal/mol]') @@ -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=40,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)