C plus plus implementation of Bethe-Heitler, 5D, polarized, γ → e+e- pair conversion event generator

We present the C++ implementation of an event generator for linearly-polarized gamma-ray conversions to gamma -> e(+)e(-) pairs. The generator uses the full 5D Bethe-Heisler probability density function for the generation of they gamma -> e(+)e(-) pair and for the recoil of the target nucleus. The same code is used for triplet generation. We present a verification of the model by comparing sample distributions with theoretical predictions.