1 Search for high-mass dilepton resonances with the ATLAS detector, Michigan State University Experimental Particle Physics Seminar University of Pennsylvania, 02.28.2012
2 Overview 1. Why are we looking for high-mass dilepton resonances? 2. 1fb-1 analysis (Phys. Rev. Lett. 107, 272002 (2011)) - event selection - backgrounds - high energy electrons and muons - signal search and limit setting 3. 5 fb-1 analysis (in preparation) - updates
Why are we looking for Physics beyond the SM? The Standard Model of Particle Physics is a very successful theory, but cannot be the end of the story... For example, it doesn't - have a dark matter candidate - explain, why gravity is so weak compared to the other fundamental forces Also: Before we find the Higgs, we cannot be sure of how electroweak symmetry is broken 3
Search for dilepton resonances (ee/μμ) q/g q/g 4 l +? l -
5 Search for dilepton resonances Dilepton resonances have been the window to a better understanding of elementary particles and forces before... Here be dragons
6 Search for dilepton resonances Dilepton resonances could be a signature of - new heavy gauge boson in the E6 model (Grand Unified Theory model) ----> spin-1 - excited Kaluza-Klein mode of the Randall Sundrum graviton ----> spin-2...and many others (resonance search is fairly model independent) ---------------------------------------Benchmark model: Sequential Standard Model Z' (same couplings as Z boson), not theoretically motivated
7 Additional symmetries Dilepton resonances could be a signature of - new heavy gauge boson in the E6 model (Grand Unified Theory model) (Phys. Rev. D 34 (1986), arxiv:0801.1345v3) GUT theories: Unification of electroweak and strong forces at high energies ----> 1 overall symmetry, which breaks down at lower energies 2 additional U(1) groups lead to Z' Several motivated choices of θe6
8 Randall Sundrum graviton Dilepton resonances could be a signature of - excited Kaluza-Klein mode of the Randall Sundrum graviton (arxiv:hep-ph/9905221v1) 1 finite warped extra dimension, 2 branes Standard brane 2 parameters: - m1 (first excitation) - k (curvature) Planck brane G(x) 1 extra dimension 0.1 Only gravitons propagate to extra brane ----> wave functions are suppressed away from extra brane ----> gravity is weak! Finite extra dimension ----> excitation like in harmonic oscillator possible ----> Kaluza-Klein tower of massive graviton states
9-1 1 fb analysis
10 Event Selection electrons - ATLAS data quality (stable beam, functioning subdetectors etc.) - pick two standard electrons, central (η < 2.47), ET > 25 GeV - require track and EM shower shape cuts ( Medium ) - require hit in innermost layer of detector - leading electron should be isolated (calorimeter based) ---> energy in a cone of 0.2 around the electron < 7 GeV - form invariant mass, require mass > 70 GeV Main background: Z/γ* Acceptance (Z', 1.5 TeV): 67%
11 Event Selection muons - ATLAS data quality (stable beam, functioning subdetectors etc.) - pick two standard combined muons, pt > 25 GeV - require 3 hits in 3 muon spectrometer layers, no overlap barrel-endcap, veto misaligned chambers - distance from primary vertex needs to be small - require tracks to be isolated - muons must have opposite charge - form invariant mass, require mass > 70 GeV Main background: Z/γ* Acceptance (Z', 1.5 TeV): 42%
12 Backgrounds Backgrounds with two prompt electrons/muons: - Drell Yan - WW,WZ,ZZ - ttbar (dileptonic decay) Backgrounds with QCD jets, which can fake prompt leptons - W+jets - QCD multijet production ----> what is the fake rate at high energies? Drell Yan is dominant background by far, all background except for QCD multijet taken from simulated samples
QCD multijet background electrons 13 Jets can fake electrons. How large is the fake rate? 1. Baseline Method Reverse Identification - dijet shape from reverted electron identification cuts - extrapolation to high invariant masses by fitting with empirical function - normalization by 2-component template fit 2. Cross-check and systematic uncertainties Isolation fit method Fake rate estimate - use calorimeter isolation distributions - measure probability for jet-like - fit signal/background templates from data for 1st and 2nd electron - system of equations to avoid double counting objects to pass Z' selection (η, ET) - apply fake rate on normalization sample (Z' selection on leading, jet selection on second electron)
14 QCD multijet background muons QCD multijet background much smaller for muons 1. Shape: Anti-track-isolated data (0.1-1.0) 2. Normalization: Ratio of isolated (0.00-0.05) /anti-isolated (0.1-1.0) dimuon events in QCD (heavy flavor) simulated samples QCD from simulation
High pt leptons 15 Looking for resonances at high invariant masses ---> need to understand properties of highly energetic objects in ATLAS Very small control sample, handles: calibration runs, cosmics, Tag-and-Probe around Z pole ---> extrapolation, simulation
16 High pt leptons resolution Electrons (Resolution 1.1 1.8% at 1 TeV): - energy measurement from electromagnetic calorimeter - resolution at high energies dominated by constant term - improvement at high energies of linearity and resolution shown in calibration runs Muons (Resolution > 15% at 1 TeV): - pt measurement from hits in inner detector and muon spectrometer - require stringent cuts on number of hits, veto misaligned areas - measured (as a function of pt) using cosmics, magnet off runs, overlap regions, inner detector vs. muon spectrometer comparisons, Z peak
17 Efficiencies and scale factors Determine trigger, reconstruction and identification efficiencies in data with Tag-and-Probe l+ - allows to get relatively unbiased control sample by applying strict cuts on Tag and test efficiency on Probe (p.ex. Z ---> l+l-) l- Z - electrons: need to subtract QCD jet background - leptons from W/Z decay: no estimate above ~200 GeV - extrapolation by observing of trends, simulation Electrons: No decrease of selection efficiency expected at high energies (careful with isolation cut) Muons: highly energetic muons occasionally radiate so much bremsstrahlung, that their tracks might be too distorted for reconstruction
18 Invariant mass distributions Electrons Muons The MC and QCD estimate are normalized to data in the mass range 70-110 GeV (Normalization factor: 99% for both electron and muon channel.)
19 Systematic Uncertainties Uncertainties on yield at invariant mass of 1.5 TeV: - only mass-dependent uncertainties are considered - no theoretical uncertainties on signal (by convention) (except for 5% Z boson cross-section uncertainty, which replaces the luminosity uncertainty) - uncertainties below 2% negligible (Pileup, energy calibration, momentum/energy resolution, electron trigger, reconstruction, identification efficiency, QCD multijet estimate)
20 Signal templates Two templates for every tested signal mass: Z' (spin-1) - limits will be set on Z' (SSM) and motivated Z' (E6) - use shape of Sequential Standard Model Z' - neglect interference with DY - reweight flat sample (with Breit-Wigner, Parton Luminosity) G* (spin-2): 1 - limits will be set for different couplings (0.01-0.1) - use shape with largest width (0.1) - fully simulated (in 5fb-1 analysis: flat sample) 0.1 Phys.Rev. D63 (2001) 075004
21 Statistical Method Likelihoods - compare invariant mass distribution in data with SM background and signal templates (for invariant masses above 130 GeV) - binned Poisson likelihood (invariant mass bins k): - Systematic uncertainties are considered through nuisance parameters (for which we assume Gaussian probability functions) - Reduced likelihood: Integral over all nuisance parameters - Convert Nsig to cross-section
Search for a signal 22 - use 2D maximum likelihood fit to find most probable M Z' and σz' - get p-value by comparing to background-only pseudo-experiments ---> what is the percentage of pseudo experiments, which show an excess at least as significant as the one seen in data - electrons: p = 54%, muons: p = 24% (evidence: 0.1%, discovery 0.00003%) ---> No significant excess found ---> Setting limits
Statistical Method Setting limits 23 - Bayes theorem: - no prior knowledge on cross-section: - 95% confidence level (C.L.) limits can be set by finding upper edge of integral, such that: ----> with 95% C.L. we say, that σsigb is below (σsigb)95
24 Limits Spin 1 - Sequential Standard Model Z' as baseline model (same couplings as the Z boson) - different E6 models Limits on SSM Z' [TeV]: ee mumu combined Observed 1.70 1.61 1.83 Expected 1.70 1.61 1.83 Combined limits on E6 models [TeV]: Z' (ψ) Z' (N) Z' (η) Z' (I) Z' (S) Z' (χ) 1.49 1.60 1.52 1.54 1.56 1.64 Previous limits on SSM Z' [TeV]: Tevatron (1.071), ATLAS (1.042), CMS (1.140) LEP (indirect, 1.79) -1 CMS 1 fb [TeV]: 1.94
25 Limits Spin 2 - Randall-Sundrum Graviton (different couplings) Limits on RS graviton (k/mpl = 0.1): ee mumu combined Observed 1.51 1.45 1.63 Expected 1.50 1.44 1.63 Combined limits for additional couplings: Coupling 0.01 0.03 0.05 0.71 1.03 1.33 Previous limits on Gravitons (k/mpl = 0.1) [TeV]: CMS, CDF, D0 (all below 1.08) Combination with diphoton channel (2.2fb -1): 1.95 TeV (see next slide)
Limits G* - Combination with diphoton channel26 - diphoton branching fraction ~2 times dielectron/dimuon - limits on Gravitons (k/mpl = 0.1) [TeV]: 1.85 (diphoton), 1.95 (diphoton+dilepton)
Technicolor search (ATLAS-CONF-2011-125) 27 Alternative mechanism of electroweak symmetry breaking (no scalar Higgs) Low Scale Technicolor model (Lane, Eichten): techni-isospin good symmetry QCD-like technihadron spectrum: near-degenerate ρt, ωt narrow spin-1 resonance, same acceptance as Z' ---> same cross-section limits as for the Z' analysis Low Scale technicolor interpretation of CDF Wjj excess (Phys. Rev. Lett. 106 (2011) 171801)
Non-resonant search: Contact Interactions arxiv:1112.4462v2 28 (same event selection as for Z' search) - looking for broad excess at high dilepton masses - models: quark/lepton compositeness, exchange particles at masses inaccessible to LHC energies 2 4 - cross-section: XS ~ XS(DY) η FI/Λ + FC/Λ ---> prior chosen flat in 1/Λ2 and alternatively 1/Λ4 p-values in signal search > 5% ---> set limits η= +/- 1 (pos./neg. interference) Λ: scale
29-1 5 fb analysis (full 2011 dataset)
30 Status of the analysis Two parts: ---> conference result showing limits on Z', G* for Moriond Electroweak (approved by collaboration last week, awaiting final sign off) ---> more extensive paper in a couple of weeks Improvements: ---> better background estimates ---> 2 station muons (already in conference note) ---> limits on more models (will be in paper)
31 2-station muons Low acceptance in muon channel due to strict cuts on the hits in the muon spectrometer (inner, middle, outer layer) ----> decent momentum resolution ----> but three layers not available for every region of the detector Inclusion of 2-station muons - 1 leg 3 stations, 1 leg 2 stations - combined muons (inner detector and spectrometer) - require 5 hits in each of the inner and outer layer - veto regions with misalignments - significance cut (ID-MS): 3σ (instead of 5 like for 3-station muons) ---> acceptance increase by ~4% (absolute) ---> resolution ~25% (dominated by inner detector)
32 Models - setting limits on more models: Z', G*, arxiv:1004.2432v2 Z*, Technicolor (MWTC, LSTC), Torsion, Kaluza-Klein Z'/γ* - model independent limits: Minimal Z' models - template production from DY for many models, using matrix elements or generator level distributions - limits on couplings instead of cross-sections, including interference effects (for Kaluza-Klein, Minimal models)
33 5 fb-1 analysis: Minimal Z' models - class of models with same coupling structure, parametrized by y' and cos(θ) (arxiv:0909.1320v2) - known models in this parametrization: - make templates by reweighting DY illustration with matrix element ---> includes interference and width - set limit on coupling y'
34 Conclusion - search for high-mass dielectron and dimuon resonances with the ATLAS detector - no significant excess found in 1 fb-1 - lower mass limits set on Z' models, RS graviton, technimesons in Low Scale Technicolor - 5 fb-1 analysis with improvements is underway
35 BACKUP
Event yield tables 36
Tight muon selection 37
Tight-loose muon selection 38
Event Display (mee = 993 GeV ) 39
Event Display (mμμ = 959 GeV ) 40
Branching fractions graviton decay 41 arxiv:hep-ph/0211205v1
Minimal Z' Models (arxiv:0909.1320v2) 42 Class of models with following couplings: ul,dl ur dr νl,el νr er 1/6 gy+1/3 gbl 2/3 gy+1/3gbl -1/3 gy+1/3 gbl -1/2 gy- gbl -gbl -gy-gbl gy = γ' cos θ, gbl = γ' sin θ ----> two parameters (γ', θ) Known models covered by this parametrization: