The word on CRESST

So I've heard back, and indeed, that was a partial exposure. I still don't know how much, what cuts were used, when things were run or anything, really. So, for the moment, I'm trusting my guesses, because that's all I have.

The moment of truth

So like I said, I was very excited for the CRESST results. I learned this morning that their talk at WONDER 10 had been posted online. I quickly downloaded the results and what did I see? Well, first the exposure and operating detectors.

Wow! OK, so 9 detectors and if we assume 90% efficiency (what they typically do), then ~ 300 kg day. That's a lot. That's about 8 times what the commissioning run had.

OK, so what did they see?

Whoa! 4 events in the low energy region (lower blue shaded band). And I was expecting (if you recall) about 15.

4 != Poisson fluctuation from 15. (In case you were about to check.) So was I totally wrong?

Well, wait. It says "several detectors added".

Several != all

I looked up several and google told me it was at least 3. Well, I already knew that 2 had zero (at least part way through). So is this 4 detectors? 8? 5?

Fortunately for me, they claim to rule out inelastic dark matter (and are not the first to do so with CRESST data) by about a factor of 4. Now, I can handle a factor of 4. Data differences, fitting sine functions and of course halo models can buy you those kinds of factors, easy. My general rule is: ruled out by O(1), no biggee. Ruled out by O(10), I am uncomfortable. Ruled out by O(100), I am very uncomfortable. So this is discomfort, but not extreme.

But I can calculate what kind of signal they expect for their parameters. So, delta = 120 keV, m=200 GeV, sigma = 4 x 10^-41 cm^2. What do I get? 17.7 (taking v0=220 km/s, vesc=600 km/s). What is the 90% lower limit on events then (i.e., for how many events would I say 17.7 is the 90% upper limit)? Answer: 12 events. 


So I'm thinking that we still have some data to see. Maybe I'm off, maybe that 12 is really 10 or 11, but it's not 4 (which is what I thought when I saw this). 12 is a lot. It's about what you'd expect scaling up the commissioning run. So, if that's the case, I think I'm going to have to wait and see what the data look like...

NW

[Note that if they are using Yellin-type analyses, like maximum gap, then they could set stronger limits with the same number of events, if the distribution if somewhat off from what is expected. So they may have yet more events than this.]

What can you learn from what they don't tell you?

So I'm a "phenomenologist." Loosely speaking, this means that I care (practically) about the results of experiments, past, present and future. Of late a great deal of my emphasis has been on dark matter experiments.

One particular experiment I'm interested in is the CRESST experiment. CRESST uses a CaWO4 crystal and measures both the phonon signal as well as a light yield, and the ratio of the two helps determine nuclear recoils from electronic recoils (and thus, hopefully, dark matter scatters from background). It's really a fantastic experiment, capable of better than 1 keV energy resolution.

The only problem is a persistent set of background events at the detector. They may be cracks, and there are plausible suggestions by the collaboration for the source of them, but they remain. It's not impossible that they could be real dark matter events if dark matter is inelastic, for which the tungsten in the target is particularly sensitive. It's speculative, but who knows?

Anyhow, I've been eager to see the updated results from the experiment for a while now. They didn't publish their last run (which was an R&D run) and I've been looking forward to the results from this current run. Unfortunately, those are not yet out.

However, they did show a "selected" set of results at the UCLA DM conference. It sounds like they have 8 or 9 working detectors and showed 2, both which were clean, i.e., had no events in them in the range up to 40 keV.

So what does one make of this?

Well, it occurs to me we can play a bit of a game. We can make a weird sort of prediction, based on psychology, of what they might see. Since clean detectors are necessarily background free (since there are no events) it's good to show them. So, if I were running an experiment, I'd show all the clean detectors I had. We can turn this around, and use the fact that we saw precisely two detectors with no events to extrapolate that these were the only event free detectors.

So, this is a big speculation, but there's no harm in idle speculation. Let's suppose they had 8 or 9 experiments, two of which had no events. Well, if the events come from a Poisson distribution with expectation N, then the probability of seeing 0 events in any given detector is exp[-N].

So, let's say there are 9 operating detectors. Then the likelihood that you have two detectors with nothing in them is 36 exp[-2N](1-exp[-N])^7. What does this look like as a function of N?

This function has a peak at 1.5 and an expectation of 1.8. This is consistent with their commissioning run.

So if they have 9 operating detectors, that's my expectation.

In essence, I'm just taking advantage of the fact that in the limit of many equivalent detectors, the fraction expected to have 0 events is just Exp[-N]. So, if that's 2/9, then N=1.5. Nine may not be many, and they may not be equivalent, but probably not too terrible of assumptions.

But we can still go a bit further. Given that we are assuming that we've seen the clean detectors, the others must have at least one event. If I take the peak expectation value of 1.5, the Poisson expectation for detectors that have at least one event is 1.93. Thus, for the seven remaining detectors I expect 13.5 events in the energy range up to 40 keV.

What does this mean? OK, well, nothing. There are lots of reasons they might have only shown two: maybe they'd only finished analyzing two, maybe they just thought that was the right number. Still, it's fun what one can do when one combines statistics+human psychology+unjustified assumptions. Just like phenomenology.

NW

[Using the expectation value of 1.8 rather than the peak of 1.5, you get an expectation of 15.1 events. If only 8 detectors are working, the peak (expectation value) is at 1.4 (1.7), and the expectation for the remaining 6 is 11 (12.5) events. And then, of course, these are just mean values. Also I should add that I'm very tired at the moment so I'm very likely making stupid mistakes with my Poisson probabilities.]

the black hole coffee bar at the perimeter institute

Second Post

OK, so it took some time for the second post. More soon.

First Post

Ah blogging, the narcissistic indulgence of the new millenium....