wiki:Korea2017-05-CMS-EXO-16-052

Version 48 (modified by Ian James Watson, 7 years ago) ( diff )

--

The first MadAnalysis 5 workshop on LHC recasting @ Korea

Analysis 5: CMS-EXO16-052

Students Jeon Dajeong [mail:djeon@cern.ch] , Kim Seulgi [mail:pikachu0617@naver.com], Lee Sang Man [mail:daniel.s.lee4@gmail.com], Watson Ian James [mail:ian.james.watson@cern.ch]
Supervisor(s) Sam Bein, Jory Sonneveld [mail:jory.sonneveld@cern.ch, samuel.bein@gmail.com]

1. Analysis summary

2l+MET

2. Analysis global information

Repository for work: ​https://github.com/daniellee1004/DilepMeT

​https://github.com/cms-sw/genproductions/tree/master/python/ThirteenTeV

​https://gitlab.cern.ch/jsonneve/ma5recasts/blob/master/cms_sus_16_033/plots/plot_saf_file.py

Low-luminosity version of this analysis:

​https://inspirehep.net/record/1618045/

it is already included in the PAD!

3. Selection description

3.1 Trigger selection

PS: This part cannot usually be described by the fast-simulation of detector.

3.2 Object definition

Electrons

  • pT > 20 GeV, |eta| < 2.4
  • Medium ID
  • PFiso04 < 0.077 (barrel electrons) or PFiso04 < 0.068) (endcap electrons)

Muons

  • pT > 20 GeV, |eta| < 2.4
  • Tight ID (Loose ID for vetoing 3 leptons events)
  • PFiso04 < 0.15

To get isolation:

double Riso = PHYSICS->Isol->eflow->relIsolation(el,event.rec(),0.4,0.,IsolationEFlow::ALL_COMPONENTS);

Jets

  • pT > 20 GeV, |eta| < 5 (PFloose ID)
  • Identified leptons in {delta}R < 0.4 disregarded
  • We categorize the selected events based on the number of selected jets with pT > 30 GeV for the corresponding requirements

b-jets

  • pT > 20 GeV, |eta| < 2.4
  • CSVv2>0.8484 (medium)

Taus

MET

  • Type-1 corrections are applied ?

3.3 Cut-flow chart

Event Selection:

  • Two same flavored lepton with opposite electric charge ( that pass lepton selection?)
  • el[0] pt > 25 and el[1] pt >20
  • mu[0] pt > 20 and mu[1] pt > 20
  • |dilep m - 91.1876( +- 0.0021)| < 15
  • dilep pt > 60
  • events having more than one jet with pT > 30 GeV - remove
  • b-tagged jet Veto : at least one bjet pt > 20, |eta| < 2.4, CSVv2 < 0.8484 - remove
  • if el[2 onwards] pt > 10 -remove
  • if mu[2 onwards] pt > 5 -remove
  • tau veto: loose id, tau pt > 18. (hadronically decaying tau leptons)

Optimization:

  • For each possible set of selections the full analysis, including the estimation of backgrounds from control samples in data and the systematic uncertainties, is repeated.

After optimization:

  • MET > 100
  • Delta phi ( dilep.pt , ptmiss) > 2.6 rad
  • |MET - dilep.pt| / dilep.pt < 0.4

If event has one jet:

  • Delta phi (jet.pt , ptmiss ) > 0.5
  • Delta R(ll) < 1.8

3.4 Definition of exotic observables

Like SUSY transverse observable. Notice that classical MT2 & MT2W are encoded within the MadAnalysis framework.

3. Detector simulation

Can we use the last release of the Delphes package with the CMS/ATLAS MA5tune-card? or do need to improve the simulation?

4. Validation material

4.1 Monte-Carlo samples

  1. LO/NLO? LO
  2. Model used:

1) dark matter production in a simplified model with a spin-1 mediator A

