Summer Research Program 2022
August 26th
Sebastián Ordoñez
Universidad Nacional de Colombia, COL
Supervisors: Agnieszka Leszczyńska and Alan Coleman



Method to improve the speedup of the IceTop simulations due to snow accumulation on tanks
Introduction
- IceTop simulation challenge: Different snow coverage of tanks for each year.
- My project: Test a method through which might be possible to improve the general speedup of the simulations.

Snow heights on top of IT tanks measured in Feb. 2012.
-
Snow height on top of the tanks increases on average by 20 cm every year.
- Significant impact on the measured signal and systematic uncertainties.
Introduction

Snow accumulation over the IceTop detector, measured in three different years: 2010, 2011, and
2012.
Introduction

Zenith angle (\(\theta\))
\(\theta\)
Snow height
Injected particles
( \(\mu^{\pm}, e^{\pm}, \gamma, \) etc )
Tank height: 1.3m
Tank radius: 0.93m
Zenith angle (\(\theta\))
\(\rho = 0.38 g/cm^{3}\)
Introduction
Introduction


- Top view of the tank from g4-tankresponse
- Lateral view of the tank from g4-tankresponse
Simulation
Simulation
-
This study uses the g4-tankresponse project:
- Calculates the response of an IceTop tank to incoming particles using Geant4.
- g4-tankresponse/resources/examples/simple_injection.py
- GCD file: GeoCalibDetectorStatus_2012.56063_V1_OctSnow.i3.gz
- Frame: MCTopPESeriesMap


Simulation
Simulation

Simulation
Simulation

Simulation
Visualization
Results
Phase space: preliminary results
- Particles: MuMinus (\(\mu^{-}\)), MuPlus (\(\mu^{+}\)), EMinus (\(e^{-}\)),EPlus (\(e^{-}\)), Gamma (\(\gamma\))
- Zenith angles: 0°,5°, ... , 75°
- Energies: [0.001, 1000] GeV > spaced evenly on a log scale
- Snow heights: [0.0, 0.5, ... , 10.0 ] m
Results: Muons
\(\mu^{-}\)

\(\mu^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 0°

Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 10°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 20°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 30°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 40°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 50°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 60°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 70°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 1 GeV
- Zenith angle: 0°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 1 GeV
- Zenith angle: 10°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 1 GeV
- Zenith angle: 20°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 1 GeV
- Zenith angle: 30°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 1 GeV
- Zenith angle: 40°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 1 GeV
- Zenith angle: 50°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 1 GeV
- Zenith angle: 60°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 1 GeV
- Zenith angle: 70°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 10 GeV
- Zenith angle: 0°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 10 GeV
- Zenith angle: 10°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 10 GeV
- Zenith angle: 20°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 10 GeV
- Zenith angle: 30°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 10 GeV
- Zenith angle: 40°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 10 GeV
- Zenith angle: 50°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 10 GeV
- Zenith angle: 60°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 10 GeV
- Zenith angle: 70°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 100 GeV
- Zenith angle: 0°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 100 GeV
- Zenith angle: 10°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 100 GeV
- Zenith angle: 20°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 100 GeV
- Zenith angle: 30°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 100 GeV
- Zenith angle: 40°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 100 GeV
- Zenith angle: 50°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 100 GeV
- Zenith angle: 60°


Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)
- Energy: 100 GeV
- Zenith angle: 70°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 0°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 10°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 20°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 30°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 40°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 50°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 60°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 0.1 GeV
- Zenith angle: 70°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 1 GeV
- Zenith angle: 0°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 1 GeV
- Zenith angle: 10°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 1 GeV
- Zenith angle: 20°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 1 GeV
- Zenith angle: 30°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 1 GeV
- Zenith angle: 40°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 1 GeV
- Zenith angle: 50°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 1 GeV
- Zenith angle: 60°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 1 GeV
- Zenith angle: 70°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 10 GeV
- Zenith angle: 0°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 10 GeV
- Zenith angle: 10°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 10 GeV
- Zenith angle: 20°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 10 GeV
- Zenith angle: 30°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 10 GeV
- Zenith angle: 40°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 10 GeV
- Zenith angle: 50°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 10 GeV
- Zenith angle: 60°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 10 GeV
- Zenith angle: 70°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 100 GeV
- Zenith angle: 0°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 100 GeV
- Zenith angle: 10°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 100 GeV
- Zenith angle: 20°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 100 GeV
- Zenith angle: 30°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 100 GeV
- Zenith angle: 40°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 100 GeV
- Zenith angle: 50°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 100 GeV
- Zenith angle: 60°


Results: Electrons
\(e^{-}\)
\(e^{+}\)
- Energy: 100 GeV
- Zenith angle: 70°


Results: Photons
\(\gamma\)
- Energy: 0.1 GeV
- Zenith angle: 0°

Results: Photons
\(\gamma\)
- Energy: 0.1 GeV
- Zenith angle: 10°

Results: Photons
\(\gamma\)
- Energy: 0.1 GeV
- Zenith angle: 20°

