Banner
    The Plot Of The Week - ATLAS Dilepton Resonance Search
    By Tommaso Dorigo | June 19th 2011 01:49 AM | 11 comments | Print | E-mail | Track Comments
    About Tommaso

    I am an experimental particle physicist working with the CMS experiment at CERN and the CDF experiment at Fermilab. In my spare time I play chess...

    View Tommaso's Profile
    The ATLAS Collaboration, one of the two high-energy physics experiments at the CERN Large Hadron Collider, has just produced updated results of their ongoing search for new heavy particles decaying into lepton pairs. They are now using up to 236 inverse picobarns of 7 TeV collisions, which is seven times more data than previous searches based on 2010 datasets. A seven-fold increase in data size grants a significant increase in sensitivity, so it is worth taking a look at what they see.

    The technical way to call the sought particles is "Z' resonances", since they can be thought of as heavier brothers of the Z boson, although they may well have different properties; they are usually the result of adding one U(1) group to the group structure of the Standard Model (which is a SU(2)xU(1) combination, plus a SU(3) for strong interactions). U(1) is a one-dimensional unitary group of transformations, and it introduces just a new particle in the theory, a new boson which mediates a new interaction. A Z' might not explain a lot about the way things really work at a fundamental level, but many like such "minimal" extensions of the Standard Model, mainly because they are easy to cook up and not in too striking disagreement with current data.

    The search for Z' resonances is pursued with momentum in CMS and ATLAS. I know what is going on in CMS being a member of that collaboration, but I only get to know where ATLAS stands when I get to read their public papers. That's why I am interesting in checking their high-mass spectrum, to see if they see anything at odds with predictions. In two words, they don't. Check it out in the dimuon invariant mass distribution below:


    The data are black points (a very stable standard by now, unfortunately about the only world-wide agreement among HEP experiments -a painful concept for those like me, who are working to reach at least lab-wide agreements on limit-setting procedures and other highly less trivial things), and backgrounds are shown with different vivid colours. Of course by far the most important process contributing to the data (which feature two muons and little else in the event) is Drell-Yan production -the process whereby a Z boson or a virtual photon is produced by quark-antiquark annihilation. This is shown in light blue. Other processes include the production of two bosons together (ZZ, or ZW), in yellow; and top pair production (in red). The data, needless to say, closely follows the sum of backgrounds, while a Z' resonance would produce a curve such as the empty histograms shown on the right (which correspond to the signal one would see if a Z' existed with a mass of one to 1.5 TeV).

    A similar plot is produced for events containing electron-positron pairs, see below:

    By comparing the two figures one notices a few things. First, the dielectron spectrum includes a non-negligible amount of QCD background (light yellow), which is basically absent in dimuon events. This is due to the non negligible fraction of electron candidates which are actually hadronic jets: you cannot get rid of that background to electron candidates, really. It arises when a jet fragments into a leading charged pion, which leaves a nice stiff track in the central tracker, and then upon entering the electromagnetic calorimeter undergoes a nasty process called "charge exchange", whereby it becomes a neutral pion. The neutral pion immediately decays into two photons, which proceed to yield an electromagnetic cascade which mimics awfully well the signal of a single electron.

    A second thing to note is that the resolution of a dielectron resonance of high mass in ATLAS is surprisingly good compared with the one for dimuons: the Z' signals are much narrower in the latter figure. That is one distinguishing feature of the ATLAS detector (for CMS things are different, due to the characteristics of their muon system and tracker).

    In the end, ATLAS can exclude new Z' bosons of a variety of models (all similar, but slightly different in the resulting new particle cross section predictions). Limits are set in the 1 TeV - 1.2 TeV range. For more details, see the public page of ATLAS results here.

    Comments

    They don't confirm CMS excess, then?

    which excess?!?!

    D

    Everybody's heard about this one already. CMS has an excess on dileptons at roughly 1 TeV.

    How much does an error bar represent? Is that +/- 1 sigma or +/- 2 sigma?

    dorigo
    There are no error bands in these plots Anon. The data have statistical uncertainties as error bars (basically sqrt(N)), and the Monte Carlo simulations overlaid have no visible uncertainty drawn -but you should consider that 10% to 20% uncertainties are possible for some of the involved processes. Not the Drell-Yan, though, which is calculated very precisely.

    Cheers,
    T.
    Ok, so the error bars represent "statistical uncertainties" over the number of observed events? What's the exact physical meaning of this?

    Take the last graph for example: in the 400-600 GeV region there are several consecutive bins with a data point corresponding to 1 event: their bars seem to go from 0.1 to 4. What do "0.1 events" and "4 events" mean here, what makes that specific interval relevant and how do you define it if it's not a X% confidence interval? A single event was counted, what's uncertain about it?

    Maybe this has something to do with the chosen bin width? (in the sense that the smaller the bin, the higher the probability that an event is wrongly assigned to that bin due to uncertainties in the event reconstruction)

    dorigo
    Dear Anon,

    0.1 (but actually 0.2) is probably a silly truncation. In general, you can define a central 68% interval for a Poisson distribution, and do things more tidily; but in most cases the data counts are just assigned a uncertainty equal to sqrt(N). I am not sure what ATLAS did, though, and you correctly noticed that they seem to be doing something else. In any case, I doubt that those bars have much information to convey: as you note, the counts are what they are, and they should not carry a uncertainty; but it has become a standard to plot them.

    I strongly doubt that the horizontal location and resolution have anything to do with the error bars here, that would really be the craft of a twisted mind, plus it would be completely incorrect.

    Best,
    T.
    Why does it look like the background is overestimated between 400 and 500 GeV? I can see it's within errorbars, yes, but it still looks odd.

    Between 900 GeV and 1 TeV there is one event when only 0.1 was expected in the dielectron distribution. At the same time, in the 300-600 GeV range more events were expected. For dimuons, there is also an excess between 500-700 GeV. Isn't all of this odd? Has perhaps CMS produced similar plots?
    thanks, BB

    dorigo
    No, it does not strike me as odd in the least. These are distributions in great agreement with expected backgrounds. As for CMS, they are surely going to produce similar results soon.

    Cheers,
    T.
    Hi Tommaso,

    One thing that always confuses me a bit in these figures is the region where the number of events per bin expected is of order one or smaller. The problem is that, due to the logarithmic scale, the bins where 0 events are observed have no data points, since 0 does not fit on a log graph.

    Lets look e.g. at the muon plot in the 500-700 GeV range. There one seems to see an excess in four consecutive bins. Only after looking closer one notes that between these four bins there is one that has zero observed events, and also the bins before and after have a downward fluctuation. Taking those "invisible" data points into account makes it much more obvious that the excess events are just upward fluctuations.

    Ok, now, assuming that those bins that have no data points are really zero measurements. Shouldn't the error be +-1? Then at least the error bars of these bins should be shown.

    Another funny feature is that the Diboson background seems to peak at 800 GeV. This is probably monte carlo statistics, right?

    Cheers,
    P.

    Add a comment

    The content of this field is kept private and will not be shown publicly.
    • Allowed HTML tags: <sup> <sub> <a> <em> <strong> <center> <cite><TH><ul> <ol> <li> <dl> <dt> <dd> <img> <br> <p> <blockquote> <strike> <object> <param> <embed> <del> <pre> <b> <i> <table> <tbody> <div> <tr> <td> <h1> <h2> <h3> <h4> <h5> <h6> <hr> <iframe><u><font>
    • Web page addresses and e-mail addresses turn into links automatically.
    CAPTCHA
    If you register, you will never be bothered to prove you are human again. And you get a real editor toolbar to use instead of this HTML thing that wards off spam bots.