From 9786b409df11ed268d0790809d8fc08ea8f5e3e4 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Renata=20Kopecn=C3=A1?= <kopecna@physi.uni-heidelberg.de>
Date: Thu, 10 Feb 2022 14:08:30 +0100
Subject: [PATCH] Update 'Event generator'

---
 Event-generator.md | 40 +++++++++++++++++++++++++++-------------
 1 file changed, 27 insertions(+), 13 deletions(-)

diff --git a/Event-generator.md b/Event-generator.md
index 9528b34..e55e7a6 100644
--- a/Event-generator.md
+++ b/Event-generator.md
@@ -1,4 +1,10 @@
-https://lhcb.github.io/starterkit-lessons/second-analysis-steps/simulation-running-gauss.html
+## Table of Contents
+
+* [Building the environment](#building-the-environment)
+* [Run without detector simulation](#run-without-detector-simulation)
+* [Convert .sim files to .root files](#convert-sim-files-to-root-files)
+
+
 # Building the environment
 
 ```
@@ -6,27 +12,35 @@ SetupProject DaVinci/v41r2
 SetupProject Gauss/v49r7
 ```
 # Run without detector simulation
-Output is .sim file
+In order to create the generator level sample, it is useful to install [Gauss](https://lhcb.github.io/starterkit-lessons/second-analysis-steps/simulation-intro.html).  I used v49r12. 
+
+In Gauss-Job_StandAlone.py, the eventy type number is defined. According to the Monte Carlo type (either signal MC, reference MC or PHSP MC), corresponding job option is necesarry. One can choose from
+
+* `JobStandAlone_KpPi0_PHSPMC.py`
+* `JobStandAlone_KpPi0_ReferenceMC.py`
+* `JobStandAlone_KpPi0_SignalMC.py`
+* `JobStandAloneKSpiSignalMC.py`
+* `JobStandAloneKSpiTightCut.py`
 
 
-for PHSP MC:
+Then, order to get eg. the PHSP MC, launch Gauss with:
+
 ```
-gaudirun.py Gauss-Job_StandAlone.py JobStandAloneKSpiTightCut.py
+./run gaudirun.py Gauss-Job_StandAlone.py JobStandAlone_KpPi0_PHSPMC.py 
 ```
 
-for Signal MC:
-```
-gaudirun.py Gauss-Job_StandAlone.py JobStandAloneKSpiSignalMC.py $APPCONFIGOPTS/Gauss/Beam6500GeV-md100-2016-nu1.6.py $DECFILESROOT/options/12115102.py $LBPYTHIA8ROOT/options/Pythia8.py
-gaudirun.py Gauss-Job_StandAlone.py JobStandAloneKSpiSignalMC.py
-```
 
 For more details see
-[https://wiki.physik.uzh.ch/lhcb/gauss:gauss](https://wiki.physik.uzh.ch/lhcb/gauss:gauss)
+[https://wiki.physik.uzh.ch/lhcb/gauss:gauss](https://wiki.physik.uzh.ch/lhcb/gauss:gauss) and [the Starterkit lesson](https://lhcb.github.io/starterkit-lessons/second-analysis-steps/simulation-running-gauss.html).
+
+
+Output is `.sim` file. Therefore we need to convert it to a ROOT file.
 
 
 # Convert .sim files to .root files
 
-The .sim input file is **hardcoded** in the tupleResultKSpi_StandAlone_TightCut.py!
+The .sim input file is **hardcoded** in the getTupleFromGauss.py. 
+
+```
+lb-run DaVinci/v44r3 gaudirun.py getTupleFromGauss.py 
 ```
-lb-run DaVinci/v41r2 gaudirun.py tupleResultKSpi_StandAlone_TightCut.py
-```
\ No newline at end of file