An experimental magnetic field detector has been constructed outside of the town of Lakeview. Although the purpose of the detector is largely unknown to the townsfolk of Lakeview you and your team have been recruited to develop a magnetic field detector that is able to detect Hawkions. Hawkions are like muons, but are slow-moving and have long lifetimes. They are a newly discovered top secret particles that the data from the probe from Artemis 13 discovered is bombarding the Lakeview area. You have a somewhat constructed model, in which the Hawkions follow a straight line trajectory, but it looks like there are some pieces of code that the team wasn't sure what to do with. You will need to select a few locations to model the magnetic field due to the Hawkions and produce arrows that represent the Hawkion's magnetic field. Best hurry, the government needs more information about the Hawkions particle before it is too late.

## Scene setup
scene.background = color.white
 
## Parameters and Initial Conditions
velocity = vector(1,0,0)
 
## Objects
charge = sphere(pos=vector(-2,0,0), radius=0.1, color=color.blue)
xaxis = cylinder(pos=vector(-3,0,0), axis=vector(6,0,0), radius = 0.01, color=color.black)
yaxis = cylinder(pos=vector(0,-3,0), axis=vector(0,6,0), radius = 0.01, color=color.black)
zaxis = cylinder(pos=vector(0,0,-3), axis=vector(0,0,6), radius = 0.01, color=color.black) 
 
## Calculation Loop
 
t = 0
dt = 0.01
 
while t < 5:
 
    rate(100)
 
    charge.pos = charge.pos + velocity*dt
 
    t = t + dt
 
 
## Not sure what to do with these
 
##p = sphere(pos=vector(-1,-1,0), radius = 0.1, color=color.cyan) 
##Barrow = arrow(color=color.red)
##Barrow.pos = p.pos
##Barrow.axis = vector(0,0,0) 

Since your team was so successful at creating the magnetic field detector for the Hawkions, the Lakeviewians have now hired you to analyze other characteristics of the Hawkions. The Hawkions appear to have a variety of masses, so you start by trying to sort the Hawkions by mass for further analysis (using a mass spectrometer). Before building the spectrometer, your team is asked to develop a model of the spectrometer. The higher ups have specified that they would like to use a constant magnetic field of 2 T and to isolate the effects of the B-field within the magnetic housing box, but they are unsure what direction the magnetic field should be relative to the catching plate. Below is the beginning of a model that was developed by another team that mysteriously disappeared.

GlowScript 3.0 VPython
 
 
## Model parameters
B = vec(0,0,0)
 
## Objects
mag_housing =  box(pos=vec(1,0,0), length=2, height=2, width=2, opacity=0.2)
catchingplate = box(pos=vec(4,0,0), length=0.1, height=10, width=10, opacity=1)
xaxis = cylinder(pos=vector(-3,0,0), axis=vector(6,0,0), radius = 0.05, color=color.white)
yaxis = cylinder(pos=vector(0,-3,0), axis=vector(0,6,0), radius = 0.05, color=color.white)
zaxis = cylinder(pos=vector(0,0,-3), axis=vector(0,0,6), radius = 0.05, color=color.white)
 
 
## Set up particles
i = 0
N = 20
particleList = []
 
while i < N:
 
    particleList.append(sphere(pos=vec(-2-4*random(),0,0), radius = 0.1, m = 20*(random()+1), v = vec(2,0,0), q = 10, color=color.red))
    i = i + 1
 
## Calculation Loop
t = 0
dt = 0.01
 
while t < 30:
 
    rate(300)
 
    for thisParticle in particleList:
 
        thisParticle.pos = thisParticle.pos + thisParticle.v*dt
 
 
        if thisParticle.pos.x > catchingplate.pos.x:
 
            thisParticle.v = vector(0,0,0)
 
    t = t + dt