Simulation analysis

This a compilation of scripts to analyse cascade-simulation results.

These scripts are written in Python 2.7. They also requires:

• Numpy/Scipy
• Matplotlib

Instruction to install it could be found here: http://www.scipy.org/install.html

Principle

There is interdependancy between modules. Don't copy one and not the others otherwise they will be errors.

In the same time, the function in each module are written in way being used independantly, just import it and python will do the rest. For example, to use the spectrum, use just add to your script: "from Modules.Spectrum import spectrum".

Few modules are available:

Modules	functions
"Read.py"	extract data from the files simulations
"Constants.py"	same constants as in "cascade-simulation"
"Analytic.py"	plot analytic comparisons on graph or estimate a value
"Integrand.py"	some integrands used in "cascade-simulation"
"Spectrum.py"	extract and plot the spectrum
"timing.py"	extract and plot delta_t/t
"Angle.py"	extract and plot theta2 in histogram or versus energy
"Map.py"	plot map of the arrival photon, points or histo2d
"Generation.py"	plot historgram of the arrival particles generation

N.B.: "draw..." function call a table of files because it allows to plots more than one figure on the same graph, rather than the rest of the function can call only one file at a time.

For example: considering two files "file1" and "file2"

• you can plot both together: drawSpectrum([file1,file2]) or only one drawSpectrum([file1]). Always use the brackets
• but you can only extract one spectrum at a time: spectrum(file1). No brackets.