During a scan over parameter the flux at source as a function of energy are not saved

[Mattia Di Mauro]

Hi everyone,

I am using the model ScalarHiggsPortal_NLO_UFO.
I would like to calculate the source spectra, as a function of energy, for different dark matter masses.
In order to do so I do:
import model ScalarHiggsPortal_NLO_UFO
define darkmatter n1
generate relic_density
add indirect_detection
output folder
launch folder
set precise
set main:numberofevents = 1000
set msdm = scan:[ 5. , 6.00227468, 7.20546027, 8.64983036, 10.38373155, 12.46520181, 14.96391304, 17.96350329, 21.5643762 , 25.88706187]

This runs Pythia for all the selected masses and makes a file scan_run_01.txt that contains the flux of neutrinos, gamma rays, positrons and antiproton.
However, it does not generate the files:
antiprotons_spectrum_pythia8.dat, gammas_spectrum_pythia8.dat, neutrinos_spectrum_pythia8.dat and positrons_spectrum_pythia8.dat that are usually printed when running MadDM for a single mass.

It would be nice to save these files also when running scan over the mass.

Best
Mattia

[Chiara Arina (carina)]

Dear Mattia,

indeed by default we do not save the source spectra, as in case of massive scans this can fill in a lot of space. If you want to store in the folder the source spectra you should modify the maddm_card.dat. The line:

off = save_output

should be replaced with:

spectra = save_output

Hope this helps

best

Chiara