Implement q-Gaussian charge distribution
The probability density function is given here.
Apparently neither numpy
nor scipy
offer an implementation so this needs to be manually.
Beam is assumed to be q-Gaussian in all dimensions (with different q, \beta
).