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

### Task

Implement the search analysis described in **CMS-SUS-16-047** and use it to predict the binned differential signal yield in `p_T^miss` for the benchmark point `T6gg_1750_1650`, in the **high-H_T^gamma signal region**, normalized to 35.9 fb⁻¹ at √s = 13 TeV.

You should:

1. Generate `T6gg_1750_1650` 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 high-H_T^gamma signal-region cuts that define this analysis. Apply them to your generated events.
3. Histogram the surviving events in `p_T^miss` using the bin edges already present in the `results/*.yaml` template (do not modify them).

### Definitions

- `p_T^miss` — the magnitude of the negative vector sum of the transverse momenta of all candidate objects in the event.
- `T6gg_1750_1650` — pair-produced squarks at `m(q̃) = 1750 GeV`, each decaying via the `T6gg` simplified-model topology to a neutralino at `1650 GeV` plus a quark jet, with two photons in the final state.
- High-H_T^gamma signal region — the higher-H_T^gamma branch of the analysis's two H_T^gamma-categorised signal regions, where `H_T^gamma` is the scalar sum of jet transverse momenta and the leading-photon `pT`.

### 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 p_T^miss > 600 GeV.
               Do not modify bin edges.
              "
description: "Event counts for binned p_T^miss in the high-H_T^gamma signal region"
target: "T6gg, squark 1750 GeV, neutralino 1650 GeV"
cm_energy_gev: 13000
luminosity_fb: 35.9
---
dependent_variables:
- header: {name: T6gg_1750_1650}
  values:
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
  - {value: null}
independent_variables:
- header: {name: PTMISS, units: GEV}
  values:
  - {low: 0.0, high: 10.0}
  - {low: 10.0, high: 20.0}
  - {low: 20.0, high: 30.0}
  - {low: 30.0, high: 40.0}
  - {low: 40.0, high: 50.0}
  - {low: 50.0, high: 60.0}
  - {low: 60.0, high: 70.0}
  - {low: 70.0, high: 80.0}
  - {low: 80.0, high: 90.0}
  - {low: 90.0, high: 100.0}
  - {low: 100.0, high: 110.0}
  - {low: 110.0, high: 120.0}
  - {low: 120.0, high: 130.0}
  - {low: 130.0, high: 140.0}
  - {low: 140.0, high: 150.0}
  - {low: 150.0, high: 160.0}
  - {low: 160.0, high: 170.0}
  - {low: 170.0, high: 180.0}
  - {low: 180.0, high: 190.0}
  - {low: 190.0, high: 200.0}
  - {low: 200.0, high: 250.0}
  - {low: 250.0, high: 300.0}
  - {low: 300.0, high: 350.0}
  - {low: 350.0, high: 450.0}
  - {low: 450.0, high: 600.0}
  - {low: 600.0, high: 800.0}