Use of event-level IceCube data in SUSY scans

Similar documents
The Sun as a Particle Physics Laboratory

Beyond the Standard Model or Bust

Global SUSY Fits with IceCube

Beyond the Standard Model or Bust

Statistical challenges in astroparticle physics

pmssm Dark Matter Searches On Ice! Randy Cotta (Stanford/SLAC) In collaboration with: K.T.K. Howe (Stanford) J.L. Hewett (SLAC) T.G.

SUSY scans and dark matter

Search for Supersymmetry at CMS in all-hadronic final states

- A Bayesian approach

Searching for Supersymmetry at the LHC David Stuart, University of California, Santa Barbara. CMS SUSY Search, D. Stuart, June 2011, Lisbon!

Prospects and Blind Spots for Neutralino Dark Matter

Collider Searches for Dark Matter

Beyond the Standard Model global fits: then, now and tomorrow

Where is SUSY? Institut für Experimentelle Kernphysik

SUSY Phenomenology & Experimental searches

Neutrinos and DM (Galactic)

Earth WIMP search with IceCube. Jan Kunnen for the IceCube Collaboration

Searches for Dark Matter Annihilations in the Sun and Earth with IceCube and DeepCore. Matthias Danninger for the IceCube collaboration

ATLAS Missing Energy Signatures and DM Effective Field Theories

Radiative natural SUSY with mixed axion-higgsino CDM

Dark matter, supersymmetry and global fits

IceCube 79 Solar WIMP Search. Matthias Danninger

An update on scalar singlet dark matter

Indirect constraints on the (C)MSSM parameter space

WIMP diffusion in the Solar System and the neutrino signal from the Sun and the Earth

PHY326/426 Dark Matter and the Universe. Dr. Vitaly Kudryavtsev F9b, Tel.:

Carsten Rott. mps. ohio-state. edu. (for the IceCube Collaboration)

BSM Higgs Searches at ATLAS

Slepton, Charginos and Neutralinos at the LHC

Dark Matter Experiments and Searches

IceCube & DeepCore Overview and Dark Matter Searches. Matthias Danninger for the IceCube collaboration

The search for particle dark matter

Razor Phenomenology ALL. Life is not just cut and count! By the Razor Guys

Search for R-parity violating Supersymmetry. III Phys. Inst. A, RWTH Aachen

Indirect Dark Matter Detection

Search for New Physics in the Multijet and Missing Transverse Momentum Final State at CMS

Exotica in CMS. ICNFP 2015 Kolymbari, Crete 25 August Claudia-Elisabeth Wulz CMS Collaboration Institute of High Energy Physics, Vienna

Search for SUSY at CMS in lepton or photon final states. PANIC 2011 K. Theofilatos (on behalf of CMS)

Indirect detection of Dark Matter with the ANTARES Neutrino Telescope

On behalf of the ATLAS and CMS Collaborations

ICHEP 2008 Philadelphia

Sneutrino dark matter and its LHC phenomenology

Preparations for SUSY searches at the LHC

Search for New Physics through Rare B Decays

Dynamics of solar system bound WIMPs

Natural SUSY and the LHC

Probing the early Universe and inflation with indirect detection

Dark Matter WIMP and SuperWIMP

Neutralino Dark Matter as the Source of the WMAP Haze

Ingredients to this analysis

Dark Matter Direct Detection in the NMSSM

Contributions by M. Peskin, E. Baltz, B. Sadoulet, T. Wizansky

Searches for New Physics in quarkonium decays at BaBar/Belle

(Neutrino indirect) detection of neutralino dark matter in (non) universals SUSY GUT models

Gustav Wikström. for the IceCube collaboration

Sho IWAMOTO. 15 Sep Osaka University. Based on [ ] in collaboration with M. Abdullah, J. L. Feng, and B. Lillard (UC Irvine)

A profile likelihood analysis of the CMSSM with Genetic Algorithms

Search for long-lived particles at CMS. TeVPA Brian Francis for the CMS Collaboration

Testing a DM explanation of the positron excess with the Inverse Compton scattering

LHC searches for dark matter.! Uli Haisch

Positron Fraction from Dark Matter Annihilation in the CMSSM

