Xenon1T-2018

This code is a first attempt at reproducing (at least approximately) the Xenon1T-2018 upper limits on the DM-proton scattering cross section, announced by the Xenon1T collaboration on 28th May. The Xenon1T-2018 paper can be found here: arXiv:1805.12562.

The python code is in the form of a jupyter notebook: Xenon1T-limit.ipynb. You can run the notebook in your browser by clicking this badge: . It's also available as an executable python script: Xenon1T-limit.py.

Feel free to use this code however the hell you want. For questions, comments or bug reports, please contact Bradley J Kavanagh (bradkav@gmail.com).

Notes & Caveats

• The code roughly reproduces the Xenon1T-2018 limit for standard spin-independent interactions (within a factor of about 2). However, it is straightforward to use it to recast results for other types of interactions (and I might add/calculate those limits myself at some point).
• At the moment, the code does a very rough Poisson upper limit, including no background uncertainties and using no spectral information from the observed events themselves. This should be improved as more information becomes available.
• With the release of the paper on arXiv, things have changed a little bit. Because we don't use event-by-event information, our best bet at getting a reasonable limit is to minimise the number of events (because not all the events look like WIMP recoils). So we use the 0.65t exposure in the signal reference region when we reproduce the limit.

Requirements

• python3
• Standard python libraries like numpy and scipy
• WIMpy, for calculating signal spectra (or your favourite DM signal calculator)

Version history

20th June 2018: Fixed a typo in the Poisson upper limit. Now uses the 0.9t info for the limit. Energy info to be added soon.
1st June 2018: Updated to include the correct efficiency and updated numbers from the paper - arXiv:1805.12562. Now uses the 0.65t information to set the limit.
28th May 2018 (2): Updated to include Brazilian bands (updated folder structure a bit)
28th May 2018 (1): Initial version, using only a simple Poisson upper limit, from crappy screenshots of the Xenon1T webcast.

License

This project is licensed under the MIT License - see the LICENSE.md file for details