diff --git a/README.rst b/README.rst index 68cc90b100fb212213975f1a6147f44d3246ccf0..d98778dd5afed0a34df80e6a50928885f5f164ef 100644 --- a/README.rst +++ b/README.rst @@ -64,6 +64,27 @@ This specific command runs all jobcards within the ``jobcards/examples`` folder and stores the output inside the folder ``output``. The folder structure is applied from the ``jobcards``\ folder. +Theory +------ + +In order to retrieve correct results and provide correct KM3NeT weights (w2) +the treatment of the GiBUU weights is an important step. A brief description +of the GiBUU weights and how to calculate actual cross sections is given on the +`GiBUU Homepage<https://gibuu.hepforge.org/trac/wiki/perWeight>`__ and +a more detailed description of the calculation can be found in the `PhD Thesis +of Tina Leitner<https://inspirehep.net/literature/849921>`__ in Chapter 8.3. +As it is mentioned in the description of the output flux file in the +`documentation<https://gibuu.hepforge.org/Documentation/code/init/neutrino/initNeutrino_f90.html#robo1685>`__ this is not taken somehow into account inside the weights. +Following the description the GiBUU event weight can be converted to a binned +cross section via + +.. math:: + \frac{d\sigma}{E} = \frac{\sum_{i\in I_\text{bin}} w_i}{\Delta E}\cdot\frac{1}{E\Phi}, + +where :math:`\Phi`__ is the simulated flux +As the weights are given for each run individually the weight also has to be divided +by the number of runs. + Tutorial -------- The python framework is build around the GiBUU workflow, i.e. a jobcard is