Saturday, 9 January 2016

Weekend Plot: The king is dead (long live the king)

The new diphoton king has been discussed at length in the blogoshpere, but the late diboson king also deserves a word or two. Recall that last summer ATLAS announced a 3 sigma excess in the dijet invariant mass distribution where each jet resembles a fast moving W or Z boson decaying to a pair of quarks. This excess can be interpreted as a 2 TeV resonance decaying to a pair of W or Z bosons. For example, it could be a heavy cousin of the W boson, W' in short, decaying to a W and a Z boson. Merely a month ago this paper argued that the excess remains statistically significant after combining several different CMS and ATLAS diboson resonance run-1 analyses in hadronic and leptonic channels of W and Z decay. However, the hammer came down seconds before the diphoton excess announced: diboson resonance searches based on the LHC 13 TeV collisions data do not show anything interesting around 2 TeV. This is a serious problem for any new physics interpretation of the excess since, for this mass scale,  the statistical power of the run-2 and run-1 data is comparable.  The tension is summarized in this plot:
The green bars show the 1 and 2 sigma best fit cross section to the diboson excess. The one on the left takes into account  only the hadronic channel in ATLAS, where the excess is most significant; the one on the right is bases on  the combined run-1 data. The red lines are the limits from run-2 searches in ATLAS and CMS, scaled to 8 TeV cross sections assuming W' is produced in quark-antiquark collisions. Clearly, the best fit region for the 8 TeV data is excluded by the new 13 TeV data. I display results for the W' hypothesis, however conclusions are similar (or more pessimistic) for other hypotheses leading to WW and/or ZZ final states. All in all,  the ATLAS diboson excess is not formally buried yet, but at this point any a reversal of fortune would be a miracle.

26 comments:

Anonymous said...

Hi Jester,
Can you explain the difference between the four red lines?

Jester said...

The limits come from ATLAS searches in three different channels: llqq (weakest, ATLAS-CONF-2015-071), lvqq (2nd weakest, ATLAS-CONF-2015-075), vvqq (2nd strongest, ATLAS-CONF-2015-068), and from the CMS search who combined different channels (strongest, CMS-PAS-EXO-15-002).

Anonymous said...

If you compare the exclusion limits presented in arXiv-1512.03371 with the ones presented by ATLAS and CMS in Run-2, you will see that the Run-2 first results can not exclude the combined Run-1 excess (but do disfavour the ATLAS hadronic result from Run-1).

Jester said...

That's exactly what I did in the plot (best fit from 1512.03371 rather than exclusion limits), and the conclusion is different.

Anonymous said...

Assuming a W', run 1 finds sigma8 = 5.3 ± 2.0 fb; run 2 finds sigma13 = 1 ± 13 fb; theory predicts sigma13/sigma8 ≈ 7. So run 2 reduces the statistical significance of the anomaly from ≈2.6 sigma to ≈1.9 sigma

Jester said...

My estimate is that the significance is reduced to ~1 sigma (not sure how you get sigma13 = 1 ± 13 fb, im getting something like sigma13 = -7 ± 8 fb since most of the 13 TeV searches mentioned above have deficit). In any case, if this was a real thing one would expect the significance to grow after doubling the dataset, instead of being reduced.

Anonymous said...

what d you think jester of the proposal in http://syymmetries.blogspot.com.au, to organize sort of a website/blog/... where all the different explanations of a given anomaly, lets say diphoton or B-decays, can be discussed in a organized way. ¿? maybe something along the lines of 1512.04357 but on the web.

Jester said...

I'm a bit skeptical this may work on internet. 1512.04357 started with a critical mass of people meeting in one workshop, and it was clear how the credit is to be assigned (individual papers first, then a common arxiv document, LH proceedings in the future). I'm not sure if this kind of structure can be reproduced without a central hierarchy, but I will follow this experiment with interest.

Unknown said...

Sorry to bother you with a question from the clueless commons, but the disappearance of the signal at higher energy provides a natural segue; how did particle physicists figure out that a signal going away at higher collision energies meant that it wasn't real at the lower energies?

Anonymous said...

Because at higher energies you are supposed to see everything you see at lower energies, even more "brightly", and then also the new stuff that appears only at higher energies.

If you fail to see something that you thought you saw at lower energies, it's because it wasn't there to begin with, and low-sample statistics has been playing tricks on you.

Unknown said...

"Because at higher energies you are supposed to see everything you see at lower energies, even more "brightly", and then also the new stuff that appears only at higher energies."

That's what I want to understand; how was that determined?

Anonymous said...

Jester - the 1 +/- 13 fb comes from Table 2 in 1512.07225. I'm not sure where that number comes from (nor your much less significant number). While one does expect the signal to increase with more data, it is always possible that the Run I data fluctuated up and the Run II data fluctuated down. We'll know for sure by summer.

