Physics89L

Notes and materials for Week 7

Topics covered: Data selection, efficiency, leakage, and rare event searches

Data analysis topics

Data selection, efficiency, and leakage in the context of dark matter (rare event) searches

This week we’re going to study techniques for selecting data of interest from a larger dataset. Oftentimes, experiments searching for particular signals need a way to sift through less interesting (a.k.a. background) data in order to see their desired signal. This is particularly important for sensitive experiments looking for rare events, where the background is dominant over any new physics.

To sensibly reduce data to see rarer and more interesting signals, we use two concepts: signal acceptance and background leakage. The former tells us how the “interesting” data is affected by our data reduction, and the latter tells us how the “less interesting” data is affected by our data reduction. Both are necessary to properly do a search for rare signals.

Data selection

Scientific context and resources

In these two notebooks, we’ll be studying a simple example of how a dark matter direct detection search works.

Dark matter is a type of matter that comprises about 85 percent of the matter density of the universe, but whose properties are still relatively unknown to us. There are a few different ways we can try to understand these properties. One of these is called “direct detection,” where we wait for dark matter particles to bump into very sensitive detectors and leave a signal. The LUX-ZEPLIN (LZ) experiment is one such detector. It uses a tub filled with liquid xenon to look for scatters of dark matter on xenon atoms:

LZ Detector

When particles scatter off of xenon atoms, they produce two signals: an initial burst of light (an S1), and a secondary burst of light (S2), created from electrons freed in the initial interaction. This light is detected by arrays of photomultiplier tubes (PMTs) in the top and bottom of the detector:

PMT Array

The S1 and S2 signal structures look slightly different, which helps us identify them (signals shown below - S1 is green, S2 is blue).

Event Waveform

Different types of interactions can create different sizes of S2, however. Dark matter scatters will generally produce less S2 light than will a background gamma ray or beta ray with a similar deposited energy. Having tools to suppress background, and equally importantly, discriminate background from our desired dark matter signal, allows us to mount ever more-sensitive searches for dark matter.

The LZ experiment will be turning on very soon, so stay tuned for word on our search!

Python functions reference

We will be using one simple new function to interpolate and subsequently evaluate an empirical/numerical relation.

Function Name What it does
scipy.interpolate.interp1d Given a set of data xvals and yvals, where yvals = f(xvals), compute y* = f(x) for some x within the range of xvals. See the SciPy interpolation documentation for more information