Here we demonstrate an example of photon-induced events in Herwig.
==Using the Les Houches Event Handler==
To start, create a new work directory:
mkdir photon_induced
cd photon_induced
Next, download below input/config/analysis files in your work directory.
1. MadGraph configure script: this is a script for generating a sample of unweighted yy > ll events with MadGraph.
{F1993095}
wget https://phab-files.hepforge.org/file/download/fmev7v6idvr75ggbdth3/PHID-FILE-i3aaqziv52hgqrsztjuo/MG_make.sh
2. MadGraph parameter and run cards: parameter space and run statistics for the MadGraph run.
{F1993058} {F1993056}
wget https://phab-files.hepforge.org/file/download/xjefkopi6whfilhore3q/PHID-FILE-kfvne4b34gu6thn2syt5/run_card.dat
wget https://phab-files.hepforge.org/file/download/6ox3c6if7ekcawjw33am/PHID-FILE-ynhnj545mdxfbtqytoqn/param_card.dat
3. Rivet analysis and plot setup: Rivet analysis, to be run alongside the Herwig setup.
{F1993057} {F1993059}
wget https://phab-files.hepforge.org/file/download/drv32ihccvchwfm6qrjx/PHID-FILE-74aubf5x2iktkajdepci/RivetAnalysis_yy.cc
wget https://phab-files.hepforge.org/file/download/mjsx7nxjbmu3curex3ew/PHID-FILE-p3xj34mwaixlritunlmp/RivetAnalysis_yy.plot
4. Herwig Input files: the main Herwig input file.
{F1993096}
wget https://phab-files.hepforge.org/file/download/ccjmisiezgugnatpxbqf/PHID-FILE-66o5cobqy5u3fydkyr7u/LesHouches.in
Once the files are downloaded, please have a look at **MG_make.sh**. It can be run by:
chmod +x MG_make.sh
[sudo] ./MG_make.sh
Once it is done, you should be able to see a new directory named "yyll" that contains the requited unweighted sample in **yyll/Events/run_01/unweighted_events.lhe.gz**.
We are now set to sun Herwig using these samples. Run
Herwig read LesHouches.in
It should generate a run file called **LesHouches.run**. Now execute
Herwig run LesHouches.run -N 10000
The last part asks Herwig to run this simulation for 10000 events only. After this is done, you should be able to see several files generated in your work directory. This will generate a number of files in your work directory:
- **LesHouches.log**: contains the details of each produced event, for the first 10 events (can be changed to include any number of events). The details include incoming/outgoing/intermediate particles and their kinematics. This is useful information to check what is actually going on under the hood.
- **LesHouches.out**: Statistics for the event handler, including the number of valid events (subtracting vetoed/scrapped events) and a list of sub-process categories with their event count and cross-section.
- **LesHouches.tex**: A list of publications, relevant to your Herwig setup, which you should cite if you are making use of this specific simulation.
- **LesHouches.yoda**: Specific results, according to the Rivet analysis you called.
To produce the plots, execute
rivet-mkhtml LesHouches.yoda -o plots
where the **-o** flag allows you to choose an output directory. Have a look at the generated plot by double-clicking on **plots/index.html**. If you are not happy with the statistics, you can improve them by running more events. This might take longer to run, but you improve your run-time by multi-threading your simulation:
Herwig run LesHouches.run -N 100000 -j 5
Here, the **-j** flag tells Herwig to split this simulation into 5 parallel positions, each running 20000 events. The results will be saved in 5 different sets:
LesHouches-1.yoda
LesHouches-2.yoda
LesHouches-3.yoda
LesHouches-4.yoda
LesHouches-5.yoda
You can combine these by
yodamerge LesHouches-*.yoda -o LesHouches-ALL.yoda
==Using Matchbox==
Another way of performing this simulation is by using Herwig's Matchbox module, which handles the matrix element generation internally. To try this, download the following input file in your work directory
{F1993789}
wget https://phab-files.hepforge.org/file/download/ipilkrnjdvoj7ynycuh4/PHID-FILE-flz6rphk35jvgu454st3/yyll-Matchbox.in
Have a look inside this file and try to understand the structure and the commands. You can run it by
Herwig read yyll-Matchbox.in
Herwig run yyll-Matchbox.run -N 1000000 -j5
yodamerge yyll-Matchbox-*.yoda -o yyll-Matchbox-ALL.yoda
Now, try to plot LesHouches and Matchbox results together:
rivet-mkhtml LesHouches-ALL.yoda:"Les Houches" Matchbox-ALL.yoda:"Matchbox" -o ALL-plots
**:"Les Houches"** and **:"Matchbox"** are optional tags for the plots. You can either include them or not. Have a look at the results by running **ALL-plots/index.html**.
== Exercise ==
1- Change the configurations in **yyll-Matchbox.in* file:
read Matchbox/LO-DefaultShower.in
# read Matchbox/LO-DipoleShower.in
# read Matchbox/LO-NoShower.in
Uncomment each line separately to change the parton shower setting from Default (Angularly-ordered) to Dipole and then turn it off and produce yoda files. Plot them together and check the difference.
2 - Herwig-7.3 supports a full QCD+QED+EW parton shower. You can control the parton shower meta by adding the command
set /Herwig/Shower/ShowerHandler:Interactions ALL # options are QCD QED QEDQCD EWOnly ALL
Alternate between these options in both Les Houches and Matchbox runs and plot them against each other and observe the differences.