add a bgk operator with fixed (hardcoded) plasma background
Simple change, adds a new event type for a BGK operator. Instead of a rotation after a charge exchange collision, this samples a new velocity from the distribution of the ion velocities.
Note, this operator would need to access a grid of velocity and temperature values from a grid, which it does not do now. Those parameters are hardcoded for the moment.