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. A short description of the plots follow:

• Top left: pT of jets with pT>15. In red is shown jets matched to a parton level W (by the usual distance \Delta R).
• Top right: Same color coding. pTsqrt(y).
• Bottom left: Same color coding. Mass of jet as given from jet algorithm.
• Bottom right: In red (signal) events with two decays of the type q_L -> q + chi+-1 -> qW chi01. Invariant mass of w+jet after cuts demanding two jets with pT>150, 1.2<y<2, 70<m_j<110. The additional jet is picked from the two hardes jets of the event.

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