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)=1            !Z decay into nu_e nu_ebar
MDME(184,1)=1            !Z decay into mu- mu+
MDME(185,1)=1            !Z decay into nu_mu nu_mubar
MDME(186,1)=1            !Z decay into tau- tau+
MDME(187,1)=1            !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
MDME(190,1)=0            !W decay into dbar u
MDME(191,1)=0            !W decay into dbar c
MDME(192,1)=0            !W decay into dbar t
MDME(194,1)=0            !W decay into sbar u
MDME(195,1)=0            !W decay into sbar c
MDME(196,1)=0            !W decay into sbar t
MDME(198,1)=0            !W decay into bbar u
MDME(199,1)=0            !W decay into bbar c
MDME(200,1)=0            !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
MSUB(23)=1               !ZW
MSUB(70)=1               !rW->ZW
MSUB(73)=1               !ZW->ZW(longitudinal)