MSEL=0
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
MDME(174,1)=0            !Z decay into d dbar
MDME(175,1)=0            !Z decay into u ubar
MDME(176,1)=0            !Z decay into s sbar
MDME(177,1)=0            !Z decay into c cbar
MDME(178,1)=0            !Z decay into b bbar
MDME(179,1)=0            !Z decay into t tbar
MDME(180,1)=0            !Z decay into b' b'bar
MDME(181,1)=0            !Z decay into t' t'bar
MDME(182,1)=1            !Z decay into e- e+
MDME(183,1)=0            !Z decay into nu_e nu_ebar
MDME(184,1)=1            !Z decay into mu- mu+
MDME(185,1)=0            !Z decay into nu_mu nu_mubar
MDME(186,1)=1            !Z decay into tau- tau+
MDME(187,1)=0            !Z decay into nu_tau nu_taubar
MDME(188,1)=0            !Z decay into tau'- tau'+
MDME(189,1)=0            !Z decay into nu'_tau nu'_tau
MSUB(15)=1               !qq gW
MSUB(30)=1               !qg qZ
CKIN(1)=60.      !lower limit for CM energy to ease computation
CKIN(3)=15.      !(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