Week 3 Guest Speaker: Marjon Moulai

WEEK 3: TALK TITLE

Week 3 Documents and REsources:

1. Week 3 slides from Dr. Shirey

Assignment #3:

• Loops learning module assigned (it’s labeled #3 but we’ll do it second)

• Complete the Week 3 assignments of CodeHS.com by Sunday, February 6 (to give Dr. Shirey time to process your answers and steer next steps.)

• Get started on the design challenge Phase 2: Design Exploration (Google Doc copy, Microsoft Word download, PDF download)

NEUTRINO FLUX EXAMPLE

(PYTHON SIMULATION)

• Flux example using Python: Use a python script to calculate the geometric acceptance of IceTop and Pierre Auger surface array. The Python text file is below.

• Need-to-knows: Geometry, trigonometric functions, introductory level calculus, online python compiler such as https://www.tutorialspoint.com/execute_python_online.php

• Do: Open the compiler and use the Python script below to calculate the total neutrino flux. Change the Zenith to adjust the angle between IceCube's surface and incoming neutrinos.

• Reflect: Besides the geometric acceptance, what other factors need to be considered in measuring cosmic ray flux using surface array?

# 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+++++++++++++++++++++++++++')