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:
| 1 Info from SLHAinterface:
| 1 Warning in SLHAinterface:
| 1 Warning in SLHAinterface:
| 1 Warning in SLHAinterface:
| 1 Warning in SLHAinterface:
| |
*------- 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
Question information
- Language:
- English Edit question
- Status:
- Open
- For:
- MadDM Edit question
- Assignee:
- Chiara Arina Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Jason Hiroshi Arakawa for more information if necessary.