From a33770dcef3f09a62f3bce7df62eba001423e116 Mon Sep 17 00:00:00 2001
From: Johannes Schumann <johannes.schumann@fau.de>
Date: Wed, 3 Nov 2021 13:33:35 +0100
Subject: [PATCH] First draft for the weight section in the Readme

---
 README.rst | 21 +++++++++++++++++++++
 1 file changed, 21 insertions(+)

diff --git a/README.rst b/README.rst
index 68cc90b..d98778d 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 
-- 
GitLab