Binning

This page explains the mathematics behind the binning functionality found in bins. The binning functions are purely mathematical in nature though they do have application in other parts of pyne that deal more directly physics with physics.

Pointwise Linear Collapse

The pointwise_linear_collapse() function takes an array of pointwise data \(y\) that has the independent variable \(x\) and collapses this into \(G\) groups as defined by the bin boundaries \(x_g\) where \(g\) indexes \(G\). Both \(x\) and \(x_g\) must be monotonic and in the same direction. Say that there are \(N\) points in \(x\) and \(y\). Let \(n\) index \(N\). Then for all points internal to a group, the collapsed value \(y_g\) is:

\[y_g = \frac{1}{x_{g+1} - x_g} \sum_{n|x_g \le x_n}^{x_n \le x_{g+1}} \frac{y_{n+1} + y_n}{2} * (x_{n+1} - x_n)\]

The term \((y_{n+1} + y_n)/2\) is the center (average) value of a linear interpolation between the two points. Therefore, \(y_g\) is the \(x\)-weighted average of \(y\) over the entire group.

In the event that the line between \(y_n\) and \(y_{n+1}\) crosses either the lower or upper bin boundary (or both) then their values are adjusted via a linear interpolation to the value at the bin boundary.

For a lower boundary crossing, the following substitutions are made to the equation above:

\[x_n \to x_g\]
\[y_n \to \frac{y_{n+1} - y_n}{x_{n+1} - x_n} (x_g - x_n) + y_n\]

For an upper boundary crossing:

\[x_{n+1} \to x_{g+1}\]
\[y_{n+1} \to \frac{y_{n+1} - y_n}{x_{n+1} - x_n} (x_{g+1} - x_n) + y_n\]