Dear Mara,
it seems the issue is that you use the sample_discrete_momentum_energy() in the source routine which is for discrete energy lines. For a contious spectrum you need to use the sample_spectrum_momentum_energy() function.
Cheers,
David
1 Like