MSEL=6   
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(190,1)=1            !W decay into dbar u
MDME(191,1)=1            !W decay into dbar c
MDME(192,1)=1            !W decay into dbar t
MDME(194,1)=1            !W decay into sbar u
MDME(195,1)=1            !W decay into sbar c
MDME(196,1)=1            !W decay into sbar t
MDME(198,1)=1            !W decay into bbar u
MDME(199,1)=1            !W decay into bbar c
MDME(200,1)=1            !W decay into bbar t
MDME(206,1)=1            !W decay into e+ nu_e
MDME(207,1)=1            !W decay into mu+ nu_mu
MDME(208,1)=1            !W decay into tau+ nu_tau