docs: def _sigma_ep(self, gam, eps):
          """
          Using sigma_1 only applies to the ultrarelativistic regime.
          Eph > 10 MeV
          ToDo: add complete e-p cross-section
          """
          with warnings.catch_warnings():
              warnings.simplefilter("ignore")
              return self._sigma_1(gam, eps)