histograms

[Haim Edri]

Hi

I would like to create density histogram for MD simulation. I can create histogram for the number of particle distribution from a list that contains the particles' positions, but cannot simply convert N(x,t) to rho(x,t)  - the local density, since "ghistogram" counts how many particles are in each bin (the position interval).

Is there a way to convert the values of  "y axis" of ghistogram? (ideally, I would just divide it by constant - the area of the bins)

Thanks

Haim Edri

 (Weizmann Inst. Rehovot, Israel)

[Bruce Sherwood]

I can't think of a way to do this in Classic VPython. In GlowScript VPython and Jupyter VPython you can replace the heights of a gvbars object named "g" by rewriting g.data, which is what is done in this program:

http://www.glowscript.org/#/user/GlowScriptDemos/folder/Examples/program/HardSphereGas-VPython

There is no ghistogram option, but the ability to rewrite existing gvbars provides generality.

[bash...@ncsu.edu]

Note that the Classic VPython graph module is written in pure Python and is found at PythonXY/Lib/site-packages/visual_common/graph.py. You could make a copy of this and modify it to suit your needs, and include that module as part of your package.