Making Herwig Samples

[lawrenceleejr]

Using Herwig steering file LEP.in with modifications to create a sample of ee to Z+Gamma events for 1 TeV CM collisions.

Running in docker with:

docker run --rm -ti -v $PWD:$PWD -w $PWD -e DISPLAY=$DISPLAY -v $XSOCK:$XSOCK -v $XAUTH:$XAUTH   -e XAUTHORITY=$XAUTH --net host hepstore/rivet-herwig:latest

Then running Herwig with:

Herwig read LEP.in

Then going into a delphes container and running it.

docker run --rm -ti -v $PWD:$PWD -w $PWD ghcr.io/scipp-atlas/mario-mapyde/delphes
/usr/local/share/delphes/delphes/DelphesHepMC2 /usr/local
/share/delphes/delphes/cards/delphes_card_CMS.tcl LEP-1000.root LEP-1000.hepmc

LEP.in:

# -*- ThePEG-repository -*-

##################################################
# Example generator based on LEP parameters
# usage: Herwig read LEP.in
##################################################

read snippets/EECollider.in

##################################################
# Change settings for the ee->Z->qq matrix element
# to avoid producing top quarks
#
# 'set' lines like this can be omitted if the
# default value is already okay.
#
# Any repository setting can be modified here
##################################################

##################################################
# Selected the hard process
##################################################
# leading-order processes
##################################################
cd /Herwig/MatrixElements
# default e+e- > q qbar (5 flavours d,u,s,c,b)
insert SubProcess:MatrixElements 0 MEee2gZ2qq
# e+e- > l+l-
# insert SubProcess:MatrixElements[0] MEee2gZ2ll
# e+e- > W+W- ZZ
# insert SubProcess:MatrixElements[0] MEee2VV
# higgs+Z
# insert SubProcess:MatrixElements[0] MEee2ZH
# higgs+e+e-/nu_enu_ebar via VBF
# insert SubProcess:MatrixElements[0] MEee2HiggsVBF
##################################################
# NLO powheg processes
##################################################
#set /Herwig/Shower/ShowerHandler:HardEmission POWHEG
#insert SubProcess:MatrixElements 0 PowhegMEee2gZ2qq
#insert SubProcess:MatrixElements 0 PowhegMEee2gZ2ll

##################################################
# LEP physics parameters (override defaults)
##################################################
cd /Herwig/Generators
set EventGenerator:EventHandler:LuminosityFunction:Energy 91.2

cd /Herwig/Generators
##################################################
## prepare for Rivet analysis or HepMC output
## when running with parton shower
##################################################
#read snippets/Rivet.in
#insert /Herwig/Analysis/Rivet:Analyses 0 XXX_2015_ABC123
read snippets/HepMC.in
#set /Herwig/Analysis/HepMC:PrintEvent -1
set /Herwig/Analysis/HepMC:PrintEvent 100000

###################################################
# Save run for later usage with 'Herwig run'
##################################################
saverun LEP EventGenerator

##################################################
# uncomment this section for an example batch run
# of two repeats with different parameters
#
# Note that a separate call of 'Herwig run'
# is not required
##################################################
set EventGenerator:NumberOfEvents 100000
#run LEP-Zpole EventGenerator
set EventGenerator:EventHandler:LuminosityFunction:Energy 1000.0
run LEP-1000 EventGenerator

[lawrenceleejr]

On SERF212a, can find the data here: /userdata/FragmentationStudy/data/Herwig/

[lawrenceleejr]

Calling this done because I think this is good enough for now.