Indirect Detection with m_DM < 100 GeV not computing spectrum

Hi,

I am currently trying to use MadDM to do indirect detection calculations for WIMPs, generating spectra, etc. MadDM works well and calculates spectra when the masses are above 100 GeV, but as soon as the dark matter masses are below 100 GeV, a Pythia error shows up.

 -------- End LHA initialization information --------
mDM: 9.1000e+01 outdir: ./
 PYTHIA Error in Pythia::check: energy-momentum not conserved
 PYTHIA Error in Pythia::next: check of event revealed problems
 PYTHIA Abort from Pythia::next: parton+hadronLevel failed; giving up
Event generation aborted prematurely, owing to error!

*------- PYTHIA Error and Warning Messages Statistics ----------------------------------------------------------*
 | |
 | times message |
 | |
 | 10 Abort from Pythia::next: parton+hadronLevel failed; giving up |
 | 100 Error in Pythia::check: energy-momentum not conserved |
 | 100 Error in Pythia::next: check of event revealed problems |
 | 1 Error in SigmaTotal::calc: too low energy |
 | 1 Info from SLHAinterface::initSLHA: No MODSEL found, keeping internal SUSY switched off |
 | 1 Info from SLHAinterface::initSLHA: importing MASS entries |
 | 1 Warning in SLHAinterface::initSLHA: ignoring DECAY tables |
 | 1 Warning in SLHAinterface::initSLHA: ignoring MASS entries |
 | 1 Warning in SLHAinterface::initSLHA: ignoring empty DECAY tables |
 | 1 Warning in SLHAinterface::initSLHA: massless particle forced stable |
 | |
 *------- End PYTHIA Error and Warning Messages Statistics ------------------------------------------------------*

Specifically, it seems to be the "energy-momentum not conserved" and "too low energy" that are causing the issues.

The gamma ray spectrum data file reads (only copying the first few lines, but the rest also reads nan for the second column):

 -8.9550e+00 nan
 -8.8650e+00 nan
 -8.7750e+00 nan
 -8.6850e+00 nan
 -8.5950e+00 nan
 -8.5050e+00 nan
 -8.4150e+00 nan
 -8.3250e+00 nan
 -8.2350e+00 nan
 -8.1450e+00 nan
 -8.0550e+00 nan
 -7.9650e+00 nan

I generated a test indirect detection run where there's a single final state through a Higgs portal: phi0 phi0 > e+ e-. However, these errors were showing up in the full runs I was trying to do.
The commands I entered were:

import model <model>
define darkmatter phi0
generate indirect_detection e+ e-
output epem_test
launch
set precise

Reshuffling and Madevent give the same error. Any number of events that I have tried also gives the same error.

I tried looking in the cards to see if there's any precision settings I could change, or the energies of the beams. But, I had no luck with that.

Is there a known fix to this? I would very much appreciate any help. I can also attach any files that would be helpful.

Thank you very much,
Jason Arakawa