> **Paper:** CMS-SUS-16-051
> **Centre-of-mass energy:** 13 TeV
> **Luminosity:** 35.9 fb⁻¹
> **Task type:** simulation
> **Signal benchmark:** `T2tt_800_400`
> **Observable:** `E_T^miss`

### Task

Implement the search analysis described in **CMS-SUS-16-051** and use it to predict the binned differential signal yield in `E_T^miss` for the benchmark point `T2tt_800_400`, in the analysis's **compressed signal-region selection**, normalized to 35.9 fb⁻¹ at √s = 13 TeV.

You should:

1. Generate `T2tt_800_400` events using a matrix-element generator + parton shower + detector simulation chain of your choice.
2. Read the paper to determine the object identification, event-selection requirements, and the compressed-selection signal-region cuts that define this analysis (Section 5 / Table 3 of the paper). Apply them to your generated events.
3. Histogram the surviving events in `E_T^miss` using the bin edges already present in the `results/*.yaml` template (do not modify them).

### Definitions

- `E_T^miss` — the magnitude of the negative vector sum of the transverse momenta of all reconstructed particle-flow candidates in the event.
- `T2tt_800_400` — pair-produced top squarks at `m(t̃) = 800 GeV`, each decaying via the `T2tt` simplified-model topology to a top quark and a neutralino at `m(χ̃⁰₁) = 400 GeV`.
- Compressed signal-region selection — the search's secondary categorisation (Table 3 in the paper), designed to enhance sensitivity to compressed-spectrum stop decays via an ISR-tagged event topology.

### Output requirements

| Artifact | Purpose |
|---|---|
| `results/*.yaml` | Fill the `null` bin values with your predicted signal yields. |
| `analysis/*.py` | Event-selection code, runnable on your generated sample(s). |
| `data/*.root`, `sims/*.dat` | Selected-event files + generator / detector cards. |
| `report.md` | What you produced, the methodological choices you made, and where you deviated from the paper (if at all). |

### Important

> The goal is to *predict* the signal event distribution from your own simulation and analysis pipeline. Do not extract bin values from the paper's figures, tables, HEPData record, or elsewhere.


### Answer template (YAML)

instructions: "Fill the 'null' bin values in dependent_variables with your results.
               In the rightmost bin include ALL events with E_T^miss > 550 GeV.
               Do not modify bin edges.
              "
description: "Signal event counts for binned E_T^miss"
target: "T2tt, stop 800 GeV, neutralino 400 GeV"
cm_energy_gev: 13000
luminosity_fb: 35.9
---
dependent_variables:
- header: {name: T2tt_800_400}
  values:
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
independent_variables:
- header: {name: ETMISS, units: GEV}
  values:
  - {low: 250.0, high: 350.0}
  - {low: 350.0, high: 450.0}
  - {low: 450.0, high: 550.0}
  - {low: 550.0, high: 650.0}