Blame view

README.md 1.94 KB
5bee6696   Thomas Fitoussi   Create README.md
1
2
3
4
5
6
7
8
9
10
# 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
8ffbee93   Thomas Fitoussi   Analytic expressi...
11
12
13
14
15
16
17
18
19
20
21

## 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:
1db105d4   Thomas Fitoussi   Update README.md
22

8ffbee93   Thomas Fitoussi   Analytic expressi...
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
| 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 theta^2 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.