MSEL=1         !(D=1) to select between full user control (0, then use MSUB)
PMAS(5,10)=4.2          !mass of b quark
PMAS(6,1)=175.         !mass of top quark
PMAS(23,1)=91.187      !mass of Z
PMAS(24,1)=80.4        !mass of W
PMAS(25,1)=1000        !mass of Higgs
MSTJ(11)=3      !Choice of the fragmentation function
MSTJ(22)=2      !Decay those unstable particles
MSTP(2)=1       !which order running alphaS
MSTP(33)=0      !(D=0) inclusion of K factors in hard cross sections 
MSTP(51)=7      !structure function chosen
MSTP(81)=1      !multiple parton interactions 1 is Pythia default
MSTP(82)=4      !Defines the multi-parton model
MSTU(21)=1      !Check on possible errors during program execution
PARP(82)=1.9    !pt cutoff for multiparton interactions
PARP(89)=1000.  !sqrts for which PARP82 is set
PARP(84)=0.4    !Multiple interactions: matter distribution parameter
PARP(90)=0.16   !Multiple interactions: rescaling power
CKIN(1)=60.     !lower limit for CM energy to ease computation
CKIN(3)=350.    !(D=0. GeV)  (CKIN 3,CKIN 4) = range of allowed values for pt hat 
CKIN(4)=-1.     !(D=-1. GeV) maximum pt hat for hard interactions GeV/c