DarkSUSY. Joakim Edsjö With Torsten Bringmann, Paolo Gondolo, Lars Bergström, Piero Ullio and Gintaras Duda. APS Meeting

Supersymmetry at the LHC

Searches for low-mass Higgs and dark bosons at BaBar

GALACTIC CENTER GEV GAMMA- RAY EXCESS FROM DARK MATTER WITH GAUGED LEPTON NUMBERS. Jongkuk Kim (SKKU) Based on Physics Letters B.

Lecture 14. Dark Matter. Part IV Indirect Detection Methods

Dark Matter Implications for SUSY

Overall CMS SUSY search strategy

Early SUSY searches at the LHC

Highlights of top-quark measurements in hadronic nal states at ATLAS

Background and sensitivity predictions for XENON1T

Preparations for SUSY searches at the LHC

Search for neutralino dark matter with the AMANDA neutrino telescope

The search for missing energy events at the LHC and implications for dark matter search (ATLAS and CMS)

Sho IWAMOTO. 7 Nov HEP phenomenology joint Cavendish DAMTP U. Cambridge

Accurate Neutralino Relic Density. Paolo Gondolo

Supersymmetry Without Prejudice at the LHC

Stau Pair Production At The ILC Tohoku University Senior Presentation Tran Vuong Tung

Supersymmetry at the LHC

Emerging the 7 TeV LHC & the LC

A Bayesian approach to supersymmetry phenomenology

Combined Search for Neutrinos from Dark Matter Annihilation in the Galactic Center using IceCube and ANTARES

Searches for electroweak production of supersymmetric gauginos and sleptons with the ATLAS detector

Neutralino Dark Matter and the 125 GeV Higgs boson measured at the LHC

SUSY at Accelerators (other than the LHC)

PoS(LHCPP2013)016. Dark Matter searches at LHC. R. Bellan Università degli Studi di Torino and INFN

Signals from Dark Matter Indirect Detection

Dark Matter in the Universe

SUSY AND COSMOLOGY. Jonathan Feng UC Irvine. SLAC Summer Institute 5-6 August 2003

Supersymmetry at the ILC

Yukawa and Gauge-Yukawa Unification

Inclusive searches in ATLAS: How can we discover SUSY in 2009

Probing SUSY Dark Matter at the LHC

SUSY at Accelerators (other than the LHC)

FERMION PORTAL DARK MATTER

Model Fitting in Particle Physics

Constraining Galactic dark matter in the Fermi-LAT sky with photon counts statistics

LHC Impact on DM searches

Searches for SUSY in events with di-leptons at CMS

Dark Matter on the Smallest Scales Annika Peter, 7/20/09

Transcription:

Use of event-level IceCube data in SUSY scans Pat Scott Department of Physics, McGill University July 3, 202 Based on: PS, Chris Savage, Joakim Edsjö and The IceCube Collaboration (esp. Matthias Danninger and Klas Hultqvist), arxiv:207.soon Slides available from www.physics.mcgill.ca/ patscott

Outline How to find DM with neutrino telescopes How to find DM with neutrino telescopes 2

Introduction The short version:

Introduction The short version: Halo WIMPs crash into the Sun

Introduction The short version: Halo WIMPs crash into the Sun 2 Some lose enough energy in the scatter to be gravitationally bound

Introduction The short version: Halo WIMPs crash into the Sun 2 Some lose enough energy in the scatter to be gravitationally bound 3 Scatter some more, sink to the core

Introduction The short version: Halo WIMPs crash into the Sun 2 Some lose enough energy in the scatter to be gravitationally bound 3 Scatter some more, sink to the core 4 Annihilate with each other, producing neutrinos

Introduction The short version: Halo WIMPs crash into the Sun 2 Some lose enough energy in the scatter to be gravitationally bound 3 Scatter some more, sink to the core 4 Annihilate with each other, producing neutrinos 5 Propagate+oscillate their way to the Earth, convert into muons in ice/water

