A fast and accessible program for computing WIMP-nucleus scattering event rates

Astronomical observations indirectly suggest a large fraction of the matter in the universe is nonbaryonic and weakly interacting. This has prompted numerous and ongoing experimental efforts to directly detect the weakly interacting massive particles (WIMPS). Recent work using an effective field framework has shown a larger number of couplings between nucleons and WIMPs than previously appreciated. While a Mathematica script to compute the nuclear target response exists, exploration of these possibilities is hampered by the size of the parameter space and the relatively high cost of computation in Mathematica. To address the large number of targets and couplings, we have developed a fast Fortran 90 code with a user-friendly Python wrapper, that can swiftly and efficiently explore the parameter space at rates between 100 and 1000 times faster than before. I cover the basics of implementation and demonstrate a practical use-case.