WimpSim

Welcome to the WimpSim WIMP Monte Carlo page for neutrino telescopes. This code calculates the annihilation of WIMPs inside the Sun or the Earth, collects all the neutrinos that emerge and let these propagate out of the Sun/Earth to the detector including neutrino interactions and neutrino oscillations in a fully consistent three-flavour way. This is done in an event-based scenario which is particulary suitable for interfacing to a neutrino telescope Monte Carlo.

On these pages, both the code to actually perform the simulations and results from an extensive set of simulations for various annihilation channels and masses are available.

If you use the results presented here, please refer to

M. Blennow, J. Edsjö and T. Ohlsson, [arXiv: 0709.3898]

If you also use the WimpSim code itself, please also refer to it as

J. Edsjö, WimpSim Neutrino Monte Carlo, http://www.physto.se/~edsjo/wimpsim/

Of course, you should also refer to the codes WimpSim uses, in particular nusigma (by J. Edsjö) and Pythia. See Torbjörn Sjöstrand's site for how they want you to refer to Pythia.

WimpSim Code

The code for WimpSim is available here.

Results

Results (figures and data files) for an extensive set of simulations run in September, 2007 are available here. These are the results referred to in the above publication. Note that these results are not re-run with the latest version of WimpSim.

Conversion script

The WimpSim results for the standard setup with default oscillation parameters have been incorporated into DarkSUSY, but are also available as a conversion script that can convert between various yields with different thresholds and angular cuts. This is particulary ueseful when comparing results from different neutrino telescopes with different thresholds. I am currently (Dec. 17, 2009) rerunning with the latest version of WimpSim. As soon as that is done, the conversion script and DarkSUSY will be accordingly updated.

Questions

If you have questions, send me an e-mail: