(r24) CMSSusyAnalysisResults < CMS

CMS Web>NewPhysics>CmsSusy>CMSSusyAnalysisResults (revision 24) (raw view)~~EditAttach~~
%TOC%

---++ Information
---++++ Analysers:

   * Janos Karancsi - PhD student (ATOMKI & UNIDEB, Debrecen, HU)
   * Viktor veszpremi - PhD (Wigner, Budapest, HU)
   * Petar Maksimovic - Prof. PhD (Johns Hopkins University, Baltimore, US)
   * Brandon Chiarito - Student (Johns Hopkins University, Baltimore, US)

---++++ Links:

For a lot of results we share a common elog with the JHU group: [[http://hep.pha.jhu.edu:8080/susy/][JHU elog]]

---+++++ Twikis:
   * [[https://twiki.cern.ch/twiki/bin/viewauth/CMS/StandardModelCrossSectionsat13TeV][SM cross sections]]
   * MC production: Run II MC Requests [[https://twiki.cern.ch/twiki/bin/viewauth/CMS/MCRequestsFor13TeVProduction][[All requests]]] [[https://twiki.cern.ch/twiki/bin/view/CMS/SMPGENThirteenTeVRequests][[Main requests]]]  --> Spreadsheets: [[https://docs.google.com/spreadsheets/d/14cbc9kulE9Rsu0mDkQHnDR7zrYz07bMyonT48rfmGCo/edit?pli=1#gid=0][[1]]], SUSY Slides:  [[https://indico.cern.ch/event/304262/session/5/contribution/9/material/slides/0.pdf][MC plans]]

---+++++ Useful info:

Papers, Notes:
|  *Links*  |  *Date*  |  *Author*  |  *Title*  |  *Comment*  |
| [[http://arxiv.org/pdf/physics/0308063v2][arXiv:physics/0308063]] | 2003 | G. Punzi | Sensitivity of searches for new signals and its optimization ||
| [[http://cms-physics.web.cern.ch/cms-physics/public/JME-09-001-pas.pdf][CMS PAS JME-09-001]] | 2009 | CMS | A Cambridge-Aachen (C-A) based Jet Algorithm for boosted top-jet tagging ||
| [[http://arxiv.org/pdf/1011.2268v3][arXiv:1011.2268]] | 2010 | J. Thaler, K. V. Tilburg | Identifying Boosted Objects with N-subjettiness | tau variables (N-sj.), Mistag rates, tagging efficiencies |
| [[http://arxiv.org/pdf/1006.2833v2][arXiv:1006.2833]] | 2010 | T. Plehn et. al | Stop Reconstruction with Tagged Tops ||
| [[http://arxiv.org/pdf/1007.2221.pdf][]] | 2010 | K. Rehermann, B. Tweedie | Efficient Identification of Boosted Semileptonic Top Quarks at the LHC | Mini Isolation, new default for SUSY group. |
| [[http://cds.cern.ch/record/1545285/files/B2G-12-005-pas.pdf][CMS-PAS-B2G-12-006]] ([[https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsB2G12005][Twiki]]) | 2013 | CMS | Search for Anomalous Top Quark Pair Production in the Boosted All-Hadronic Final State | Mistag rate of top tagging algo, Multijet bkg estimation |
| [[http://cds.cern.ch/record/1543467/files/B2G-12-006-pas.pdf][CMS-PAS-B2G-12-006]] ([[https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsB2G12006][Twiki]]) | 2013 | CMS | Search for ttbar resonances in semileptonic final state | 2D cut for semi-leptonic top |
| [[https://cds.cern.ch/record/1647419/files/JME-13-007-pas.pdf][CMS-PAS-JME-13-007]] | 2014 february | CMS | Boosted Top Jet Tagging at CMS | |
| [[http://cms.cern.ch/iCMS/analysisadmin/cadi?ancode=SUS-14-007][SUS-14-007]] ([[https://indico.cern.ch/event/373813/contribution/3/material/slides/0.pdf][Approval]]) | 2015 | Razor group | Search for SUSY with boosted W's using the razor variables | |


Presentations:
   * [[https://indico.cern.ch/event/376578/contribution/0/7/material/slides/0.pdf][Boosted top trigger - AK8PFJet360_TrimMass30]] - ,March 26, 2015
   * [[https://indico.cern.ch/event/364108/contribution/3/material/slides/0.pdf][Top Tagging at 13 TeV]] - JMAR Meeting, March 11, 2015
   * [[https://indico.cern.ch/event/331025/contribution/1/material/slides/0.pdf][Update on Differential tt Cross Section with Boosted Tops]] - Top Cross Section Meeting, July 24, 2014
   * [[https://indico.cern.ch/event/315003/contribution/1/material/slides/0.pdf][Jet Substructure Hands-on Advanced Tutorial Session]] - LPC, June 18-19, 2014
   * [[https://indico.cern.ch/event/301310/session/5/contribution/53/material/slides/0.pdf][Physics performance and potentials of boosted topologies]] - CMS Physics week, February 24, 2014
   * [[https://indico.cern.ch/event/239415/contribution/0/material/slides/0.pdf][QCD-Multijet Background Modelling Using The ABCD Method]] - ttH meeting, February 28, 2013
   * [[http://wwwiexp.desy.de/users/auterman/talks/20080311_HADSUSY_backgrounds.pdf][Background Estimation Methods (ABCD - slide 8)]]

---++ Important results
---+++ Links to latest results 

   * [[https://indico.cern.ch/event/377190/contribution/0/material/slides/1.pdf][15/04/15 - Razor meeting]]
   * [[http://hep.pha.jhu.edu:8080/susy/234][15/03/29 - JHU elog #234]] - Background estimation (NTop sideband, tau32 inversion)
   * [[http://hep.pha.jhu.edu:8080/susy/232][15/03/16 - JHU elog #232]] - Background Estimation (DPhi/NTop sideband - tau32/mass inversion)
   * [[http://hep.pha.jhu.edu:8080/susy/231][15/03/09 - JHU elog #231]] - 
   * [[https://indico.cern.ch/event/377183/contribution/1/material/slides/0.pdf][15/02/25 - Razor Meeting]]
   * [[http://hep.pha.jhu.edu:8080/susy/230][15/02/24 - JHU elog #230]] - Distributions (no cuts - cross section reweighted)
   * [[http://hep.pha.jhu.edu:8080/susy/229][15/02/10 - JHU elog #229]] - Ntuples (Feb 13 Feb 18)
   * [[http://hep.pha.jhu.edu:8080/susy/228][15/01/20 - JHU elog #228]] - Distributions (2 cuts - HTall, DPhi cuts, R)
   * [[http://hep.pha.jhu.edu:8080/susy/226][15/01/17 - JHU elog #226]] - Distributions (1 cut - HTall, DPhi cuts)
   * [[http://hep.pha.jhu.edu:8080/susy/225][15/01/14 - JHU elog #225]] - Distributions (1/2D, no cuts)
   * [[http://hep.pha.jhu.edu:8080/susy/224][15/01/12 - JHU elog #224]] - first plots (from Brandon's ntuple, my new one)
   * [[http://hep.pha.jhu.edu:8080/susy/220][14/12/23 - JHU elog #220]] - warmup

---+++ Object/Event selections
---++++ Object definitions
All objects come from MiniAOD slimmed Pat collections: [[https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookMiniAOD#High_level_physics_objects][MiniAOD objects]]

---+++++ Jets
 Anti-kt algorithm jets with cone size: R=0.8, B2G baseline cuts:
   * pT > 100 GeV/c
   * | eta | < 4

---+++++ Electrons
 gedGsfElectrons with cuts (B2G baseline: pT > 30 GeV/c, |eta | < 2.5):
   * pT > 35 GeV/c
   * | eta | < 2.5

---+++++ Muons
 muons with cuts (B2G baseline: pT > 30 GeV/c, |eta | < 2.5)
   * pT > 45 GeV/c
   * | eta | < 2.1
   * tight ID

---+++++ MET
 type1 PF MET

---+++++ Event Selection
Exactly 2 AK8 Jets with cuts (loose hadronic top selection):
   * pT > 400 GeV/c
   * tau3 / tau2 < 0.75
   * Pruned Mass > 140 Gev/c^2
Signal selection cuts (See below)
   * R (Razor variable) - [[https://twiki.cern.ch/twiki/bin/view/CMSPublic/RazorLikelihoodHowTo][Recipe]]
      * Calculated for jets with | eta | < 3
   * H_T,all = Sum ( AK8Jet pT ) + Sum ( Lepton pT ) + MET
   * DPhi = | DeltaPhi ( top_had1, top_had2 ) |

---+++ Signal Selection

---++++ Cut flow chart 
---+++++ Best cuts using Smin maximisation script - New
   * 2 Hadronic tops

Maximum Smin, for all combination of cuts:
| *Best cut* |  *T5ttttDeg_4bodydec_mGo1300*  |  *T5ttttDeg_3bodydec_mGo1300*  |  *TTJets*  |  *WJets_HT*  |  *ZJets_HT*  |  *QCD_HT*  |  *T_tW*  |  *All bg*  |  *Smin = Eff./(5/2+sqrt(B)) *  |                                                           
| No cut |  23.47 |  24.17 |  4213.88 |  104.94 |  7.60 |  14906.68 |  72.73 |  19305.84 |  0.00707 |                         
| HTall>2050 |  15.74 |  15.24 |  283.71 |  30.66 |  2.09 |  1306.98 |  15.65 |  1639.08 |  0.01467 |                         
| DPhi<2.3 |  10.80 |  10.70 |  107.54 |  3.18 |  0.43 |  258.67 |  4.74 |  374.57 |  0.02025 |                               
| R>0.45 |  10.35 |  10.11 |  2.29 |  0.76 |  0.15 |  0.00 |  0.36 |  3.55 |  0.09538 |                                       
| DPhi<2.7 && HTall>2050 |  11.78 |  11.38 |  82.04 |  5.30 |  0.57 |  258.67 |  5.10 |  351.69 |  0.02216 |                  
| R>0.45 && HTall>1200 |  10.34 |  10.11 |  2.29 |  0.76 |  0.15 |  0.00 |  0.36 |  3.55 |  0.09538 |                         
| R>0.36 && DPhi<2.75 |  13.46 |  13.33 |  2.29 |  1.11 |  0.30 |  0.00 |  0.00 |  3.70 |  0.12461 |                          
| R>0.36 && DPhi<2.75 && HTall>1100 |  13.46 |  13.33 |  2.29 |  1.11 |  0.30 |  0.00 |  0.00 |  3.70 |  0.12461 | 

Using cuts (chosen by hand):
   * R > 0.4
   * DPhi < 2.8

| *Cut* |  *T5ttttDeg_4bodydec_mGo1300*  |  *T5ttttDeg_3bodydec_mGo1300*  |  *TTJets*  |  *WJets_HT*  |  *ZJets_HT*  |  *QCD_HT*  |  *T_tW*  |  *All Background*  |  *Smin = Eff./(5/2+sqrt(B)) *  |                                                        
| No Cut |  23.47 |  24.17 |  4213.88 |  104.94 |  7.60 |  14906.68 |  72.73 |  19305.84 |  0.00707 |                         
| DPhi < 2.8 |  17.94 |  18.34 |  580.83 |  14.60 |  1.72 |  1493.78 |  14.89 |  2105.82 |  0.01568 |                         
| R > 0.4 |  12.77 |  12.50 |  9.81 |  1.46 |  0.38 |  0.00 |  0.73 |  12.38 |  0.08590 |                                     
| R > 0.4 + DPhi < 2.8 |  12.03 |  11.98 |  2.29 |  0.86 |  0.27 |  0.00 |  0.00 |  3.41 |  0.11400 |                                                                                          

---+++ Background Estimation
---++++ Signal selection cuts:
   * 2 Hadronic tops
   * R > 0.4
   * DPhi < 2.8

---++++ Introducing the method
ABCD Method:
   * Define NTop<2 sideband, where the sideband  is defined as 2 loose top-like jets:
      * Has at least 2 AK8 jets with pt > 400, prunedMass > 140
      * Sideband cut: NTopHad < 2 (at least one tau32 cut is inverted)
   * Measure event yields in the R<0.4 sideband both for A) NTop<2 and C) NTop==2
   * Fit NTop<2 sideband (A, B) with an exponential decay curve: exp(const + x * slope)
      * Extract fit parameters and their error: const, slope
   * Scale function by the event yield ratio C/A to get Signal region prediction: D

Recipe to plot very similar distributions by hand (for leading and subleading jets):
<pre>
root -l
TChain c("c");
c.Add("/data/gridout/jkarancs/SusyAnalysis/B2G/TTreeNtuple/Feb28_edm_Feb20/WJetsToLNu_HT-600toInf/*.root/B2GTTreeMaker/B2GTree");
c.Draw("evt_R","jetAK8_Pt[0]>400&&jetAK8_Pt[1]>400&&jetAK8_prunedMass[0]>100&&jetAK8_prunedMass[1]>100&&evt_NTopHad<2");
c.Draw("evt_R","evt_NTopHad==2","SAME")
c.Draw("evt_R","evt_NTopHad==2&&evt_TTHadDPhi<2.8","SAME")
</pre>

---++++ Plots - Extrapolation of R shape using an exponential fit in the NTop<2 SideBand

| *TTbar* | *QCD* | *Comment* |
|<a href="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_TTBar.png"><img src="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_TTBar.png" height="300"> |<a href="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_QCD_HT.png"><img src="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_QCD_HT.png" height="300"> | R shape fit with exponential in the Ntop<2 sideband and NTop==2 signal bands. Range [0.15,1.0] for sideband, [0.15,0.4] for signal band. |
| *W+Jets* | *Z+Jets --> NuNu* | *Comment* |
|<a href="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_WJets_HT.png"><img src="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_WJets_HT.png" height="300"> |<a href="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_ZJets_HT.png"><img src="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_ZJets_HT.png" height="300"> | R shape fit with exponential in the Ntop<2 sideband and NTop==2 signal bands. Range [0.15,1.0] for sideband, [0.15,0.4] for signal band, except for WJets for which low range is 0.2. |
| *G + Jets* | *Single top (tW channel)* | *Comment* |
|<a href="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_GJets_HT.png"><img src="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_GJets_HT.png" height="300"> |<a href="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_T_tW.png"><img src="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_T_tW.png" height="300"> | R shape fit with exponential in the Ntop<2 sideband and NTop==2 signal bands. Range [0.15,1.0] for sideband, [0.15,0.4] for signal band. |
| *DY + Jets* | *Comment* |
|<a href="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_DYJets_HT.png"><img src="%ATTACHURLPATH%/R_NTopHad_DPhiBelow2p8_DYJets_HT.png" height="300"> | R shape fit with exponential in the Ntop<2 sideband and NTop==2 signal bands. Range [0.15,1.0] for sideband, [0.15,0.4] for signal band. |

---++++ Table
To see how much background is there with [[http://www.grid.kfki.hu/twiki/bin/view/CMS/CMSSusySimulation][given]] statistics:
   * No lepton veto applied

| *Sample* | *A (R<0.4, SB)* | *B (R fit, SB) pred.* | *B (R>0.4, SB)* | *C (R<0.4, Sig.B.)* | *D (R fit, Sig.B.) pred.* | *D = B*C/A pred.* | *D = B (R fit, SB) * C/A pred.* | *D (R>0.4, Sig.B.) obs.* | *Statistics in D* |
| TTJets |  291.89 +- 9.77 |  1.26 +- 0.34 |  2.29 +- 0.86 |  578.22 +- 13.75 |  1.68 +- 0.40 |  4.53 +- 1.72 |  2.49 +- 0.68 |  2.29 +- 0.86 |  7 |
| WJets_HT |  37.42 +- 1.37 |  1.90 +- 0.28 |  2.12 +- 0.33 |  13.59 +- 0.83 |  0.66 +- 0.19 |  0.77 +- 0.13 |  0.69 +- 0.11 |  0.86 +- 0.21 |  17 |
| ZJets_HT |  3.76 +- 0.19 |  0.75 +- 0.07 |  0.92 +- 0.09 |  1.38 +- 0.11 |  0.25 +- 0.05 |  0.34 +- 0.05 |  0.28 +- 0.04 |  0.27 +- 0.05 |  29 |
| QCD_HT |  4260.50 +- 170.30 |  0.08 +- 0.11 |  0.00 +- 0.00 |  1493.78 +- 100.84 |  0.03 +- 0.07 |  0.00 +- 0.00 |  0.03 +- 0.04 |  0.00 +- 0.00 |  0 |
| T_tW |  12.72 +- 2.12 |  0.03 +- 0.05 |  0.00 +- 0.00 |  14.89 +- 2.30 |  0.06 +- 0.11 |  0.00 +- 0.00 |  0.03 +- 0.06 |  0.00 +- 0.00 |  0 |
| SUM |  4606.28 +- 170.60 |  4.02 +- 0.46 |  5.33 +- 0.93 |  2101.85 +- 101.80 |  2.68 +- 0.47 |  5.64 +- 2.99 |  3.51 +- 0.70 |  3.41 +- 0.89 |  |

---++++ Main backgrounds identified on tested samples - TTBar, QCD, W+Jets
http://hep.pha.jhu.edu:8080/susy/230

---++++ QCD background estimation technique
http://hep.pha.jhu.edu:8080/susy/115

---+++ Estimation of Uncertainties


-- Main.JanosKarancsi - 2015-02-26
Attachments
Topic attachments
I	Attachment	History	Action	Size	Date	Who
png	R_NTopHad_DPhiBelow2p8_DYJets_HT.png	r1	manage	14.1 K	2015-03-14 - 14:53	JanosKarancsi
png	R_NTopHad_DPhiBelow2p8_GJets_HT.png	r1	manage	14.1 K	2015-03-14 - 14:53	JanosKarancsi
png	R_NTopHad_DPhiBelow2p8_QCD_HT.png	r1	manage	15.0 K	2015-03-14 - 14:54	JanosKarancsi
png	R_NTopHad_DPhiBelow2p8_TTBar.png	r1	manage	14.0 K	2015-03-14 - 14:54	JanosKarancsi
png	R_NTopHad_DPhiBelow2p8_T_tW.png	r1	manage	14.7 K	2015-03-14 - 14:54	JanosKarancsi
png	R_NTopHad_DPhiBelow2p8_WJets_HT.png	r1	manage	14.3 K	2015-03-14 - 14:54	JanosKarancsi
png	R_NTopHad_DPhiBelow2p8_ZJets_HT.png	r1	manage	14.1 K	2015-03-14 - 14:54	JanosKarancsi
png	R_NTopSideBand_WJets_Ht600.png	r1	manage	26.4 K	2015-03-11 - 11:26	JanosKarancsi
png	R_NTopSideBand_ZJetsToNuNu_Ht600.png	r1	manage	19.7 K	2015-03-11 - 11:26	JanosKarancsi
Topic revision: r24 - 2015-05-26 - JanosKarancsi
CMS (CERN)
Copyright &© by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback