muon monte carlo: a versatile tool for lepton propagation through matter dmitry chirkin, lbnl,...
TRANSCRIPT
![Page 1: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/1.jpg)
Muon Monte Carlo: a versatile tool for lepton propagation through matter
Dmitry Chirkin, LBNL, Berkeley, USA
October 31, 2006, Dortmund University
![Page 2: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/2.jpg)
Introduction
Muon propagation: why do we need it?
Muon/neutrino detectors?
![Page 3: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/3.jpg)
Particles observed by neutrino detectors
![Page 4: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/4.jpg)
Muon Monte Carlo
A tool for muon propagation simulation
![Page 5: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/5.jpg)
Structure of the program
![Page 6: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/6.jpg)
Simulation of muon propagation
Starting with Ei Ending with Ef
f(E)dx
P(E)dx
Continuous losses
Stochastic losses
1.0
vcut=0.05
Ecut=500 MeV
![Page 7: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/7.jpg)
Method of propagation
Distribution of the final energy of the muons that crossed 300 m of Frejus rock starting with 100 TeV
0.05
10-4
0.05
0.01
10-3
10-4
![Page 8: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/8.jpg)
Method of propagation
Distribution of the final energy of the muons that crossed 300 m of Frejus rock starting with 100 TeV
0.05
10-4
0.05
0.01
10-3
10-4
0.05
0.05
0.0110-3
10-4
![Page 9: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/9.jpg)
Muon cross sections
Ionization losses + knock-on electrons
Bremsstrahlung
Photonuclear
Electron pair production
Decay
10 TeV muon
![Page 10: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/10.jpg)
Bremsstrahlung
muons electrons
![Page 11: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/11.jpg)
Photonuclear interaction
Photon-nucleon Photonuclear
![Page 12: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/12.jpg)
Q2
1 GeV2soft hard
photoproduction DIS
GVDM CKMT ALLM
Bezrukov-Bugaev Butkevich-Mikheyev Abramowicz Levin Levy Maor
1991 19972002BB 1981
BB + Hard 03Bugaev Shlepin
ZEUS 94
Kokoulin 99
Nuclear effects
Dutta Smirnov
Muon propagator (MMC) settings: ph-nu settings
![Page 13: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/13.jpg)
Mass effects
Delta-correction to ionization (included into the ionization cross section)
LPM suppression of the bremsstrahlung and direct electron pair production
Dialectric suppression of the bremsstrahlung cross section
![Page 14: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/14.jpg)
Moliere scattering
![Page 15: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/15.jpg)
Electron, tau, and monopolemuon electron
tau monopole
![Page 16: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/16.jpg)
Neutrino propagation
Neutrino cross sections
Also: oscillations
Earth density profile is implemented
![Page 17: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/17.jpg)
Interpolation errors
Distribution of the final energy of the muons that crossed 300 m of Frejus rock starting with 100 TeV
Comparison: parameterized vs. non-parameterized
Interpolation precision: (epa-enp)/epa
Interpolation order: g=2,…6
vcut=0.01vcut=10-4vcut=10-4
Elow=10 TeV
![Page 18: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/18.jpg)
Algorithm errors: average propagation
Deviation from average energy loss (with vcut=1)Propagating 4 106 muons through 100 m of Frejus rock
![Page 19: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/19.jpg)
Algorithm errors: survival probability
106 muons with energy 9 TeV propagated through 10 km of water
![Page 20: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/20.jpg)
Comparison with other codes: MUM (MUons + Medium)
MUM code by E. Bugaev, I. Sokalski, S. Klimushin
![Page 21: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/21.jpg)
Spectra of the secondaries
MMC MUMLOHLIP
![Page 22: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/22.jpg)
Number and energy of secondaries
![Page 23: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/23.jpg)
Implementation for muon/neutrino detector
3 propagation regions:
• before the detector: propagation with fixed vcut
• inside the detector: propagation with fixed vcut or Ecut
• after the detector: fast propagation with vcut=1.0
![Page 24: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/24.jpg)
Parameterization of atmospheric lepton fluxes withCORSIKA
Primaries with Z=1,…,26: Poli-gonato composition model
Run CORSIKA
Parameterize simultated fluxes with
With corrections for zenith angle, muon energy loss and decay
![Page 25: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/25.jpg)
Parameterization of the atmosphere
![Page 26: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/26.jpg)
Muon energy losses
![Page 27: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/27.jpg)
Atmospheric lepton fluxes
muons muon neutrinos electron neutrinos
![Page 28: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/28.jpg)
Integrated fluxes
![Page 29: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/29.jpg)
Quality of the fits
fit quality stability of the result
![Page 30: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/30.jpg)
Things to remember• mmc was written in 2000 and has been updated a few times with new cross sections and features
• mmc has been used by AMANDA and now IceCube, also in data analysis of Frejus
• mmc is available at http://icecube.berkeley.edu/~dima/work/MUONPR
• mmc stands for Muon Monte Carlo and propagates muons
• perhaps more appropriate name is ALMC: All Lepton Monte Carlo, since it propagates muons, taus, electrons, all neutrinos
• mmc can also stand for monopole monte carlo
![Page 31: Muon Monte Carlo: a versatile tool for lepton propagation through matter Dmitry Chirkin, LBNL, Berkeley, USA October 31, 2006, Dortmund University](https://reader036.vdocument.in/reader036/viewer/2022062500/5697bfa91a28abf838c9a1be/html5/thumbnails/31.jpg)
Applet demonstration