The code here is designed to calculate the properties of the electric field producd from electromagnetic showers in Ice. I have so far implemented the 1-D parameterized form of calculating the field as well as a simple analytic approximation around the Cherenkov angle. The goal is to implement the full 3-D parameterization where one uses the individual tracks in an event.

##########################################

##########################################

Physics -- This code houses some useful physics constants.

Dim1 -- This piece has all the methods needed to perform the 1D analysis of a shower shop. It takes as input a vector which is supposed to be representative of the shower profile. the entries should be N(e-p) at each step in radiation length. It can also implement the simple analytic approximation around the Cherenkov Angle

Analysis -- This class is where the analysis is actually done. It outputs the necessary histograms and figures that are useful to study the properties. The user could use this class and add methods, or use it as an example on how to use Dim1 (and later Dim3).

TreeAnalyzer -- In order to do the event by event analysis, I figured I would need convenient access to the Geant4 data. So I wrote a very basic flat tree that has basic variables needed to analyze the shower profile.

##########################################

##########################################

After compiling, the executable 'analysis' will appear. There are three options the user can use to run the program:

-a -- This will make the Rx|E| vs Theta using the average profile produced from TrackAna.py -b -- This will make Rx|E|vs Theta using the simple expression valid around the Cherenkov angle. -c -- Make Rx|E| vs Theta using event-by-event profiles that are calculated from a simple TTree. In IceBlock, the code is makeNtuple.py It contains very simple structure.