Results: Photons
\(\gamma\)
- Energy: 0.1 GeV
- Zenith angle: 30°

Results: Photons
\(\gamma\)
- Energy: 0.1 GeV
- Zenith angle: 40°

Results: Photons
\(\gamma\)
- Energy: 0.1 GeV
- Zenith angle: 50°

Results: Photons
\(\gamma\)
- Energy: 0.1 GeV
- Zenith angle: 60°

Results: Photons
\(\gamma\)
- Energy: 0.1 GeV
- Zenith angle: 70°

Results: Photons
\(\gamma\)
- Energy: 1 GeV
- Zenith angle: 0°

Results: Photons
\(\gamma\)
- Energy: 1 GeV
- Zenith angle: 10°

Results: Photons
\(\gamma\)
- Energy: 1 GeV
- Zenith angle: 20°

Results: Photons
\(\gamma\)
- Energy: 1 GeV
- Zenith angle: 30°

Results: Photons
\(\gamma\)
- Energy: 1 GeV
- Zenith angle: 40°

Results: Photons
\(\gamma\)
- Energy: 1 GeV
- Zenith angle: 50°

Results: Photons
\(\gamma\)
- Energy: 1 GeV
- Zenith angle: 60°

Results: Photons
\(\gamma\)
- Energy: 1 GeV
- Zenith angle: 70°

Results: Photons
\(\gamma\)
- Energy: 10 GeV
- Zenith angle: 0°

Results: Photons
\(\gamma\)
- Energy: 10 GeV
- Zenith angle: 10°

Results: Photons
\(\gamma\)
- Energy: 10 GeV
- Zenith angle: 20°

Results: Photons
\(\gamma\)
- Energy: 10 GeV
- Zenith angle: 30°

Results: Photons
\(\gamma\)
- Energy: 10 GeV
- Zenith angle: 40°

Results: Photons
\(\gamma\)
- Energy: 10 GeV
- Zenith angle: 50°

Results: Photons
\(\gamma\)
- Energy: 10 GeV
- Zenith angle: 60°

Results: Photons
\(\gamma\)
- Energy: 10 GeV
- Zenith angle: 70°

Results: Photons
\(\gamma\)
- Energy: 100 GeV
- Zenith angle: 0°

Results: Photons
\(\gamma\)
- Energy: 100 GeV
- Zenith angle: 10°

Results: Photons
\(\gamma\)
- Energy: 100 GeV
- Zenith angle: 20°

Results: Photons
\(\gamma\)
- Energy: 100 GeV
- Zenith angle: 30°

Results: Photons
\(\gamma\)
- Energy: 100 GeV
- Zenith angle: 40°

Results: Photons
\(\gamma\)
- Energy: 100 GeV
- Zenith angle: 50°

Results: Photons
\(\gamma\)
- Energy: 100 GeV
- Zenith angle: 60°

Results: Photons
\(\gamma\)
- Energy: 100 GeV
- Zenith angle: 60°

Results: Muons
\(\mu^{-}\)
\(\mu^{+}\)


Results: Update
\(\mu^{-}\)
\(\mu^{+}\)


Results: Update
\(\mu^{-}\)
\(\mu^{+}\)


Results: Update
\(\mu^{-}\)
\(\mu^{+}\)


Results: Update
\(\mu^{-}\)
\(\mu^{+}\)


Results: Update
\(\mu^{-}\)
\(\mu^{+}\)


Results: Update
\(\mu^{-}\)
\(\mu^{+}\)


Results: Update
\(\mu^{-}\)
\(\mu^{+}\)


Results: Update
\(\mu^{-}\)
\(\mu^{+}\)


Results: Muons
\(\mu^{-}\)
Fixed snow height: 0.379m
Results: Update
\(\mu^{+}\)
Results: Muons
Different snow heights

Fixed zenith angle
nEvents: 10223
\(\mu^{-}\)
Results: Muons
\(\mu^{-}\)

Results: Muons
\(\mu^{-}\)

Results: Muons
\(\mu^{-}\)

Results: Muons
\(\mu^{-}\)

Results: Muons
\(\mu^{-}\)

Results: Muons
\(\mu^{-}\)

Results: Muons
\(\mu^{-}\)

Results: Muons
\(\mu^{-}\)

Results: Muons
\(\mu^{-}\)

Results: Muons
\(\mu^{-}\)

Results: Muons
\(\mu^{-}\)

Results: Muons
\(\mu^{-}\)

Results: Muons
Fixed zenith angle: 0°
Results: Interpolation
Phase space simulations
Fixed steps: Energy, zenith angle, snow heights
Multidimensional linear interpolation
"Lookup tables"
Sampling distributions
Any energy, zenith angle and snow height
Results: Interpolation


Multidimensional linear interpolation
Preliminary results
Results: Interpolation

- Dark blue: interpolation results
- Light blue: simulations
Preliminary results
Third presentation: IceCube internship
By Sebastian Ordoñez
Third presentation: IceCube internship
- 404