PYTHIA 8 error for mDM < 100 GeV

Hi,

I am using MadDM for studying DM indirect detection. In my model, there are 3 kinds of CP-even neutral Higgs bosons, and DM can annihilate into a pair of them, as well as other final states.

I was trying to generate gamma-ray spectrum with the 'precise' mode. I found that the gamma-ray spectrum can be successfully generated for mDM > 100 GeV. However, for mDM < 100 GeV, PYTHIA 8 failed to generate the spectrum. The error information in pythia8.log is as follows.

 -------- LHA initialization information ------------

  beam kind energy pdfgrp pdfset
     A 52 97.198 0 247000
     B 52 97.198 0 247000

  Event weighting strategy = -4

  Processes, with strategy-dependent cross section info
  number xsec (pb) xerr (pb) xmax (pb)
    1995 5.4686e+02 8.3103e-01 5.4686e+02

 -------- End LHA initialization information --------
mDM: 9.7198e+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
 PYTHIA Warning in TauDecays::decay: unknown correlated tau production, assuming from unpolarized photon
Event generation aborted prematurely, owing to error!

 *------- PYTHIA Event and Cross Section Statistics -------------------------------------------------------------*
 | |
 | Subprocess Code | Number of events | sigma +- delta |
 | | Tried Selected Accepted | (estimated) (mb) |
 | | | |
 |-----------------------------------------------------------------------------------------------------------------|
 | | | |
 | Les Houches User Process(es) 9999 | 10 10 0 | 0.000e+00 0.000e+00 |
 | ... whereof user classification code 1995 | 10 10 0 | |
 | | | |
 | sum | 10 10 0 | 0.000e+00 0.000e+00 |
 | |
 *------- End PYTHIA Event and Cross Section Statistics ----------------------------------------------------------*

 *------- 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: attempting to set properties |
 | 1 Info from SLHAinterface::initSLHA: No MODSEL found, keeping internal SUSY switched off |
 | 1 Info from SLHAinterface::initSLHA: importing DECAY tables |
 | 1 Info from SLHAinterface::initSLHA: importing MASS entries |
 | 2 Warning in SLHAinterface::initSLHA: ignoring DECAY tables |
 | 1 Warning in SLHAinterface::initSLHA: ignoring MASS entries |
 | 4 Warning in SLHAinterface::initSLHA: ignoring QNUMBERS |
 | 1 Warning in SLHAinterface::initSLHA: massless particle forced stable |
 | 1 Warning in SLHAinterface::initSLHA: using QNUMBERS for id codes < 1000000 may clash with SM. |
 | 10 Warning in TauDecays::decay: unknown correlated tau production, assuming from unpolarized photon |
 | |
 *------- End PYTHIA Error and Warning Messages Statistics ------------------------------------------------------*

Any idea for fixing this?

Thanks for the help!

Best regards,
Zhao-Huan