# General settings for Pythia8 wrapper program for inclusive top
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 setup
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: t-channel single top
Top:qq2tq(t:W) = on


## Set matrix-element level cuts
# Not applied to ttbar since cuts would pertain to hard (ie ttbar) system
#PhaseSpace:pTHatMin = 0   # min pT
#PhaseSpace:pTHatMax = 13000   # max pT
#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:preferLHAPDF=2
Tune:pp = 0
# The "LHC" thermal tune from Table 1 in Fischer and Sjostrand
# Note hadron rescattering was changed, so here switching on new model
# with default parameters. An alternatiev without rescattering is also
# provided, to show differences. 
StringPT:thermalModel = on
StringPT:temperature = 0.21
StringFlav:BtoMratio = 0.357
StringFlav:StrangeSuppression = 0.5
StringPT:expNSP = 0.13
ColourReconnection:range = 1.1
HadronLevel:Rescatter = on
MultiPartonInteractions:pT0Ref = 2.5