kopecna/EWP-BplusToKstMuMu-AngAna
kopecna/EWP-BplusToKstMuMu-AngAna
2
1
Fork 0
Code Issues Pull Requests repo.project_board Releases Wiki Activity

Update 'Event generator'

Renata Kopecná 2022-02-10 14:08:30 +01:00
parent 6e6ae21382
commit 9786b409df
1 changed files with 27 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

40
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
```