Identifying massive hadronic decays in SUSY decay chains.

This page is intended to help discussion between Are, Erckan, John and Jon and isn't intended to be used by anyone else.

Project plan

We'll start by assuming that in a SUSY pair production, one decay chain has allowed us to remove all non-SUSY backgrounds, and what we are trying to do is interrogate/unpick the second decay chain to learn more about the SUSY spectrum. In particular we are looking at high ET jets and seeing they come from hadronic decays of W or Z.

John: I suggest that we use the benchmarks postulated in hep-ph/0508198, particularly scenarios \alpha, \beta, \gamma and \delta. The first three have neutralino dark matter but non-universal Higgs mases (NUHM), with parameters chosen to get varying branching ratios for \chi_2^0 \to Z,H \chi_1^0, as seen in Fig. 2. I do not know offhand about the \chi^\pm \to W \chi decays.

The point \delta has gravitino dark matter, but the neutralino is the next-lightest sparticle and its decays are invisible. Hence, all four points have 'standard' missing energy signatures.

Since there are issues in treating NUHM models using ISAJET, we might want to start with point \delta. It is in any case the most interesting, as it has relatively large masses and hence is quite challenging using conventional techniques (see Fig. 9).

After our discussion at CERN on 19 July, we decided the immediate way ahead was as follows:

• Are to produce the PYTHIA config files from these and upload to this page.
• Are to run the [http://hepforge.cedar.ac.uk/hzsteer HZSteer] analysis and show us the key distributions (mass, y, ET, eta etc... as shown previously).
• Erckan to run similar jobs with Athena ATLFast and compare to Are's.

After this, we should have a good idea how much mileage there is in this kind of study, and hopefully be in a position to start looking at more challenging SUSY scenarios and treating non-SUSY background events.

The SLHA file for John's \delta point has been generated with ISAJET 7.69, and SUSY decays added with SDecay 1.1. I checked the decays with the ones given by ISAJET and modulo small numerical differences they look the same. The file is attached. Also attached is the NLO cross sections for \delta as calculated with Prospino2.

Plots for \delta have been produced and attached (NEW version as of 10.08.06). As a first attempt we are looking for events were two W have been produced in chargino decays. The decay chain we are looking for is ~q_L -> q + ~chargino_1 -> q + W + ~neutralino_1.

A short description of the plots follow:

• Top left: Plot of pT of jets. Cut on pT>15 and |\eta| < 2.5. In red is shown jets matched to a parton level W (by the usual distance \Delta R < 0.1).
• Top right: Plot of jet mass m_W as given from jet algorithm. Same color coding. Cut on jet pT: 150 < pT < 600.
• Bottom left: Plot of jet log(pT*sqrt(y)). Same color coding. Cut on jet pT: 150 < pT < 600. Cut on jet mass: 80 < m_W < 90.
• Bottom right: Plot of invariant mass of w+jet system. We require two W's identified in each event from the cuts given above pluss a cut on separation scale: 1.3 < log(pT*sqrt(y)) < 2.0. The W candidates are combined with the two hardest jets in the event. In red (signal) events with two decays of the type mentioned above.

Note that the three first plots have been normalized, and the cuts accumulated through the plots, which was not done in the first version. From the distributions of jet pT and mass, cuts look to be potentially very effective. However, in the log(pT*sqrt(y)) distribution the W-tagged jets have almost the same distribution as other jets, with only a slight excess just below log(pT*sqrt(y))=2. As a result the cut on separation scale is not very effective, and in the invariant mass distribution we see that while the signal plot (red) has roughly the expected endpoints, one can show analytically that for the \delta masses we should have 247 < m_{qW} < 1272, it is burried in SUSY background.

As you can see there is no longer a top peak as in the first version, so this is no longer a problem. There are still too many background jets, and rather surprisingly the separation scale cut does not work very well. I

We have discrimination between the W's and other jets as seen from the first three plots, but the invariant mass plot shows there are too many W candidates around after cuts. Note the a part of the problem is the copious production of top quarks in this model, which can be seen as peaks in the invariant mass plot. I'm still working on the numbers for how efficient the cuts are and how much of the problem is caused by the real W's produced in other decays, and how much is being faked by other jets.

-- JonButterworth - 19 Jul 2006

-- ErkcanOzcan

• delta_slha.dat: SLHA file for the \delta benchmark

• delta.xsec: Cross sections at NLO for \delta calculated with Prospino2

• delta.ps: Som start-up plots for \delta