Introduction The short version: Halo WIMPs crash into the Sun 2 Some lose enough energy in the scatter to be gravitationally bound 3 Scatter some more, sink to the core 4 Annihilate with each other, producing neutrinos 5 Propagate+oscillate their way to the Earth, convert into muons in ice/water 6 Look for Čerenkov radiation from the muons in IceCube, ANTARES, etc

What can the muon signal tell me? Roughly: Number how much annihilation is going on in the Sun = info on σ SD, σ SI and σv Spectrum sensitive to WIMP mass m χ and branching fractions BF into different annihilation channels X Direction how likely it is that they come from the Sun In model-independent analyses a lot of this information is either discarded or not given with final limits Goal: Use as much of this information on σ SD, σ SI, σv, m χ and BF (X) as possible to directly constrain specific points and regions in WIMP model parameter spaces (+LHC+DD+... )

What can the muon signal tell me? The focus here is supersymmetry (SUSY) but this talk is on a framework, applicable to any model. Goal: Use as much of this information on σ SD, σ SI, σv, m χ and BF (X) as possible to directly constrain specific points and regions in WIMP model parameter spaces (+LHC+DD+... )

What can the muon signal tell me? The focus here is supersymmetry (SUSY) but this talk is on a framework, applicable to any model. All the methods discussed here are available in DarkSUSY 5.: www.darksusy.org Goal: Use as much of this information on σ SD, σ SI, σv, m χ and BF (X) as possible to directly constrain specific points and regions in WIMP model parameter spaces (+LHC+DD+... )

What can the muon signal tell me? The focus here is supersymmetry (SUSY) but this talk is on a framework, applicable to any model. All the methods discussed here are available in DarkSUSY 5.: www.darksusy.org All IceCube data used are available at http://icecube.wisc.edu/science/data/ic22-solar-wimp (and in DarkSUSY, for convenience) Goal: Use as much of this information on σ SD, σ SI, σv, m χ and BF (X) as possible to directly constrain specific points and regions in WIMP model parameter spaces (+LHC+DD+... )

Outline How to find DM with neutrino telescopes How to find DM with neutrino telescopes 2

SUSY Scanning with IceCube Simple Likelihood Simplest way to do anything is to make it a counting problem... Compare observed number of events n and predicted number θ for each model, taking into account error σ ɛ on acceptance: L num(n θ BG + θ sig ) = 2πσɛ 0 (θ BG + ɛθ sig ) n e (θ BG+ɛθ sig ) n! [ ɛ exp 2 ( ) ] ln ɛ 2 dɛ. σ ɛ () Nuisance parameter ɛ takes into account systematic errors on effective area, from theory, etc. σ ɛ 20% for IceCube.

SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal θbg 0 E disp(n i E i ) dpsignal de i (E i, Ξ) de i (3) dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i

SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts L = L num(n ) with Number of lit channels (energy estimator) L spec,i (N i, Ξ) = θbg dp BG (N i )+ θsignal E disp(n i E dn i and L ang,i ( cos φ i ) = θbg n L spec,i L ang,i (2) i= 0 i ) dpsignal de i (E i, Ξ) de i (3) dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i

SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts L = L num(n ) with Number of lit channels (energy estimator) L spec,i (N i, Ξ) = θbg dp BG (N i )+ θsignal E disp(n i E dn i and L ang,i ( cos φ i ) = θbg n L spec,i L ang,i (2) i= 0 i ) dpsignal de SUSY parameters i (E i, Ξ) de i (3) dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i

SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal θbg Predicted signal spectrum (from theory) 0 E disp(n i E i ) dpsignal de i (E i, Ξ) de i (3) dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i

SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal θbg Predicted signal spectrum (from theory) E disp(n i E (E de i i, Ξ) de i (3) Instrument response function 0 i ) dpsignal dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i

SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal Observed BG distribution θbg Predicted signal spectrum (from theory) 0 E disp(n i E (E de i i, Ξ) de i (3) Instrument response function i ) dpsignal dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i

SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal θbg Event arrival angle 0 E disp(n i E i ) dpsignal de i (E i, Ξ) de i (3) dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i

SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal θbg 0 E disp(n i E i ) dpsignal de i (E i, Ξ) de i (3) dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i Predicted signal direction (δ function at Sun)

SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal θbg 0 E disp(n i E (E de i i, Ξ) de i (3) Instrument response function i ) dpsignal dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i Predicted signal direction (δ function at Sun)

SUSY Scanning with IceCube Full Likelihood Full unbinned likelihood with number (L num ), spectral (L spec ) and angular (L ang ) parts with L spec,i (N i, Ξ) = and L ang,i ( cos φ i ) = L = L num(n ) n L spec,i L ang,i (2) i= θbg dp BG dn i (N i )+ θsignal Observed BG distribution θbg 0 E disp(n i E (E de i i, Ξ) de i (3) Instrument response function i ) dpsignal dp BG ( cos φ i ) + θsignal PSF( cos φ i ) (4) d cos φ i Predicted signal direction (δ function at Sun)

SUSY Scanning with IceCube Statistics 0 Why simple IN/OUT analyses are not enough:... Only partial goodness of fit, no measure of convergence, no idea how to generalise to regions or whole space. Frequency/density of models in IN/OUT scans means essentially nothing. More information comes from a global statistical fit. = parameter estimation exercise Composite likelihood made up of observations from all over: dark matter relic density from WMAP precision electroweak tests at LEP LEP limits on sparticle masses B-factory data (rare decays, b sγ) muon anomalous magnetic moment LHC searches, direct detection (only roughly implemented for now)

Example: SUSY Scanning with IceCube Global Fits CMSSM, IceCube-22 events m 0 m /2 and m χ 0 nuclear scattering cross-sections Scott, Danninger, Savage, Edsjö & Hultqvist 202 42 Scott, Danninger, Savage, Edsjö & Hultqvist 202 38 Scott, Danninger, Savage, Edsjö & Hultqvist 202 IC22 global fit IC22 global fit m0 (TeV) 3 2 Mean IC22 global fit log 0 ( σsi,p/cm 2 ) 43 44 45 46 log 0 ( σsd,p/cm 2 ) 39 40 4 42 43 Mean Mean 0.0 0.5.5 2.0 2.5 m (TeV) 2 0.0 m χ 0 (TeV) 0.0 m χ 0 (TeV) Contours indicate σ and 2σ credible regions Grey contours correspond to fit without IceCube data Shading+contours indicate relative probability only, not overall goodness of fit

Example: SUSY Scanning with IceCube Global Fits Base Observables Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) 42 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) 38 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) Old data Old data m0 (TeV) 3 2 Mean Old data log 0 ( σsi,p/cm 2 ) 43 44 45 46 log 0 ( σsd,p/cm 2 ) 39 40 4 42 43 Mean Mean 0.0 0.5.5 2.0 2.5 m (TeV) 2 0.0 m χ 0 (TeV) 0.0 m χ 0 (TeV)

Example: SUSY Scanning with IceCube Global Fits Base Observables + XENON-00 Grey contours correspond to Base Observables only Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) 42 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) 38 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) XENON00 XENON00 m0 (TeV) 3 2 Mean XENON00 log 0 ( σsi,p/cm 2 ) 43 44 45 46 log 0 ( σsd,p/cm 2 ) 39 40 4 42 43 Mean Mean 0.0 0.5.5 2.0 2.5 m (TeV) 2 0.0 m χ 0 (TeV) 0.0 m χ 0 (TeV)

Example: SUSY Scanning with IceCube Global Fits Base Observables + XENON-00 + CMS 5 fb Grey contours correspond to Base Observables only m0 (TeV) 3 2 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) Mean XENON00 +CMS 5 fb log 0 ( σsi,p/cm 2 ) 42 43 44 45 46 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) XENON00 +CMS 5 fb Mean log 0 ( σsd,p/cm 2 ) 38 39 40 4 42 43 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) XENON00 +CMS 5 fb Mean 0.0 0.5.5 2.0 2.5 m (TeV) 2 0.0 m χ 0 (TeV) 0.0 m χ 0 (TeV)

