Searching for Lepton Flavour Violation with the
Mu3e Experiment
Niklaus Berger
Institut für Kernphysik, Johannes-Gutenberg Universität Mainz PSI 2016
October 2016
Niklaus Berger – PSI, October 2016 – Slide 2
LFV Muon Decays: Experimental Situation
μ
+→ e
+γ μ
-N → e
-N μ
+→ e
+e
-e
+LFV Muon Decays: Experimental Situation
μ
+→ e
+γ μ
-N → e
-N μ
+→ e
+e
-e
+MEG (PSI) SINDRUM II (PSI) SINDRUM (PSI)
B(μ
+→ e
+γ) < 4.2 ∙ 10
-13(2016) B(μ
-Au → e
-Au) < 7 ∙ 10
-13(2006) B(μ
+→ e
+e
-e
+) < 1.0 ∙ 10
-12(1988)
Niklaus Berger – PSI, October 2016 – Slide 4
LFV Muon Decays: Experimental Situation
μ
+→ e
+γ μ
-N → e
-N μ
+→ e
+e
-e
+MEG (PSI) SINDRUM II (PSI) SINDRUM (PSI)
B(μ
+→ e
+γ) < 4.2 ∙ 10
-13(2016) B(μ
-Au → e
-Au) < 7 ∙ 10
-13(2006) B(μ
+→ e
+e
-e
+) < 1.0 ∙ 10
-12(1988)
upgrading Mu2e/Comet Mu3e
LFV Muon Decays: Experimental signatures
μ
+→ e
+γ μ
-N → e
-N μ
+→ e
+e
-e
+Kinematics
• 2-body decay
• Monoenergetic e+, γ
• Back-to-back
Kinematics
• Quasi 2-body decay
• Monoenergetic e-
• Single particle detected
Kinematics
• 3-body decay
• Invariant mass constraint
• Σ pi = 0
Niklaus Berger – PSI, October 2016 – Slide 6
LFV Muon Decays: Experimental signatures
μ
+→ e
+γ μ
-N → e
-N μ
+→ e
+e
-e
+Kinematics
• 2-body decay
• Monoenergetic e+, γ
• Back-to-back Background
• Accidental background
Kinematics
• Quasi 2-body decay
• Monoenergetic e-
• Single particle detected Background
• Decay in orbit
• Antiprotons, pions, cosmics
Kinematics
• 3-body decay
• Invariant mass constraint
• Σ pi = 0 Background
• Radiative decay
• Accidental background
LFV Muon Decays: Experimental signatures
μ
+→ e
+γ μ
-N → e
-N μ
+→ e
+e
-e
+Kinematics
• 2-body decay
• Monoenergetic e+, γ
• Back-to-back Background
• Accidental background
Kinematics
• Quasi 2-body decay
• Monoenergetic e-
• Single particle detected Background
• Decay in orbit
• Antiprotons, pions
Kinematics
• 3-body decay
• Invariant mass constraint
• Σ pi = 0 Background
• Radiative decay
• Accidental background
Con tinuous Be am
Con tinuous Be am Pul sed Be
am
Niklaus Berger – PSI, October 2016 – Slide 8
e +
e + e -
• μ+ → e+e-e+
• Two positrons, one electron
• From same vertex
• Same time
• Sum of 4-momenta corresponds to muon at rest
• Maximum momentum: ½ mμ = 53 MeV/c
The signal
• Combination of positrons from ordinary muon decay with electrons from:
- photon conversion, - Bhabha scattering, - Mis-reconstruction
• Need very good timing, vertex and momentum resolution
Accidental Background
e
+e
+e
-Niklaus Berger – PSI, October 2016 – Slide 10
• Allowed radiative decay with internal conversion:
μ
+→ e
+e
-e
+νν
• Only distinguishing feature:
Missing momentum carried by neutrinos
Internal conversion background
• Need excellent
momentum resolution
• New: NLO available from Matteo Fael and Signer et al. - now 10-20% easier
Branching Ratio
10-12
10-16 10-18 10-14
e+e-e+ mass (MeV/c2)
105 106 104
103 102
101
Internal conversion background
Signal
Building the
Mu3e Experiment
aiming for a branching ratio sensitivity of 10
-16Niklaus Berger – PSI, October 2016 – Slide 12
• Apply magnetic field (e.g. 1 Tesla)
• Measure curvature of particles in field
• Limited by detector resolution and scattering in detector
Momentum measurement
2 Billion Muon Decays/s
50 ns, 1 Tesla field
Niklaus Berger – PSI, October 2016 – Slide 14
High voltage monolithic active pixel sensors - Ivan Perić
• Use a high voltage commercial process (automotive industry)
Fast and thin sensors: HV-MAPS
P-substrate
N-well E field
High voltage monolithic active pixel sensors - Ivan Perić
• Use a high voltage commercial process (automotive industry)
• Small active region, fast charge collection via drift
Fast and thin sensors: HV-MAPS
P-substrate N-well
Particle
E field
Niklaus Berger – PSI, October 2016 – Slide 16
High voltage monolithic active pixel sensors - Ivan Perić
• Use a high voltage commercial process (automotive industry)
• Small active region, fast charge collection via drift
Fast and thin sensors: HV-MAPS
P-substrate N-well
Particle E field
• Implement logic directly in N-well in the pixel - smart diode array
• Can be thinned down to < 50 μm
(I.Perić, P. Fischer et al., NIM A 582 (2007) 876 )
Submitting large (2×1 cm2) prototype in next weeks
Performance
row-axis [mm]
0 0.5 1 1.5 2 2.5 3
column-axis [mm]
0 0.5 1 1.5 2 2.5 3
efficiency_pixeluv
Entries 900390
Mean x 1.557
Mean y 1.803
RMS x 0.922
RMS y 0.8324
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
efficiency_pixeluv
Entries 900390
Mean x 1.557
Mean y 1.803
RMS x 0.922
RMS y 0.8324
Mupix7, 735 mV threshold, HV = -85 V
Threshold [V]
0.68 0.7 0.72 0.74 0.76 0.78
Efficiency
0.9 0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1
Efficiency Noise
99 %
Noiserate per pixel [1/s]
1 10 102
103
Niklaus Berger – PSI, October 2016 – Slide 18
• 50 μm silicon
• 25 μm Kapton™ flexprint with aluminium traces
• 25 μm Kapton™ frame as support
• About 1‰ of a radiation length per layer
Mechanics
Niklaus Berger – PSI, October 2016 – Slide 20
• Add no material:
Cool with gaseous Helium (low scattering, high mobility)
• ~ 250 mW/cm2 - total ~3 kW
• Simulations: Need ~ several m/s flow
Cooling
• Full scale heatable prototype built
• 36 cm active length
• Vibrations studied using Michelson-Interferometer
• Can keep temperature below 70°C
• 1 T magnetic field
• Resolution dominated by multiple scattering
• Momentum resolution to first order:
Σ
P/P ~ θ
MS/Ω
• Precision requires large lever arm (large bending angle Ω) and low multiple scattering θMS
Momentum measurement
Ω MS
θ
MSB
Niklaus Berger – PSI, October 2016 – Slide 22
Precision vs. Acceptance
50 MeV/c 25 MeV/c 12 MeV/c B→
33 cm
Precision vs. Acceptance
50 MeV/c 25 MeV/c 12 MeV/c B→
Niklaus Berger – PSI, October 2016 – Slide 24
Precision vs. Acceptance
50 MeV/c 25 MeV/c 12 MeV/c B→
Precision vs. Acceptance
50 MeV/c 25 MeV/c 12 MeV/c B→
Niklaus Berger – PSI, October 2016 – Slide 26
Precision vs. Acceptance
50 MeV/c 25 MeV/c 12 MeV/c B→
Ω ~ π MS
θMS
B
Detector Design
muon beam
target
Niklaus Berger – PSI, October 2016 – Slide 28
Detector Design
muon beam
target
Detector Design
muon beam
target
inner pixel layers
Niklaus Berger – PSI, October 2016 – Slide 30
Detector Design
outer pixel layers
muon beam
target
inner pixel layers
Detector Design
outer pixel layers
muon beam
target inner pixel layers recurl pixel
layers
scintillating fibres
Niklaus Berger – PSI, October 2016 – Slide 32
2
] [MeV/c m
rec96 98 100 102 104 106 108 110
2
Events per 0.2 MeV/c
−3
10
−2
10
−1
10 1 10 10
2at 10
-12→ eee µ
at 10
-13→ eee µ
at 10
-14→ eee µ
at 10
-15→ eee µ
ν ν
→ eee µ
muons/s muon stops at 10
810
15Performance Simulations: Mass reconstruction
Work in progress
Poster by Alexandr Kozlinskiy on track reconstruction
Need suppression of accidental background:
Timing
Niklaus Berger – PSI, October 2016 – Slide 34
Detector Design
scintillating fibres
outer pixel layers
muon beam
target
inner pixel layers
Detector Design
outer pixel layers
muon beam
target inner pixel layers recurl pixel
layers
recurl pixel layers
scintillating fibres
Scintillating tiles
Niklaus Berger – PSI, October 2016 – Slide 36
• 3 layers of 250 μm scintillating fibres
• Read-out by silicon photomultipliers (SiPMs) and custom ASIC (MuTRiG)
• Timing resolution O(0.5 - 1 ns)
(See posters by Giada Rutar, Angela Papa)
Timing Detector: Scintillating Fibres
Single photon Efficiency > 98%
(≥ 2 photons)
Timing Detector: Scintillating tiles
• ~ 0.5 cm3 scintillating tiles
• Read-out by silicon photomultipliers (SiPMs) and custom ASIC (STiC)
Scin ator Tiles
SiPM Readout
Electronics
Time Difference [ps]
-7500 -500 -250 0 250 500 750
2000 4000 6000 8000 10000
σ = 79.2 ps
• Test beam with tiles, SiPMs and readout ASIC
• Timing resolution ~ 80 ps
Niklaus Berger – PSI, October 2016 – Slide 38
Sensitivity
Phase IA: Starting 2018 2∙107 μ/s
Target Inner pixel layers
Outer pixel layers μ Beam
Target Inner pixel layers
Scintillating fibres
Outer pixel layers Recurl pixel layers
Scintillator tiles
μ Beam
Sensitivity
Phase IB: 2019+
1∙108 μ/s
Niklaus Berger – PSI, October 2016 – Slide 40
Sensitivity
Phase II: 2021+
New Beam Line 2∙109 μ/s
Target Inner pixel layers
Scintillating fibres
Outer pixel layers Recurl pixel layers
Scintillator tiles
μ Beam
• Mu3e aims for μ → eee at the 10-16 level
• First large scale use of HV-MAPS
• Build detector layers thinner than a hair
• Timing at the 100 ps level
• Reconstruct 2 billion tracks/s in 1 Tbit/s on ~50 GPUs
• Start data taking in 2018
• 2 billion muons/s not before 2021
Conclusion
1940 1960 1980 2000 2020
Year
90%–CL bound
10–14 10–12 10–10 10–8 10–6 10–4 10–2 100
μ eγ
μ 3e
μN eN
τ μγ
τ 3μ
10–16
SINDRUM SINDRUM II MEG
MEG plan Mu3e Phase I
Mu3e Phase II
Niklaus Berger – PSI, October 2016 – Slide 42
Backup Material
New physics in μ
+→ e
+e
-e
+Tree diagrams
• Higgs triplet model
• Extra heavy vector bosons (Z’)
• Extra dimensions (Kaluza-Klein tower) Loop diagrams
• Supersymmetry
• Little Higgs models
• Seesaw models
• GUT models (leptoquarks)
• and much more...
Niklaus Berger – PSI, October 2016 – Slide 44
New physics in μ
+→ e
+e
-e
+Muon decays at the 10
-16level sensitive to new physics
at O (1000 TeV) scale for O (1) couplings!
MUPIX electronics
Niklaus Berger – PSI, October 2016 – Slide 46
A general effective Lagrangian
Tensor terms (dipole)
L
μ → eee= 2 G
F( m
μA
Rμ
Rσ
μνe
LF
μν+ m
μA
Lμ
Lσ
μνe
RF
μν+ g
1(μ
Re
L) (e
Re
L) + g
2(μ
Le
R) (e
Le
R)
+ g
3(μ
Rγ
μe
R) (e
Rγ
μe
R) + g
4(μ
Lγ
μe
L) (e
Lγ
μe
L)
+ g
5(μ
Rγ
μe
R) (e
Lγ
μe
L) + g
6(μ
Lγ
μe
L) (e
Rγ
μe
R) + H. C. )
e.g. supersymmetry
Four-fermion terms scalar
vector
e.g. Z’
(Y. Kuno, Y. Okada,
Rev.Mod.Phys. 73 (2001) 151)
Comparison with μ
+→ e
+γ
L
LFV= A m
μ Rμ
Rσ
μνe
LF
μν+ (μ
Lγ
μe
L) (e
Lγ
μe
L) (κ+1)Λ
2κ (κ+1)Λ
2• One loop term and one contact term
• Ratio κ between them
• Common mass scale Λ
• Allows for sensitivity comparisons between μ → eee and μ → eγ
• In case of dominating dipole couplings (κ = 0):
B(μ → eee) = 0.006 (essentially αem) B(μ → eγ)
Niklaus Berger – PSI, October 2016 – Slide 48
Detector Design
outer pixel layers
muon beam
target inner pixel layers recurl pixel
layers
recurl pixel layers
scintillating fibres
Scintillating tiles
Muons from PSI
DC muon beams at PSI:
• πE5 beamline: ~ 108 muons/s
(MEG experiment, Mu3e phase I)
• Surface muons, p = 29.7 MeV/c Stopped in < 1 mm of plastic
Niklaus Berger – PSI, October 2016 – Slide 50
Muons from PSI
DC muon beams at PSI:
• πE5 beamline: ~ 108 muons/s
(MEG experiment, Mu3e phase I)
• Surface muons, p = 29.7 MeV/c Stopped in < 1 mm of plastic
• The μ → eee experiment (final stage) requires 2 × 109 muons/s focused and collimated on a ~2 cm spot
Muons from PSI
DC muon beams at PSI:
• πE5 beamline: ~ 108 muons/s
(MEG experiment, Mu3e phase I)
• Surface muons, p = 29.7 MeV/c Stopped in < 1 mm of plastic
• The μ → eee experiment (final stage) requires 2 × 109 muons/s focused and collimated on a ~2 cm spot
• More than ~ 1011 muons/s are produced;
bring magnetic elements closer to cap- ture them:
High intensity muon beamline (HiMB) study currently ongoing
Niklaus Berger – PSI, October 2016 – Slide 52
HV-MAPS
3 m m
HV-MAPS
3 m m
Pixels with amplifier
40 x 32 pixels
80 x 103 μm pixel size
Niklaus Berger – PSI, October 2016 – Slide 54
HV-MAPS
3 m m
Pixels with amplifier
40 x 32 pixels
80 x 103 μm pixel size
Comparator and digital pixel logic
Tests done at
• CERN 250 GeV pions
• DESY 5 GeV electrons
• PSI 250 MeV pions
• Mainz 1.5 GeV electrons
• Thanks for all the beam time and support!
Beam tests
Niklaus Berger – PSI, October 2016 – Slide 56
Introduction
Y
• X
Introduction
Y
• X
Niklaus Berger – PSI, October 2016 – Slide 58
Introduction
Y
• X
Introduction
Y
• X
Niklaus Berger – PSI, October 2016 – Slide 62
Timing Detector: Scintillating tiles
• Test beam with tiles, SiPMs and readout ASIC
• Timing resolution ~ 80 ps
Time Difference [ps]
-7500 -500 -250 0 250 500 750
2000 4000 6000 8000 10000
σ = 79.2 ps
Front
Back
3.5 cm
Data Acquisition
Niklaus Berger – PSI, October 2016 – Slide 64
• 280 Million pixels (+ fibres and tiles)
• No trigger
• ~ 1 Tbit/s
• FPGA-based switching network
• O(50) PCs with GPUs
Data Acquisition
2928 Pixel Sensors
up to 36 1.25 Gbit/s links
FPGA FPGA FPGA
...
86 FPGAs
1 6 Gbit/s link each
GPU
PC GPU
PC
GPU 12 PCs PC
12 10 Gbit/s links per
8 Inputs each
~ 3072 Fibre Readout Channels
FPGA FPGA
...
48 FPGAs
~ 3500 Tiles
FPGA FPGA
...
48 FPGAs
Data Collection
Server
Mass Storage Gbit Ethernet
Switching
Board Switching
Board Switching
Board
Front-end(inside magnet)
Switching Board
• 280 Million pixels (+ fibres and tiles)
• No trigger
• ~ 1 Tbit/s
• Need to find and fit billions of tracks/s
Online reconstruction
Niklaus Berger – PSI, October 2016 – Slide 66
• PCs with Graphics Processing Units (GPUs)
• Online track and event reconstruction
• 109 3D track fits/s achieved
• Data reduction by factor ~1000
• Data to tape < 100 Mbyte/s
Online filter farm
Introduction
Y
• X
Niklaus Berger – PSI, October 2016 – Slide 68
Cooling tests
Global helium stream
Local helium stream
• Limited by detector resolution and scattering in detector