INSTRUCTIONS ON USING THE SIMULATION:

Execute prgrm.py from the terminal and follow the prompts.
For the number of divisions, a larger number will give a more precise answer,
however, the simulation will take longer, the larger the product of these
numbers are.

The program will return the total amount of photoelectrons created by the
photocathode, as well as creating a log file in the appropriately named folder.
The resulting log file will be of the following format:

photon incidence x-coordinate	photon incidence y-coordinate	photon incidence time	photon emission x-coordinate   photon emission z-coordinate  photon emission time   omega angle   zeta angle   wavelength   photon distance   particle distance   probability

photon incidence x-coordinate and y-coordinate give where a particular photon created by cherenkov radiation would hit the
photocathode, in micrometers, with the incident particle always entering the glass at x=y=0

photon emission x-coordinate and z-coordinate give where a particular photon was created. y-coordinate set to 0.

photon incidence time is the calculated interval between the incident particle entering the glass and the photon 
hitting the photocathode, in picoseconds

photon emission time is the calculated interval between the incident particle entering the glass and the photon 
being created, in picoseconds

wavelength is the wavelength of photon in nanometers

for omega and zeta angle, see the coordinate system sketch

photon distance is the distance the photon travelled from its creation to the photocathode, in micrometers

particle distance is the distance the particle travelled between entering the glass and creating this particular
photon, in micrometers

probability is the percentage chance of this particular photon creating a photoelectron


INSTRUCTIONS ON MODIFYING LIST OF WINDOWS:

Adding a new photocathode follows from these instructions easily
1) Open sim.py
2) Under the current list of windows, add
 
new_window = {"name"    : "Name of your new window",
                "d"         :  0.011,
                "b"         : Sellmeier coefficients part 1
                "c"         : and part 2,
                "file"      : location of your window datafile} 
                
3) Below it, find window_list and add new_window to it:

window_list = [new_window, fused_silica, fused_silica2, fused_silica5, fused_silica8, fused_silica11, MgF2]
photocathode_list = [QE58]