Alex K said...

@Unknown

It is not a priori true that higher collision energy in a hadron collider will give you a better "view" of your signal of choice. One always has to take into account the so-called background processes, which give you final states which are the same as or will be indistinguishable from the signal you are looking for. If these background processes become much stronger to a greater extent than your signal at larger collision energy, you are worse off at a high energy collider. However, the production of very heavy multi-TeV particles that could explain such an excess (in this plot a new 1.9 TeV boson) is a typical example where going from 4 TeV proton energy to 6.5 TeV proton energy gives you a sizable increase in signal rate, and so these processes tend to profit from a higher energy LHC. In this case, to such an extent that, as Jester says, the 13 TeV data have comparable power to the 7+8 TeV data despite there being much less integrated luminosity available.

Jester said...

Anon: I see, thanks. 1512.07225 refers to only one ATLAS search at 13 TeV (the one that has a small excess), and ignores the other two that have deficits. I'm getting a similar result if I omit these two, while naively combining the CMS and 3 ATLAS searches I'm getting sigma13 = -3 ± 10 fb (the number I quoted in an earlier comment had a bug). In any case, we agree the significance drops a little below 2 sigma after the 13 TeV data are taken in account

mfb said...

A local significance of 2 sigma? If I had a dollar for every 2 sigma excess...

Anonymous said...

Jester, you said " I'm getting a similar result if I omit these two, while naively combining the CMS and 3 ATLAS searches I'm getting sigma13 = -3 ± 10 fb (the number I quoted in an earlier comment had a bug). In any case, we agree the significance drops a little below 2 sigma after the 13 TeV data are taken in account". But is sigma13 is 7-8 times as big as sigma8, then sigma8 would have to be (with sigma13 < 7 fb) less than 1 fb. That seems like no significant signal at all. Or did you scale sigma13 to be the same as sigma8?

Jester said...

Sigma13 quoted above is not scaled (the corresponding cross section scaled down to LHC 8 TeV is 6 times smaller). There is a clear tension between the run-2 and run-1 diboson data at 2 TeV. Both datasets have a similar statistical power after scaling the cross sections to a common reference. So when you combine the 3 sigma excess in run-1 with the half-a-sigma deficit in run-2 there remains a bit of an excess: 1.7 sigma according to my calculation, and 1.9 sigma according to 1512.07225. My interpretation is that this residual excess is not interesting at all.

David Brown said...

If the W± bosons both have heavy W' cousins, then how likely is it that the Z boson has a heavy Z' cousin?

Anonymous said...

this article seems to kill the other king 1601.03153 ;)

Anonymous said...

...and the other king is dead...1601.03153

Jester said...

failed gunpowder plot ;)

Anonymous said...

I find it quite amusing that it takes 100+ theory papers, with varying level of quality, until somebody comes along writing the first really relevant thing about the 750 GeV bump (I'm also talking about 1601.03153). This study seems to be very reasonable, at and least I am convinced that the "true" significance is much lower than what ATLAS and CMS claimed...

Ervin Goldfain said...

Anon @ 17:13,

I concur with you, yet it is still surprising that both collaborations see the bump roughly in the same place. So, the anomaly continues to beg for a good explanation, whatever that might turn out to be...

Anonymous said...

not criticizing the general message in 1601.03153 of background modeling being important - however I don't understand the use of the BIC in this context. BIC is an approximation to the Bayes-Factor - but I can hardly imagine anyone being able to properly define a BF in this context (BF + having no or very weak prior knowledge is a can of worms) -> BIC is meaningless here.

Also their background only fit seems to be over-predicting for >800 GeV. In the end this is a non-trivial model selection problem. It all boils down to how much ATLAS is able to simulate backgrounds. Is it really true that ATLAS decided for a1=0 after seeing the data? I'd imagine this was justified from their own MC studies?

Anonymous said...

@anon the BIC is an approximation in special cases to -2 log Bayesian evidence, not a Bayes-factor. A Bayesian evidence (and indeed a Bayes-factor) is independent of priors for models (but is a functional of prior densities for parameters in models).

I agree that the use of BIC is strange, however, as it's a combination of Bayesian model selection for the background shape followed by strictly frequentist hypothesis testing of the background only model (with a LLR test-statistic).

Obviously getting the background predictions/uncertainties is important, and I'd be astounded if ATLAS had made such a trivial error as using an unduly restrictive functional form for the background shape. If they have, it's a major embaressment ranking alongside FTL neutrinos, especially given the 150 hep-ph papers.

Anonymous said...

@anon

by prior knowledge I was referring to priors on model parameters, which the BF as ratio of the marginal densities is highly sensitive to (not priors on models themselves). I agree with the rest.