SusySubjet < Main

---+ 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 [[http://www-spires.dur.ac.uk/cgi-bin/spiface/hep/www?rawcmd=FIND+EPRINT+HEP-PH+AND+EPRINT+0508198][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. 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 GeV < pT < 600 GeV. 
   * Bottom left: Plot of jet log(pT*sqrt(y)). Same color coding. Cut on jet pT: 150 GeV < pT < 600 GeV. Cut on jet mass: 80 GeV < m_W < 90 GeV.
   * 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 have been added 

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.
 


-- Main.JonButterworth - 19 Jul 2006

-- Main.ErkcanOzcan 
   * [[%ATTACHURL%/delta_slha.dat][delta_slha.dat]]: SLHA file for the \delta benchmark

   * [[%ATTACHURL%/delta.xsec][delta.xsec]]: Cross sections at NLO for \delta calculated with Prospino2

   * [[%ATTACHURL%/delta.ps][delta.ps]]: Som start-up plots for \delta
Attachments
Topic attachments
I	Attachment	Action	Size	Date	Who	Comment
ps	delta.ps	manage	28.4 K	2006-07-20 - 21:19	AreRaklev	Som start-up plots for \delta
xsec	delta.xsec	manage	10.2 K	2006-07-20 - 19:30	AreRaklev	Cross sections at NLO for \delta calculated with Prospino2
dat	delta_slha.dat	manage	45.8 K	2006-07-20 - 13:02	AreRaklev	SLHA file for the \delta benchmark
This topic: Main > TWikiUsers > JonButterworth > SusySubjet
Topic revision: r16 - 2006-08-10 - AreRaklev