= 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 == * CMS contact person: Nicholas Charles Smith * Center of mass energy: 13 TeV * Integrated luminosity: 35.9 fb−1 * Link to arXiv: none * Link to the collaboration page: http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/EXO-16-052/ * Link to CDS: http://cds.cern.ch/record/2264690?ln=en * Link to Inspire: http://inspirehep.net/record/1599663 * Other useful links... 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 * pT > 18 GeV, |eta| < 2.3 * DecayModeFinding ? * Loose isolation < 5 GeV 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 3. MC generator program (name + version): MAD-GRAPH5_aMC@NLO 2.2.2 4. Shower program (name + version): PYTHIA 8.205 5. 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 === * http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/EXO-16-052/CMS-PAS-EXO-16-052_Table_004.png == 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_1j_qcd_Qcut100.txt || 0.1251 || 0.1134 || 509.470135635 || 0.165344589911 || || ./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 ==