# settings of Pythia 8 wrapper program for Z(->ee,mumu) + 1 jet
Main:numberOfEvents = 100000          ! number of events to generate
Next:numberShowEvent = 0           ! suppress full listing of first events

# random seed
Random:setSeed = on
Random:seed = 48

# Beam parameter settings.
Beams:idA = 2212                ! first beam, p = 2212, pbar = -2212
Beams:idB = 2212                ! second beam, p = 2212, pbar = -2212
Beams:eCM = 13000               ! CM energy of collision

# Process setup: Drell-Yan + 1 jet 
# (note: Z+photon not included in these runs; must be revised if future analyses are included that do not veto Z+photon)
# Switch on gamma*/Z + jet
WeakBosonAndParton:qqbar2gmZg = on 
WeakBosonAndParton:qg2gmZq = on    
WeakBosonAndParton:fgm2gmZf = on   ! Photon-initiated (only if PDFs contain photons)
# Switch off all Z decays
23:onMode = off
# Switch back on those to electrons and muons
23:onIfAny = 11 13

# Set cuts
# Use this for hard leading-jets in a certain pT window (at ME level)
PhaseSpace:pTHatMin = 40   # min pT
PhaseSpace:pTHatMax = 13000   # max pT

# Note: mHat cuts would apply to Z+jet invariant mass, not lepton pair; not used
# PhaseSpace:mHatMin = 0   # min mHat
# PhaseSpace:mHatMax = 13000   # max mHat

# Makes particles with c*tau > 10 mm stable:
ParticleDecays:limitTau0 = On
ParticleDecays:tau0Max = 10.0

# Tune setup:
Tune:pp = 0

# The "Ropes" tune by Christian Bierlich, details:
# https://gitlab.com/Pythia8/releases/-/issues/80#note_667759779

# Parameter of the MPI model to keep total multiplicity reasonable
MultiPartonInteractions:pT0Ref = 2.15

# Parameters related to Junction formation/QCD based CR
BeamRemnants:remnantMode = 1
BeamRemnants:saturation = 5
ColourReconnection:mode = 1
ColourReconnection:allowDoubleJunRem = off
ColourReconnection:m0 = 0.3
ColourReconnection:allowJunctions = on
ColourReconnection:junctionCorrection = 1.2
ColourReconnection:timeDilationMode = 2
ColourReconnection:timeDilationPar = 0.18

# Enable rope hadronization
Ropewalk:RopeHadronization = on

# Also enable string shoving, but don't actually do anything.
# This is just to allow strings to free stream until hadronization
# where the overlaps between strings are calculated.
Ropewalk:doShoving = on
Ropewalk:tInit = 1.5 # Propagation time
Ropewalk:deltat = 0.05
Ropewalk:tShove = 0.1
Ropewalk:gAmplitude = 0. # Set shoving strength to 0 explicitly

# Do the ropes.
Ropewalk:doFlavour = on

# Parameters of the rope model
Ropewalk:r0 = 0.5 # in units of fm
Ropewalk:m0 = 0.2 # in units of GeV
Ropewalk:beta = 0.1

# Enabling setting of vertex information is necessary
# to calculate string overlaps.
PartonVertex:setVertex = on
PartonVertex:protonRadius = 0.7
PartonVertex:emissionWidth = 0.1