I graduated from UNC in 2012 and went on to pursue a PhD in hadronic physics at MIT. After graduating, I decided to look more at neutrino physics, and, as part of that, I spent a winter at the South Pole in 2020. Since then I have been working as a post-doc on IceCube primarily investigating sterile neutrinos with a focus on the relevant statistics.
EMAIL: jmhardin@mit.edu
This talk will describe sterile neutrinos and the way we discover new things. I will also briefly mention how IceCube runs at the South Pole.
Guest Speaker John Hardin (20 minutes + 10 min Q&A)
Dr. Shirey's Week 4 slides.
Your feedback so far.
Break (5 minutes)
Next steps in the design challenge:
Indicate your favorite design idea on your Phase 2: Design Exploration document (Google Doc copy, Microsoft Word download, PDF download).
Email or share your Phase 2 Ideas doc to Dr. Shirey (katey@edukatey.com) by Sunday, February 13.
Our Python experts will review, gather resources, comment back to you.
Complete the Week 4 assignments of CodeHS.com by Sunday, February 13 (to give Dr. Shirey time to process your answers and steer next steps.)
Last year, Dr. Xinhua Bai spoke to IceCube Afterschool Interns. His lecture was on cosmic ray experiments, particle flux and detector geometric acceptance (lecture recording here). He also provided the Python code below which calculates potential particle flux based on zenith angle. Need-to-knows: Geometry, trigonometric functions, introductory-level calculus, online python compiler such as https://www.tutorialspoint.com/execute_python_online.php
Flux sandbox to fork: https://codehs.com/sandbox/id/flux-sandbox-7ZR5Qf
# https://www.tutorialspoint.com/execute_python_online.php
# https://repl.it/languages/python3
import os
import sys
import random
import numpy as np
import math
import matplotlib.pyplot as plt
pi = 3.1415926
NLoop = 15
zenith = 0.0 # zenith angle in degree
zenithRad = 0.0 # zenith angle in rad
area = 3000.0 # surface array area = 3000.0 km^2
#area = 1.0 # surface array area = 3000.0 km^2
Acceptance = 0.0
AcceptanceOld = 0.0
print ('+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++')
print ('+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++')
print (' Maximum Angle Integrated Acceptance Increment')
print (' Theta (degree) [Km^2 sr] [Km^2 sr]')
for n in range (NLoop): # loop to NLoop*5 degrees
zenith = (n+1)*5.0
#print '\n zenith angle = ', zenith, ' (degree)'
zenithRad = pi*zenith/180.0
Acceptance = area*pi*math.sin(zenithRad)*math.sin(zenithRad)
# print zenith, Acceptance, Acceptance-AcceptanceOld # messy output
print (" %8.2f %8.2f %8.2f" % (zenith, Acceptance, Acceptance-AcceptanceOld))
AcceptanceOld = Acceptance
print ('+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++')
print ('++++++++++++++++++++++++++++1111+++++++++++++++++++++++++++')