First studies towards topquark pair differential cross section measurement in the dilepton channel at s = 13 TeV with the CMS detector Mykola Savitskyi, Maria Aldaya, Carmen Diez Pardos et al. DESYCMS Top Group, Hamburg 1
The top quark is special Heaviest elementary particle known to date: mt 173 GeV Mass close to scale of electroweak symmetry breaking (EWSB): yt important role in EWSB? Decays before hadronizing: 5 1025 s << QCD unique window on bare quark Sensitive to physics phenomena beyond the Standard Model (BSM): new physics may preferentially couple/decay to top Major source of background for many Higgs and BSM searches 2
Why measure differentially? Precise understanding of top quark distributions is crucial: Precision tests of perturbative QCD for top quark production at different phase space regions Tune and test theory predictions and models: potential to reduce signal modelling systematics Enhance sensitivity to BSM physics LHC is a top factory Several millions of top quark pair (tt) events produced already in Run I ( s = 7, 8 TeV) Run II: expect much larger data sets at s = 13(14) TeV important to measure top quark distributions with very high precision 3
Top Pair Production and Decay Production cross sections: Weak interaction: Top decay Energy, TeV 8 13 σinclnnlo+nnll(tt), pb 252.89 831.76 t W+b ~ 100 % Strong interaction: Top pairs gg fusion qq annihilation LHC Energy, TeV 8 13 gg tt ~82% ~90% qq tt ~18% ~10% Samples are classified according to Wdecay: dilepton (ℓℓ), lepton+jets (ℓ+jets), all jets 4
Selection in Dilepton Channel 1. Trigger conditions 2. Lepton pair selection: opposite charge e and μ isolation criteria pt > 20 GeV, η < 2.4 invariant mass: mll > 20 GeV 3. Exclusion of Z region: in ee and μμ: 76 GeV < mll < 106 GeV 4. Presence of two jets (antikt R=0.4 ) with pt > 30 GeV and η < 2.4 5. Missing ET > 40 GeV in ee and μμ 6. At least one b tagged jet 7. Meaningful solution for kinematic event reconstruction Largest backgrounds: other tt, Single top, DrellYan events 5
Kinematic Reconstruction Measured input: 2 jets, 2 leptons, MET Unknowns: pν, pν 6 Constraints: > mt, mt 2 > mw(+), mw() 2 > (pν+ pν)t = MET 2 Reconstructing each event 100 times and smearing inputs by their resolution: > top mass fixed to 172.5 GeV > W mass at RECO level smeared accordingly to W mass distribution > Jet and lepton energies are corrected for detector effects Consider weighted average of solutions for all smeared points: 100 1 top ptop = w ( p x, y, z i x, y, z )i w i=0 6
Selection in Dilepton Channel Channel: combined (ee+eμ+μμ) Pseudodata: poissonsmeared sum of all signal and background simulation samples Normalized to: Lint = 5.0 fb1 Observable: Lepton η Signal: MadGraph+Pythia8 Selection steps: 1. Trigger conditions 2. Lepton pair selection: opposite charge e and μ isolation criteria pt > 20 GeV, η < 2.4 invariant mass: mll > 20 GeV 7
Selection in Dilepton Channel Channel: combined (ee+eμ+μμ) Pseudodata: poissonsmeared sum of all signal and background simulation samples Normalized to: Lint = 5.0 fb1 Observable: Lepton η Signal: MadGraph+Pythia8 Selection steps: 1. Trigger conditions 2. Lepton pair selection: opposite charge e and μ isolation criteria pt > 20 GeV, η < 2.4 invariant mass: mll > 20 GeV 3. Exclusion of Z region: in ee and μμ: 76 GeV < mll < 106 GeV 7
Selection in Dilepton Channel Channel: combined (ee+eμ+μμ) Pseudodata: poissonsmeared sum of all signal and background simulation samples Normalized to: Lint = 5.0 fb1 Observable: Lepton η Signal: MadGraph+Pythia8 Selection steps: 1. Trigger conditions 2. Lepton pair selection: opposite charge e and μ isolation criteria pt > 20 GeV, η < 2.4 invariant mass: mll > 20 GeV 3. Exclusion of Z region: in ee and μμ: 76 GeV < mll < 106 GeV 4. Presence of two jets: with pt > 30 GeV and within η < 2.4 7
Selection in Dilepton Channel Channel: combined (ee+eμ+μμ) Pseudodata: poissonsmeared sum of all signal and background simulation samples Normalized to: Lint = 5.0 fb1 Observable: Lepton η Signal: MadGraph+Pythia8 Selection steps: 1. Trigger conditions 2. Lepton pair selection: opposite charge e and μ isolation criteria pt > 20 GeV, η < 2.4 invariant mass: mll > 20 GeV 3. Exclusion of Z region: in ee and μμ: 76 GeV < mll < 106 GeV 4. Presence of two jets: with pt > 30 GeV and within η < 2.4 5. Missing ET > 40 GeV in ee and μμ 7
Selection in Dilepton Channel Channel: combined (ee+eμ+μμ) Pseudodata: poissonsmeared sum of all signal and background simulation samples Normalized to: Lint = 5.0 fb1 Observable: Lepton η Signal: MadGraph+Pythia8 Selection steps: 1. Trigger conditions 2. Lepton pair selection: opposite charge e and μ isolation criteria pt > 20 GeV, η < 2.4 invariant mass: mll > 20 GeV 3. Exclusion of Z region: in ee and μμ: 76 GeV < mll < 106 GeV 4. Presence of two jets: with pt > 30 GeV and within η < 2.4 5. Missing ET > 40 GeV in ee and μμ 6. At least one b tagged jet 7
Selection in Dilepton Channel Channel: combined (ee+eμ+μμ) Pseudodata: poissonsmeared sum of all signal and background simulation samples Normalized to: Lint = 5.0 fb1 Observable: Lepton η Signal: MadGraph+Pythia8 Selection steps: 1. Trigger conditions 2. Lepton pair selection: opposite charge e and μ isolation criteria pt > 20 GeV, η < 2.4 invariant mass: mll > 20 GeV 3. Exclusion of Z region: in ee and μμ: 76 GeV < mll < 106 GeV 4. Presence of two jets: with pt > 30 GeV and within η < 2.4 5. Missing ET > 40 GeV in ee and μμ 6. At least one b tagged jet 7. Meaningful solution for kinematic event reconstruction 7
Selection in Dilepton Channel Channel: combined (ee+eμ+μμ) Pseudodata: poissonsmeared sum of all signal and background simulation samples Normalized to: Lint = 5.0 fb1 Observable: Lepton η Signal: MadGraph+Pythia8 Table given after all cuts: Process Fraction, % tt Signal 78.7 tt Other 14.1 tw 3.3 W+Jets 0.2 DY ee/μμ 2.4 DY ττ 0.6 tt+z/w/γ 0.7 7
Control Distributions Channel: combined (ee+eμ+μμ), Lint = 5.0 fb1 Plots given after full selection Pseudodata: poissonsmeared sum of all signal and background simulation samples 8
Differential Cross Section For each observable X the normalized differential cross section in the ith bin is defined as: i i ( 1 dσ 1 N events σ dx = σ Δ i L X ) N ievents number of events after background subtraction, efficiency, acceptance and bintobin migration correction σ total tt cross section in same phase space L integrated luminosity ΔiX bin width Phase space definition: Top quarks, tt system (obtained via kinematic reconstruction of event) extrapolated to full phase space after corrections for detector and hadronization effects Leptons or bjets measured in visible phase space (leptons: pt > 20 GeV, η < 2.4; jets: pt > 30 GeV, η < 2.4) after correction for detector effects Bintobin migrations are reduced by binning optimization and corrected by unfolding 9
Binning and Migrations Migration effects studied by: rec & gen Ni pi= rec Ni & gen N rec i si = εi = N gen i & sel N rec i N all generated i purity: sensitive to migrations to ith bin stability: sensitive to migrations out of ith bin efficiency in ith bin Binning criteria: stability or purity ~0.5 Response matrices are constructed from signal MC phase space rebin unfolding eμ, Private Work 10
Measured Cross Sections: leptons Normalization allows to: reduce systematic uncertainties perform shape comparisons of different theory models to data Systematic uncertainty: include JES and JER uncertainties for now 11
Measured Cross Sections: tops Normalization allows to: reduce systematic uncertainties perform shape comparisons of different theory models to data Systematic uncertainty: include JES and JER uncertainties for now 12
Measured Cross Sections: ttpair Sensitive to high order contributions Sensitive to resonances Normalization allows to: reduce systematic uncertainties perform shape comparisons of different theory models to data Systematic uncertainty: include JES and JER uncertainties for now 13
Conclusions & Outlook Top quark pair differential cross section measurements: Essential for constraining the SM Ideal probe for looking for new physics beyond the SM Needed for tuning of PDF sets and modern art Monte Carlo event generators First studies towards measurements at 13 TeV were presented: Not latest status & results currently Top secret at CMS and work in progress! Optimize binning with final configuration of the data analysis Evaluate systematic uncertainties Compare to different theory predictions 14
Conclusions & Outlook Top quark pair differential cross section measurements: Essential for constraining the SM Ideal probe for looking for new physics beyond the SM Needed for tuning of PDF sets and modern art Monte Carlo event generators First studies towards measurements at 13 TeV were presented: Not latest status & results currently Top secret at CMS and work in progress! Optimize binning with final configuration of the data analysis Evaluate systematic uncertainties Compare to different theory predictions Thank you for your attention! :) 14
Backup
LHC Collider Run 1: Lint 25 fb1 in 20102012 at 7/8 TeV Run 2: Phase 0 in 20152017 design energy: s=13 TeV~14 TeV nominal luminosity: L~ 1 1034cm2s1 bunch spacing: 25 ns pile up: ~49 Lint per year ~ 45 fb1 Planned to collect: Lint ~ 75 100 fb1 15 s=13 TeV in 2015 p p
The CMS Detector General purpose 4 detectors: Tracker: Detection and momentum measurement for charged particles Calorimeter: Identification and energy measurement of jets and electrons Muon system: Identification and momentum measurement of muons Phase 0 (20152017) upgrade: 16 Complete muon coverage Colder tracker Photodetectors in HCAL New beampipe and infrastructure updates
Visible Phase Space Definition Object definition at generator level: > particles after radiation and hadronization > leptons: from Wdecay > jets: antikt (with cone of ΔR=0.5) algorithm > bjets: identified by Bhadrons Directly measured quantities: leptons and bjets Visible particle LVL Phase Space Correct for: Detector Effects Measure at visible Phase Space: Leptons: pt > 20 GeV, η < 2.4 Jets: pt > 30 GeV, η < 2.4 17
Unfolding Unfolding techniques correct migrations between bins Response matrix (A): represents binbybin correlations Unfolding problem is transformed to unfolding minimization problem: regularization 2 2 A x )T COV 1 χ =( N ( N A x ) τ K ( x ) N N: BG corrected data x: unfolded result Nonphysical fluctuations removed by means of the regularization: continuous regularization parameter selected at minimum of average global correlation eμ, Private Work eμ, Private Work In this measurement regularized unfolding is used 18