2) dark matter production in a simplified model with a spin-0 mediator phi

  1. MC generator program (name + version): MAD-GRAPH5_aMC@NLO 2.2.2
  2. Shower program (name + version): PYTHIA 8.205
  3. Tune of Pythia: CUETP8M1

+ The coupling g(kai) is set to one.

For 0.25 < g(q) < 1.0,

the width ranges 1-5% of the mediator mass for g(q) = 0.25,

30-50% of the mediator mass for g(q) = 1.00

Events for ADD extra-dimension scenario

LO, EFT implementaion in PYTHIA 8

M(D) = 1, 2 and 3 TeV, each with n = 2, 3, 4, 5, 6, 7

The events for the unparticle model

LO, PYTHIA 8, tune CUETP8M1

assuming cutoff scale ㅅ(U) =15TeV

ZH & gluon-gluon fusion = WZ & ZZ processes

4.1.1 Running Madgraph

import DMsimp_s_spin0
generate p p > xd xd~ z [QCD]
add process p p > xd xd~ z j [QCD]
output DM_spin0_noreallythistime
launch

Turn on madspin and the showering. Then fix:

madspin_card.dat -> set "z > l+ l-"

run_card.dat -> can set number of events, but also change shower to PYTHIA8, ickkw to 3 (FxFx merging)

shower_card.dat -> Qcut to 15 (50?) and njmax to 1 (since were generating 0 + 1 extra partons)

param_card.dat -> set the masses of Y0 (or Y1 for spin 1), the mediator, and Xd, the dark matter particle (the paper uses MY0=500, MXD=150GeV)

we will also have to set the pythia settings to match those of CMS (UE tune, etc.).

4.2 Reference plots / cut-flows

4. Validation results

Sample Cross Section (pb) Eff. Expected
./Samples/DMsimp_s_spin1-axial_0j.txt 0.021318904968 0.1272 97.3523531583
./Samples/DMsimp_s_spin1-axial_ee_0j.txt 0.010659482475 0.0978 37.4256561594
./Samples/DMsimp_s_spin1-axial_mm_0j.txt 0.010659482475 0.1481 56.6742298283
./Samples/DMsimp_s_spin1-vector_0j.txt 0.0393553945622 0.1227 173.357758169
./Samples/DMsimp_s_spin1-vector_0j_qcd.txt 0.06788 0.1251 304.8551892
./Samples/DMsimp_s_spin1-vector_1j.txt 0.0666819164253 0.0752 180.037839919
./Samples/DMsimp_s_spin1-vector_1j_qcd.txt 0.1251 0.1125 505.082719261
./Samples/DMsimp_s_spin1-vector_ee_0j.txt 0.0196775973111 0.1003 70.8545020699
./Samples/DMsimp_s_spin1-vector_mm_0j.txt 0.0196775973111 0.1572 111.050126873

From the paper:

  • vector 500 GeV med, 150 GeV DM: 40.9+-1.6 ee, 57.9+-2.3 mm, 98.8+-3.9 total
  • axial 500 GeV med, 150 GeV DM: 27.5 ee, 38.0+-1.5 mm, 65.5+-2.6 total
Sample Cross Section (pb) Eff. Expected Exclusion (sigma/sigma_0)
Run 30/10/2017
./Samples/DMsimp_s_spin0P_0j_1000_1.txt 5.28328517293e-06 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0P_0j_100_1.txt 5.29581021579e-05 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0P_0j_10_1.txt 9.81497810099e-05 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0P_0j_500_1.txt 2.10019490453e-05 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0P_0j_50_1.txt 8.02167613779e-05 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0S_0j_1000_1.txt 6.46834888936e-06 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0S_0j_100_1.txt 9.65746285084e-05 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0S_0j_10_1.txt 0.00028727155715 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0S_0j_500_1.txt 2.66324370678e-05 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0S_0j_50_1.txt 0.000177391817776 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0_default.txt 6.50211256127e-06 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin1-axial_0j.txt 0.021318904968 0.1272 97.3523531583 0.864208980016
./Samples/DMsimp_s_spin1-axial_ee_0j.txt 0.010659482475 0.0978 37.4256561594 2.24554648455
./Samples/DMsimp_s_spin1-axial_mm_0j.txt 0.010659482475 0.1481 56.6742298283 1.48558733966
./Samples/DMsimp_s_spin1-vector_0j.txt 0.0393553945622 0.1227 173.357758169 0.484974465556
./Samples/DMsimp_s_spin1-vector_0j_10_1.txt 0.299915034487 0.0076 81.8288180094 1.02663162079
./Samples/DMsimp_s_spin1-vector_0j_50_1.txt 0.217641714014 0.0263 205.490777121 0.410468667038
./Samples/DMsimp_s_spin1-vector_0j_100_1.txt 0.168503454019 0.0490 296.414425965 0.283063751571
./Samples/DMsimp_s_spin1-vector_0j_500_1.txt 0.0419776079768 0.1291 194.553199914 0.430585232638
./Samples/DMsimp_s_spin1-vector_0j_1000_1.txt 0.0112635212815 0.1509 61.0179864735 1.37606381584
./Samples/DMsimp_s_spin1-vector_0j_qcd.txt 0.06788 0.1251 304.8551892 0.275613171719
./Samples/DMsimp_s_spin1-vector_1j.txt 0.0666819164253 0.0752 180.037839919 0.469328247868
./Samples/DMsimp_s_spin1-vector_1j_qcd.txt 0.1251 0.1125 505.082719261 0.166834705882
./Samples/DMsimp_s_spin1-vector_1j_qcd_Qcut10.txt 0.1251 -1.0000 -4491.09 ERR
./Samples/DMsimp_s_spin1-vector_1j_qcd_Qcut100.txt 0.1251 -1.0000 -4491.09 ERR
./Samples/DMsimp_s_spin1-vector_ee_0j.txt 0.0196775973111 0.1003 70.8545020699 1.1868190876
./Samples/DMsimp_s_spin1-vector_mm_0j.txt 0.0196775973111 0.1572 111.050126873 0.759155804105
Run 30/11/2017
./Samples/DMsimp_s_spin0P_0j_1000_1.txt 5.28486997514e-06 0.0064 0.00121425172549 69140.5077591
./Samples/DMsimp_s_spin0P_0j_100_1.txt 5.29739882874e-05 0.0181 0.0344219678493 2447.54969572
./Samples/DMsimp_s_spin0P_0j_100_1_QCD.txt 0 -1.0000 -0.0 99999.9996275
./Samples/DMsimp_s_spin0P_0j_10_1.txt 9.8179223129e-05 0.0024 0.00845912186479 9942.6415468
./Samples/DMsimp_s_spin0P_0j_10_1_QCD.txt 0.0001969 -1.0000 -7.06871 ERR
./Samples/DMsimp_s_spin0P_0j_500_1.txt 2.10082494163e-05 0.0126 0.00950287154097 8841.97930621
./Samples/DMsimp_s_spin0P_0j_50_1.txt 8.0240824427e-05 0.0077 0.0221809710964 3788.30687204
./Samples/DMsimp_s_spin0P_0j_50_1_QCD.txt 0.0001702 -1.0000 -6.11018 ERR
./Samples/DMsimp_s_spin0S_0j_1000_1.txt 6.47028916275e-06 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0S_0j_100_1.txt 9.660359773e-05 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0S_0j_10_1.txt 0.0002873577256 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0S_0j_500_1.txt 2.66404256188e-05 0.0000 0.0 99999.9996275
./Samples/DMsimp_s_spin0S_0j_50_1.txt 0.000177445033 0.0000 0.0 99999.9996275

5. Reinterpretation (optional)

6. Questions / issues

Note: See TracWiki for help on using the wiki.