Example: SUSY Scanning with IceCube Global Fits Base Observables + XENON-00 + CMS 5 fb + IC22 00 Grey contours correspond to Base Observables only m0 (TeV) 3 2 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) Mean IC22 00 +XENON00 +CMS 5 fb log 0 ( σsi,p/cm 2 ) 42 43 44 45 46 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) IC22 00 +XENON00 +CMS 5 fb Mean log 0 ( σsd,p/cm 2 ) 38 39 40 4 42 43 Scott, Danninger, Savage, Edsjö, Hultqvist & The IceCube Collab. (202) Mean IC22 00 +XENON00 +CMS 5 fb 0.0 0.5.5 2.0 2.5 m (TeV) 2 0.0 m χ 0 (TeV) 0.0 m χ 0 (TeV) CMSSM, IceCube-22 with 00 boosted effective area (kinda like IceCube-86+DeepCore)

Example: Model Recovery m0 (TeV) 3 2 Scott, Danninger, Savage, Edsjö & Hultqvist 202 Mean IC22 00 Predicted signal events 40 20 00 80 60 40 20 IC22 00 Scott, Danninger, Savage, Edsjö & Hultqvist 202 Mean True value CMSSM, IceCube-22 00 signal reconstruction 60 signal events, 500 GeV, χχ W + W 0.0 0.5.5 2.0 m (TeV) 2 0 200 400 600 800 m χ 0 (GeV) 42 Scott, Danninger, Savage, Edsjö & Hultqvist 202 38 Scott, Danninger, Savage, Edsjö & Hultqvist 202 IC22 00 IC22 00 log 0 ( σsi,p/cm 2 ) 43 44 45 Mean log 0 ( σsd,p/cm 2 ) 39 40 4 42 43 Mean Grey contours correspond to reconstruction without energy information 0 200 400 600 800 m χ 0 (GeV) 0 200 400 600 800 m χ 0 (GeV)

Closing remarks A framework for directly comparing event-level neutrino telescope data to individual points in theory parameter spaces is in place These iclike extensions are available in DarkSUSY 5. IceCube event data has been released in a form digestible by DarkSUSY 5. Direct SUSY analyses of IC79 data are on the way Many models exist that only IC86 will be sensitive to The extensions can be used equally well for non-susy BSM scenarios too

Outline Backup Slides 3 Backup Slides

Advertisement Backup Slides If you liked any of the plots in this talk... Flat priors pippi v Best fit Mean True value 2.5 3.0 ) 3.5 log 0 (m /TeV 2 log 0 (m0/tev) 3.5 3.0 2.5 2.0 Best fit Mean True value pippi v Flat priors 2.5 3.0 ) 3.5 log 0 (m /TeV 2 A0 (TeV) 4 2 0 2 Flat priors Prof. likelihood Best fit Mean pippi v 0 20 30 40 50 60 tan β Profile likelihood ratio Λ = L/Lmax Pippi parse it, plot it arxiv:206.2245 http://github.com/patscott/pippi Generic pdflatex sample parser, post-processor & plotter

CMS 5 fb analyses Backup Slides [GeV] m /2 800 700 600 500 400 = LSP τ CMS Preliminary m(q ~ ) = 500 m(q ~ ) = 000 SS Dilepton Razor L int m(q ~ ) = 2000 - = 4.98 fb, m(g ~ ) = 500 MT2 m( ~ q) = 2500 s = 7 TeV tan(β)=0 A 0 = 0 GeV µ > 0 m t = 73.2 GeV ~ ± LEP2 l ± LEP2 χ m(g ~ ) = 000 300 200 00 m(g ~ ) = 500 OS Dilepton Multi-Lepton Lepton Jets+MHT Non-Convergent RGE's 500 000 500 2000 2500 3000 [GeV] m 0 No EWSB

Backup Slides XENON-00 00-day analysis WIMP-Nucleon Cross Section [cm 2 ] 0 0 0 0 0 0-39 -40-4 -42-43 -44 DAMA/Na CoGeNT DAMA/I CDMS (20) CDMS (200) XENON0 (S2 only, 20) EDELWEISS (20) XENON00 (20) observed limit (90% CL) Expected limit of this run: ± σ expected ± 2 σ expected XENON00 (200) Trotta et al. 0-45 Buchmueller et al. 6 7 8 90 20 30 40 50 00 200 300 400 000 2 WIMP Mass [GeV/c ]