In this section we will simulate the process 'P + P -> Z+X' at 7 TeV center-of-mass energy and
study the properties of the Z-boson with help of the [[https://rivet.hepforge.org/analyses/MC_ZINC.html | MC_ZINC ]] Rivet analysis.
There are different aspects of event generation, each of which plays a crucial role in the correct simulation of Z-bosons at the LHC:
* Hard matrix element
* Parton Shower
* Initial and final state radiation
* Multi parton interactions
* Hadronizaion
During the course of this exercise we will see the different contributions to the Z shape observables, like the Z pT spectrum.
Please download the following input file and move it to your work directory:
Input file:{F543842}
Alternatively, you can use the following command-line codes:
```
```
The input file can be used for Z production and Z+Jets production. For now, all parts except for the hard matrix element are commented out. We will uncomment the necessary parts as we progress with the exercise.
Run the input file without any modifications
```
Herwig read LHC-Z.in
Herwig run LHC-Z.run -N 10000
```
It is recommended to change the name of the newly created .yoda file after every new run to something meaningful s.t. you don't have to repeat runs in the end when we compare the runs with each other.
For example:
mv LHC-Z.yoda LHC-noShower.yoda
Now plot the results with
```
rivet-mkhtml LHC-noShower.yoda
```
Have a look at the produced histograms and try to explain why there is no Z pT and why the Z pT in the peak region basically is just a delta function. You should also have a look at the .log file and understand the first few events.
To have a more realistic simulation of the Z-boson we are going to uncomment some sections in the input file and include the LO Parton shower without initial and final state radiation.
```
##################################################
## ShowerHandler(s)
##################################################
# Initial and final state radiation
set /Herwig/Shower/ShowerHandler:DoFSR No
set /Herwig/Shower/ShowerHandler:DoISR No
cd /Herwig/EventHandlers
#set EventHandler:CascadeHandler NULL
set EventHandler:CascadeHandler:MPIHandler NULL # switches off multi parton interactions
set EventHandler:DecayHandler NULL # switches off particle decays
set EventHandler:HadronizationHandler NULL # switches off hadronization
```
Now run Herwig again and compare your results with the first run. Why do the pT spectra change but the reconstructed mass of the Z boson looks roughly the same?
Activate initial and final state radiation and repeat the run.
```
set /Herwig/Shower/ShowerHandler:DoFSR Yes
set /Herwig/Shower/ShowerHandler:DoISR Yes
```
For the next task we are making a jump to a next-to-leading order parton shower
Have a look at the parton shower selection section from the input file
```
##################################################
## Matching and shower selection
##################################################
#read Matchbox/MCatNLO-DefaultShower.in
# read Matchbox/Powheg-DefaultShower.in
## use for strict LO/NLO comparisons
# read Matchbox/MCatLO-DefaultShower.in
## use for improved LO showering
# read Matchbox/LO-DefaultShower.in
# read Matchbox/MCatNLO-DipoleShower.in
# read Matchbox/Powheg-DipoleShower.in
## use for strict LO/NLO comparisons
# read Matchbox/MCatLO-DipoleShower.in
## use for improved LO showering
# read Matchbox/LO-DipoleShower.in
# read Matchbox/NLO-NoShower.in
read Matchbox/LO-NoShower.in
```
In the input file we choose the "read Matchbox/MCatLO-DefaultShower.in" option activating an MC@LO angular ordered parton shower. Now please choose the "read Matchbox/MCatNLO-DefaultShower.in" option, to choose the MC@NLO matched shower and run the event again.
By now you should have a sheer amount of .yoda files created which you can compare with
rivet-mkhtml *.yoda
For the final task in this section activate multi parton interactions, decays and hadronization by uncommenting the relevant lines in the input file and compare the leading order partonshower with the next-to-leading order parton shower simulation.
By now you have finished all the obligatory parts of the Herwig tutorial.
You can choose to [[https://herwig.hepforge.org/tutorials/gettingstarted/matrixelements.html#some-scale-variations | dive deeper]] into parton shower variations, play around with different cuts and scales
or continue with [[https://phab.hepforge.org/w/mcnet2019/thirdrun/ | Third run: minimum bias events and soft physics]] (recommended).
==A note on matrix elements==
As you can see, the process selection section has been completely deactivated. This is because for a process like `P + P -> Z + X`, Herwig needs to use the external matrix element (ME) providers. However, the use of these tools in the docker/VM would be problematic due to some linking issues. In an installed (and properly linked) version of Herwig, one can easily use these external ME providers by selecting the hard-process in the `Process selection` and choosing the ME providers from the `Matrix element library selection`.
```
##################################################
## Process selection
##################################################
## Model assumptions
# read Matchbox/StandardModelLike.in
# read Matchbox/DiagonalCKM.in
## Set the order of the couplings
# cd /Herwig/MatrixElements/Matchbox
# set Factory:OrderInAlphaS 1
# set Factory:OrderInAlphaEW 1
## Select the process
# do Factory:Process p p -> Z0 j
##################################################
## Matrix element library selection
##################################################
## Select a generic tree/loop combination or a
## specialized NLO package
# read Matchbox/MadGraph-GoSam.in
# read Matchbox/MadGraph-MadGraph.in
# read Matchbox/MadGraph-NJet.in
# read Matchbox/MadGraph-OpenLoops.in
```