ScriptPlotInteractionForcesHistogram

<source lang="python"> import pylab # the matlab-like interface of matplotlib import numpy

1. run uniax.py to get this file
2. or just get some data from TriaxialTest, whatever

O.load('/tmp/uniax-tension.xml.bz2')

1. axis normal to the plane in which we do the histogram

axis=0 # x, i.e. plot the yz plane ax1,ax2=(axis+1)%3,(axis+2)%3 ## get the other two indices, i.e. 1 and 2 in this case angles,forces=[],[] for i in O.interactions: if not i.isReal: continue norm=i.geom['normal'] angle=atan(norm[ax2]/norm[ax1]) force=i.phys['normalForce'].Length() angles.append(angle) forces.append(force)

1. easier: plain histogram

pylab.hist(angles,weights=forces,bins=20) pylab.show()

1. polar histogram
   1. prepare data

values,bins=numpy.histogram(angles,weights=forces,bins=20)

1. 1. prepare polar plot

pylab.subplot(111,polar=True);

1. 1. plot bar chart, with the histogram data
      1. bins has one edge extra, remove it: [:-1]

pylab.bar(left=bins[:-1],height=values,width=.7*pi/20); pylab.show() </source>